Ab initio molecular dynamics of liquid hydrogen chloride

Anharmonic effects in the quantum cluster equilibrium method

The well-established quantum cluster equilibrium (QCE) model provides a statistical thermodynamic framework to apply high-level ab initio calculations of finite cluster structures to macroscopic liquid phases using the partition function. So far, the harmonic approximation has been applied throughout the calculations. In this article, we apply an important correction in the evaluation of the one-particle partition function and account for anharmonicity. Therefore, we implemented an analytical approximation to the Morse partition function and the derivatives of its logarithm with respect to temperature, which are required for the evaluation of thermodynamic quantities. This anharmonic QCE approach has been applied to liquid hydrogen chloride and cluster distributions, and the molar volume, the volumetric thermal expansion coefficient, and the isobaric heat capacity have been calculated. An improved description for all properties is observed if anharmonic effects are considered.

Aqueous solutions: state of the art in ab initio molecular dynamics

The simulation of liquids by ab initio molecular dynamics (AIMD) has been a subject of intense activity over the last two decades. The significant increase in computational resources as well as the development of new and efficient algorithms has elevated this method to the status of a standard quantum mechanical tool that is used by both experimentalists and theoreticians. As AIMD computes the electronic structure from first principles, it is free of ad hoc parametrizations and has thus been applied to a large variety of physical and chemical problems. In particular, AIMD has provided microscopic insight into the structural and dynamical properties of aqueous solutions which are often challenging to probe experimentally. In this review, after a brief theoretical description of the Born-Oppenheimer and Car-Parrinello molecular dynamics formalisms, we show how AIMD has enhanced our understanding of the properties of liquid water and its constituent ions: the proton and the hydroxide ion. Thereafter, a broad overview of the application of AIMD to other aqueous systems, such as solvated organic molecules and inorganic ions, is presented. We also briefly describe the latest theoretical developments made in AIMD, such as methods for enhanced sampling and the inclusion of nuclear quantum effects.

Statistical thermodynamics of biomembranes

An overview of the major issues involved in the statistical thermodynamic treatment of phospholipid membranes at the atomistic level is summarized: thermodynamic ensembles, initial configuration (or the physical system being modeled), force field representation as well as the representation of long-range interactions. This is followed by a description of the various ways that the simulated ensembles can be analyzed: area of the lipid, mass density profiles, radial distribution functions (RDFs), water orientation profile, deuterium order parameter, free energy profiles and void (pore) formation; with particular focus on the results obtained from our recent molecular dynamic (MD) simulations of phospholipids interacting with dimethylsulfoxide (Me(2)SO), a commonly used cryoprotective agent (CPA).

Kohn-Sham-like approach toward a classical density-functional theory of inhomogeneous polar molecular liquids: an application to liquid hydrogen chloride

The Gordian knot of density-functional theories for classical molecular liquids remains finding an accurate free-energy functional in terms of the densities of the atomic sites of the molecules. Following Kohn and Sham, we show how to solve this problem by considering noninteracting molecules in a set of effective potentials. This shift in perspective leads to an accurate and computationally tractable description in terms of simple three-dimensional functions. We also treat both the linear- and saturation- dielectric responses of polar systems, presenting liquid hydrogen chloride as a case study.

A molecular site-site integral equation that yields the dielectric constant

Our recent derivation [K. M. Dyer et al., J. Chem. Phys. 127, 194506 (2007)] of a diagrammatically proper, site-site, integral equation theory using molecular angular expansions is extended to polar fluids. With the addition of atomic site charges we take advantage of the formal long-ranged potential field cancellations before renormalization to generate a set of numerically stable equations. Results for calculations in a minimal (spherical) angular basis set are presented for the radial distribution function, the first dipolar (110) projection, and the dielectric constant for two model diatomic systems. All results, when compared to experiment and simulation, are a significant quantitative and qualitative improvement over previous site-site theories. More importantly, the dielectric constant is not trivial and close to simulation and experiment.