A study of the phase behavior of a simple model of chiral molecules and enantiomeric mixtures

Pre-Smectic Ordering and the Unwinding Helix in Monte Carlo Simulations of Cholesteric Liquid-Crystals

In this paper, molecular chirality is studied for liquid-crystal fluids represented by hard rods with the addition of an attractive chiral dispersion term. Chiral forces between molecular pairs are assumed to be long-ranged and are described in terms of the pseudotensor of Goossens [W. J. A. Goossens, Mol. Cryst. Liq. Cryst. 1971, 12, 237-244]. Following Varga and Jackson [S. Varga and G. Jackson, Chem. Phys. Lett. 2003, 377, 6-12], this is combined with a hard-spherocylinder core. We investigate the relationship between molecular chirality and the helical pitch of the system, which occurs in the absence of full three-dimensional periodic boundary conditions. The dependence of the wavenumber of this pitch on the thermodynamic variables, temperature, and density is measured. We also explore the use of a novel surface boundary interaction model. As a result of this approach, we are able to lower the temperature of the system without the occurrence of nematic droplets, which would interfere with the formation of a uniaxial pitch. Regarding the theoretical predictions of Wensink and Jackson [H. H. Wensink and G. Jackson, J. Chem. Phys. 2009, 130, 234911], on the one hand, we have qualitative agreement with the observed non-monotonic density dependence of the wavenumber. Initially increasing with density, the wavenumber reaches a maximum, before falling as the density moves toward the point of phase transition from cholesteric to smectic. However, further analysis for shorter rods, in the presence of novel boundary conditions, reveals some disagreement with the theory, at least in this case; the unwinding of the cholesteric helix in the cholesteric phase occurs simultaneously with subtle increases in smectic ordering. These pre-smectic fluctuations have not been accounted for so far in theories on cholesterics but turn out to play a key role in controlling the pitch of cholesteric phases of rod-shaped mesogens with a small to moderate aspect ratio.

Crystal Prediction via Genetic Algorithms in a Model Chiral System

Chiral crystals and their constituent molecules play a prominent role in theories about the origin of biological homochirality and in drug discovery, design, and stability. Although the prediction and identification of stable chiral crystal structures is crucial for numerous technologies, including separation processes and polymorph selection and control, predictive ability is often complicated by a combination of many-body interactions and molecular complexity and handedness. In this work, we address these challenges by applying genetic algorithms to predict the ground-state crystal lattices formed by a chiral tetramer molecular model, which we have previously shown to exhibit complex fluid-phase behavior. Using this approach, we explore the relative stability and structures of the model's conglomerate and racemic crystals, and present a structural phase diagram for the stable Bravais crystal types in the zero-temperature limit.

Interconversion-controlled liquid-liquid phase separation in a molecular chiral model

Liquid-liquid phase separation of fluids exhibiting interconversion between alternative states has been proposed as an underlying mechanism for fluid polyamorphism and may be of relevance to the protein function and intracellular organization. However, molecular-level insight into the interplay between competing forces that can drive or restrict phase separation in interconverting fluids remains elusive. Here, we utilize an off-lattice model of enantiomers with tunable chiral interconversion and interaction properties to elucidate the physics underlying the stabilization and tunability of phase separation in fluids with interconverting states. We show that introducing an imbalance in the intermolecular forces between two enantiomers results in nonequilibrium, arrested phase separation into microdomains. We also find that in the equilibrium case, when all interaction forces are conservative, the growth of the phase domain is restricted only by the system size. In this case, we observe phase amplification, in which one of the two alternative phases grows at the expense of the other. These findings provide novel insights on how the interplay between dynamics and thermodynamics defines the equilibrium and steady-state morphologies of phase transitions in fluids with interconverting molecular or supramolecular states.

Thermodynamic mechanism for solution phase chiral amplification via a lattice model

We formulate a two-dimensional lattice model to study the equilibrium phase behavior of a ternary mixture composed of two enantiomeric forms of a chiral molecule and a nonchiral liquid solvent. Numerical solution of the model invoking the mean-field approximation of statistical mechanics allows the calculation of a ternary phase diagram. A prominent feature of the phase diagram is the appearance of a mirror-image pair of triple points involving coexistence of a liquid phase enriched in one of the enantiomers with two solid phases: a racemic and an enantiopure crystal. Thus, over broad ranges of initial composition, including liquid mixtures containing almost equal amounts of the two enantiomers, the equilibrium state of the system produces liquid-phase chiral amplification. The calculations predict that chiral amplification is favored at low temperatures, and by strengthening those molecular interactions that stabilize the racemic crystal. The phase behavior that we obtain is qualitatively identical to that reported in a recent experimental study of solutions of amino acids, aimed at probing the liquid-phase control of asymmetric catalysis. Those results, and the present calculations that provide molecular-level insight into their underlying causes, suggest a possible thermodynamic scenario for the liquid-phase emergence of chiral imbalance in a prebiotic and presumably nearly racemic world.

Calculation of Free Energies and Chemical Potentials for Gas Hydrates Using Monte Carlo Simulations

We describe a method for calculating free energies and chemical potentials for molecular models of gas hydrate systems using Monte Carlo simulations. The method has two components: (i) thermodynamic integration to obtain the water and guest molecule chemical potentials as functions of the hydrate occupancy; (ii) calculation of the free energy of the zero-occupancy hydrate system using thermodynamic integration from an Einstein crystal reference state. The approach is applicable to any classical molecular model of a hydrate. We illustrate the methodology with an application to the structure-I methane hydrate using two molecular models. Results from the method are also used to assess approximations in the van der Waals-Platteeuw theory and some of its extensions. It is shown that the success of the van der Waals-Platteeuw theory is in part due to a cancellation of the error arising from the assumption of a fixed configuration of water molecules in the hydrate framework with that arising from the neglect of methane-methane interactions.

Chiral separation on a model adsorbent with periodic surface heterogeneity

Optimization of enantioselectivity in heterogeneous catalysis and chiral chromatography is a challenging task for the production of enantiopure chemicals. Enantioselective adsorbents usually consist of a surface with chiral receptors being either chiral molecules linked to the surface or chiral pockets formed by molecular templating of the surface. In both cases, the enantioselectivity is controlled mainly by the strength of the receptor-enantiomer interaction, such that one-to-one correspondence is usually preserved. The authors use Monte Carlo calculations to show that this steric requirement is not a necessary condition for the effective separation of chiral molecules. In particular, they propose a way in which a chiral surface can be constructed by a suitable spatial distribution of active sites for which the classical concept of a chiral receptor is no longer useful. Their calculations indicate that the effectiveness of the separation is affected mainly by the difference in shape of the adsorption energy distribution functions corresponding to the enantiomers.