@article{21408,
  abstract     = {Rational design strategies for self-assembly require a detailed understanding of both the equilibrium state and the assembly kinetics. While the former is starting to be well understood, the latter remains a major theoretical challenge, especially in programmable systems and the so-called semi-addressable regime, where binding is often nondeterministic and the formation of off-target structures negatively influences the assembly. Here, we show that it is possible to simultaneously sculpt the assembly outcome and the assembly kinetics through the underexplored design space of binding energies and particle concentrations. By formulating the assembly process as a complex reaction network, we calculate and optimize the tradeoff between assembly speed and quality and show that parameter optimization can speed up assembly by many orders of magnitude without lowering the yield of the target structure. Although the exact speedup varies from design to design, we find the largest speedups for nondeterministic systems where unoptimized assembly is the slowest, sometimes even making them assemble faster than optimized, fully addressable designs. Therefore, these results not only solve a key challenge in semi-addressable self-assembly but further emphasize the utility of semi-addressability, where designs have the potential to be faster as well as cheaper (fewer particle species) and better (higher yield). More broadly, our results highlight the importance of parameter optimization in programmable self-assembly and provide practical tools for simultaneous optimization of kinetics and yield in a wide range of systems.},
  author       = {Hübl, Maximilian and Goodrich, Carl Peter},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {8},
  publisher    = {AIP Publishing},
  title        = {{Simultaneous optimization of assembly time and yield in programmable self-assembly}},
  doi          = {10.1063/5.0304731},
  volume       = {164},
  year         = {2026},
}

@article{21748,
  abstract     = {Cells are defined by lipid membranes that differ in their structure across the tree of life. While the membranes of most bacteria and eukaryotes consist of single-headed bilayer lipids, the membranes of archaea are composed of mixtures of single-headed bilayer lipids and double-headed bolalipids. Archaeal bolalipids can adopt straight or u-shaped conformations, enabling them—together with bilayer lipids—to control whether membranes form bilayer or monolayer structures. Yet, the physical principles governing archaeal membranes remain largely unexplored, especially how membrane structure couples to externally imposed curvature during membrane remodeling. Here, we perform coarse-grained molecular dynamics simulations of toroidal vesicles to systematically probe the effects of all relevant combinations of mean and Gaussian curvatures on shape stability and lipid organization. We find that soft bilayer membranes can sustain all curvatures induced, whereas rigid bolalipid monolayer membranes either transition to different vesicle shapes or rupture. Bilayer-mimicking u-shaped bolalipids and bilayer lipids are spatially accumulated in regions of high mean membrane curvature independent of Gaussian curvature. Our work identifies curvature–composition coupling as a physical signature of archaeal membrane remodeling.},
  author       = {Frey, Felix F and Santana de Freitas Amaral, Miguel and Šarić, Anđela},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {14},
  publisher    = {AIP Publishing},
  title        = {{Cracking donuts and sorting lipids: Geometry controls archaeal membrane stability and lipid organization}},
  doi          = {10.1063/5.0325170},
  volume       = {164},
  year         = {2026},
}

@article{21840,
  abstract     = {The transport properties of nanofluidic channels are usually studied under constant (DC) voltage or pressure driving. However, the frequency response under sinusoidal (AC) drivings offers rich insights into the time-dependent transport mechanisms. Inspired by recent electrochemical approaches, we investigate the couplings between ionic and electronic transport under AC driving. We show that conduction electrons of the channel walls participate in ionic current via capacitive electrochemical coupling, defining a critical frequency and length scale where electron-dominated conductivity emerges. We further analyze how electron–ion coupling modifies electro-osmotic flows and demonstrate that fluctuation-induced momentum transfer between the electrolyte and wall electrons produces distinct AC transport signatures, depending on the charge carrier polarity. Altogether, we establish a frequency-dependent transport matrix that couples ionic, electronic, and hydrodynamic flows. These findings establish AC nanofluidic transport as a powerful probe of interfacial phenomena under confinement and suggest new directions for engineering nanofluidic functionalities through electron–electrolyte coupling.},
  author       = {Coquinot, Baptiste and Lizée, Mathieu and Bocquet, Lydéric and Kavokine, Nikita},
  issn         = {1089-7690},
  journal      = {The Journal of Chemical Physics},
  number       = {13},
  publisher    = {AIP Publishing},
  title        = {{Electron–electrolyte coupling in AC transport through nanofluidic channels}},
  doi          = {10.1063/5.0313352},
  volume       = {164},
  year         = {2026},
}

