Atomistic Simulation of a Nematogen Using a Force Field Derived from Quantum Chemical Calculations

Predicting Spontaneous Orientational Self-Assembly: In Silico Design of Materials with Quantum Mechanically Derived Force Fields

De novo design of self-assembled materials hinges upon our ability to relate macroscopic properties to individual building blocks, thus characterizing in such supramolecular architectures a wide range of observables at varied time/length scales. This work demonstrates that quantum mechanical derived force fields (QMD-FFs) do satisfy this requisite and, most importantly, do so in a predictive manner. To this end, a specific FF, built solely based on the knowledge of the target molecular structure, is employed to reproduce the spontaneous transition to an ordered liquid crystal phase. The simulations deliver a multiscale portrait of such self-assembly processes, where conformational changes within the individual building blocks are intertwined with a 200 ns ensemble reorganization. The extensive characterization provided not only is in quantitative agreement with the experiment but also connects the time/length scales at which it was performed. Realizing QMD-FF predictive power and unmatched accuracy stands as an important leap forward for the bottom-up design of advanced supramolecular materials.

Automated Parameterization of Quantum Mechanically Derived Force Fields for Soft Materials and Complex Fluids: Development and Validation

The reliability of molecular dynamics (MD) simulations in predicting macroscopic properties of complex fluids and soft materials, such as liquid crystals, colloidal suspensions, or polymers, relies on the accuracy of the adopted force field (FF). We present an automated protocol to derive specific and accurate FFs, fully based on ab initio quantum mechanical (QM) data. The integration of the Joyce and Picky procedures, recently proposed by our group to provide an accurate description of simple liquids, is here extended to larger molecules, capable of exhibiting more complex fluid phases. While the standard Joyce protocol is employed to parameterize the intramolecular FF term, a new automated procedure is here proposed to handle the computational cost of the QM calculations required for the parameterization of the intermolecular FF term. The latter is thus obtained by integrating the old Picky procedure with a fragmentation reconstruction method (FRM) that allows for a reliable, yet computationally feasible sampling of the intermolecular energy surface at the QM level. The whole FF parameterization protocol is tested on a benchmark liquid crystal, and the performances of the resulting quantum mechanically derived (QMD) FF were compared with those delivered by a general-purpose, transferable one, and by the third, "hybrid" FF, where only the bonded terms were refined against QM data. Lengthy atomistic MD simulations are carried out with each FF on extended 5CB systems in both isotropic and nematic phases, eventually validating the proposed protocol by comparing the resulting macroscopic properties with other computational models and with experiments. The QMD-FF yields the best performances, reproducing both phases in the correct range of temperatures and well describing their structure, dynamics, and thermodynamic properties, thus providing a clear protocol that may be explored to predict such properties on other complex fluids or soft materials.

Atomistic Mechanism of Anisotropic Heat Conduction in the Liquid Crystal 4-Heptyl-4'-cyanobiphenyl: All-Atom Molecular Dynamics

All-atom molecular dynamics simulations were performed on 4-heptyl-4'-cyanobiphenyl (7CB) to study the mechanism of heat conduction in this nematic liquid crystal atomistically. To describe 7CB properly, the AMBER-type force field was optimized for the dihedral parameter of biphenyl and the Lennard-Jones parameters. The molecular dynamics calculation using the optimized force field well reproduced the experimental values of the isotropic-nematic phase transition temperature, density, and anisotropy of the thermal conductivity. Furthermore, the contributions of convection, intramolecular interaction, and intermolecular interaction to the thermal conductivity were determined by performing thermal conductivity decomposition analysis. According to the analysis, the contributions of convection, bond stretching, and bond bending interactions were higher in the direction parallel to the nematic director than that perpendicular to the director, which is the origin of the anisotropy in the nematic phase. This result indicates that the anisotropy is caused by well-aligned covalent bonds and high mobility parallel to the director. This quantitative description of the mechanism of heat conduction of 7CB is foreseen to provide new insights toward designing highly thermally conductive liquid-crystalline materials.

Hierarchical bounding structures for efficient virial computations: Towards a realistic molecular description of cholesterics

We detail the application of bounding volume hierarchies to accelerate second-virial evaluations for arbitrary complex particles interacting through hard and soft finite-range potentials. This procedure, based on the construction of neighbour lists through the combined use of recursive atom-decomposition techniques and binary overlap search schemes, is shown to scale sub-logarithmically with particle resolution in the case of molecular systems with high aspect ratios. Its implementation within an efficient numerical and theoretical framework based on classical density functional theory enables us to investigate the cholesteric self-assembly of a wide range of experimentally relevant particle models. We illustrate the method through the determination of the cholesteric behavior of hard, structurally resolved twisted cuboids, and report quantitative evidence of the long-predicted phase handedness inversion with increasing particle thread angles near the phenomenological threshold value of 45°. Our results further highlight the complex relationship between microscopic structure and helical twisting power in such model systems, which may be attributed to subtle geometric variations of their chiral excluded-volume manifold.

