Communication: molecular dynamics and (1)H NMR of n-hexane in liquid crystals.

The NMR spectrum of n-hexane orientationally ordered in the nematic liquid crystal ZLI-1132 is analysed using covariance matrix adaptation evolution strategy (CMA-ES). The spectrum contains over 150 000 transitions, with many sharp features appearing above a broad, underlying background signal that results from the plethora of overlapping transitions from the n-hexane as well as from the liquid crystal. The CMA-ES requires initial search ranges for NMR spectral parameters, notably the direct dipolar couplings. Several sets of such ranges were utilized, including three from MD simulations and others from the modified chord model that is specifically designed to predict hydrocarbon-chain dipolar couplings. In the end, only inaccurate dipolar couplings from an earlier study utilizing proton-proton double quantum 2D-NMR techniques on partially deuterated n-hexane provided the necessary estimates. The precise set of dipolar couplings obtained can now be used to investigate conformational averaging of n-hexane in a nematic environment.

A model-free temperature-dependent conformational study of n-pentane in nematic liquid crystals.

The proton NMR spectra of n-pentane orientationally ordered in two nematic liquid-crystal solvents are studied over a wide temperature range and analysed using covariance matrix adaptation evolutionary strategy. Since alkanes possess small electrostatic moments, their anisotropic intermolecular interactions are dominated by short-range size-and-shape effects. As we assumed for n-butane, the anisotropic energy parameters of each n-pentane conformer are taken to be proportional to those of ethane and propane, independent of temperature. The observed temperature dependence of the n-pentane dipolar couplings allows a model-free separation between conformer degrees of order and conformer probabilities, which cannot be achieved at a single temperature. In this way for n-pentane 13 anisotropic energy parameters (two for trans trans, tt, five for trans gauche, tg, and three for each of gauche+ gauche+, pp, and gauche+ gauche-, pm), the isotropic trans-gauche energy difference Etg and its temperature coefficient Etg(') are obtained. The value obtained for the extra energy associated with the proximity of the two methyl groups in the gauche+ gauche- conformers (the pentane effect) is sensitive to minute details of other assumptions and is thus fixed in the calculations. Conformer populations are affected by the environment. In particular, anisotropic interactions increase the trans probability in the ordered phase.

Conformational problem of alkanes in liquid crystals by NMR spectroscopy: a mini-review.

Recent discoveries of the role of alkane flexibility in determining liquid-crystal behaviour are surveyed. With the impetus for understanding the alkane conformational problem established, recent model dependent (1)H NMR work on the topic will be reviewed where progress is made but the need to circumvent models eventually becomes evident. A closer look at the rigid basic units of alkanes will provide the way forward where it is shown that the orientational ordering and anisotropic potentials of these molecules dissolved in liquid crystals scale with each other. Once this relationship is established, a series of works using anisotropic and isotropic (1)H NMR spectroscopy to study alkane conformational statistics will be covered, wherein the influence of the gas, isotropic condensed and anisotropic condensed phases will be described.

NMR of short-chain hydrocarbons in nematic and smectic a liquid crystals.

The proton NMR spectra of ethane, propane, and n-butane codissolved and orientationally ordered in two liquid-crystal solvents that exhibit both nematic (N) and smectic A (SmA) phases (one of which also has a reentrant nematic (RN) phase) are analyzed using CMA-ES (covariance-matrix adaptation evolution strategy). To deal with problems arising from the broad liquid-crystal background signal, a smoothed experimental spectrum is fitted to a smoothed calculated spectrum. The ethane and propane dipolar couplings and anisotropic energy parameters scale with each other, and the n-butane couplings are assumed to behave likewise. This restriction on the relative values of the n-butane energy parameters facilitates a fit to the temperature dependence of the n-butane dipolar couplings, in which the six n-butane energy parameters (three for trans and three for gauche), the isotropic trans-gauche energy difference E(tg), and its temperature coefficient E(tg)', and the methyl CCH angle decrease are obtained. Unlike earlier studies, the fit does not employ a model for the anisotropic intermolecular potential. The nematic potential in the SmA phase of the liquid-crystal mixture 6OCB/8OCB is estimated by interpolated values obtained from the ethane spectra for the N and RN regions. Analysis of results for the SmA phase involves adding a nematic-smectic coupling prefactor and a smectic prefactor. Spectra obtained in the liquid-crystal 8OCB are calculated with no adjustable parameters, scaling the potentials using the ethane values, and the agreement between experiment and calculation is outstanding. Conformer populations are affected by the environment, and the isotropic solute-solvent interactions result in increased gauche populations, whereas anisotropic interactions increase the trans probability. The trans probability in the SmA phase is slightly lower than expected from the nematic results.

Efficient analysis of highly complex nuclear magnetic resonance spectra of flexible solutes in ordered liquids by using molecular dynamics.

The NMR spectra of n-pentane as solute in the liquid crystal 5CB are measured at several temperatures in the nematic phase. Atomistic molecular dynamics simulations of this system are carried out to predict the dipolar couplings of the orientationally ordered pentane, and the spectra predicted from these simulations are compared with the NMR experimental ones. The simulation predictions provide an excellent starting point for analysis of the experimental NMR spectra using the covariance matrix adaptation evolutionary strategy. This shows both the power of atomistic simulations for aiding spectral analysis and the success of atomistic molecular dynamics in modeling these anisotropic systems.

Nuclear magnetic resonance study of alkane conformational statistics.

NMR spectra of ethane, propane, and n-butane as solutes in the nematic liquid crystals 4-n-pentyl-4(')-cyanobiphenyl (5CB) and Merck ZLI 1132 (1132) are investigated over a wide temperature range. The ratios of dipolar couplings of ethane to propane are constant over the entire temperature range. Assuming that this constancy applies to the butane conformers facilitates the separation of probability from order parameter. This separation allows the investigation of conformational distribution without the need of invoking any model for the anisotropic intermolecular potential. The results give an order matrix that is consistent with that predicted from model potentials that describe the orientational potential in terms of short-range size and shape effects. The isotropic intermolecular potential contribution to the trans-gauche energy difference E(tg) is found to be temperature dependent with the values and variation in agreement with that found when the same results are analyzed using the chord model for anisotropic interactions [A. C. J. Weber and E. E. Burnell, Chem. Phys. Lett. 506, 196 (2011)]. The fit obtained for 9 spectra in 5CB (63 dipolar couplings) has an RMS difference between experimental and calculated dipolar couplings of 2.7 Hz, while that for the 16 spectra in 1132 (112 couplings) is 6.2 Hz; this excellent fit with nine adjustable parameters suggests that the assumption of equal temperature dependencies of the order parameters for ethane, propane, and each conformer of butane is correct. Also the fit parameters (E(tg) and the methyl angle increase) obtained for 1132 and 5CB agree. The results indicate that the chord model, which was designed to treat hydrocarbon chains, is indeed the model of choice for these chains. The temperature variation of E(tg) provides a challenge for theoreticians. Finally, even better fits to the experimental dipolar couplings are obtained when the energy in the Boltzmann factor is used for scaling ethane to butane results. However, in this case the values obtained for E(tg) differ between 1132 and 5CB.

Scope and limitations of accurate structure determination of solutes dissolved in liquid crystals.

The determination of accurate structures of relatively small molecules dissolved in liquid-crystal solvents is no trivial matter. Extensive vibrational corrections to the observed dipolar couplings are required. Vibrational force fields are often available, but the usual harmonic corrections are strictly limited to small-amplitude internal motions. Moreover, the need to also apply anharmonic corrections and to include the elusive vibration-reorientation interaction is problematic and can only be fulfilled for a very limited set of small molecules. In this paper we discuss the implications for the accuracy of the structure of larger molecules for which this information is not available. We discuss the examples of azulene and biphenylene, and set realistic limits on their proton structures, derived from (1)H dipolar couplings extracted from NMR spectra obtained in different liquid-crystal solvents and analyzed with sophisticated evolutionary algorithms.

The smectic effect on solute order parameters rationalized by double Maier-Saupe Kobayashi-McMillan theory.

From the dipolar couplings obtained by NMR spectroscopy we have calculated the order parameters of a wide variety of solutes in the nematic and smectic A phases of the liquid crystals 8CB and 8OCB. These measurements are then rationalized with the previously tested two Maier-Saupe Kobayashi-McMillan interaction potential from which smectic order parameters are calculated.

NMR determination of smectic ordering of probe molecules.

The NMR spectra of the three solutes ortho-, meta-, and para-dichlorobenzene in the nematic and smectic A phases of the liquid crystals 8CB and 8OCB are analyzed to yield two orientational order parameters for each solute. Extrapolation of the asymmetry in the energy parameters that describe the orientational ordering in the nematic phase are used to provide estimates of the strength of the nematic potential in the smectic A phase. The experimentally determined asymmetry of the orientational order parameters in the smectic A phase is then used in conjunction with Kobayashi-McMillan theory applied to solutes to give information about the smectic A layering and the nematic/smectic A coupling. In both smectic A solvents, the solute smectic coupling constant, tau, is negative (with the origin fixed at the center of the smectic layer) for all solutes. The signs and relative values of tau indicate that the ortho and para solutes favor the interlayer region while the meta solute is more evenly distributed throughout the layers.