@article{21381,
  abstract     = {The lack of long-range electrostatics is a key limitation of modern machine learning interatomic potentials (MLIPs), hindering reliable applications to interfaces, charge-transfer reactions, polar and ionic materials, and biomolecules. In this Perspective, we distill two design principles behind the Latent Ewald Summation framework, which can capture long-range interactions, charges, and electrical response just by learning from standard energy and force training data: (i) use a Coulomb functional form with environment-dependent charges to capture electrostatic interactions, and (ii) avoid explicit training on ambiguous density functional theory partial charges. When both principles are satisfied, substantial flexibility remains: essentially any short-range MLIP can be augmented; charge equilibration schemes can be added when desired; dipoles and Born effective charges can be inferred or fine-tuned; and charge/spin-state embeddings or tensorial targets can be further incorporated. We also discuss current limitations and open challenges. Together, these minimal, physics-guided design rules suggest that incorporating long-range electrostatics into MLIPs is simpler and perhaps more broadly applicable than is commonly assumed.},
  author       = {Kim, Dongjin and Cheng, Bingqing},
  issn         = {1089-7690},
  journal      = {The Journal of Chemical Physics},
  number       = {6},
  publisher    = {AIP Publishing},
  title        = {{Long-range electrostatics for machine learning interatomic potentials is easier than we thought}},
  doi          = {10.1063/5.0316886},
  volume       = {164},
  year         = {2026},
}

@article{19276,
  abstract     = {Impurity motion in a many-body environment has been a central issue in the field of low-temperature physics for decades. In bosonic quantum fluids, the onset of a drag force experienced by point-like objects is due to collective environment excitations, driven by the exchange of linear momentum between the impurity and the many-body bath. In this work we consider a rotating impurity, with the aim of exploring how angular momentum is exchanged with the surrounding bosonic environment. In order to elucidate these issues, we employ a quasiparticle approach based on the angulon theory, which allows us to effectively deal with the non-trivial algebra of quantized angular momentum in the presence of a many-body environment. We uncover how impurity dressing by environmental excitations can establish an exchange channel, whose effectiveness crucially depends on the initial state of the impurity. Remarkably, we find that there is a critical value of initial angular momentum, above which this channel effectively freezes.},
  author       = {Cappellaro, Alberto and Bighin, Giacomo and Cherepanov, Igor and Lemeshko, Mikhail},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {7},
  publisher    = {AIP Publishing},
  title        = {{Environment-limited transfer of angular momentum in Bose liquids}},
  doi          = {10.1063/5.0253451},
  volume       = {162},
  year         = {2025},
}

@article{19279,
  abstract     = {Recent experimental advances in nanofluidics have allowed to explore ion transport across molecular-scale pores, in particular, for iontronic applications. Two-dimensional nanochannels—in which a single molecular layer of electrolyte is confined between solid walls—constitute a unique platform to investigate fluid and ion transport in extreme confinement, highlighting unconventional transport properties. In this work, we study ionic association in 2D nanochannels, and its consequences on non-linear ionic transport, using both molecular dynamics simulations and analytical theory. We show that under sufficient confinement, ions assemble into pairs or larger clusters in a process analogous to a Kosterlitz–Thouless transition, here modified by the dielectric confinement. We further show that the breaking of pairs results in an electric-field dependent conduction, a mechanism usually known as the second Wien effect. However the 2D nature of the system results in non-universal, temperature-dependent, scaling of the conductivity with electric field, leading to ionic coulomb blockade in some regimes. A 2D generalization of the Onsager theory fully accounts for the non-linear transport. These results suggest ways to exploit electrostatic interactions between ions to build new nanofluidic devices.},
  author       = {Toquer, Damien and Bocquet, Lydéric and Robin, Paul},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {6},
  publisher    = {AIP Publishing},
  title        = {{Ionic association and Wien effect in 2D confined electrolytes}},
  doi          = {10.1063/5.0241949},
  volume       = {162},
  year         = {2025},
}