Validating an optimized GAFF force field for liquid crystals: TNI predictions for bent-core mesogens and the first atomistic predictions of a dark conglomerate phase

Atomistic simulations of bent core mesogens provide excellent T_NI predictions and show the formation of a dark conglomerate phase.

Noncovalent Interactions in the Catechol Dimer

Noncovalent interactions play a significant role in a wide variety of biological processes and bio-inspired species. It is, therefore, important to have at hand suitable computational methods for their investigation. In this paper, we report on the contribution of dispersion and hydrogen bonds in both stacked and T-shaped catechol dimers, with the aim of delineating the respective role of these classes of interactions in determining the most stable structure. By using second-order Møller⁻Plesset (MP2) calculations with a small basis set, specifically optimized for these species, we have explored a number of significant sections of the interaction potential energy surface and found the most stable structures for the dimer, in good agreement with the highly accurate, but computationally more expensive coupled cluster single and double excitation and the perturbative triples (CCSD(T))/CBS) method.

Tuning Coulombic interactions to stabilize nematic and smectic ionic liquid crystal phases in mixtures of charged soft ellipsoids and spheres

We have investigated the effect of electrostatic interactions in mixtures of soft ellipsoids and spheres based on the well-known Gay-Berne (GB) and Lennard-Jones (LJ) potential, respectively. These model systems, in their original version, that is without any electrostatic charge, have been thoroughly investigated in the literature both as pure components and mixtures. In particular, mixtures of particles of different shapes, such as spheres and ellipsoids, tend to phase separate because of the excluded volume effects. Common ionic liquid crystals, based on imidazolium or other quaternary ammonium salts, are usually composed of roughly elongated (although flexible) cations and roughly spherical anions, that is, particles with a similar shape such as the GB and LJ models. Therefore, in this work, we present the results of molecular dynamics simulations of mixtures of positively charged GB and negatively charged LJ particles as models of ionic liquid crystals. Interestingly, by modulating the charge of the particles it is possible to stabilize isotropic, nematic, smectic and crystalline ionic phases. The relative weight of Coulomb (a radial, therefore isotropic interaction) and van der Waals (an anisotropic interaction) contributions is a key parameter to tune the stability of various mesophases.

Comparison of the accuracy of periodic reaction field methods in molecular dynamics simulations of a model liquid crystal system

A periodic reaction field (PRF) method is a technique to estimate long-range interactions. The method has the potential to effectively reduce the computational cost while maintaining adequate accuracy. We performed molecular dynamics (MD) simulations of a model liquid-crystal system to assess the accuracy of some variations of the PRF method in low-charge-density systems. All the methods had adequate accuracy compared with the results of the particle mesh Ewald (PME) method, except for a few simulation conditions. Furthermore, in all of the simulation conditions, one of the PRF methods had the same accuracy as the PME method.

Parametrization and Validation of Intramolecular Force Fields Derived from DFT Calculations

The energy and its first and second geometrical derivatives obtained by DFT calculations for a number of conformations of a single molecule are used to parametrize intramolecular force fields, suitable for computer simulations. A systematic procedure is proposed to adequately treat either fully atomistic or more simplified force fields, as within the united atom approach or other coarse grained models. The proposed method is tested and validated by performing molecular dynamics simulations on several different molecules, comparing the results with literature force fields and relevant experimental data. Particular emphasis is given to the united atom approach for flexible molecules characterized by "soft" torsional potentials which are known to retain a high degree of chemical specificity.

Optimization of the GAFF force field to describe liquid crystal molecules: the path to a dramatic improvement in transition temperature predictions

Systematic optimization of the General Amber Force Field (GAFF) for mesogenic fragments leads to a dramatic improvement in the modelling of liquid crystal clearing points.

Non-Newtonian rheological properties of shearing nematic liquid crystal model systems based on the Gay–Berne potential

The logarithm of the viscosity of a nematic liquid crystal is a linear function of the square root of the shear rate in the non-Newtonian regime.

Comparison of the dielectric and NMR results for liquid crystals: dynamic aspects

Critical analysis of the results of studies of molecular rotational dynamics in liquid crystalline substances with the aid of the dielectric spectroscopy (DS) and nuclear magnetic resonance (NMR) is given. Both methods are known to be sensitive to different aspects of molecular rotations: the polarization vector and the relaxation time τ(DS) in the case of DS, a tensor describing a nuclear interaction and the correlation time τ(NMR) for NMR method. Furthermore, both methods provide correlation functions with different rank values. A common basis for the comparison between τ(DS) and τ(NMR) is postulated. Several examples of the temperature dependence of the correlation times coming from the two spectroscopic methods are presented. Qualitative agreements of the correlation times were achieved in most cases.

