Liquid dibromomethane under pressure: a computational study.

Molecular dynamics simulations have been performed on liquid dibromomethane at thermodynamic states corresponding to temperature in the range 268-328 K and pressure varying from 1 bar to 3000 bar. The interaction model is a simple effective two-body pair potential with atom-atom Coulomb and Lennard-Jones interactions and molecules are rigid. Thermodynamic properties have been studied, including the isobaric thermal expansion coefficient, the isothermal compressibility, the heat capacities and the speed of sound. The simulation results exhibit a crossing of the isotherms of the isobaric thermal expansion coefficient at about 800 bar in very good agreement with the prediction of an isothermal fluctuation equation of state predicting such a crossing in the pressure range 650-900 bar, though experimental results up to 1000 bar do not find any crossing.

Translational and Rotational Diffusion in Liquid Water at Very High Pressure: A Simulation Study.

Translational and rotational diffusion coefficients of liquid water have been computed from molecular dynamics simulation with a recent polarizable potential at 298, 400, and 550 K at very high pressure. At 298 K, the model reproduces the initial increase and the occurrence of a maximum for the translational and rotational diffusion coefficients when the pressure increases. At 400 and 550 K, translational and rotational diffusion coefficients are found to monotonically decrease when pressure increases in the gigapascal range, with the translational coefficient decreasing faster than the rotational one. At 400 K, such an evolution of the rotational diffusion coefficient contrasts with quasielastic neutron scattering results predicting a near independence of the rotational diffusion with a pressure increase above ≃0.5 GPa. An interpretation is proposed to explain this discrepancy. The pressure dependence of the translation-rotation coupling is analyzed. The anisotropy of rotational diffusion is investigated by computing the rotational diffusion tensor in a molecular system of axes and the reorientational correlation times of rank 1 and rank 2 of the inertia axes and of the OH bond vector. Deviation of the simulation data with respect to the predictions of the isotropic Debye model of rotational diffusion are quantified and can be used to estimate experimental rotational diffusion coefficients from experimental reorientational correlation times.

Dipole and quadrupole polarizabilities of the water molecule as a function of geometry.

Dipolar, dipole-quadrupole and quadrupole-quadrupole static polarizabilities of the water molecule have been determined by ab initio calculations at coupled cluster level of theory with single, double and perturbative triple excitations CCSD(T) with an aug-cc-pVTZ basis set using a finite field and field-gradient method. The geometry dependence of polarizability tensor components has been explored and modeled by power series expansion in bond length and angle variations up to sum of powers equal to 4. The results provide a very detailed description of the static polarizability of water up to quadrupolar rank which can be used for the test and development of novel accurate polarizable interaction potentials for modeling aqueous solutions. © 2016 Wiley Periodicals, Inc.

Geometry-dependent distributed polarizability models for the water molecule.

Geometry-dependent distributed polarizability models have been constructed by fits to ab initio calculations at the coupled cluster level of theory with up to noniterative triple excitations in an augmented triple-zeta quality basis set for the water molecule in the field of a point charge. The investigated models include (i) charge-flow polarizabilities between chemically bonded atoms, (ii) isotropic or anisotropic dipolar polarizabilities on oxygen atom or on all atoms, and (iii) combinations of models (i) and (ii). For each model, the polarizability parameters have been optimized to reproduce the induction energy of a water molecule polarized by a point charge successively occupying a grid of points surrounding the molecule. The quality of the models is ascertained by examining their ability to reproduce these induction energies as well as the molecular dipolar and quadrupolar polarizabilities. The geometry dependence of the distributed polarizability models has been explored by changing bond lengths and HOH angle to generate 125 molecular structures (reduced to 75 symmetry-unique ones). For each considered model, the distributed polarizability components have been fitted as a function of the geometry by a Taylor expansion in monomer coordinate displacements up to the sum of powers equal to 4.

Distributed polarizability models for imidazolium-based ionic liquids.