@article{19880,
  abstract     = {We investigate quantum transport in a two-dimensional electron system coupled to a chiral molecular potential, demonstrating how molecular chirality and orientation affect charge and spin transport properties. We propose a minimal model for realizing true chiral symmetry breaking on a magnetized surface, with a crucial role played by the tilt angle of the molecular dipole with respect to the surface. For non-zero tilting, we show that the Hall response exhibits clear signatures of chirality-induced effects, in both charge- and spin-resolved observables. Concerning the former, tilted enantiomers produce asymmetric Hall conductances and, even more remarkably, the persistence of this feature in the absence of spin–orbit coupling (SOC) signals how the enantiospecific charge response results from electron scattering off the molecular potential. Concerning spin-resolved observables where SOC plays a relevant role, we reveal that chiral symmetry breaking is crucial in enabling spin-flipping processes.},
  author       = {Al Hyder, Ragheed and Lemeshko, Mikhail and Cappellaro, Alberto},
  issn         = {1089-7690},
  journal      = {The Journal of Chemical Physics},
  number       = {23},
  publisher    = {AIP Publishing},
  title        = {{Quantum transport in the presence of a chiral molecular potential}},
  doi          = {10.1063/5.0271155},
  volume       = {162},
  year         = {2025},
}

@article{19595,
  abstract     = {We investigate the locality of magnetic response in polycyclic aromatic molecules using a novel deep-learning approach. Our method employs graph neural networks (GNNs) with a graph-of-rings representation to predict nucleus independent chemical shifts (NICS) in the space around the molecule. We train a series of models, each time reducing the size of the largest molecules used in training. The accuracy of prediction remains high (MAE < 0.5 ppm), even when training the model only on molecules with up to four rings, thus providing strong evidence for the locality of magnetic response. To overcome the known problem of generalization of GNNs, we implement a k-hop expansion strategy and succeed in achieving accurate predictions for molecules with up to 15 rings (almost 4 times the size of the largest training example). Our findings have implications for understanding the magnetic response in complex molecules and demonstrate a promising approach to overcoming GNN scalability limitations. Furthermore, the trained models enable rapid characterization, without the need for more expensive DFT calculations.},
  author       = {Davidson, Yair and Philipp, Aviad and Chakraborty, Sabyasachi and Bronstein, Alexander and Gershoni-Poranne, Renana},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {14},
  publisher    = {AIP Publishing},
  title        = {{How local is “local”? Deep learning reveals locality of the induced magnetic field of polycyclic aromatic hydrocarbons}},
  doi          = {10.1063/5.0257558},
  volume       = {162},
  year         = {2025},
}

@article{17278,
  abstract     = {An azeotrope is a constant boiling point mixture, and its behavior is important for fluid separation processes. Predicting azeotropes from atomistic simulations is difficult due to the complexities and convergence problems of Monte Carlo and free-energy perturbation techniques. Here, we present a methodology for predicting the azeotropes of binary mixtures, which computes the compositional dependence of chemical potentials from molecular dynamics simulations using the S0 method and employs experimental boiling point and vaporization enthalpy data. Using this methodology, we reproduce the azeotropes, or lack thereof, in five case studies, including ethanol/water, ethanol/isooctane, methanol/water, hydrazine/water, and acetone/chloroform mixtures. We find that it is crucial to use the experimental boiling point and vaporization enthalpy for reliable azeotrope predictions, as empirical force fields are not accurate enough for these quantities. Finally, we use regular solution models to rationalize the azeotropes and reveal that they tend to form when the mixture components have similar boiling points and strong interactions.},
  author       = {Wang, Xiaoyu and Cheng, Bingqing},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {3},
  publisher    = {AIP Publishing},
  title        = {{Integrating molecular dynamics simulations and experimental data for azeotrope predictions in binary mixtures}},
  doi          = {10.1063/5.0217232},
  volume       = {161},
  year         = {2024},
}

@article{15024,
  abstract     = {Electrostatic correlations between ions dissolved in water are known to impact their transport properties in numerous ways, from conductivity to ion selectivity. The effects of these correlations on the solvent itself remain, however, much less clear. In particular, the addition of salt has been consistently reported to affect the solution’s viscosity, but most modeling attempts fail to reproduce experimental data even at moderate salt concentrations. Here, we use an approach based on stochastic density functional theory, which accurately captures charge fluctuations and correlations. We derive a simple analytical expression for the viscosity correction in concentrated electrolytes, by directly linking it to the liquid’s structure factor. Our prediction compares quantitatively to experimental data at all temperatures and all salt concentrations up to the saturation limit. This universal link between the microscopic structure and viscosity allows us to shed light on the nanoscale dynamics of water and ions under highly concentrated and correlated conditions.},
  author       = {Robin, Paul},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {6},
  publisher    = {AIP Publishing},
  title        = {{Correlation-induced viscous dissipation in concentrated electrolytes}},
  doi          = {10.1063/5.0188215},
  volume       = {160},
  year         = {2024},
}