Joyce and Ulysses: integrated and user-friendly tools for the parameterization of intramolecular force fields from quantum mechanical data

The Joyce program is augmented with several new features, including the user friendly Ulysses GUI, the possibility of complete excited state parameterization and a more flexible treatment of the force field electrostatic terms. A first validation is achieved by successfully comparing results obtained with Joyce2.0 to literature ones, obtained for the same set of benchmark molecules. The parameterization protocol is also applied to two other larger molecules, namely nicotine and a coumarin based dye. In the former case, the parameterized force field is employed in molecular dynamics simulations of solvated nicotine, and the solute conformational distribution at room temperature is discussed. Force fields parameterized with Joyce2.0, for both the dye's ground and first excited electronic states, are validated through the calculation of absorption and emission vertical energies with molecular mechanics optimized structures. Finally, the newly implemented procedure to handle polarizable force fields is discussed and applied to the pyrimidine molecule as a test case.

An automated approach for the parameterization of accurate intermolecular force-fields: pyridine as a case study

An automated protocol is proposed and validated, which integrates accurate quantum mechanical calculations with classical numerical simulations. Intermolecular force fields, (FF) suitable for molecular dynamics (MD) and Monte Carlo simulations, are parameterized through a novel iterative approach, fully based on quantum mechanical data, which has been automated and coded into the PICKY software, here presented. The whole procedure is tested and validated for pyridine, whose bulk phase, described through MD simulations performed with the specifically parameterized FF, is characterized by computing several of its thermodynamic, structural, and transport properties, comparing them with their experimental counterparts.

Liquid Crystal Properties of the n-Alkyl-cyanobiphenyl Series from Atomistic Simulations with Ab Initio Derived Force Fields

Lengthy molecular dynamics (MD) simulations were performed at constant atmospheric pressure and different temperatures for the series of the 4-n-alkyl-4'-cyanobiphenyls (nCB) with n = 6, 7, and 8. The accurate atomistic force field (Bizzarri, M.; Cacelli, I.; Prampolini, G; Tani, A. J. Phys. Chem. A 2004, 108, 10336), successfully employed to reproduce thermodynamic and transport properties of the 5CB molecule, has here been extended to higher homologues. Nematic and isotropic phases were found for all members of the series, and also, a smectic phase was (tentatively) identified for 8CB at 1 atm and 300 K. Transition temperatures reproduce the experimental values within +/-10 K. Also, structural properties as second and fourth rank orientational order parameters are in good agreement with the corresponding experimental quantities. This means that the well-known odd-even effect, observed for many properties along the nCB series, is well reproduced, despite the narrow range of oscillations, e.g., in clearing temperatures. A detailed analysis of the correlation between molecular properties and odd-even effects is presented.

Atomistic simulation of a model liquid crystal

We present atomistic molecular dynamics computer simulations of the bulk phases of a model liquid crystal system based on 8CB. The model differs from real 8CB because it employs a united-atom description to eliminate all hydrogen atoms, and neglects all long-range electrostatic interactions. Despite this simplification, the pressure-temperature phase diagram shows an order-disorder transition, in which isotropic, smectic, and nematiclike behaviors are observed. A detailed analysis of the inter- and intramolecular structures of the ordered phases is given, together with an examination of finite size effects and the equilibration times of the system. It is shown that, whereas a system may appear to be thermodynamically and mechanically equilibrated after a period of 10-15 ns, it is possible for an imprint of the starting configuration to persist for much longer time scales. In the present case, however, such an imprint does not appear to affect the observed phase behavior.

Modeling a Liquid Crystal Dynamics by Atomistic Simulation with an Ab Initio Derived Force Field

Atomistic molecular dynamics (MD) simulations of 4-n-pentyl 4'-cyano-biphenyl (5CB) have been performed, adopting a specific ab initio derived force field. Two state points in the nematic phase and three in the isotropic phase, as determined in a previous work, have been considered. At each state point, at least 10 ns have been produced, allowing us to accurately calculate single-molecule properties. In the isotropic phase, the values of the translational diffusion coefficient, and even more so the activation energy for the process, agree well with experimental data. Qualitatively, also the dynamic anisotropy of the nematic phase is correctly accounted for. Rotational diffusion coefficients, which describe spinning and tumbling motions, fall well within the range of experimental values. The reorientational dynamics of our model 5CB covers diverse time regimes. The longest one is strongly temperature dependent and characterized by a relaxation time in accord with experimental dielectric relaxation data. Shear viscosity and Landau-de Gennes relaxation times, typically collective variables, reproduce the experimental results very well in the isotropic phase. In the nematic phase, despite a large statistical uncertainty due to the extremely slow relaxation of the correlation functions involved, our simulation yields the correct relative order of the three experimental Miesowicz viscosities.