Search

Monte Carlo simulations in the grand canonical ensemble were used to obtain critical parameters and conditions leading to microphase separation for block copolymers with solvophilic and solvophobic segments. Solvent selectivity was systematically varied to distinguish between systems that undergo macrophase separation to ones that microphase separate in the dilute phase, prior to macrophase separating. Finite-size scaling was used to obtain the critical parameters. Interestingly, corrections to scaling increase significantly for systems that form finite aggregates. The threshold value of solvent selectivity for aggregation was determined for symmetric diblock chains of varying length. The results indicate that long diblock copolymers form micelles in the dilute phase prior to macrophase separation, even in marginally selective solvents. The dependence of critical temperature on solvent selectivity was obtained for triblock, multiblock, and alternating chains. For highly selective solvents, strong structuring of both dilute and dense phases makes it harder to reach equilibrium.

Histogram-reweighting grand canonical Monte Carlo simulations are used to obtain the critical properties of lattice chains composed of solvophilic and solvophobic monomers. The model is a modification of one proposed by Larson \emph{et al.} [J. Chem. Phys. 83, 2411 (1985)], lowering the ``contrast'' between beads of different type to prevent aggregation into finite-size micelles that would mask true phase separation between bulk high- and low-density phases. Oligomeric chains of length between 5 and 24 beads are studied. Mixed-field finite-size scaling methods are used to obtain the critical properties with typical relative accuracies of better than $10^{-4}$ for the critical temperature and $10^{-3}$ for the critical volume fraction. Diblock chains are found to have lower critical temperatures and volume fractions relative to the corresponding homopolymers. Addition of solvophilic blocks of increasing length to a fixed-length solvophobic segment results in a decrease of both critical temperature and critical volume fraction, with an eventual slow asymptotic approach to the long-chain limiting behavior. Moving a single solvophobic or solvophilic bead along a chain leads to a minimum or maximum in the critical temperature, with no change in critical volume fraction. Chains of identical length and composition have a significant spread in their critical properties, depending on their precise sequence. The present study has implications on understanding biomolecular phase separation and for developing design rules for synthetic polymers with specific phase separation properties. It also provides data potentially useful for the further development of theoretical models for polymer and surfactant phase behavior.

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters., Comment: has accepted by Physical Review Letters

The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water at deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. We apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 1 nm from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30,000 atoms. Future work will be directed towards the calculation of nucleation free energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces on ice nucleation., Comment: 12 pages, 6 figures

This article focuses on phase and aggregation behavior for linear chains composed of blocks of hydrophilic and hydrophobic segments. Phase and conformational transitions of patterned chains are relevant for understanding liquid-liquid separation of biomolecular condensates, which play a prominent role in cellular biophysics, but also for surfactant and polymer applications. Previous studies of simple models for multiblock chains have shown that, depending on the sequence pattern and chain length, such systems can fall into one of two categories: displaying either phase separation or aggregation into finite-size clusters. The key new result of the present study is that both formation of finite-size aggregates and phase separation can be observed for certain chain architectures at appropriate conditions of temperature and concentration. For such systems, a bulk dense liquid condenses from a dilute phase that already contains multi-chain finite-size aggregates. The computational approach involves several distinct steps using histogram-reweighting grand canonical Monte Carlo simulations, which are described here in some level of detail., Comment: 9 pages, 11 figures

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here we utilize the Deep Potential methodology -- a machine learning approach -- to study this ubiquitous phase transition, starting from the phase diagram in the liquid-vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K, and evaluate the Deep Potential model performance against experimental results and the semi-empirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model, due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 $\pm$ 0.008) nm at 296.4 K. Lastly, we identify that water molecules display a preferential orientation in the liquid-vapor interface, in which H atoms tend to point towards the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behaviour is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid-vapor coexistence and water cavitation.