Quantum chemical calculations are used to derive distributed polarizability models sufficiently accurate and compact to be used in classical molecular dynamics simulations of imidazolium-based room temperature ionic liquids. Two distributed polarizability models are fitted to reproduce the induction energy of three imidazolium cations (1,3-dimethyl-, 1-ethyl-3-methyl-, and 1-butyl-3-methylimidazolium) and four anions (tetrafluoroborate, hexafluorophosphate, nitrate, and thiocyanate) polarized by a point charge located successively on a grid of surrounding points. The first model includes charge-flow polarizabilities between first-neighbor atoms and isotropic dipolar polarizability on all atoms (except H), while the second model includes anisotropic dipolar polarizabilities on all atoms (except H). For the imidazolium cations, particular attention is given to the transferability of the distributed polarizability sets. The molecular polarizability and its anisotropy rebuilt by the distributed models are found to be in good agreement with the exact ab initio values for the three cations and 23 additional conformers of 1-ethyl-3-methyl-, 1-butyl-3-methyl-, 1-pentyl-3-methyl-, and 1-hexyl-3-methylimidazolium cations.

Molecular dynamics simulations of quinoline in the liquid phase.

Molecular dynamics simulations of liquid quinoline have been performed at experimental densities corresponding to the temperature range 276-320 K. The intermolecular potential is a simple effective two-body potential between rigid molecules having 17 atomic Lennard-Jones and electrostatic Coulomb interaction sites. The vaporization enthalpy is overestimated by 8-9% with respect to the experimental value. The translational diffusion coefficient exhibits a small non-Arrhenius behavior with a change in temperatures near 290 and 303 K. The rotational diffusion tensor is rotated around the z axis perpendicular to the molecular plane by an angle of 4-6° with respect to the frame of reference defined by the principal axes of inertia. The rotational diffusion tensor presents a significant anisotropy with D(rot,y)/D(rot,x) ≃ 0.6-0.5 and D(rot,z)/D(rot,x) ≃ 1.6-1.3 between 276 and 320 K when the x axis is defined as the long molecular axis and the y axis is situated nearly along the central C-C bond. The rotational diffusion coefficients, the reorientational correlation times of the C-H vectors, and the T1(13)C NMR relaxation times present a non-Arrhenius break around 288-290 K in agreement with several experimental results. In addition, a non-Arrhenius break can also be observed at 303 K for these properties. It has been found that the structure evolves smoothly in the studied temperature range. Center of mass-center of mass and atom-atom radial distribution functions show a monotonous evolution with temperature. Various types of first-neighbor dimers have been defined, and their population analysis has revealed a continuous monotonous evolution with temperature. Thus, the non-Arrhenius behavior observed for translational and rotational diffusion is correlated with the monotonous evolution of the population of first-neighbor dimers at a microscopic level and not with a sharp structural transition.

Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: a sequential molecular dynamics/quantum mechanical approach.

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.

Hit and miss of classical nucleation theory as revealed by a molecular simulation study of crystal nucleation in supercooled sulfur hexafluoride.

Classical nucleation theory pictures the homogeneous nucleation of a crystal as the formation of a spherical crystalline embryo, possessing the properties of the macroscopic crystal, inside a parent supercooled liquid. In this work we study crystal nucleation in moderately supercooled sulfur hexafluoride by umbrella sampling simulations. The nucleation free energy evolves from 5.2kBT at T=170 K to 39.1kBT at T=195 K. The corresponding critical nucleus size ranges from 40 molecules at T=170 K to 266 molecules at T=195 K. Both nucleation free energy and critical nucleus size are shown to evolve with temperature according to the equations derived from the classical nucleation theory. Inspecting the obtained nuclei we show, however, that they present quite anisotropic shapes in opposition to the spherical assumption of the theory. Moreover, even though the critical nuclei possess the structure of the stable bcc plastic phase, the only mechanically stable crystal phase for SF6 in the temperature range investigated, they are shown to be less ordered than the corresponding macroscopic crystal. Their crystalline order is nevertheless shown to increase regularly with their size. This is confirmed by a study of a nucleus growth from a critical size to a size of the order of 10(4) molecules. Similarly to the fact that it does not affect the temperature dependence of the nucleation free energy and of the critical nucleus size, the ordering of the nucleus with size does not affect the growth rate of the nucleus.

Long-range electrostatic interactions in hybrid quantum and molecular mechanical dynamics using a lattice summation approach.