@article{18174,
  abstract     = {We investigate the phase ordering (pattern formation) of systems of two-dimensional core–shell particles using Monte Carlo (MC) computer simulations and classical density functional theory (DFT). The particles interact via a pair potential having a hard core and a repulsive square shoulder. Our simulations show that on cooling, the liquid state structure becomes increasingly characterized by long wavelength density modulations and on further cooling forms a variety of other phases, including clustered, striped, and other patterned phases. In DFT, the hard core part of the potential is treated using either fundamental measure theory or a simple local density approximation, whereas the soft shoulder is treated using the random phase approximation. The different DFTs are benchmarked using large-scale grand-canonical-MC and Gibbs-ensemble-MC simulations, demonstrating their predictive capabilities and shortcomings. We find that having the liquid state static structure factor S(k) for wavenumber k is sufficient to identify the Fourier modes governing both the liquid and solid phases. This allows us to identify from easier-to-obtain liquid state data the wavenumbers relevant to the periodic phases and to predict roughly where in the phase diagram these patterned phases arise.},
  author       = {Wassermair, Michael and Kahl, Gerhard and Roth, Roland and Archer, Andrew J.},
  issn         = {1089-7690},
  journal      = {The Journal of chemical physics},
  number       = {12},
  publisher    = {AIP Publishing},
  title        = {{Fingerprints of ordered self-assembled structures in the liquid phase of a hard-core, square-shoulder system}},
  doi          = {10.1063/5.0226954},
  volume       = {161},
  year         = {2024},
}

@article{12705,
  abstract     = {The elasticity of disordered and polydisperse polymer networks is a fundamental problem of soft matter physics that is still open. Here, we self-assemble polymer networks via simulations of a mixture of bivalent and tri- or tetravalent patchy particles, which result in an exponential strand length distribution analogous to that of experimental randomly cross-linked systems. After assembly, the network connectivity and topology are frozen and the resulting system is characterized. We find that the fractal structure of the network depends on the number density at which the assembly has been carried out, but that systems with the same mean valence and same assembly density have the same structural properties. Moreover, we compute the long-time limit of the mean-squared displacement, also known as the (squared) localization length, of the cross-links and of the middle monomers of the strands, showing that the dynamics of long strands is well described by the tube model. Finally, we find a relation connecting these two localization lengths at high density and connect the cross-link localization length to the shear modulus of the system.},
  author       = {Sorichetti, Valerio and Ninarello, Andrea and Ruiz-Franco, José and Hugouvieux, Virginie and Zaccarelli, Emanuela and Micheletti, Cristian and Kob, Walter and Rovigatti, Lorenzo},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {7},
  publisher    = {American Institute of Physics},
  title        = {{Structure and elasticity of model disordered, polydisperse, and defect-free polymer networks}},
  doi          = {10.1063/5.0134271},
  volume       = {158},
  year         = {2023},
}

@article{12831,
  abstract     = {The angulon, a quasiparticle formed by a quantum rotor dressed by the excitations of a many-body bath, can be used to describe an impurity rotating in a fluid or solid environment. Here, we propose a coherent state ansatz in the co-rotating frame, which provides a comprehensive theoretical description of angulons. We reveal the quasiparticle properties, such as energies, quasiparticle weights, and spectral functions, and show that our ansatz yields a persistent decrease in the impurity’s rotational constant due to many-body dressing, which is consistent with experimental observations. From our study, a picture of the angulon emerges as an effective spin interacting with a magnetic field that is self-consistently generated by the molecule’s rotation. Moreover, we discuss rotational spectroscopy, which focuses on the response of rotating molecules to a laser perturbation in the linear response regime. Importantly, we take into account initial-state interactions that have been neglected in prior studies and reveal their impact on the excitation spectrum. To examine the angulon instability regime, we use a single-excitation ansatz and obtain results consistent with experiments, in which a broadening of spectral lines is observed while phonon wings remain highly suppressed due to initial-state interactions.},
  author       = {Zeng, Zhongda and Yakaboylu, Enderalp and Lemeshko, Mikhail and Shi, Tao and Schmidt, Richard},
  issn         = {1089-7690},
  journal      = {The Journal of Chemical Physics},
  number       = {13},
  publisher    = {American Institute of Physics},
  title        = {{Variational theory of angulons and their rotational spectroscopy}},
  doi          = {10.1063/5.0135893},
  volume       = {158},
  year         = {2023},
}