A longstanding question in water research is the possibility that supercooled liquid water can undergo a liquid-liquid phase transition (LLT) into high- and low-density liquids. We used several complementary molecular simulation techniques to evaluate the possibility of an LLT in an ab initio neural network model of water trained on density functional theory calculations with the SCAN exchange correlation functional. We conclusively show the existence of a first-order LLT and an associated critical point in the SCAN description of water, representing the first definitive computational evidence for an LLT in water from first principles., Comment: 14 pages, 3 figures. Supplemental material contains 11 pages, 8 figures

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine learning model trained on density-functional theory (DFT) energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation) which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates., Comment: 20 pages, 5 figures

Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water's structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited., Comment: 22 pages, 4 figures

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make possible molecular simulations with the accuracy of quantum mechanical density functional theory, at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei+core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the deep potential long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples, the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization, and compared the MB-pol results to experimental data. We also compared rigid vs. fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system size or the rigid vs. flexible treatment of the intramolecular degrees of freedom. These results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye towards studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ice Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic., Comment: 20 pages, 9 figures

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., Comment: The following article has been accepted by the Journal of Chemical Physics AIP. After it is published, it will be found at https://aip.scitation.org/journal/jcp

Phase separation of intrinsically disordered proteins is important for the formation of membraneless organelles, or biomolecular condensates, which play key roles in the regulation of biochemical processes within cells. In this work, we investigated the phase separation of different sequences of a coarse-grained model for intrinsically disordered proteins and discovered a surprisingly rich phase behavior. We studied both the fraction of total hydrophobic parts and the distribution of hydrophobic parts. Not surprisingly, sequences with larger hydrophobic fractions showed conventional liquid-liquid phase separation. The location of the critical point was systematically influenced by the terminal beads of the sequence, due to changes in interfacial composition and tension. For sequences with lower hydrophobicity, we observed not only conventional liquid-liquid phase separation, but also reentrant phase behavior, in which the liquid phase density decreases at lower temperatures. For some sequences, we observed formation of open phases consisting of aggregates, rather than a normal liquid. These aggregates had overall lower densities than the conventional liquid phases, and exhibited complex geometries with large interconnected string-like or membrane-like clusters. Our findings suggest that minor alterations in the ordering of residues may lead to large changes in the phase behavior of the protein, a fact of significant potential relevance for biology.

We present an algorithm for neighbor search in molecular simulations on graphics processing units (GPUs) based on bounding volume hierarchies (BVHs). The BVH is compressed into a low-precision, quantized representation to increase the BVH traversal speed compared to a previous implementation. We find that neighbor search using the quantized BVH is roughly two to four times faster than current state-of-the-art methods using uniform grids (cell lists) for a suite of benchmarks for common molecular simulation models. Based on the benchmark results, we recommend using the BVH instead of a single cell list for neighbor list generation in molecular simulations on GPUs., Comment: 28 pages, 8 figures

In this work, we use large-scale molecular dynamics simulations coupled to free energy calculations to identify for the first time a limit of stability (spinodal) and a change in the nucleation mechanism in aqueous NaCl solutions. This is a system of considerable atmospheric, geological and technical significance. We find that the supersaturated metastable NaCl solution reaches its limit of stability at sufficiently high salt concentrations, as indicated by the composition dependence of the salt chemical potential, indicating the transition to a phase separation by spinodal decomposition. However, the metastability limit of the NaCl solution does not correspond to spinodal decomposition with respect to crystallization. We find that beyond this spinodal, a liquid/amorphous separation occurs in the aqueous solution, whereby the ions first form disordered clusters. We term these clusters as "amorphous salt". We also identify a transition from one- to two-step crystallization mechanism driven by a spinodal. In particular, crystallization from aqueous NaCl solution beyond the spinodal is a two-step process, in which the ions first phase-separate into disordered amorphous salt clusters, followed by the crystallization of ions in the amorphous salt phase. In contrast, in the aqueous NaCl solution at concentrations lower than the spinodal, crystallization occurs via a one-step process, as the ions aggregate directly into crystalline nuclei. The change of mechanism with increasing supersaturation underscores the importance of an accurate determination of the driving force for phase separation. The study has broader implications on the mechanism for nucleation of crystals from solutions at high supersaturations.

