1
|
Kříž K, van Maaren PJ, van der Spoel D. Impact of Combination Rules, Level of Theory, and Potential Function on the Modeling of Gas- and Condensed-Phase Properties of Noble Gases. J Chem Theory Comput 2024; 20:2362-2376. [PMID: 38477573 DOI: 10.1021/acs.jctc.3c01257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2024]
Abstract
The systems of noble gases are particularly instructive for molecular modeling due to the elemental nature of their interactions. They do not normally form bonds nor possess a (permanent) dipole moment, and the only forces determining their bonding/clustering stems from van der Waals forces─dispersion and Pauli repulsion, which can be modeled by empirical potential functions. Combination rules, that is, formulas to derive parameters for pair potentials of heterodimers from parameters of corresponding homodimers, have been studied at length for the Lennard-Jones 12-6 potentials but not in great detail for other, more accurate, potentials. In this work, we examine the usefulness of nine empirical potentials in their ability to reproduce quantum mechanical (QM) benchmark dissociation curves of noble gas dimers (He, Ne, Ar, Kr, and Xe homo- and heterodimers), and we systematically study the efficacy of different permutations of combination relations for each parameter of the potentials. Our QM benchmark comprises dissociation curves computed by several different coupled cluster implementations as well as symmetry-adapted perturbation theory. The two-parameter Lennard-Jones potentials were decisively outperformed by more elaborate potentials that sport a 25-30 times lower root-mean-square error (RMSE) when fitted to QM dissociation curves. Very good fits to the QM dissociation curves can be achieved with relatively inexpensive four- or even three-parameter potentials, for instance, the damped 14-7 potential (Halgren, J. Am. Chem. Soc. 1992, 114, 7827-7843), a four-parameter Buckingham potential (Werhahn et al., Chem. Phys. Lett. 2015, 619, 133-138), or the three-parameter Morse potential (Morse, Phys. Rev. 1929, 34, 57-64). Potentials for heterodimers that are generated from combination rules have an RMSE that is up to 20 times higher than potentials that are directly fitted to the QM dissociation curves. This means that the RMSE, in particular, for light atoms, is comparable in magnitude to the well-depth of the potential. Based on a systematic permutation of combination rules, we present one or more combination rules for each potential tested that yield a relatively low RMSE. Two new combination rules are introduced that perform well, one for the van der Waals radius σij as ( 1 2 ( σ i 3 + σ j 3 ) ) 1 / 3 and one for the well-depth ϵij as ( 1 2 ( ϵ i - 2 + ϵ j - 2 ) ) - 1 / 2 . The QM data and the fitted potentials were evaluated in the gas phase against experimental second virial coefficients for homo- and heterodimers, the latter of which allowed evaluation of the combination rules. The fitted models were used to perform condensed phase molecular dynamics simulations to verify the melting points, liquid densities at the melting point, and the enthalpies of vaporization produced by the models for pure substances. Subtle differences in the benchmark potentials, in particular, the well-depth, due to the level of theory used were found here to have a profound effect on the macroscopic properties of noble gases: second virial coefficients or the bulk properties in simulations. By explicitly including three-body dispersion in molecular simulations employing the best pair potential, we were able to obtain accurate melting points as well as satisfactory densities and enthalpies of vaporization.
Collapse
Affiliation(s)
- Kristian Kříž
- Department of Cell and Molecular Biology, Uppsala University, Box 596, Uppsala SE-75124, Sweden
| | - Paul J van Maaren
- Department of Cell and Molecular Biology, Uppsala University, Box 596, Uppsala SE-75124, Sweden
| | - David van der Spoel
- Department of Cell and Molecular Biology, Uppsala University, Box 596, Uppsala SE-75124, Sweden
| |
Collapse
|
2
|
Bop CT, Quintas-Sánchez E, Sur S, Robin M, Lique F, Dawes R. Inelastic scattering in isotopologues of O 2-Ar: the effects of mass, symmetry, and density of states. Phys Chem Chem Phys 2021; 23:5945-5955. [PMID: 33666616 DOI: 10.1039/d1cp00326g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The two species considered here, O2 (oxygen molecule) and Ar (argon-atom), are both abundant components of Earth's atmosphere and hence familiar collision partners in this medium. O2 is quite reactive and extensively involved in atmospheric chemistry, including Chapman's cycle of the formation and destruction of ozone; while Ar, like N2, typically plays the nevertheless crucial role of inert collider. Inert species can provide stabilization to metastable encounter-complexes through the energy transfer associated with inelastic collisions. The interplay of collision frequency and energy transfer efficiency, with state lifetimes and species concentrations, contributes to the rich and varied chemistry and dynamics found in diverse environments ranging from planetary atmospheres to the interstellar and circumstellar media. The nature and density of bound and resonance states, coupled electronic states, symmetry, and nuclear spin-statistics can all play a role. Here, we systematically investigate some of those factors by looking at the O2-Ar system, comparing rigorous quantum-scattering calculations for the 16O16O-40Ar, 18O16O-40Ar, and 18O18O-40Ar isotope combinations. A new accurate potential energy surface was constructed for this purpose holding the O2 bond distance at its vibrationally averaged distance.
Collapse
Affiliation(s)
- Cheikh T Bop
- Laboratoire Ondes et Milieux Complexes, UMR 6294, Centre National de la Recherche Scientifique-Universite du Havre, F-76063 Le Havre, France. and Universite du Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | | | - Sangeeta Sur
- Missouri University of Science and Technology, Rolla, MO 65409-0010, USA.
| | - Mathurin Robin
- Laboratoire Ondes et Milieux Complexes, UMR 6294, Centre National de la Recherche Scientifique-Universite du Havre, F-76063 Le Havre, France.
| | - François Lique
- Laboratoire Ondes et Milieux Complexes, UMR 6294, Centre National de la Recherche Scientifique-Universite du Havre, F-76063 Le Havre, France. and Universite du Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | - Richard Dawes
- Missouri University of Science and Technology, Rolla, MO 65409-0010, USA.
| |
Collapse
|
3
|
Walz MM, van der Spoel D. Systematically improved melting point prediction: a detailed physical simulation model is required. Chem Commun (Camb) 2019; 55:12044-12047. [DOI: 10.1039/c9cc06177k] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Our detailed physical simulation model allows for an unprecedented and systematically improved prediction of melting points of alkali halides.
Collapse
Affiliation(s)
- Marie-Madeleine Walz
- Department of Cell and Molecular Biology
- Uppsala University
- SE-75124 Uppsala
- Sweden
| | - David van der Spoel
- Department of Cell and Molecular Biology
- Uppsala University
- SE-75124 Uppsala
- Sweden
| |
Collapse
|
4
|
Walz MM, Ghahremanpour MM, van Maaren PJ, van der Spoel D. Phase-Transferable Force Field for Alkali Halides. J Chem Theory Comput 2018; 14:5933-5948. [DOI: 10.1021/acs.jctc.8b00507] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Marie-Madeleine Walz
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Mohammad M. Ghahremanpour
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - Paul J. van Maaren
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| | - David van der Spoel
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, SE-75124 Uppsala, Sweden
| |
Collapse
|
5
|
Wei L, Li P, Tang K. Iterative combining rules for the van der Waals potentials of mixed rare gas systems. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.02.066] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
6
|
Abstract
Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.
Collapse
Affiliation(s)
| | - Kenneth M. Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, and Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
7
|
Ivlev DV, Kiselev MG. Helium diffusion in aqueous sodium chloride solution at high pressures. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2012. [DOI: 10.1134/s003602441206012x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
8
|
|
9
|
Peng Z, Ewig CS, Hwang MJ, Waldman M, Hagler AT. Derivation of Class II Force Fields. 4. van der Waals Parameters of Alkali Metal Cations and Halide Anions. J Phys Chem A 1997. [DOI: 10.1021/jp964080y] [Citation(s) in RCA: 152] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Zhengwei Peng
- Molecular Simulations Inc., 9685 Scranton Road, San Diego, California 92121
| | - Carl S. Ewig
- Molecular Simulations Inc., 9685 Scranton Road, San Diego, California 92121
| | - Ming-Jing Hwang
- Molecular Simulations Inc., 9685 Scranton Road, San Diego, California 92121
| | - Marvin Waldman
- Molecular Simulations Inc., 9685 Scranton Road, San Diego, California 92121
| | - Arnold T. Hagler
- Molecular Simulations Inc., 9685 Scranton Road, San Diego, California 92121
| |
Collapse
|
10
|
|
11
|
Mirsky K. Carbon monoxide molecules in an argon matrix: empirical evaluation of the Ar·Ar, C·Ar and O·Ar potential parameters. Chem Phys 1980. [DOI: 10.1016/0301-0104(80)85220-7] [Citation(s) in RCA: 67] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
12
|
|
13
|
Shashkov A, Abramenko T, Nesterov N, Joshi R, Afshar R, Saxena S. Thermal conductivity of argon, krypton and their mixtures at low temperatures (90-270 k). Chem Phys 1978. [DOI: 10.1016/0301-0104(78)85088-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
14
|
Nain VPS, Aziz RA, Jain PC, Saxena SC. Interatomic potentials and transport properties for neon, argon, and krypton. J Chem Phys 1976. [DOI: 10.1063/1.433497] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
|
16
|
Taylor WL, Weissman S. Experimental thermal diffusion factors for20Ne–22Ne and their application as a test of the neon interatomic potential. J Chem Phys 1974. [DOI: 10.1063/1.1681589] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
17
|
Kong CL. Combining rules for intermolecular potential parameters. II. Rules for the Lennard‐Jones (12–6) potential and the Morse potential. J Chem Phys 1973. [DOI: 10.1063/1.1680358] [Citation(s) in RCA: 201] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
Kong CL. Atomic distortion and the repulsive interactions of the noble gas atoms. J Chem Phys 1973. [DOI: 10.1063/1.1680124] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
19
|
Kistemaker P, Hanna M, Tom A, De Vries A. Rotational relaxation in mixtures of methane with helium, argon and xenon. ACTA ACUST UNITED AC 1972. [DOI: 10.1016/0031-8914(72)90115-2] [Citation(s) in RCA: 20] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
20
|
Sitarski M, Popielawski J. The kinetic theory of diffusion and thermal diffusion in adsorbed gases. Mol Phys 1972. [DOI: 10.1080/00268977200100361] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
21
|
|