Molecular dynamics studies of smectic B liquid crystals of soft parallel spherocylinders with sixfold bond orientational order

Entropy and heat capacity calculations of simulated crystal-hexatic smectic-B-smectic-A liquid-crystal phase transitions

Entropy and heat capacity are calculated in phase sequence of crystal, hexatic smectic- B , and smectic- A liquid crystals through constant pressure and temperature molecular-dynamics simulations of parallel soft spherocylinders. The transition from crystal to hexatic smectic- B phase is continuous while the transition to smectic- A phase is first order. The dependence of the phase sequence against the molecular shape anisotropy is investigated and there exists a triple point at a rather small anisotropy. Hopping diffusion of molecules is observed in the hexatic smectic- B phase.

Stability of Nematic and Smectic Phases in Rod-Like Mesogens with Orientation−Dependent Attractive Interactions

The stability of isotropic (I), nematic (N), smectic A (Sm A), and hexatic (Hex) liquid crystalline phases is studied for a fluid of molecules with a rod-like shape and dispersive interactions dependent on orientation. The fluid is modeled with the spherocylindrical Gay-Berne-Kihara interaction potential proposed in a recent work, with parameters favoring parallel pair orientations. The liquid crystal phase diagram is characterized for different molecular aspect ratios by means of Monte Carlo simulations in the isobaric-isothermal ensemble. Three types of triple points are observed, namely, I-Sm A-Hex, I-N-Sm A, and N-Sm A-Hex, leading to island-shape domains for the smectic A phase. The resulting phase diagrams are compared with those derived previously for prolate fluids of ellipsoidal and spherocylindrical symmetry. It is concluded that the stability of the layered phases with respect to the nematic phase is enhanced in the spherocylindrical fluids due to geometrical constraints. Furthermore, the anisotropy of the dispersive interactions induces a stronger dependence of the overall phase diagram on temperature and aids in the energetic stabilization of the hexatic crystalline phase with respect to the fluid smectic A phase.

Sixfold bond orientational properties of a model liquid crystal in the dimensional crossover of B phases: A computer simulation study

A wide range of NPT simulations of a bead necklace liquid crystal model in the crystal B, smectic B, smectic A, and nematic phases have been performed. Systems with up to 21 600 molecules have been studied to observe the behavior of slowly decaying spatial correlation functions. The pair correlation function and its in-plane restriction are consistent with a crystalline phase made of independent two-dimensional crystalline layers. Smectic B phase is studied by the bond orientational pair correlation functions g(6) and its extension g(6ext). The first reaches a constant value, which seems to rule out a classical hexatic phase. The latter shows a power-law decay within the layers: its typical decay exponent (eta(6ext)) is evaluated. Relationships between multiple harmonics of the C(6n) order parameter have been evaluated through the whole range of existence of B phases (crystalline and smectic): the extension to the crystalline phase holds and provides an excellent fit of the simulation data.

Solvent-free simulations of fluid membrane bilayers

A molecular level model for lipid bilayers is presented. Lipids are represented by rigid, asymmetric, soft spherocylinders in implicit solvent. A simple three parameter potential between pairs of lipids gives rise to a rich assortment of phases including (but not limited to) micelles, fluid bilayers, and gel-like bilayers. Monte Carlo simulations have been carried out to verify self-assembly, characterize the phases corresponding to different potential parametrizations, and to quantify the physical properties associated with those parameter sets corresponding to fluid bilayer behavior. The studied fluid bilayers have compressibility moduli in agreement with experimental systems, but display bending moduli at least three times larger than typical biological membranes without cholesterol.

The role of molecular shape in bilayer elasticity and phase behavior

A previously developed molecular level model for lipid bilayers [G. Brannigan and F. L. H. Brown, J. Chem. Phys. 120, 1059 (2004)] is extended to allow for variations in lipid length and simulations under constant surface tension conditions. The dependence of membrane elasticity on bilayer thickness is obtained by adjusting lipid length at constant temperature and surface tension. Additionally, bilayer fluidity at various lipid lengths is quantified by analysis of a length versus temperature phase diagram at vanishing tension. Regions of solid, gel-like (hexatic) and fluid bilayer behavior are established by identification of phase boundaries. The main melting transition is found to be density driven; the melting temperature scales inversely with lipid length since thermal expansion increases with lipid aspect ratio.