Can we approach the gas–liquid critical point using slab simulations of two coexisting phases?

The Intrinsic Fragility of the Liquid-Vapor Interface: A Stress Network Perspective

The evolution of the liquid-vapor interface of a Lennard-Jones fluid is examined with molecular dynamics simulations using the intrinsic sampling method. Results suggest clear damping of the intrinsic profiles with increasing temperature. Investigating the surface stress distribution, we have identified a linear variation of the space-filling nature (fractal dimension) of the stress clusters at the intrinsic surface with increasing surface tension or, equivalently, with decreasing temperature. A percolation analysis of these stress networks indicates that the stress field is more disjointed at higher temperatures. This leads to more fragile (or poorly connected) interfaces which result in a reduction in surface tension.

Interfacial properties of binary mixtures of Lennard-Jones chains in planar interfaces by molecular dynamics simulation

Binary mixtures of fully flexible linear tangent chains composed of bonded Lennard-Jones interaction sites (monomers) were studied using the molecular dynamics simulation in the NVT ensemble. Their interfacial properties were investigated in planar interfaces by direct simulation of an explicit liquid film in equilibrium with its vapor. A method for the calculation of long-range interactions in inhomogeneous fluids was implemented to take into account the potential truncation effects. Surface tension and the pressure tensor were calculated via the classical Irving-Kirkwood method; vapor pressure, orthobaric densities, density profiles, and Gibbs relative adsorption of the volatile component with respect to the heavy component were also obtained. The properties were studied as a function of the temperature, molar concentration of the heavy component, and the asymmetry of the mixture. According to the results of this work, the temperature loses influence on the surface tension, vapor pressure, and Gibbs relative adsorption curves as the molecular length of the heavy component increases. This suggests that the universal behavior observed in pure fluids of Lennard-Jones chains also holds for binary mixtures. The contribution of the long-range interactions turned out to account for about 60%, 20%, and 10% of the surface tension, vapor pressure, and orthobaric density final values, respectively. This contribution was even larger at high temperatures and for large molecules. Strong enrichment of the volatile component at the interface was observed in the asymmetric mixtures. One of these mixtures even showed a barotropic effect at elevated pressures and a class III phase behavior.

Single-Protein Collapse Determines Phase Equilibria of a Biological Condensate

Recent advances in microscopy of living cells have established membraneless organelles as critical elements of diverse biological processes. The body of experimental work suggests that formation of such organelles is driven by liquid-liquid phase separation, a physical process that has been studied extensively for both simple liquids and mixtures of polymers. Here, we combine molecular dynamics simulations with polymer theory to show that the thermodynamic behavior of one particular biomolecular condensate-fused in sarcoma (FUS)-can be quantitatively accounted for at the level of the chain collapse theory. First, we show that a particle-based molecular dynamics model can reproduce known phase separation properties of a FUS condensate, including its critical concentration and susceptibility to mutations. Next, we obtain a polymer physics representation of a FUS condensate by examining the behavior of a single FUS protein as a function of temperature. We use the chain collapse theory to determine the thermodynamic properties of the condensate and to characterize changes in the single-chain conformation at the onset of phase separation. Altogether, our findings suggest that the phase behavior of FUS condensates can be explained by the properties of individual FUS proteins and that the change in the FUS conformation is the main force driving for the phase separation.

Molecular Dynamics Simulation of Pure n-Alkanes and Their Mixtures at Elevated Temperatures Using Atomistic and Coarse-Grained Force Fields

The properties of higher n-alkanes and their mixtures is a topic of significant interest for the oil and chemical industry. However, the experimental data at high temperatures are scarce. The present study focuses on simulating n-dodecane, n-octacosane, their binary mixture at a n-dodecane mole fraction of 0.3, and a model mixture of the commercially available hydrocarbon wax SX-70 to evaluate the performance of several force fields on the reproduction of properties such as liquid densities, surface tension, and viscosities. Molecular dynamics simulations over a broad temperature range from 323.15 to 573.15 K were employed in examining a broad set of atomistic molecular models assessed for the reproduction of experimental data. The well-established united atom TraPPE (TraPPE-UA) was compared against the all atom optimized potentials for liquid simulations (OPLS) reparametrization for long n-alkanes, L-OPLS, as well as Lipid14 and MARTINI force fields. All models qualitatively reproduce the temperature dependence of the aforementioned properties, but TraPPE-UA was found to reproduce liquid densities most accurately and consistently over the entire temperature range. TraPPE-UA and MARTINI were very successful in reproducing surface tensions, and L-OPLS was found to be the most accurate in reproducing the measured viscosities as compared to the other models. Our simulations show that these widely used force fields originating from the world of biomolecular simulations are suitable candidates in the study of n-alkane properties, both in the pure and mixture states.