A method based on a lattice summation technique for treating long-range electrostatic interactions in hybrid quantum mechanics/molecular mechanics simulations is presented in this article. The quantum subsystem is studied at the semiempirical level, whereas the solvent is described by a two-body potential of molecular mechanics. Molecular dynamics simulations of a (quantum) chloride ion in (classical) water have been performed to test this technique. It is observed that the application of the lattice summations to solvent-solvent interactions as well as on solute-solvent ones has a significant effect on solvation energy and diffusion coefficient. Moreover, two schemes for the computation of the long-range contribution to the electrostatic interaction energy are investigated. The first one replaces the exact charge distribution of the quantum solute by a Mulliken charge distribution. The long-range electrostatic interactions are then calculated for this charge distribution that interacts with the solvent molecule charges. The second one is more accurate and involves a modified Fock operator containing long-range electron-charge interactions. It is shown here that both schemes lead to similar results, the method using Mulliken charges for the evaluation of long-range interactions being, however, much more computationally efficient.

Molecular simulation of the homogeneous crystal nucleation of carbon dioxide.

We report on a molecular simulation study of the homogeneous nucleation of CO2 in the supercooled liquid at low pressure (P = 5 MPa) and for degrees of supercooling ranging from 32% to 60%. In all cases, regardless of the degree of supercooling, the structure of the crystal nuclei is that of the Pa3 phase, the thermodynamically stable phase. For the more moderate degree of supercooling of 32%, the nucleation is an activated process and requires a method to sample states of high free energy. In this work, we apply a series of bias potentials, which promote the ordering of the centers of mass of the molecules and allow us to gradually grow crystal nuclei. The reliability of the results so obtained is assessed by studying the evolution of the nuclei in the absence of any bias potential, and by determining their probability of growth. We estimate that the size of the critical nucleus, for which the probability of growth is 0.5, is approximately 240 molecules. Throughout the nucleation process, the crystal nuclei clearly exhibit a Pa3 structure, in apparent contradiction with Ostwald's rule of stages. The other polymorphs have a much larger free energy. This makes their formation highly unlikely and accounts for the fact that the nucleation of CO2 proceeds directly in the stable Pa3 structure.

Atomistic simulation of the homogeneous nucleation and of the growth of N2 crystallites.

We report on a computer simulation study of the early stages of the crystallization of molecular nitrogen. First, we study how homogeneous nucleation takes place in supercooled liquid N(2) for a moderate degree of supercooling. Using the umbrella sampling technique, we determine the free energy barrier of formation for a critical nucleus of N(2). We show that, in accord with Ostwald's rule of stages, the structure of the critical nucleus is predominantly that of a metastable polymorph (alpha-N(2) for the state point investigated). We then monitor the evolution of several critical nuclei through a series of unbiased molecular dynamics trajectories. The growth of N(2) crystallites is accompanied by a structural evolution toward the stable polymorph beta-N(2). The microscopic mechanism underlying this evolution qualitatively differs from that reported previously. We do not observe any dissolution or reorganization of the alpha-like core of the nucleus. On the contrary, we show that alpha-like and beta-like blocks coexist in postcritical nuclei. We relate the structural evolution to a greater adsorption rate of beta-like molecules on the surface and show that this transition actually starts well within the precritical regime. We also carefully investigate the effect of the system size on the height of the free energy barrier of nucleation and on the structure and size of the critical nucleus.

Reorganization and growth of metastable alpha-N2 critical nuclei into stable beta-N2 crystals.

We report on results on the crystal nucleation and growth of nitrogen. Using molecular dynamics simulations, we show that while nucleation proceeds into the metastable alpha-phase (i.e., the crystalline phase associated with the lowest free energy barrier of formation), growth of the crystallite proceeds through a reorganization of the nucleus into the thermodynamically stable beta-phase.

OPEP: a tool for the optimal partitioning of electric properties.

OPEP is a suite of FORTRAN programs targeted at the optimal partitioning of molecular electric properties. It includes an interactive module for the construction of Cartesian grids of points, on which either the molecular electrostatic potential or the induction energy is mapped. The generation of distributed multipoles and polarizabilities is achieved using either the formalism of the normal equations of the least-squares problem, which restates the fitting procedure in terms of simple matrix operations, or a statistical approach, which provides a pictorial description of the distributed models of multipoles and polarizabilities, thereby allowing the pinpointing of pathological cases. Molecular symmetry is accounted for by means of local atomic frames, which are generated in an automated fashion. A Tcl/Tk graphical user interface wraps the suite of programs, thereby making OPEP a user-friendly package for building models of distributed multipoles and polarizabilities. OPEP is an open-source suite of programs distributed free of charge under the GNU general public license (GPL) at http://www.lctn.uhp-nancy.fr/Opep.