Atomistic analysis of nematic phase transition in 4-cyano-4$^{\prime}$-$n$-alkyl biphenyl liquid crystals: Sampling for the first-order phase transition and the free-energy decomposition

Shunsuke Ogita,Yoshiki Ishii,Go Watanabe,Hitoshi Washizu,Kang Kim,Nobuyuki Matubayasi
2024-10-07
Abstract:Molecular dynamics (MD) simulations were conducted using the generalized replica exchange method (gREM) on the 4-cyano-4$^{\prime}$-$n$-alkylbiphenyl ($n$CB) 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 $n$CB molecule for both nematic and isotropic phases, revealing a change in the temperature dependence across the NI phase transition. Further decomposition into intermolecular interaction energetic and entropic terms shows quantitatively the balance between these contributions at the NI phase transition temperature.
Soft Condensed Matter
What problem does this paper attempt to address?
The problem this paper attempts to address is how to study the free energy changes and their decomposition during the nematic-isotropic (NI) phase transition of the 4-cyano-4'-n-alkylbiphenyl (nCB) system using molecular dynamics (MD) simulations combined with the generalized replica exchange method (gREM). Specifically, the authors focus on the following aspects: 1. **Enhanced Sampling**: Since the NI phase transition is a first-order phase transition with significant energy gaps, effective sampling near the transition temperature is very challenging. The authors use the gREM method to improve sampling efficiency in this region. 2. **Free Energy Calculation**: Based on the energy representation theory (ER theory), the authors calculate the insertion free energy of nCB molecules in the nematic and isotropic phases to characterize the thermodynamic stability of the NI phase transition. 3. **Free Energy Decomposition**: The insertion free energy is further decomposed into intermolecular interaction energy and entropy contributions to quantitatively analyze the equilibrium relationship of these factors at the NI phase transition temperature. Through these studies, the authors hope to achieve a more precise molecular-level understanding, surpassing traditional mean-field models, and reveal the microscopic mechanisms of the NI phase transition.