1
|
von Domaros M, Perlt E. Anharmonic effects in the quantum cluster equilibrium method. J Chem Phys 2017; 146:124114. [PMID: 28388115 DOI: 10.1063/1.4978958] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Michael von Domaros
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 6, D-53115 Bonn, Germany
| | - Eva Perlt
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 6, D-53115 Bonn, Germany
| |
Collapse
|
2
|
English NJ, Waldron CJ. Perspectives on external electric fields in molecular simulation: progress, prospects and challenges. Phys Chem Chem Phys 2016; 17:12407-40. [PMID: 25903011 DOI: 10.1039/c5cp00629e] [Citation(s) in RCA: 143] [Impact Index Per Article: 17.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
In this review, the application of a wide variety of external electric fields in molecular simulation shall be discussed, including time-varying and electromagnetic, as well as the utility and potential impact and prospects for exploitation of such simulations for real-world and industrial end use. In particular, non-equilibrium molecular dynamics will be discussed, as well as challenges in addressing adequate thermostatting and scaling field amplitudes to more experimentally relevant levels. Attention shall be devoted to recent progress and advances in external fields in ab initio molecular simulation and dynamics, as well as elusive challenges thereof (and, to some extent, for molecular dynamics from empirical potentials), such as timescales required to observe low-frequency and intensity field effects. The challenge of deterministic molecular dynamics in external fields in sampling phase space shall be discussed, along with prospects for application of fields in enhanced-sampling simulations. Finally, the application of external electric fields to a wide variety of aqueous, nanoscale and biological systems will be discussed, often motivated by the possibility of exploitation in real-world applications, which serve to underpin our molecular-level understanding of field effects in terms of microscopic mechanisms, and possibly with a view to control thereof.
Collapse
Affiliation(s)
- Niall J English
- School of Chemical and Bioprocess Engineering, University College Dublin, Belfield, Dublin 4, Ireland.
| | | |
Collapse
|
3
|
Riniker S, Kunz APE, van Gunsteren WF. On the Calculation of the Dielectric Permittivity and Relaxation of Molecular Models in the Liquid Phase. J Chem Theory Comput 2011; 7:1469-75. [DOI: 10.1021/ct100610v] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sereina Riniker
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093 Zürich, Switzerland
| | - Anna-Pitschna E. Kunz
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093 Zürich, Switzerland
| | - Wilfred F. van Gunsteren
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093 Zürich, Switzerland
| |
Collapse
|
4
|
Kunz APE, Eichenberger AP, van Gunsteren WF. A simple, efficient polarizable molecular model for liquid carbon tetrachloride. Mol Phys 2011. [DOI: 10.1080/00268976.2010.533208] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
5
|
Determination of the Structure of Fluid Deuterium Chloride by Neutron Diffraction at High Pressure. ACTA ACUST UNITED AC 2010. [DOI: 10.1002/bbpc.19910950507] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
6
|
Horn HW, Swope WC, Pitera JW, Madura JD, Dick TJ, Hura GL, Head-Gordon T. Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. J Chem Phys 2007; 120:9665-78. [PMID: 15267980 DOI: 10.1063/1.1683075] [Citation(s) in RCA: 1410] [Impact Index Per Article: 82.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A re-parameterization of the standard TIP4P water model for use with Ewald techniques is introduced, providing an overall global improvement in water properties relative to several popular nonpolarizable and polarizable water potentials. Using high precision simulations, and careful application of standard analytical corrections, we show that the new TIP4P-Ew potential has a density maximum at approximately 1 degrees C, and reproduces experimental bulk-densities and the enthalpy of vaporization, DeltaH(vap), from -37.5 to 127 degrees C at 1 atm with an absolute average error of less than 1%. Structural properties are in very good agreement with x-ray scattering intensities at temperatures between 0 and 77 degrees C and dynamical properties such as self-diffusion coefficient are in excellent agreement with experiment. The parameterization approach used can be easily generalized to rehabilitate any water force field using available experimental data over a range of thermodynamic points.
Collapse
Affiliation(s)
- Hans W Horn
- IBM Almaden Research Center, San Jose, California 95120, USA.
| | | | | | | | | | | | | |
Collapse
|
7
|
Lomba E, Lombardero M, Abascal J. New aspects in the simulation and behaviour of polar molecular fluids. Mol Phys 2007. [DOI: 10.1080/00268978900102741] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- E. Lomba
- a Instituto de Química Física Rocasolano, C.S.I.C. , Serrano 119, E-28006 , Madrid , Spain
- b Depto. Química Física, Facultad de CC. Químicas , Universidad Complutense de Madrid , E-28040 , Madrid , Spain
| | - M. Lombardero
- a Instituto de Química Física Rocasolano, C.S.I.C. , Serrano 119, E-28006 , Madrid , Spain
- b Depto. Química Física, Facultad de CC. Químicas , Universidad Complutense de Madrid , E-28040 , Madrid , Spain
| | - J.L.F. Abascal
- c Depto. Química Física, Facultad de CC. Químicas , Universidad Complutense de Madrid , E-28040 , Madrid , Spain
| |
Collapse
|
8
|
|
9
|
Dubois V, Pasquarello A. Ab initio molecular dynamics of liquid hydrogen chloride. J Chem Phys 2006; 122:114512. [PMID: 15836234 DOI: 10.1063/1.1869972] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We carried out an ab initio molecular dynamics simulation of liquid hydrogen chloride (l-HCl) at a temperature of 313 K. Comparison with inelastic neutron scattering data shows that the simulation achieves an overall good description of the structural correlations, improving significantly upon a description based on classical interaction potentials. Despite some minor differences between theory and experiment in the H-H partial structure factor, the simulation gives a description of the hydrogen bonding in impressive agreement with experiment, for both the amount and the bond-length distribution of the bonds. In the simulation, 40% of the molecules are nonbonded, while the hydrogen-bonded chains are short, principally consisting of dimers (25%) and trimers (15%). Neighboring molecules in the simulation are found to form L-shaped arrangements, like in the isolated (HCl)(2) dimer and in crystalline phases of HCl. The time correlation of the molecular-axis orientation is found to be characterized by a very short decay time (0.13 ps), consistent with the short length of the hydrogen-bonded chains. Other dynamical properties investigated in this work include the diffusion coefficient and the vibrational density of states. We evaluated the molecular dipole of the HCl molecule in the liquid using a definition based on the coupling of rotational modes to an external electric field. The average dipole moment (1.53 D) derived in this way is found to be considerably larger than for the isolated molecule (1.11 D). Our results show that the dipole moment in [script-l]-HCl undergoes large fluctuations, both in orientation and in modulus. Upon the onset of an external field, such dipole fluctuations concur to reduce the fluctuations of the dielectric response.
Collapse
Affiliation(s)
- Vincent Dubois
- Institut de Théorie des Phénomènes Physiques (ITP), Ecole Polytechnique Fédérale de Lausanne (EPFL), Switzerland
| | | |
Collapse
|
10
|
Dubois V, Umari P, Pasquarello A. Dielectric susceptibility of dipolar molecular liquids by ab initio molecular dynamics: application to liquid HCl. Chem Phys Lett 2004. [DOI: 10.1016/j.cplett.2004.04.021] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
11
|
|
12
|
Martín C, Lomba E, Lombardero M, Lado F, Ho/ye JS. Integral equation and simulation studies of a realistic model for liquid hydrogen chloride. J Chem Phys 1994. [DOI: 10.1063/1.466586] [Citation(s) in RCA: 20] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
13
|
Zhou H, Bagchi B, Papazyan A, Maroncelli M. Solvation dynamics in a Brownian dipole lattice: A comparison between theory and computer simulation. J Chem Phys 1992. [DOI: 10.1063/1.463307] [Citation(s) in RCA: 41] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Zhou H, Bagchi B. Dielectric and orientational relaxation in a Brownian dipolar lattice. J Chem Phys 1992. [DOI: 10.1063/1.462944] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Adams D. Hard dumbbell fluid g1and g2from Monte Carlo. Mol Phys 1991. [DOI: 10.1080/00268979100101521] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
16
|
Levesque D, Weis J, Oxtoby DW. A molecular dynamics simulation of rotational and vibrational relaxation in liquid HCl. J Chem Phys 1983. [DOI: 10.1063/1.445868] [Citation(s) in RCA: 60] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|