Role of three-body interactions in formation of bulk viscosity in liquid argon

With the aim of locating the origin of discrepancy between experimental and computer simulation results on bulk viscosity of liquid argon, a molecular dynamic simulation of argon interacting via ab initio pair potential and triple-dipole three-body potential has been undertaken. Bulk viscosity, obtained using Green-Kubo formula, is different from the values obtained from modeling argon using Lennard-Jones potential, the former being closer to the experimental data. The conclusion is made that many-body inter-atomic interaction plays a significant role in formation of bulk viscosity.Comment: 4 pages, 3 figure

Correlations in liquid water for the TIP3P-Ewald, TIP4P-2005, TIP5P-Ewald, and SWM4-NDP models.

Water is one of the simplest molecules in existence, but also one of the most important in biological and engineered systems. However, understanding the structure and dynamics of liquid water remains a major scientific challenge. Molecular dynamics simulations of liquid water were performed using the water models TIP3P-Ewald, TIP4P-2005, TIP5P-Ewald, and SWM4-NDP to calculate the radial distribution functions (RDFs), the relative angular distributions, and the excess enthalpies, entropies, and free energies. In addition, lower-order approximations to the entropy were considered, identifying the fourth-order approximation as an excellent estimate of the full entropy. The second-order and third-order approximations are ~20% larger and smaller than the true entropy, respectively. All four models perform very well in predicting the radial distribution functions, with the TIP5P-Ewald model providing the best match to the experimental data. The models also perform well in predicting the excess entropy, enthalpy, and free energy of liquid water. The TIP4P-2005 and SWM4-NDP models are more accurate than the TIP3P-Ewald and TIP5P-Ewald models in this respect. However, the relative angular distribution functions of the four water models reveal notable differences. The TIP5P-Ewald model demonstrates an increased preference for water molecules to act both as tetrahedral hydrogen bond donors and acceptors, whereas the SWM4-NDP model demonstrates an increased preference for water molecules to act as planar hydrogen bond acceptors. These differences are not uncovered by analysis of the RDFs or the commonly employed tetrahedral order parameter. However, they are expected to be very important when considering water molecules around solutes and are thus a key consideration in modelling solvent entropy.Acknowledgements go to Mike Payne for careful reading of the paper; Peter Freddolino, Chris Baker, David Payne, and Bracken King for helpful discussions; Stuart Rankin for technical help; and the NVIDIA CUDA Centre of Excellence at the Cambridge HPCS for use of the CUDA-accelerated GPUs. Thanks also go to the reviewers for their helpful comments. All calculations were performed using the Darwin Supercomputer of the University of Cambridge High Performance Computing Service (http://www.hpc.cam.ac.uk/) provided by Dell Inc. using Strategic Research Infrastructure Funding from the Higher Education Funding Council for England and were funded by the Engineering and Physical Sciences Research Consul (United Kingdom) (EPSRC) under Grant No. EP/F032773/1. Thanks for financial support go to the MRC, Wellcome Trust, and EPSRC.This is the author accepted manuscript. The final version is available from AIP Publishing via http://dx.doi.org/10.1063/1.368344

Mires in the Maluti Mountains of Lesotho.

Contains fulltext : 149842.pdf (publisher's version ) (Open Access)11 p

Full optimization of Jastrow-Slater wave functions with application to the first-row atoms and homonuclear diatomic molecules

We pursue the development and application of the recently-introduced linear optimization method for determining the optimal linear and nonlinear parameters of Jastrow-Slater wave functions in a variational Monte Carlo framework. In this approach, the optimal parameters are found iteratively by diagonalizing the Hamiltonian matrix in the space spanned by the wave function and its first-order derivatives, making use of a strong zero-variance principle. We extend the method to optimize the exponents of the basis functions, simultaneously with all the other parameters, namely the Jastrow, configuration state function and orbital parameters. We show that the linear optimization method can be thought of as a so-called augmented Hessian approach, which helps explain the robustness of the method and permits us to extend it to minimize a linear combination of the energy and the energy variance. We apply the linear optimization method to obtain the complete ground-state potential energy curve of the C_2 molecule up to the dissociation limit, and discuss size consistency and broken spin-symmetry issues in quantum Monte Carlo calculations. We perform calculations of the first-row atoms and homonuclear diatomic molecules with fully optimized Jastrow-Slater wave functions, and we demonstrate that molecular well depths can be obtained with near chemical accuracy quite systematically at the diffusion Monte Carlo level for these systems.Comment: 15 pages, 3 figures, to appear in Journal of Chemical Physic

A discrete model of water with two distinct glassy phases

We investigate a minimal model for non-crystalline water, defined on a Husimi lattice. The peculiar random-regular nature of the lattice is meant to account for the formation of a random 4-coordinated hydrogen-bond network. The model turns out to be consistent with most thermodynamic anomalies observed in liquid and supercooled-liquid water. Furthermore, the model exhibits two glassy phases with different densities, which can coexist at a first-order transition. The onset of a complex free-energy landscape, characterized by an exponentially large number of metastable minima, is pointed out by the cavity method, at the level of 1-step replica symmetry breaking.Comment: expanded version: 6 pages, 7 figure

Addressing hysteresis and slow equilibration issues in cavity-based calculation of chemical potentials.

In this paper, we explore the strengths and weaknesses of a cavity-based method to calculate the excess chemical potential of a large molecular solute in a dense liquid solvent. Use of the cavity alleviates some technical problems associated with the appearance of (integrable) divergences in the integrand during alchemical particle growth. The excess chemical potential calculated using the cavity-based method should be independent of the cavity attributes. However, the performance of the method (equilibration time and the robustness) does depend on the cavity attributes. To illustrate the importance of a suitable choice of the cavity attributes, we calculate the partition coefficient of pyrene in toluene and heptane using a coarse-grained model. We find that a poor choice for the functional form of the cavity may lead to hysteresis between growth and shrinkage of the cavity. Somewhat unexpectedly, we find that, by allowing the cavity to move as a pseudo-particle within the simulation box, the decay time of fluctuations in the integrand of the thermodynamic integration can be reduced by an order of magnitude, thereby increasing the statistical accuracy of the calculation.BP ICA

Predicting vapor liquid equilibria using density functional theory: a case study of argon

Predicting vapor liquid equilibria (VLE) of molecules governed by weak van der Waals (vdW) interactions using the first principles approach is a significant challenge. Due to the poor scaling of the post Hartree-Fock wave function theory with system size/basis functions, the Kohn-Sham density functional theory (DFT) is preferred for systems with a large number of molecules. However, traditional DFT cannot adequately account for medium to long range correlations which are necessary for modeling vdW interactions. Recent developments in DFT such as dispersion corrected models and nonlocal van der Waals functionals have attempted to address this weakness with a varying degree of success. In this work, we predict the VLE of argon and assess the performance of several density functionals and the second order Møller-Plesset perturbation theory (MP2) by determining critical and structural properties via first principles Monte Carlo simulations. PBE-D3, BLYP-D3, and rVV10 functionals were used to compute vapor liquid coexistence curves, while PBE0-D3, M06-2X-D3, and MP2 were used for computing liquid density at a single state point. The performance of the PBE-D3 functional for VLE is superior to other functionals (BLYP-D3 and rVV10). At T = 85 K and P = 1 bar, MP2 performs well for the density and structural features of the first solvation shell in the liquid phase

Infrared spectra of protonated polycyclic aromatic hydrocarbon molecules: Azulene

The infrared (IR) spectrum of protonated azulene (AzuH(+), C10H9+) has been measured in the fingerprint range (600-1800 cm(-1)) by means of IR multiple photon dissociation (IRMPD) spectroscopy in a Fourier transform ion cyclotron resonance mass spectrometer equipped with an electrospray ionization source using a free electron laser. The potential energy surface of AzuH(+) has been characterized at the B3LYP/6-311G** level in order to determine the global and local minima and the corresponding transition states for interconversion. The energies of the local and global minima, the dissociation energies for the lowest-energy fragmentation pathways, and the proton affinity have been evaluated at the CBS-QB3 level. Comparison with calculated linear IR absorption spectra supports the assignment of the IRMPD spectrum to C4-protonated AzuH(+), the most stable of the six distinguishable C-protonated AzuH(+) isomers. Comparison between Azu and C4-AzuH(+) reveals the effects of protonation on the geometry, vibrational properties, and the charge distribution of these fundamental aromatic molecules. Calculations at the MP2 level indicate that this technique is not suitable to predict reliable IR spectra for this type of carbocations even for relatively large basis sets. The IRMPD spectrum of protonated azulene is compared to that of isomeric protonated naphthalene and to an astronomical spectrum of the unidentified IR emission bands

Improving the efficiency of G0W0G_0W_0 calculations with approximate spectral decompositions of dielectric matrices

Recently it was shown that the calculation of quasiparticle energies using the $G_0W_0$ approximation can be performed without computing explicitly any virtual electronic states, by expanding the Green function and screened Coulomb interaction in terms of the eigenstates of the static dielectric matrix. Avoiding the evaluation of virtual electronic states leads to improved efficiency and ease of convergence of $G_0W_0$ calculations. Here we propose a further improvement of the efficiency of these calculations, based on an approximation of density-density response functions of molecules and solids. The approximation relies on the calculation of a subset of eigenvectors of the dielectric matrix using the kinetic operator instead of the full Hamiltonian, and it does not lead to any substantial loss of accuracy for the quasiparticle energies. The computational savings introduced by this approximation depend on the system, and they become more substantial as the number of electrons increases