@article{14321,
  abstract     = {We demonstrate the possibility of a coupling between the magnetization direction of a ferromagnet and the tilting angle of adsorbed achiral molecules. To illustrate the mechanism of the coupling, we analyze a minimal Stoner model that includes Rashba spin–orbit coupling due to the electric field on the surface of the ferromagnet. The proposed mechanism allows us to study magnetic anisotropy of the system with an extended Stoner–Wohlfarth model and argue that adsorbed achiral molecules can change magnetocrystalline anisotropy of the substrate. Our research aims to motivate further experimental studies of the current-free chirality induced spin selectivity effect involving both enantiomers.},
  author       = {Al Hyder, Ragheed and Cappellaro, Alberto and Lemeshko, Mikhail and Volosniev, Artem},
  issn         = {1089-7690},
  journal      = {The Journal of Chemical Physics},
  keywords     = {Physical and Theoretical Chemistry, General Physics and Astronomy},
  number       = {10},
  publisher    = {AIP Publishing},
  title        = {{Achiral dipoles on a ferromagnet can affect its magnetization direction}},
  doi          = {10.1063/5.0165806},
  volume       = {159},
  year         = {2023},
}

@article{12912,
  abstract     = {The chemical potential of adsorbed or confined fluids provides insight into their unique thermodynamic properties and determines adsorption isotherms. However, it is often difficult to compute this quantity from atomistic simulations using existing statistical mechanical methods. We introduce a computational framework that utilizes static structure factors, thermodynamic integration, and free energy perturbation for calculating the absolute chemical potential of fluids. For demonstration, we apply the method to compute the adsorption isotherms of carbon dioxide in a metal-organic framework and water in carbon nanotubes.},
  author       = {Schmid, Rochus and Cheng, Bingqing},
  issn         = {1089-7690},
  journal      = {The Journal of Chemical Physics},
  number       = {16},
  publisher    = {AIP Publishing},
  title        = {{Computing chemical potentials of adsorbed or confined fluids}},
  doi          = {10.1063/5.0146711},
  volume       = {158},
  year         = {2023},
}

@article{14603,
  abstract     = {Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances.},
  author       = {Reinhardt, Aleks and Chew, Pin Yu and Cheng, Bingqing},
  issn         = {1089-7690},
  journal      = {Journal of Chemical Physics},
  number       = {18},
  publisher    = {AIP Publishing},
  title        = {{A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals}},
  doi          = {10.1063/5.0173341},
  volume       = {159},
  year         = {2023},
}

@article{11400,
  abstract     = {By varying the concentration of molecules in the cytoplasm or on the membrane, cells can induce the formation of condensates and liquid droplets, similar to phase separation. Their thermodynamics, much studied, depends on the mutual interactions between microscopic constituents. Here, we focus on the kinetics and size control of 2D clusters, forming on membranes. Using molecular dynamics of patchy colloids, we model a system of two species of proteins, giving origin to specific heterotypic bonds. We find that concentrations, together with valence and bond strength, control both the size and the growth time rate of the clusters. In particular, if one species is in large excess, it gradually saturates the binding sites of the other species; the system then becomes kinetically arrested and cluster coarsening slows down or stops, thus yielding effective size selection. This phenomenology is observed both in solid and fluid clusters, which feature additional generic homotypic interactions and are reminiscent of the ones observed on biological membranes.},
  author       = {Palaia, Ivan and Šarić, Anđela},
  issn         = {1089-7690},
  journal      = {The Journal of Chemical Physics},
  keywords     = {Physical and Theoretical Chemistry, General Physics and Astronomy},
  number       = {19},
  publisher    = {AIP Publishing},
  title        = {{Controlling cluster size in 2D phase-separating binary mixtures with specific interactions}},
  doi          = {10.1063/5.0087769},
  volume       = {156},
  year         = {2022},
}

