Chiral symmetry breaking in a microscopic model with asymmetric autocatalysis and inhibition

Amino acid chiral amplification using Monte Carlo dynamic

This study investigates the stability of chiral-molecule solution phases, with a specific focus on amino acids. The model framework is based on a two-dimensional square lattice model, where individual sites may be occupied by oriented chiral molecules or structureless solvent particles. Utilizing the Glauber dynamics and statistical mechanical formalism, as previously introduced and examined by Lombardo et al., we explore the influence of temperature, amino acid concentration, enantiomeric excess, and homochiral interaction strength on nucleation mechanisms, equilibrium phase behavior, and crystal composition. Our findings offer thermodynamic insights into the chiral amplification process of amino acids, contributing to a deeper understanding of the underlying processes.

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.

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.

Effect of configuration-dependent multi-body forces on interconversion kinetics of a chiral tetramer model

We describe a reformulation of the four-site molecular model for chiral phenomena introduced by Latinwo et al. ["Molecular model for chirality phenomena," J. Chem. Phys. 145, 154503 (2016)]. The reformulation includes an additional eight-body force that arises from an explicit configuration-dependent term in the potential energy function, resulting in a coarse-grained energy-conserving force field for molecular dynamics simulations of chirality phenomena. In this model, the coarse-grained interaction energy between two tetramers depends on their respective chiralities and is controlled by a parameter λ, where λ < 0 favors local configurations involving tetramers of opposite chirality and λ > 0 gives energetic preference to configurations involving tetramers of the same chirality. We compute the autocorrelation function for a quantitative chirality metric and demonstrate that the multi-body force modifies the interconversion kinetics such that λ ≠ 0 increases the effective barrier for enantiomer inversion. Our simulations reveal that for λ > 0 and temperatures below a sharply defined threshold value, this effect is dramatic, giving rise to spontaneous chiral symmetry breaking and locking molecules into their chiral identity.

2D Ising Model for Adsorption-induced Enantiopurification of Racemates

Mechanisms for the spontaneous transformation of achiral chemical systems into states of enantiomeric purity have important ramifications in modern pharmacology and potential relevance to the origins of homochirality in life on Earth. Such mechanisms for enantiopurification are needed for production of chiral pharmaceuticals and other bioactive compounds. Previously proposed chemical mechanisms leading from achiral systems to near homochirality are initiated by a symmetry-breaking step resulting in a minor excess of one enantiomer via statistical fluctuations in enantiomer concentrations. Subsequent irreversible processes then amplify the majority enantiomer concentration while simultaneously suppressing minority enantiomer production. Herein, equilibrium adsorption of amino acid enantiomer mixtures onto chiral and achiral surfaces reveals amplification of surface enantiomeric excess relative to the gas phase; i. e. enantiopurification of chiral adsorbates by adsorption. This adsorption-induced amplification of enantiomeric excess is shown to be well-describe by the 2D Ising model. More importantly, the 2D-Ising model predicts formation of homochiral monolayers from adsorption of racemic mixtures or prochiral molecules on achiral surfaces; i. e. enantiopurification with no apparent chiral driving force.

Deterministic and Stochastic Chiral Symmetry Breaking Exhibited by Racemic Aminomethylenehelicene Oligomers

Racemic mixtures of aminomethylenehelicene (P)- and (M)-pentamers exhibited deterministic and stochastic chiral symmetry breaking during hetero-double-helix formation and self-assembly in solution. Heating a 50:50 mixture of (P)- and (M)-pentamers at 90 °C, and then cooling the mixture to 70 °C resulted in hetero-double-helix formation; a Cotton effect with negative Δε at λ=315 nm appeared. Chiral self-assembly occurred when the mixture was cooled to 25 °C. A strong tendency of deterministic chiral symmetry breaking appeared at the molecular and self-assembled levels, which was indicated by the negative Δε at λ=315 nm that appeared in most cases in repeated experiments. Mixtures containing 60:40 and 40:60 (P)-/(M)-pentamers also self-assembled with the same chirality. When a homo-double-helix (P)-/(M)-pentamer and a random coil (M)-/(P)-pentamer were mixed, the chiral self-assembly formed stochastically, and heating and cooling resulted in deterministic chiral symmetry breaking.

Modeling spontaneous chiral symmetry breaking and deracemization phenomena: discrete versus continuum approaches

We derive the class of population balance equations (PBE), recently applied to model the Viedma deracemization experiment, from an underlying microreversible kinetic reaction scheme. The continuum limit establishing the relationship between the micro- and macroscopic processes and the associated particle fluxes erases the microreversible nature of the molecular interactions in the population growth rate functions and limits the scope of such PBE models to strict kinetic control. The irreversible binary agglomeration processes modeled in those PBEs contribute an additional source of kinetic control. These limitations are crucial regarding the question of the origin of biological homochirality, where the interest in any model lies precisely in its ability for absolute asymmetric synthesis and the amplification of the tiny inherent statistical chiral fluctuations about the ideal racemic composition up to observable enantiometric excess levels.