Entropic Description of Gas Hydrate Ice-Liquid Equilibrium via Enhanced Sampling of Coexisting Phases

Atomistic analysis of nematic phase transition in 4-cyano-4'-n-alkyl biphenyl liquid crystals: Sampling for the first-order phase transition and the free-energy decomposition

Molecular dynamics simulations were conducted using the generalized replica exchange method (gREM) on the 4-cyano-4'-n-alkyl biphenyl (nCB) system with n = 5, 6, 7, and 8, which exhibits a nematic-isotropic (NI) phase transition. Sampling near the phase transition temperature in systems undergoing first-order phase transitions, such as the NI phase transition, is demanding due to the substantial energy gap between the two phases. To address this, gREM, specifically designed for first-order phase transitions, was utilized to enhance sampling near the NI phase transition temperature. Free-energy calculations based on the energy representation (ER) theory were employed to characterize the NI phase transition. ER evaluates the insertion free energy of the nCB molecule for both nematic and isotropic phases, revealing a change in the temperature dependence across the NI phase transition. Further decomposition into energetic and entropic terms quantitatively shows the balance between these contributions at the NI phase transition temperature.

Simulating the nematic-isotropic phase transition of liquid crystal model via generalized replica-exchange method

The nematic-isotropic (NI) phase transition of 4-cyano-4'-pentylbiphenyl was simulated using the generalized replica-exchange method (gREM) based on molecular dynamics simulations. The effective temperature is introduced in the gREM, allowing for the enhanced sampling of configurations in the unstable region, which is intrinsic to the first-order phase transition. The sampling performance was analyzed with different system sizes and compared with that of the temperature replica-exchange method (tREM). It was observed that gREM is capable of sampling configurations at sufficient replica-exchange acceptance ratios even around the NI transition temperature. A bimodal distribution of the order parameter at the transition region was found, which is in agreement with the mean-field theory. In contrast, tREM is ineffective around the transition temperature owing to the potential energy gap between the nematic and isotropic phases.

Simulation of fluid/gel phase equilibrium in lipid vesicles

Simulation of single component dipalmitoylphosphatidylcholine (DPPC) coarse-grained DRY-MARTINI lipid vesicles of diameter 10 nm (1350 lipids), 20 nm (5100 lipids) and 40 nm (17 600 lipids) is performed using statistical temperature molecular dynamics (STMD), to study finite size effects upon the order-disorder gel/fluid transition. STMD obtains enhanced sampling using a generalized ensemble, obtaining a flat energy distribution between upper and lower cutoffs, with little computational cost over canonical molecular dynamics. A single STMD trajectory of moderate length is sufficient to sample 20+ transition events, without trapping in the gel phase, and obtain well averaged properties. Phase transitions are analyzed via the energy-dependence of the statistical temperature, T_S(U). The transition temperature decreases with decreasing diameter, in agreement with experiment, and the transition changes from first order to borderline first-second order. The size- and layer-dependence of the structure of both stable phases, and of the pathway of the phase transition, are determined. It is argued that the finite size effects are primarily caused by the disruption of the gel packing by curvature. Inhomogeneous states with faceted gel patches connected by unusual fluid seams are observed at high curvature, with visually different structure in the inner and outer layers due to the different curvatures. Thus a simple physical picture describes phase transitions in nanoscale finite systems far from the thermodynamic limit.

Revealing Transition States during the Hydration of Clay Minerals

A molecular-scale understanding of the transition between hydration states in clay minerals remains a challenging problem because of the very fast stepwise swelling process observed from X-ray diffraction (XRD) experiments. XRD profile modeling assumes the coexistence of multiple hydration states in a clay sample to fit the experimental XRD pattern obtained under humid conditions. While XRD profile modeling provides a macroscopic understanding of the heterogeneous hydration structure of clay minerals, a microscopic model of the transition between hydration states is still missing. Here, for the first time, we use molecular dynamics simulation to investigate the transition states between a dry interlayer, one-layer hydrate, and two-layer hydrate. We find that the hydrogen bonds that form across the interlayer at the clay particle edge make an important contribution to the energy barrier to interlayer hydration, especially for initial hydration.

Enhanced Sampling of Phase Transitions in Coarse-Grained Lipid Bilayers

Freezing and melting of dipalmitoylphosphatidylcholine (DPPC) bilayers are simulated in both the explicit (Wet) and implicit solvent (Dry) coarse-grained MARTINI force fields with enhanced sampling, via the isobaric, molecular dynamics version of the generalized replica exchange method (gREM). Phase transitions are described with the entropic viewpoint, based upon the statistical temperature as a function of enthalpy, T_S(H) = 1/(dS(H)/dH), where S is the configurational entropy. Bilayer thickness, area per lipid, and the second-rank order parameter (P2) are calculated vs temperature in the transition range. In a 32-lipid Wet MARTINI system, transitions in the lipid and water subsystems are strongly coupled, giving rise to considerable structure in T_S(H) and the need to specify the state of the water when reporting a lipid transition temperature. For gel lipid + liquid water → fluid lipid + liquid water, we find 292.4 K. The small system is influenced by finite-size effects, but it is argued that the entropic approach is well suited to revealing them, which will be particularly relevant for studies of finite nanosystems where there is no thermodynamic limit. In a 390-lipid Dry MARTINI system, two-dimensional analogues of the topographies of coexisting states ("subphases") seen in pure fluids are found. They are not seen in the 32-lipid Wet or Dry system, but the Dry lipids show a new type of state with gel in one leaflet and tilted gel in the other. Dry bilayer transition temperatures are 333.3 K (390 lipids) and 338 K (32 lipids), indicating that the 32-lipid system is not too small for a qualitative study of the transition. Physical arguments are given for Dry lipid system size dependence and for the difference between Wet and Dry systems.