1
|
Peeters J, Vanommeslaeghe K. A Simple Model for the Pauli Repulsion with Possible Utility in QM, MM and Chemical Education. J Chem Theory Comput 2024. [PMID: 39038213 DOI: 10.1021/acs.jctc.4c00748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/24/2024]
Abstract
The Pauli repulsion is the intermolecular force responsible for the volume and low compressibility of condensed-phase matter at normal conditions. A simple model for this force is presented, wherein per-atom electron densities are represented as spherical charge distributions that are prevented from significantly overlapping. In the example of two noble gas atoms approaching one another beyond their van der Waals radii, the distance between the centers of the electronic charge distributions becomes larger than the distance between the nuclei, giving rise to an unfavorable electrostatic interaction. For the purpose of calculating this interaction, the model is further simplified by representing the per-atom electron density as a negative point charge, loosely inspired by the classical Drude oscillator. The dispersion interaction is simplified to an R-6 term, centered on the aforementioned point charges. Despite the gross simplicity of the resulting formalism, near-quantitative agreement with high-level QM interaction energies of noble gas dimers is achieved. Accordingly, the present model is thought to have utility in force fields, in post-HF and post-DFT dispersion corrections, and in chemical education.
Collapse
Affiliation(s)
- Jordy Peeters
- Department of Analytical Chemistry, Applied Chemometrics and Molecular Modelling, Vrije Universiteit Brussel (VUB), Laarbeeklaan 103, 1090 Brussels, Belgium
| | - Kenno Vanommeslaeghe
- Department of Analytical Chemistry, Applied Chemometrics and Molecular Modelling, Vrije Universiteit Brussel (VUB), Laarbeeklaan 103, 1090 Brussels, Belgium
| |
Collapse
|
2
|
Schwerdtfeger P, Wales DJ. 100 Years of the Lennard-Jones Potential. J Chem Theory Comput 2024; 20:3379-3405. [PMID: 38669689 DOI: 10.1021/acs.jctc.4c00135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2024]
Abstract
It is now 100 years since Lennard-Jones published his first paper introducing the now famous potential that bears his name. It is therefore timely to reflect on the many achievements, as well as the limitations, of this potential in the theory of atomic and molecular interactions, where applications range from descriptions of intermolecular forces to molecules, clusters, and condensed matter.
Collapse
Affiliation(s)
- Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, Private Bag 102904, Auckland 0745, New Zealand
| | - David J Wales
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
3
|
Khabibrakhmanov A, Fedorov DV, Tkatchenko A. Universal Pairwise Interatomic van der Waals Potentials Based on Quantum Drude Oscillators. J Chem Theory Comput 2023; 19:7895-7907. [PMID: 37875419 PMCID: PMC10653113 DOI: 10.1021/acs.jctc.3c00797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Revised: 09/30/2023] [Accepted: 10/05/2023] [Indexed: 10/26/2023]
Abstract
Repulsive short-range and attractive long-range van der Waals (vdW) forces play an appreciable role in the behavior of extended molecular systems. When using empirical force fields, the most popular computational methods applied to such systems, vdW forces are typically described by Lennard-Jones-like potentials, which unfortunately have a limited predictive power. Here, we present a universal parameterization of a quantum-mechanical vdW potential, which requires only two free-atom properties─the static dipole polarizability α1 and the dipole-dipole C6 dispersion coefficient. This is achieved by deriving the functional form of the potential from the quantum Drude oscillator (QDO) model, employing scaling laws for the equilibrium distance and the binding energy, and applying the microscopic law of corresponding states. The vdW-QDO potential is shown to be accurate for vdW binding energy curves, as demonstrated by comparing to the ab initio binding curves of 21 noble-gas dimers. The functional form of the vdW-QDO potential has the correct asymptotic behavior at both zero and infinite distances. In addition, it is shown that the damped vdW-QDO potential can accurately describe vdW interactions in dimers consisting of group II elements. Finally, we demonstrate the applicability of the atom-in-molecule vdW-QDO model for predicting accurate dispersion energies for molecular systems. The present work makes an important step toward constructing universal vdW potentials, which could benefit (bio)molecular computational studies.
Collapse
Affiliation(s)
- Almaz Khabibrakhmanov
- Department of Physics and Materials
Science, University of Luxembourg, L-1511 Luxembourg
City, Luxembourg
| | - Dmitry V. Fedorov
- Department of Physics and Materials
Science, University of Luxembourg, L-1511 Luxembourg
City, Luxembourg
| | - Alexandre Tkatchenko
- Department of Physics and Materials
Science, University of Luxembourg, L-1511 Luxembourg
City, Luxembourg
| |
Collapse
|
4
|
Lang J, Przybytek M, Lesiuk M, Jeziorski B. Collision-induced three-body polarizability of helium. J Chem Phys 2023; 158:114303. [PMID: 36948830 DOI: 10.1063/5.0137879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2023] Open
Abstract
We present the first-principles determination of the three-body polarizability and the third dielectric virial coefficient of helium. Coupled-cluster and full configuration interaction methods were used to perform electronic structure calculations. The mean absolute relative uncertainty of the trace of the polarizability tensor, resulting from the incompleteness of the orbital basis set, was found to be 4.7%. Additional uncertainty due to the approximate treatment of triple and the neglect of higher excitations was estimated at 5.7%. An analytic function was developed to describe the short-range behavior of the polarizability and its asymptotics in all fragmentation channels. We calculated the third dielectric virial coefficient and its uncertainty using the classical and semiclassical Feynman-Hibbs approaches. The results of our calculations were compared with experimental data and with recent Path-Integral Monte Carlo (PIMC) calculations [Garberoglio et al., J. Chem. Phys. 155, 234103 (2021)] employing the so-called superposition approximation of the three-body polarizability. For temperatures above 200 K, we observed a significant discrepancy between the classical results obtained using superposition approximation and the ab initio computed polarizability. For temperatures from 10 K up to 200 K, the differences between PIMC and semiclassical calculations are several times smaller than the uncertainties of our results. Except at low temperatures, our results agree very well with the available experimental data but have much smaller uncertainties. The data reported in this work eliminate the main accuracy bottleneck in the optical pressure standard [Gaiser et al., Ann. Phys. 534, 2200336 (2022)] and facilitate further progress in the field of quantum metrology.
Collapse
Affiliation(s)
- J Lang
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - M Przybytek
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - M Lesiuk
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - B Jeziorski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
5
|
Kayang KW, Volkov AN, Zhilyaev PA, Sharipov F. The ab initio potential energy curves of atom pairs and transport properties of high-temperature vapors of Cu and Si and their mixtures with He, Ar, and Xe gases. Phys Chem Chem Phys 2023; 25:4872-4898. [PMID: 36692492 DOI: 10.1039/d2cp04981c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The potential energy curves (PECs) for the homonuclear He-He, Ar-Ar, Cu-Cu, and Si-Si dimers, as well as heteronuclear Cu-He, Cu-Ar, Cu-Xe, Si-He, Si-Ar, and Si-Xe dimers, are obtained in quantum Monte Carlo (QMC) calculations. It is shown that the QMC method provides the PECs with an accuracy comparable with that of the state-of-the-art coupled cluster singles and doubles with perturbative triples corrections [CCSD(T)] calculations. The QMC data are approximated by the Morse long range (MLR) and (12-6) Lennard-Jones (LJ) potentials. The MLR and LJ potentials are used to calculate the deflection angles in binary collisions of corresponding atom pairs and transport coefficients of Cu and Si vapors and their mixtures with He, Ar, and Xe gases in the range of temperature from 100 K to 10 000 K. It is shown that the use of the LJ potentials introduces significant errors in the transport coefficients of high-temperature vapors and gas mixtures. The mixtures with heavy noble gases demonstrate anomalous behavior when the viscosity and thermal conductivity can be larger than that of the corresponding pure substances. In the mixtures with helium, the thermal diffusion factor is found to be unusually large. The calculated viscosity and diffusivity are used to determine parameters of the variable hard sphere and variable soft sphere molecular models as well as parameters of the power-law approximations for the transport coefficients. The results obtained in the present work include all information required for kinetic or continuum simulations of dilute Cu and Si vapors and their mixtures with He, Ar, and Xe gases.
Collapse
Affiliation(s)
- Kevin W Kayang
- Department of Mechanical Engineering, University of Alabama, 7th Avenue, Tuscaloosa, AL 35487, USA.
| | - Alexey N Volkov
- Department of Mechanical Engineering, University of Alabama, 7th Avenue, Tuscaloosa, AL 35487, USA.
| | - Petr A Zhilyaev
- Skolkovo Institute of Science and Technology, 121205, Bolshoy Boulevard 30, bld. 1, Moscow, Russia
| | - Felix Sharipov
- Departamento de Física, Universidade Federal do Paraná, Caixa Postal 19044, Curitiba 81531-980, Brazil
| |
Collapse
|
6
|
Shanks BL, Potoff JJ, Hoepfner MP. Transferable Force Fields from Experimental Scattering Data with Machine Learning Assisted Structure Refinement. J Phys Chem Lett 2022; 13:11512-11520. [PMID: 36469859 DOI: 10.1021/acs.jpclett.2c03163] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Deriving transferable pair potentials from experimental neutron and X-ray scattering measurements has been a longstanding challenge in condensed matter physics. State-of-the-art scattering analysis techniques estimate real-space microstructure from reciprocal-space total scattering data by refining pair potentials to obtain agreement between simulated and experimental results. Prior attempts to apply these potentials with molecular simulations have revealed inaccurate predictions of thermodynamic fluid properties. In this Letter, a machine learning assisted structure-inversion method applied to neutron scattering patterns of the noble gases (Ne, Ar, Kr, and Xe) is shown to recover transferable pair potentials that accurately reproduce both microstructure and vapor-liquid equilibria from the triple to critical point. Therefore, it is concluded that a single neutron scattering measurement is sufficient to predict macroscopic thermodynamic properties over a wide range of states and provide novel insight into local atomic forces in dense monatomic systems.
Collapse
Affiliation(s)
- Brennon L Shanks
- Department of Chemical Engineering, University of Utah, Salt Lake City, UT84112-9202, United States
| | - Jeffrey J Potoff
- Department of Chemical Engineering and Materials Science, Wayne State University, Detroit, MI48202, United States
| | - Michael P Hoepfner
- Department of Chemical Engineering, University of Utah, Salt Lake City, UT84112-9202, United States
| |
Collapse
|
7
|
Graham RS, Wheatley RJ. Machine learning for non-additive intermolecular potentials: quantum chemistry to first-principles predictions. Chem Commun (Camb) 2022; 58:6898-6901. [PMID: 35642644 DOI: 10.1039/d2cc01820a] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Prediction of thermophysical properties from molecular principles requires accurate potential energy surfaces (PES). We present a widely-applicable method to produce first-principles PES from quantum chemistry calculations. Our approach accurately interpolates three-body non-additive interaction data, using the machine learning technique, Gaussian Processes (GP). The GP approach needs no bespoke modification when the number or type of molecules is changed. Our method produces highly accurate interpolation from significantly fewer training points than typical approaches, meaning ab initio calculations can be performed at higher accuracy. As an exemplar we compute the PES for all three-body cross interactions for CO2-Ar mixtures. From these we calculate the CO2-Ar virial coefficients up to 5th order. The resulting virial equation of state (EoS) is convergent for densities up to the critical density. Where convergent, the EoS makes accurate first-principles predictions for a range of thermophysical properties for CO2-Ar mixtures, including the compressibility factor, speed of sound and Joule-Thomson coefficient. Our method has great potential to make wide-ranging first-principles predictions for mixtures of comparably sized molecules. Such predictions can replace the need for expensive, laborious and repetitive experiments and inform the continuum models required for applications.
Collapse
Affiliation(s)
- Richard S Graham
- School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, UK.
| | | |
Collapse
|
8
|
Mizuse K, Sato U, Tobata Y, Ohshima Y. Rotational spectroscopy of the argon dimer by time-resolved Coulomb explosion imaging of rotational wave packets. Phys Chem Chem Phys 2022; 24:11014-11022. [PMID: 35470358 DOI: 10.1039/d2cp01113a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report time-domain rotational spectroscopy of the argon dimer, Ar2, by implementing time-resolved Coulomb explosion imaging of rotational wave packets. The rotational wave packets are created in Ar2 with a linearly polarized, nonresonant, ultrashort laser pulse, and their spatiotemporal evolution is fully characterized by measuring angular distribution of the fragmented Ar+ promptly ejected from Ar22+ generated by the more intense probe pulse. The pump-probe measurements have been carried out up to a delay time of 16 ns. The alignment parameters, derived from the observed images, exhibit periodic oscillation lasting for more than 15 ns. The pure rotational spectrum of Ar2 is obtained by Fourier transformation of the time traces of the alignment parameters. The frequency resolution in the spectrum is about 90 MHz, the highest ever achieved for Ar2. The rotational constant and the centrifugal distortion constant are determined with much improved precision than the previous experimental results: B0 = 1.72713 ± 0.00009 GHz and D0 = 0.0310 ± 0.0005 MHz. The present B0 value does not match within the quoted experimental uncertainty with that from the VUV spectroscopy, so far accepted as an experimental reference to assess theories. The present improved constants would stand as new references to calibrate state-of-the-art theoretical investigations and an indispensable experimental source for the construction of an accurate empirical intermolecular potential.
Collapse
Affiliation(s)
- Kenta Mizuse
- Department of Chemistry, School of Science, Tokyo Institute of Technology, 2-12-1-W4-9 Ookayama, Meguro, Tokyo 152-8550, Japan. .,Department of Chemistry, School of Science, Kitasato University, 1-15-1 Kitazato, Minami, Sagamihara, Kanagawa 252-0373, Japan.
| | - Urara Sato
- Department of Chemistry, School of Science, Kitasato University, 1-15-1 Kitazato, Minami, Sagamihara, Kanagawa 252-0373, Japan.
| | - Yuya Tobata
- Department of Chemistry, School of Science, Tokyo Institute of Technology, 2-12-1-W4-9 Ookayama, Meguro, Tokyo 152-8550, Japan.
| | - Yasuhiro Ohshima
- Department of Chemistry, School of Science, Tokyo Institute of Technology, 2-12-1-W4-9 Ookayama, Meguro, Tokyo 152-8550, Japan.
| |
Collapse
|
9
|
Deiters UK, Sadus RJ. First-principles determination of the solid-liquid-vapor triple point: The noble gases. Phys Rev E 2022; 105:054128. [PMID: 35706194 DOI: 10.1103/physreve.105.054128] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Accepted: 04/21/2022] [Indexed: 06/15/2023]
Abstract
We report first-principles calculations of the triple point that allow us to predict the triple point temperature of atomic fluids to an accuracy that has not been previously possible. This is achieved by proposing a molecular simulation technique that can be used for solid-liquid equilibria at arbitrarily low pressures. It is demonstrated that the triple point is significantly influenced by the choice of two-body, three-body and quantum interactions. An improved theoretical understanding of triple points is important for both science in general, and metrology in particular, as it links the Boltzmann constant and the Kelvin temperature scale to fundamental constants.
Collapse
Affiliation(s)
- Ulrich K Deiters
- Institute for Physical Chemistry, Faculty of Mathematics and Natural Sciences, University of Cologne, Greinstrasse 4-6, D-50939 Köln, Germany
| | - Richard J Sadus
- Department of Computer Science and Software Engineering, Swinburne University of Technology, Wurundjeri Country, PO Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
10
|
Alboul L, Lishchuk SV. Bulk viscosity of gaseous argon from molecular dynamics simulations. Phys Rev E 2022; 105:054135. [PMID: 35706273 DOI: 10.1103/physreve.105.054135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Accepted: 05/06/2022] [Indexed: 06/15/2023]
Abstract
The bulk viscosity of dilute argon gas is calculated using molecular dynamics simulations in the temperature range 150-500 K and is found to be proportional to density squared in the investigated range of densities 0.001-1 kg m^{-3}. A comparison of the results obtained using Lennard-Jones and Tang-Toennies models of pair interaction potential reveals that the value of the bulk viscosity coefficient is sensitive to the choice of the pair interaction model. The inclusion of the Axilrod-Teller-Muto three-body interaction in the model does not noticeably affect the values of the bulk viscosity in dilute states, contrary to the previously investigated case of dense fluids.
Collapse
Affiliation(s)
- Lyuba Alboul
- Industry & Innovation Research Institute, Sheffield Hallam University, Sheffield S1 1WB, United Kingdom
| | - Sergey V Lishchuk
- Thermodynamics and Process Engineering, Technische Universität Berlin, 10587 Berlin, Germany
| |
Collapse
|
11
|
Bader F, Tremblay JC, Paulus B. Theoretical modeling of molecules in weakly interacting environments: trifluoride anions in argon. Phys Chem Chem Phys 2022; 24:3555-3567. [PMID: 35080559 DOI: 10.1039/d1cp02338a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The properties of molecules can be affected by the presence of a host environment. Even in inert rare gas matrices such effects are observable, as for instance in matrix isolation spectroscopy. In this work we study the trifluoride anion in cryogenic argon environments. To investigate the structure and vibrational properties of the guest-host systems, a potential energy surface of compound F-3-argon structures is determined from ab initio calculations with the CCSD(T)-F12b approach. Argon environments are probed with minima hopping optimizations of extended trifluoride-argon clusters. The vibrations of F-3 within the optimized environments are examined with anharmonic vibrational analyses. Among the three identified structural surroundings for the trifluoride, two are characterized by relatively favorable guest-host and host-host interactions as well as vibrational zero-point energies. A striking dependence of the trifluoride properties on the particular argon environment reveals the delicate influence of the host atoms on the guest molecule. Very good agreement with measured data suggests that in experiment F-3 occupies a double-vacancy site.
Collapse
Affiliation(s)
- Frederik Bader
- Institut für Chemie und Biochemie, Freie Universität Berlin, D-14195, Berlin, Germany.
| | - Jean Christophe Tremblay
- Laboratoire de Physique et Chimie Théoriques, CNRS-Université de Lorraine, UMR 7019, ICPM, 1Bd Arago, 57070 Metz, France
| | - Beate Paulus
- Institut für Chemie und Biochemie, Freie Universität Berlin, D-14195, Berlin, Germany.
| |
Collapse
|
12
|
Simulations of internal energy and pρT properties of monatomic fluids from accurate ab initio potentials and their uncertainty analysis. J Supercrit Fluids 2022. [DOI: 10.1016/j.supflu.2021.105425] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
13
|
Garberoglio G, Harvey AH, Jeziorski B. Path-integral calculation of the third dielectric virial coefficient of noble gases. J Chem Phys 2021; 155:234103. [PMID: 34937356 DOI: 10.1063/5.0077684] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a rigorous framework for fully quantum calculation of the third dielectric virial coefficient Cɛ(T) of noble gases, including exchange effects. The quantum effects are taken into account with the path-integral Monte Carlo method. Calculations employing state-of-the-art pair and three-body potentials and pair polarizabilities yield results generally consistent with the few scattered experimental data available for helium, neon, and argon, but rigorous calculations with well-described uncertainties will require the development of surfaces for the three-body nonadditive polarizability and the three-body dipole moment. The framework, developed here for the first time, will enable new approaches to primary temperature and pressure metrology based on first-principles calculations of gas properties.
Collapse
Affiliation(s)
- Giovanni Garberoglio
- European Centre for Theoretical Studies in Nuclear Physics and Related Areas (FBK-ECT*), Trento I-38123, Italy
| | - Allan H Harvey
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| | - Bogumił Jeziorski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
14
|
Rock CA, Arradondo SN, Tschumper GS. Solvation of Isoelectronic Halide and Alkali Metal Ions by Argon Atoms. J Phys Chem A 2021; 125:10524-10531. [PMID: 34851634 DOI: 10.1021/acs.jpca.1c08069] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
This work systematically examines the interactions of alkali metal cations and their isoelectronic halide counterparts with up to six solvating Ar atoms (M+Arn and X-Arn, where M = Li, Na, K, and Rb; X = H, F, Cl, and Br; and n = 1-6) via full geometry optimizations with the MP2 method and robust, correlation-consistent quadruple-ζ (QZ) basis sets. 116 unique M+Arn and X-Arn stationary points have been characterized on the MP2/QZ potential energy surface. To the best of our knowledge, approximately two dozen of these stationary points have been reported here for the first time. Some of these new structures are either the lowest-energy stationary point for a particular cluster or energetically competitive with it. The CCSD(T) method was employed to perform additional single-point energy computations upon all MP2/QZ-optimized structures using the same basis set. CCSD(T)/QZ results indicate that internally solvated structures with the ion at/near the geometric center of the cluster have appreciably higher energies than those placing the ion on the periphery. While this study extends the prior investigations of M+Arn clusters found within the literature, it notably provides one of the first thorough characterizations of and comparisons to the corresponding negatively charged X-Arn clusters.
Collapse
Affiliation(s)
- Carly A Rock
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Sarah N Arradondo
- Department of Chemistry, Washington College, Chestertown, Maryland 21620-1438, United States
| | - Gregory S Tschumper
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
15
|
Sharipov F, Benites VJ. Transport coefficients of isotopic mixtures of noble gases based on ab initio potentials. Phys Chem Chem Phys 2021; 23:16664-16674. [PMID: 34337630 DOI: 10.1039/d1cp01971f] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The transport coefficients such as viscosity, thermal conductivity, diffusion and thermal diffusion of neon, argon, krypton, and xenon are computed for a wide range of temperatures taking into consideration their real isotopic compositions. A new concept of isotopic thermal diffusion factor is introduced and calculated. The Chapman-Enskog method based on the 10th order approximation with respect to the Sonine polynomial expansion is applied. Ab initio potentials of interatomic interactions are employed to compute the transport cross-sections as they are part of the coefficient expressions. To study the influence of the isotopic composition, the same transport coefficients have been calculated for the single gases having an average atomic mass. The estimated numerical error of the present results is a function of the temperature and is different for each coefficient. At the room temperature, the relative numerical error of viscosity, thermal conductivity and diffusion coefficient is on the order of 10-6. The numerical error of the thermal diffusion factor affects the fifth decimal digit. The influence of the isotopic composition on viscosity and thermal conductivity depends on the gas species. It is negligible for argon and significant (about 0.02%) for xenon, while neon and krypton are weakly affected by the isotopic composition. The diffusion coefficient for each pair of isotopes differs from the corresponding self-diffusion coefficient by about 3%. The thermal diffusion factor of each isotope differs from the thermal self-diffusion factor in the third decimal digit.
Collapse
Affiliation(s)
- Felix Sharipov
- Departamento de Física, Universidade Federal do Paraná, Caixa Postal 19044, Curitiba, 81531-990, Brazil.
| | | |
Collapse
|
16
|
Hellmann R, Gaiser C, Fellmuth B, Vasyltsova T, Bich E. Thermophysical properties of low-density neon gas from highly accurate first-principles calculations and dielectric-constant gas thermometry measurements. J Chem Phys 2021; 154:164304. [PMID: 33940840 DOI: 10.1063/5.0047999] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
New interatomic potential energy and interaction-induced polarizability curves for two ground-state neon atoms were developed and used to predict the second density, acoustic, and dielectric virial coefficients and the dilute gas shear viscosity and thermal conductivity of neon at temperatures up to 5000 K. The potential energy curve is based on supermolecular coupled-cluster (CC) calculations at very high levels up to CC with single, double, triple, quadruple, and perturbative pentuple excitations [CCSDTQ(P)]. Scalar and spin-orbit relativistic effects, the diagonal Born-Oppenheimer correction, and retardation of the dispersion interactions were taken into account. The interaction-induced polarizability curve, which in this work is only needed for the calculation of the second dielectric virial coefficient, is based on supermolecular calculations at levels up to CCSDT and includes a correction for scalar relativistic effects. In addition to these first-principles calculations, highly accurate dielectric-constant gas thermometry (DCGT) datasets measured at temperatures from 24.5 to 200 K were analyzed to obtain the difference between the second density and dielectric virial coefficients with previously unattained accuracy. The agreement of the DCGT values with the ones resulting from the first-principles calculations is, despite some small systematic deviations, very satisfactory. Apart from this combination of two virial coefficients, the calculated thermophysical property values of this work are significantly more accurate than any available experimental data.
Collapse
Affiliation(s)
- Robert Hellmann
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Christof Gaiser
- Physikalisch-Technische Bundesanstalt (PTB), Abbestraße 2-12, 10587 Berlin, Germany
| | - Bernd Fellmuth
- Physikalisch-Technische Bundesanstalt (PTB), Abbestraße 2-12, 10587 Berlin, Germany
| | - Tatjana Vasyltsova
- Bereich Maschinenbau/Verfahrens- und Umwelttechnik, Hochschule Wismar, Philipp-Müller-Str. 14, 23966 Wismar, Germany
| | - Eckard Bich
- Institut für Chemie, Universität Rostock, Albert-Einstein-Str. 3a, 18059 Rostock, Germany
| |
Collapse
|
17
|
Singh AN, Dyre JC, Pedersen UR. Solid–liquid coexistence of neon, argon, krypton, and xenon studied by simulations. J Chem Phys 2021; 154:134501. [DOI: 10.1063/5.0045398] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Aditya N. Singh
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53703, USA
| | - Jeppe C. Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P. O. Box 260, DK-4000 Roskilde, Denmark
| | - Ulf R. Pedersen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P. O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
18
|
Sheng XW, Tang KT. The development of a full range analytical interatomic potential. Phys Chem Chem Phys 2021; 23:7748-7757. [PMID: 32959835 DOI: 10.1039/d0cp04083e] [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
A chronological account is given to the development of a full range interatomic potential. Starting with a simple phenomenological model, the terms in the model are gradually modified, so that they can carry some definite physical meaning. To gain insight, a systematic, order by order interaction potential theory is developed. Conversely, this theory suggests the functional form for the potential model. At present, we have a simple interaction model that is capable of describing the van der Waals potentials of many systems from R = 0 to R→∞.
Collapse
Affiliation(s)
- X W Sheng
- Department of Physics, Anhui Normal University, Anhui, Wuhu 24100, China
| | | |
Collapse
|
19
|
Schwerdtfeger P, Burrows A, Smits OR. The Lennard-Jones Potential Revisited: Analytical Expressions for Vibrational Effects in Cubic and Hexagonal Close-Packed Lattices. J Phys Chem A 2021; 125:3037-3057. [PMID: 33787272 DOI: 10.1021/acs.jpca.1c00012] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Analytical formulas are derived for the zero-point vibrational energy and anharmonicity corrections of the cohesive energy and the mode Grüneisen parameter within the Einstein model for the cubic lattices (sc, bcc, and fcc) and for the hexagonal close-packed structure. This extends the work done by Lennard-Jones and Ingham in 1924, Corner in 1939, and Wallace in 1965. The formulas are based on the description of two-body energy contributions by an inverse power expansion (extended Lennard-Jones potential). These make use of three-dimensional lattice sums, which can be transformed to fast converging series and accurately determined by various expansion techniques. We apply these new lattice sum expressions to the rare gas solids and discuss associated critical points. The derived formulas give qualitative but nevertheless deep insight into vibrational effects in solids from the lightest (helium) to the heaviest rare gas element (oganesson), both presenting special cases because of strong quantum effects for the former and strong relativistic effects for the latter.
Collapse
Affiliation(s)
- Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study (NZIAS), Massey University Albany, Private Bag 102904, Auckland 0745, New Zealand
| | - Antony Burrows
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study (NZIAS), Massey University Albany, Private Bag 102904, Auckland 0745, New Zealand
| | - Odile R Smits
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study (NZIAS), Massey University Albany, Private Bag 102904, Auckland 0745, New Zealand
| |
Collapse
|
20
|
Sheng X, Toennies JP, Tang KT. Conformal Analytical Potential for All the Rare Gas Dimers over the Full Range of Internuclear Distances. PHYSICAL REVIEW LETTERS 2020; 125:253402. [PMID: 33416396 DOI: 10.1103/physrevlett.125.253402] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2020] [Accepted: 11/18/2020] [Indexed: 06/12/2023]
Abstract
An analytical model for the potential between two rare gas atoms at distances between R=0 to R→∞ is assumed to be conformal with the previously published potential for He_{2} [J. Chem. Phys. 142, 131102 (2015)JCPSA60021-960610.1063/1.4916740]. The potential curves of the rare gas dimers all have the same shape and only depend on the well parameters D_{e} and R_{e}. The potentials and the vibrational levels for the 11 homonuclear and heteronuclear dimers for which recent ab initio calculations are available agree, within several percent, with the ab initio results. For the other rare gas dimers, the new potential provides the first realistic estimates for the potentials.
Collapse
Affiliation(s)
- Xiaowei Sheng
- Department of Physics, Anhui Normal University, Anhui, Wuhu 241000, China
| | - J Peter Toennies
- Max-Planck-Institut für Dynamik und Selbstorganisation, Am Fassberg 17, D-37077, Göttingen, Germany
| | - K T Tang
- Department of Physics, Pacific Lutheran University, Tacoma, Washington 98447, USA
| |
Collapse
|
21
|
Abstract
The role of cohesive r-4 interactions on the existence of a vapor phase and the formation of vapor-liquid equilibria is investigated by performing molecular simulations for the n-4 potential. The cohesive r-4 interactions delay the emergence of a vapor phase until very high temperatures. The critical temperature is up to 5 times higher than normal fluids, as represented by the Lennard-Jones potential. The greatest overall influence on vapor-liquid equilibria is observed for the 5-4 potential, which is the lowest repulsive limit of the potential. Increasing n initially mitigates the influence of r-4 interactions, but the moderating influence declines for n > 12. A relationship is reported between the critical temperature and the Boyle temperature, which allows the critical temperature to be determined for a given n value. The n-4 potential could provide valuable insight into the behavior of non-conventional materials with both very low vapor pressures at elevated temperatures and highly dipolar interactions.
Collapse
Affiliation(s)
- Richard J Sadus
- Centre for Computational Innovations, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
22
|
Garberoglio G, Harvey AH. Path-Integral Calculation of the Second Dielectric and Refractivity Virial Coefficients of Helium, Neon, and Argon. JOURNAL OF RESEARCH OF THE NATIONAL INSTITUTE OF STANDARDS AND TECHNOLOGY 2020; 125:125022. [PMID: 39081565 PMCID: PMC11239192 DOI: 10.6028/jres.125.022] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Accepted: 07/02/2020] [Indexed: 08/02/2024]
Abstract
We present a method to calculate dielectric and refractivity virial coefficients using the path-integral Monte Carlo formulation of quantum statistical mechanics and validate it by comparing our results with equivalent calculations in the literature and with more traditional quantum calculations based on wavefunctions. We use state-of-the-art pair potentials and polarizabilities to calculate the second dielectric and refractivity virial coefficients of helium (both 3He and 4He), neon (both 20Ne and 22Ne), and argon. Our calculations extend to temperatures as low as 1 K for helium, 4 K for neon, and 50 K for argon. We estimate the contributions to the uncertainty of the calculated dielectric virial coefficients for helium and argon, finding that the uncertainty of the pair polarizability is by far the greatest contribution. Agreement with the limited experimental data available is generally good, but our results have smaller uncertainties, especially for helium. Our approach can be generalized in a straightforward manner to higher-order coefficients.
Collapse
Affiliation(s)
- Giovanni Garberoglio
- European Centre for Theoretical Studies in Nuclear
Physics and Related Areas (FBK-ECT*) and Trento Institute for Fundamental Physics and
Applications (TIFPA-INFN), Trento, I-38123, Italy
| | - Allan H. Harvey
- European Centre for Theoretical Studies in Nuclear
Physics and Related Areas (FBK-ECT*) and Trento Institute for Fundamental Physics and
Applications (TIFPA-INFN), Trento, I-38123, Italy
| |
Collapse
|
23
|
Yang M, Song B. Accurate second Kerr virial coefficient of rare gases from the state-of-the-art ab initio potentials and (hyper)polarizabilities. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1706005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Min Yang
- Key Laboratory of Thermo-Fluid Science and Engineering, Ministry of Education, School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, China
| | - Bo Song
- Key Laboratory of Thermo-Fluid Science and Engineering, Ministry of Education, School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, China
| |
Collapse
|
24
|
Uribe FJ, Velasco RM. Exact solutions for shock waves in dilute gases. Phys Rev E 2019; 100:023118. [PMID: 31574705 DOI: 10.1103/physreve.100.023118] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2019] [Indexed: 11/07/2022]
Abstract
In 1922 Becker found an exact solution for shock waves in gases using the Navier-Stokes-Fourier constitutive equations for a Prandtl number of value 3/4 with constant transport coefficients. His analysis has been extended to study some cases where an implicit solution can be found in an exact way. In this work we consider this problem for the so-called soft-spheres model in which the viscosity and thermal conductivity are proportional to a power of the temperature η,κ∝T^{σ}. In particular, we give implicit exact solutions for the Maxwell model (σ=1), hard spheres (σ=1/2), and when σ (the viscosity index) is a natural number.
Collapse
Affiliation(s)
- F J Uribe
- Department of Physics, Universidad Autónoma Metropolitana-Iztapalapa 09340, CDMX, México
| | - R M Velasco
- Department of Physics, Universidad Autónoma Metropolitana-Iztapalapa 09340, CDMX, México
| |
Collapse
|
25
|
Kodrycka M, Patkowski K. Platinum, gold, and silver standards of intermolecular interaction energy calculations. J Chem Phys 2019; 151:070901. [DOI: 10.1063/1.5116151] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Monika Kodrycka
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| |
Collapse
|
26
|
Bezrukov DS, Kleshchina NN, Kalinina IS, Buchachenko AA. Empirically Modified Potentials of Interaction between Rare Gases for Matrix Isolation Problems. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2019. [DOI: 10.1134/s003602441908003x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
27
|
Luo QY, Song B. Accurate internal energy of argon fluid from a state-of-the-art ab initio potential with uncertainty estimations. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.110980] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
28
|
Rutskoy BV, Bezrukov DS. Ab Initio Description of the Structure and Interaction Energy of Perhalomethane Dimers. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2019. [DOI: 10.1134/s0036024419080259] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
29
|
Deiters UK, Sadus RJ. Fully a priori prediction of the vapor-liquid equilibria of Ar, Kr, and Xe from ab initio two-body plus three-body interatomic potentials. J Chem Phys 2019; 151:034509. [DOI: 10.1063/1.5109052] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ulrich K. Deiters
- Institute of Physical Chemistry, University of Cologne, Luxemburger Str. 116, D-50939 Köln, Germany
| | - Richard J. Sadus
- Centre for Computational Innovations, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
30
|
Rivlin T, McKemmish LK, Spinlove KE, Tennyson J. Low temperature scattering with the R-matrix method: argon-argon scattering. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1615143] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Tom Rivlin
- Department of Physics and Astronomy, University College London, London, UK
| | - Laura K. McKemmish
- Department of Physics and Astronomy, University College London, London, UK
- Department of Chemistry, University of New South Wales, Sydney, Australia
| | - K. Eryn Spinlove
- Department of Physics and Astronomy, University College London, London, UK
| | - Jonathan Tennyson
- Department of Physics and Astronomy, University College London, London, UK
| |
Collapse
|
31
|
Deiters UK, Sadus RJ. Two-body interatomic potentials for He, Ne, Ar, Kr, and Xe fromab initiodata. J Chem Phys 2019; 150:134504. [DOI: 10.1063/1.5085420] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ulrich K. Deiters
- Institute of Physical Chemistry, University of Cologne, Luxemburger Str. 116, D-50939 Köln, Germany
| | - Richard J. Sadus
- Centre for Computational Innovations, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
32
|
Sadus RJ. Two-body intermolecular potentials from second virial coefficient properties. J Chem Phys 2019; 150:024503. [DOI: 10.1063/1.5080308] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Richard J. Sadus
- Centre for Computational Innovations, Swinburne University of Technology, P.O. Box 218 Hawthorn, Victoria 3122, Australia
| |
Collapse
|
33
|
Sadus RJ. Molecular simulation of orthobaric isochoric heat capacities near the critical point. Phys Rev E 2019; 99:012139. [PMID: 30780235 DOI: 10.1103/physreve.99.012139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2018] [Indexed: 06/09/2023]
Abstract
A molecular simulation strategy is investigated for detecting the divergence of the isochoric heat capacity (C_{V}) on the vapor and liquid coexistence branches of a fluid near the critical point. The procedure is applied to the empirical Lennard-Jones potential and accurate state-of-the-art ab initio two-body and two-body + three-body potentials for argon. Simulations with the Lennard-Jones potential predict the divergence of C_{V}, and the phenomenon is also observed for both two-body and two-body + three body potentials. The potentials also correctly predict the crossover between vapor and liquid C_{V} values and the subcritical liquid C_{V} minimum, which marks the commencement C_{V} divergence. The effect of three-body interactions is to delay the onset of divergence to higher subcritical temperatures.
Collapse
Affiliation(s)
- Richard J Sadus
- Centre for Computational Innovations, Swinburne University of Technology, P.O. Box 218 Hawthorn, Victoria 3122, Australia
| |
Collapse
|
34
|
Shahbaz M, Szalewicz K. Do Semilocal Density-Functional Approximations Recover Dispersion Energies at Small Intermonomer Separations? PHYSICAL REVIEW LETTERS 2018; 121:113402. [PMID: 30265106 DOI: 10.1103/physrevlett.121.113402] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2018] [Revised: 06/12/2018] [Indexed: 06/08/2023]
Abstract
The methods that add dispersion energies to interaction energies computed using density-functional theory (DFT), known as DFT+D methods, taper off the dispersion energies at distances near van der Waals minima and smaller based on an assumption that DFT starts to reproduce the dispersion energies there. We show that this assumption is not correct as the alleged contribution behaves unphysically and originates to a large extent from nonexchange-correlation terms. Thus, dispersion functions correct DFT in this region for deficiencies unrelated to dispersion interactions.
Collapse
Affiliation(s)
- Muhammad Shahbaz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
35
|
Przybytek M. Dispersion Energy of Symmetry-Adapted Perturbation Theory from the Explicitly Correlated F12 Approach. J Chem Theory Comput 2018; 14:5105-5117. [DOI: 10.1021/acs.jctc.8b00470] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Michał Przybytek
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
36
|
Sadus RJ. Second virial coefficient properties of the n-m Lennard-Jones/Mie potential. J Chem Phys 2018; 149:074504. [DOI: 10.1063/1.5041320] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Richard J. Sadus
- Computational Science Laboratory, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
37
|
Goel H, Ling S, Ellis BN, Taconi A, Slater B, Rai N. Predicting vapor liquid equilibria using density functional theory: A case study of argon. J Chem Phys 2018; 148:224501. [PMID: 29907054 DOI: 10.1063/1.5025726] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Himanshu Goel
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Sanliang Ling
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Breanna Nicole Ellis
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Anna Taconi
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| | - Ben Slater
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Neeraj Rai
- Dave C. Swalm School of Chemical Engineering, and Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, Mississippi 39762, USA
| |
Collapse
|
38
|
Hellmann R, Jäger B, Bich E. State-of-the-art ab initio potential energy curve for the xenon atom pair and related spectroscopic and thermophysical properties. J Chem Phys 2018; 147:034304. [PMID: 28734299 DOI: 10.1063/1.4994267] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new ab initio interatomic potential energy curve for two ground-state xenon atoms is presented. It is based on supermolecular calculations at the coupled-cluster level with single, double, and perturbative triple excitations [CCSD(T)] employing basis sets up to sextuple-zeta quality, which were developed as part of this work. In addition, corrections were determined for higher coupled-cluster levels up to CCSDTQ as well as for scalar and spin-orbit relativistic effects at the CCSD(T) level. A physically motivated analytical function was fitted to the calculated interaction energies and used to compute the vibrational spectrum of the dimer, the second virial coefficient, and the dilute gas transport properties. The agreement with the best available experimental data for the investigated properties is excellent; the new potential function is superior not only to previous ab initio potentials but also to the most popular empirical ones.
Collapse
Affiliation(s)
- Robert Hellmann
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| | - Benjamin Jäger
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| | - Eckard Bich
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| |
Collapse
|
39
|
Dutta NN, Patkowski K. Improving “Silver-Standard” Benchmark Interaction Energies with Bond Functions. J Chem Theory Comput 2018; 14:3053-3070. [DOI: 10.1021/acs.jctc.8b00204] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Narendra Nath Dutta
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
40
|
Myatt PT, Dham AK, Chandrasekhar P, McCourt FRW, Le Roy RJ. A new empirical potential energy function for Ar2. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1437932] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Philip T. Myatt
- Department of Chemistry, University of Waterloo, Waterloo, Ontario, Canada
| | - Ashok K. Dham
- Department of Chemistry, University of Waterloo, Waterloo, Ontario, Canada
- Department of Physics, Punjabi University, Patiala, India
| | | | | | - Robert J. Le Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario, Canada
| |
Collapse
|
41
|
Shirkov L, Sladek V. Benchmark CCSD-SAPT study of rare gas dimers with comparison to MP-SAPT and DFT-SAPT. J Chem Phys 2017; 147:174103. [DOI: 10.1063/1.4997569] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Leonid Shirkov
- Faculty of Chemistry, Adam Mickiewicz University, Umultowska 89b, 61-614 Poznań, Poland
| | - Vladimir Sladek
- Institute of Chemistry–Centre for Glycomics, Slovak Academy of Sciences, 845 38 Bratislava, Slovakia
| |
Collapse
|
42
|
Jerabek P, Smits O, Pahl E, Schwerdtfeger P. A relativistic coupled-cluster interaction potential and rovibrational constants for the xenon dimer. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1359347] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Paul Jerabek
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study and the Institute for Natural and Mathematical Sciences, Massey University, Auckland, New Zealand
| | - Odile Smits
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study and the Institute for Natural and Mathematical Sciences, Massey University, Auckland, New Zealand
| | - Elke Pahl
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study and the Institute for Natural and Mathematical Sciences, Massey University, Auckland, New Zealand
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study and the Institute for Natural and Mathematical Sciences, Massey University, Auckland, New Zealand
| |
Collapse
|
43
|
Shizgal BD. A comparison of pseudospectral methods for the solution of the Schrödinger equation; the Lennard-Jones ( n , 6) potential. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.05.009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
44
|
Vlasiuk M, Sadus RJ. Ab initio interatomic potentials and the thermodynamic properties of fluids. J Chem Phys 2017; 147:024505. [DOI: 10.1063/1.4991012] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Vlasiuk M, Sadus RJ. Predicting vapor-liquid phase equilibria with augmented ab initio interatomic potentials. J Chem Phys 2017; 146:244504. [DOI: 10.1063/1.4986917] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Maryna Vlasiuk
- Centre for Molecular Simulation, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| | - Richard J. Sadus
- Centre for Molecular Simulation, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
46
|
Jäger B, Bich E. Thermophysical properties of krypton-helium gas mixtures from ab initio pair potentials. J Chem Phys 2017; 146:214302. [PMID: 28595411 DOI: 10.1063/1.4984100] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new potential energy curve for the krypton-helium atom pair was developed using supermolecular ab initio computations for 34 interatomic distances. Values for the interaction energies at the complete basis set limit were obtained from calculations with the coupled-cluster method with single, double, and perturbative triple excitations and correlation consistent basis sets up to sextuple-zeta quality augmented with mid-bond functions. Higher-order coupled-cluster excitations up to the full quadruple level were accounted for in a scheme of successive correction terms. Core-core and core-valence correlation effects were included. Relativistic corrections were considered not only at the scalar relativistic level but also using full four-component Dirac-Coulomb and Dirac-Coulomb-Gaunt calculations. The fitted analytical pair potential function is characterized by a well depth of 31.42 K with an estimated standard uncertainty of 0.08 K. Statistical thermodynamics was applied to compute the krypton-helium cross second virial coefficients. The results show a very good agreement with the best experimental data. Kinetic theory calculations based on classical and quantum-mechanical approaches for the underlying collision dynamics were utilized to compute the transport properties of krypton-helium mixtures in the dilute-gas limit for a large temperature range. The results were analyzed with respect to the orders of approximation of kinetic theory and compared with experimental data. Especially the data for the binary diffusion coefficient confirm the predictive quality of the new potential. Furthermore, inconsistencies between two empirical pair potential functions for the krypton-helium system from the literature could be resolved.
Collapse
Affiliation(s)
- Benjamin Jäger
- Institut für Chemie, Universität Rostock, D-18059 Rostock, Germany
| | - Eckard Bich
- Institut für Chemie, Universität Rostock, D-18059 Rostock, Germany
| |
Collapse
|
47
|
Kleshchina NN, Korchagina KA, Bezrukov DS, Buchachenko AA. Modeling of Manganese Atom and Dimer Isolated in Solid Rare Gases: Structure, Stability, and Effect on Spin Coupling. J Phys Chem A 2017; 121:2429-2441. [PMID: 28263615 DOI: 10.1021/acs.jpca.6b12444] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Structures and energies of the trapping sites of manganese atom and dimer in solid Ar, Kr, and Xe are investigated within the classical model, which balances local distortion and long-range crystal order of the host and provides a means to estimate the relative site stabilities. The model is implemented with the additive pairwise potential field based on the ab initio and best empirical interatomic potential functions. In agreement with experiment, Mn single substitution (SS) and tetrahedral vacancy (TV) occupation are identified as stable for Ar and Kr, whereas the SS site is only found for Xe. Stable trapping sites of the weakly bound Mn2 dimer are shown to be the mergers of SS and/or TV atomic sites. For Ar, (SS + SS) and (TV + TV) sites are close in energy, whereas (SS + TV) site lies higher. The (SS + SS) accommodation is identified as the only stable site in Kr and Xe at low energies. The results are compared with the resonance Raman, electron spin resonance, and absorption spectroscopy data. Reproducing the numbers of stable sites, the calculations tend to underestimate the matrix effect on the dimer vibrational frequency and spin-spin coupling constant. Nonetheless, the level of agreement is found to be informative for tentative assignments of the complex features seen in Mn2 matrix isolation spectroscopy.
Collapse
Affiliation(s)
- Nadezhda N Kleshchina
- Department of Chemistry, M. V. Lomonosov Moscow State University , Moscow 119991, Russia
| | - Kseniia A Korchagina
- Department of Chemistry, M. V. Lomonosov Moscow State University , Moscow 119991, Russia
| | - Dmitry S Bezrukov
- Department of Chemistry, M. V. Lomonosov Moscow State University , Moscow 119991, Russia.,Skolkovo Institute of Science and Technology , Skolkovo Innovation Center, Building 3, Moscow 143026, Russia
| | - Alexei A Buchachenko
- Skolkovo Institute of Science and Technology , Skolkovo Innovation Center, Building 3, Moscow 143026, Russia
| |
Collapse
|
48
|
Hellmann R. Nonadditive three-body potential and third to eighth virial coefficients of carbon dioxide. J Chem Phys 2017; 146:054302. [DOI: 10.1063/1.4974995] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Robert Hellmann
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| |
Collapse
|
49
|
Schwerdtfeger P, Tonner R, Moyano GE, Pahl E. Towards J/mol Accuracy for the Cohesive Energy of Solid Argon. Angew Chem Int Ed Engl 2016; 55:12200-5. [PMID: 27593519 DOI: 10.1002/anie.201605875] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2016] [Indexed: 11/07/2022]
Abstract
The cohesive energies of argon in its cubic and hexagonal closed packed structures are computed with an unprecedented accuracy of about 5 J mol(-1) (corresponding to 0.05 % of the total cohesive energy). The same relative accuracy with respect to experimental data is also found for the face-centered cubic lattice constant deviating by ca. 0.003 Å. This level of accuracy was enabled by using high-level theoretical, wave-function-based methods within a many-body decomposition of the interaction energy. Static contributions of two-, three-, and four-body fragments of the crystal are all individually converged to sub-J mol(-1) accuracy and complemented by harmonic and anharmonic vibrational corrections. Computational chemistry is thus achieving or even surpassing experimental accuracy for the solid-state rare gases.
Collapse
Affiliation(s)
- Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, Private Bag 102904, 0632, Auckland, New Zealand.
| | - Ralf Tonner
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032, Marburg, Germany
| | - Gloria E Moyano
- Instituto de Química, Universidad de Antioquia, AA 1126, Medellín, Colombia
| | - Elke Pahl
- Centre for Theoretical Chemistry and Physics, Institute for Natural and Mathematical Sciences, Massey University Auckland, Private Bag 102904, 0632, Auckland, New Zealand
| |
Collapse
|
50
|
Hochgenaue Berechnung der Kohäsionsenergie von festem Argon mit Abweichungen im J/mol-Bereich. Angew Chem Int Ed Engl 2016. [DOI: 10.1002/ange.201605875] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|