@article{10827,
  abstract     = {Titanium dioxide has been extensively studied in the rutile or anatase phase, while its high-pressure phases are less well-understood, despite that many are thought to have interesting optical, mechanical, and electrochemical properties. First-principles methods, such as density functional theory (DFT), are often used to compute the enthalpies of TiO2 phases at 0 K, but they are expensive and, thus, impractical for long time scale and large system-size simulations at finite temperatures. On the other hand, cheap empirical potentials fail to capture the relative stabilities of various polymorphs. To model the thermodynamic behaviors of ambient and high-pressure phases of TiO2, we design an empirical model as a baseline and then train a machine learning potential based on the difference between the DFT data and the empirical model. This so-called Δ-learning potential contains long-range electrostatic interactions and predicts the 0 K enthalpies of stable TiO2 phases that are in good agreement with DFT. We construct a pressure–temperature phase diagram of TiO2 in the range 0 < P < 70 GPa and 100 < T < 1500 K. We then simulate dynamic phase transition processes by compressing anatase at different temperatures. At 300 K, we predominantly observe an anatase-to-baddeleyite transformation at about 20 GPa via a martensitic two-step mechanism with a highly ordered and collective atomic motion. At 2000 K, anatase can transform into cotunnite around 45–55 GPa in a thermally activated and probabilistic manner, accompanied by diffusive movement of oxygen atoms. The pressures computed for these transitions show good agreement with experiments. Our results shed light on how to synthesize and stabilize high-pressure TiO2 phases, and our method is generally applicable to other functional materials with multiple polymorphs.},
  author       = {Lee, Jacob G. and Pickard, Chris J. and Cheng, Bingqing},
  issn         = {1089-7690},
  journal      = {The Journal of chemical physics},
  number       = {7},
  publisher    = {AIP Publishing},
  title        = {{High-pressure phase behaviors of titanium dioxide revealed by a Δ-learning potential}},
  doi          = {10.1063/5.0079844},
  volume       = {156},
  year         = {2022},
}

@article{12249,
  abstract     = {The chemical potential of a component in a solution is defined as the free energy change as the amount of that component changes. Computing this fundamental thermodynamic property from atomistic simulations is notoriously difficult because of the convergence issues involved in free energy methods and finite size effects. This Communication presents the so-called S0 method, which can be used to obtain chemical potentials from static structure factors computed from equilibrium molecular dynamics simulations under the isothermal–isobaric ensemble. This new method is demonstrated on the systems of binary Lennard-Jones particles, urea–water mixtures, a NaCl aqueous solution, and a high-pressure carbon–hydrogen mixture. },
  author       = {Cheng, Bingqing},
  issn         = {1089-7690},
  journal      = {The Journal of Chemical Physics},
  keywords     = {Physical and Theoretical Chemistry, General Physics and Astronomy},
  number       = {12},
  publisher    = {AIP Publishing},
  title        = {{Computing chemical potentials of solutions from structure factors}},
  doi          = {10.1063/5.0107059},
  volume       = {157},
  year         = {2022},
}

@article{7956,
  abstract     = {When short-range attractions are combined with long-range repulsions in colloidal particle systems, complex microphases can emerge. Here, we study a system of isotropic particles, which can form lamellar structures or a disordered fluid phase when temperature is varied. We show that, at equilibrium, the lamellar structure crystallizes, while out of equilibrium, the system forms a variety of structures at different shear rates and temperatures above melting. The shear-induced ordering is analyzed by means of principal component analysis and artificial neural networks, which are applied to data of reduced dimensionality. Our results reveal the possibility of inducing ordering by shear, potentially providing a feasible route to the fabrication of ordered lamellar structures from isotropic particles.},
  author       = {Pȩkalski, J. and Rzadkowski, Wojciech and Panagiotopoulos, A. Z.},
  issn         = {1089-7690},
  journal      = {The Journal of chemical physics},
  number       = {20},
  publisher    = {AIP Publishing},
  title        = {{Shear-induced ordering in systems with competing interactions: A machine learning study}},
  doi          = {10.1063/5.0005194},
  volume       = {152},
  year         = {2020},
}