We simulated two particle-based fluid models, namely multiparticle collision dynamics and dissipative particle dynamics, under shear using reverse nonequilibrium simulations (RNES). In cubic periodic simulation boxes, the expected shear flow profile for a Newtonian fluid developed, consistent with the fluid viscosities. However, unexpected secondary flows along the shear gradient formed when the simulation box was elongated in the flow direction. The standard shear flow profile was obtained when the simulation box was longer in the shear-gradient dimension than the flow dimension, while the secondary flows were always present when the flow dimension was at least 25% larger than the shear-gradient dimension. The secondary flows satisfy the boundary conditions imposed by the RNES and have a lower rate of viscous dissipation in the fluid than the corresponding unidirectional flows. This work highlights a previously unappreciated limitation of RNES for generating shear flow in simulation boxes that are elongated in the flow dimension, an important consideration when applying RNES to complex fluids like polymer solutions.

Nonequilibrium molecular dynamics simulations are used to investigate the influence of hydrodynamic interactions on vertical segregation (stratification) in drying mixtures of long and short polymer chains. In agreement with previous computer simulations and theoretical modeling, the short polymers stratify on top of the long polymers at the top of the drying film when hydrodynamic interactions between polymers are neglected. However, no stratification occurs at the same drying conditions when hydrodynamic interactions are incorporated through an explicit solvent model. Our analysis demonstrates that models lacking hydrodynamic interactions do not faithfully represent stratification in drying mixtures, in agreement with recent analysis of an idealized model for diffusiophoresis, and must be incorporated into such models in future.

We present an efficient open-source implementation of the multiparticle collision dynamics (MPCD) algorithm that scales to run on hundreds of graphics processing units (GPUs). We especially focus on optimizations for modern GPU architectures and communication patterns between multiple GPUs. We show that a mixed-precision computing model can improve performance compared to a fully double-precision model while still providing good numerical accuracy. We report weak and strong scaling benchmarks of a reference MPCD solvent and a benchmark of a polymer solution with research-relevant interactions and system size. Our MPCD software enables simulations of mesoscale hydrodynamics at length and time scales that would be otherwise challenging or impossible to access.

We discuss the detailed balance condition for hybrid Monte Carlo method

We used molecular dynamics simulations and the path sampling technique known as forward flux sampling to study homogeneous nucleation of NaCl crystals from supersaturated aqueous solutions at 298 K and 1 bar. Nucleation rates were obtained for a range of salt concentrations for the Joung-Cheatham NaCl force field combined with the SPC/E water model. The calculated nucleation rates are significantly lower than available experimental measurements. The estimates for the nucleation rates in this work do not rely on classical nucleation theory, but the pathways observed in the simulations suggest that the nucleation process is better described by classical nucleation theory than an alternative interpretation based on Ostwald's step rule, in contrast to some prior simulations of related models. In addition to the size of NaCl nucleus, we find that the crystallinity of a nascent cluster plays an important role in the nucleation process. Nuclei with high crystallinity were found to have higher growth probability and longer lifetimes, possibly because they are less exposed to hydration water.

Colloidal particles, amphiphiles, and functionalized nanoparticles are examples of systems that frequently exhibit short-range attractions coupled with long-range repulsions. In this work, we observe striking differences in the dynamics of self-assembled clusters that form in a simple isotropic model of such systems when the strength of attraction is varied. We find that while attraction-dominated particles self-assemble into compact clusters with properties similar to micelles formed by amphiphilic molecules, repulsion-dominated particles self-assemble into open clusters which have much shorter life-times. There is also a significantly different dependence of the solution osmotic pressure versus composition: formation of compact clusters causes a decrease in the pressure vs. density slope, while formation of open clusters does not effect the pressure. This thermodynamic quantity turns out to be much more sensitive in picking out different clustering characteristics than the overall aggregation curves or cluster shapes. Our results have significant implications in developing design principles for stable cluster self-assembly and detection in both laboratory settings and in computer simulations., Comment: The cluster life-times for the attraction-dominated (epsilon = 1.6) case, shown as the red curve in Fig. 4, was previously at the incorrect temperature. Previously the temperature was kT = 0.67, in the new version the temperature is kT = 0.75, as the figure label states. The text is altered accordingly