1
|
Moscato P, Haque MN. New alternatives to the Lennard-Jones potential. Sci Rep 2024; 14:11169. [PMID: 38750117 PMCID: PMC11096349 DOI: 10.1038/s41598-024-60835-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2023] [Accepted: 04/28/2024] [Indexed: 05/18/2024] Open
Abstract
We present a new method for approximating two-body interatomic potentials from existing ab initio data based on representing the unknown function as an analytic continued fraction. In this study, our method was first inspired by a representation of the unknown potential as a Dirichlet polynomial, i.e., the partial sum of some terms of a Dirichlet series. Our method allows for a close and computationally efficient approximation of the ab initio data for the noble gases Xenon (Xe), Krypton (Kr), Argon (Ar), and Neon (Ne), which are proportional to r - 6 and to a very simple d e p t h = 1 truncated continued fraction with integer coefficients and depending on n - r only, where n is a natural number (with n = 13 for Xe, n = 16 for Kr, n = 17 for Ar, and n = 27 for Neon). For Helium (He), the data is well approximated with a function having only one variable n - r with n = 31 and a truncated continued fraction with d e p t h = 2 (i.e., the third convergent of the expansion). Also, for He, we have found an interesting d e p t h = 0 result, a Dirichlet polynomial of the formk 1 6 - r + k 2 48 - r + k 3 72 - r (withk 1 , k 2 , k 3 all integers), which provides a surprisingly good fit, not only in the attractive but also in the repulsive region. We also discuss lessons learned while facing the surprisingly challenging non-linear optimisation tasks in fitting these approximations and opportunities for parallelisation.
Collapse
Affiliation(s)
- Pablo Moscato
- School of Information and Physical Sciences, The University of Newcastle, Callaghan, NSW, 2308, Australia.
| | - Mohammad Nazmul Haque
- School of Information and Physical Sciences, The University of Newcastle, Callaghan, NSW, 2308, Australia
- ResTech Pty Ltd, CE Building, Design Drive, Callaghan, NSW, 2308, Australia
| |
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
|
Paterson T, Bannerman MN, Lue L. Using the Zeno line to assess and refine molecular models. J Chem Phys 2024; 160:154503. [PMID: 38624125 DOI: 10.1063/5.0192770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2023] [Accepted: 04/01/2024] [Indexed: 04/17/2024] Open
Abstract
The Zeno line is the locus of points on the temperature-density plane where the compressibility factor of the fluid is equal to one. It has been observed to be straight for a broad variety of real fluids, although the underlying reasons for this are still unclear. In this work, a detailed study of the Zeno line and its relation to the vapor-liquid coexistence curve is performed for two simple model pair-potential fluids: attractive square-well fluids with varying well-widths λ and Mie n-6 fluids with different repulsive exponents n. Interestingly, the Zeno lines of these fluids are curved, regardless of the value of λ or n. We find that for square-well fluids, λ ≈ 1.8 presents a Zeno line, which is the most linear over the largest temperature range. For Mie n-6 fluids, we find that the straightest Zeno line occurs for n between 8 and 10. Additionally, the square-well and Mie fluids with the straightest Zeno line showed the closest quantitative agreement with the vapor-liquid coexistence curve for experimental fluids that follow the principle of corresponding states (e.g., argon, xenon, krypton, methane, nitrogen, and oxygen). These results suggest that the Zeno line can provide a useful additional feature, in complement to other properties, such as the phase envelope, to evaluate molecular models.
Collapse
Affiliation(s)
- Thomas Paterson
- Department of Chemical and Process Engineering, University of Strathclyde, James Weir Building, 75 Montrose Street, Glasgow G1 1XJ, United Kingdom
| | - Marcus N Bannerman
- School of Engineering, University of Aberdeen, Aberdeen AB24 3UE, United Kingdom
| | - Leo Lue
- Department of Chemical and Process Engineering, University of Strathclyde, James Weir Building, 75 Montrose Street, Glasgow G1 1XJ, United Kingdom
| |
Collapse
|
4
|
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
|
5
|
Saric D, Bell IH, Guevara-Carrion G, Vrabec J. Influence of repulsion on entropy scaling and density scaling of monatomic fluids. J Chem Phys 2024; 160:104503. [PMID: 38456532 DOI: 10.1063/5.0196592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 02/18/2024] [Indexed: 03/09/2024] Open
Abstract
Entropy scaling is applied to the shear viscosity, self-diffusion coefficient, and thermal conductivity of simple monatomic fluids. An extensive molecular dynamics simulation series is performed to obtain these transport properties and the residual entropy of three potential model classes with variable repulsive exponents: n, 6 Mie (n = 9, 12, 15, and 18), Buckingham's exponential-six (α = 12, 14, 18, and 30), and Tang-Toennies (αT = 4.051, 4.275, and 4.600). A wide range of liquid and supercritical gas- and liquid-like states is covered with a total of 1120 state points. Comparisons to equations of state, literature data, and transport property correlations are made. Although the absolute transport property values within a given potential model class may strongly depend on the repulsive exponent, it is found that the repulsive steepness plays a negligible role when entropy scaling is applied. Hence, the plus-scaled transport properties of n, 6 Mie, exponential-six, and Tang-Toennies fluids lie basically on one master curve, which closely corresponds with entropy scaling correlations for the Lennard-Jones fluid. This trend is confirmed by literature data of n, 6 Mie, and exponential-six fluids. Furthermore, entropy scaling holds for state points where the Pearson correlation coefficient R is well below 0.9. The condition R > 0.9 for strongly correlating liquids is thus not necessary for the successful application of entropy scaling, pointing out that isomorph theory may be a part of a more general framework that is behind the success of entropy scaling. Density scaling reveals a strong influence of the repulsive exponent on this particular approach.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587 Berlin, Germany
| | - Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| | | | - Jadran Vrabec
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587 Berlin, Germany
| |
Collapse
|
6
|
Ströker P, Meier K. Vapor-liquid equilibrium and thermodynamic properties of saturated argon and krypton from Monte Carlo simulations using ab initio potentials. J Chem Phys 2024; 160:094503. [PMID: 38426525 DOI: 10.1063/5.0196466] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2024] [Accepted: 01/29/2024] [Indexed: 03/02/2024] Open
Abstract
Vapor-liquid equilibria and thermodynamic properties of saturated argon and krypton were calculated by semi-classical Monte Carlo simulations with the NpT + test particle method using ab initio potentials for the two-body and nonadditive three-body interactions. The NpT + test particle method was extended to the calculation of second-order thermodynamic properties, such as the isochoric and isobaric heat capacities or the speed of sound, of the saturated liquid and vapor by using our recently developed approach for the systematic calculation of arbitrary thermodynamic properties in the isothermal-isobaric ensemble. Generally, the results for all simulated properties agree well with experimental data and the current reference equations of state for argon and krypton. In particular, the results for the vapor pressure and for the density and speed of sound of the saturated liquid and vapor agree with the most accurate experimental data for both noble gases almost within the uncertainty of these data, a level of agreement unprecedented for many-particle simulations. This study demonstrates that the vapor-liquid equilibrium and thermodynamic properties at saturation of a pure fluid can be predicted by Monte Carlo simulations with high accuracy when the intermolecular interactions are described by state-of-the-art ab initio pair and nonadditive three-body potentials and quantum effects are accounted for.
Collapse
Affiliation(s)
- Philipp Ströker
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Karsten Meier
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| |
Collapse
|
7
|
Deiters UK, Sadus RJ. An intermolecular potential for hydrogen: Classical molecular simulation of pressure-density-temperature behavior, vapor-liquid equilibria, and critical and triple point properties. J Chem Phys 2023; 158:2890478. [PMID: 37184016 DOI: 10.1063/5.0148248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Accepted: 04/25/2023] [Indexed: 05/16/2023] Open
Abstract
An intermolecular potential is reported for molecular hydrogen that combines two-body interactions from ab initio data with three-body interactions. The accuracy of the two-body potential is validated by comparison with experimental second virial coefficient data. Experimental pressure-density-temperature data are used to validate the addition of three-body interactions, often yielding very accurate predictions. Classical Monte Carlo simulations that neglect quantum effects are reported for the vapor-liquid equilibria (VLE), critical properties, and the triple point. A comparison with experimental data indicates that the effect of quantum interactions is to narrow the VLE phase envelope and to lower the critical temperature. The three-body interactions have a considerable influence on the phase behavior, resulting in good agreement with the experimental density. The critical properties of the two-body + three-body potential for hydrogen provide an alternative set of input parameters to improve the accuracy of theoretical predictions at temperatures above 100 K. In the vicinity of the critical point, the coexistence densities do not obey the law of rectilinear diameters, which is a feature that has largely been overlooked in both experimental data and reference equations of state.
Collapse
Affiliation(s)
- Ulrich K Deiters
- Institute of Physical Chemistry, Faculty of Mathematics and Natural Sciences, University of Cologne, Greinstr 4-6, D-50939 Köln, Germany
| | - Richard J Sadus
- Department of Computer Technologies, Swinburne University of Technology, Wurundjeri Country, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
8
|
Lambo RL, Koyanagi GK, Ragyanszki A, Horbatsch M, Fournier R, Hessels EA. Calculation of the local environment of a barium monofluoride molecule in an argon matrix: a step towards using matrix-isolated BaF for determining the electron electric dipole moment. Mol Phys 2023. [DOI: 10.1080/00268976.2023.2198044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/07/2023]
Affiliation(s)
- Ricardo L. Lambo
- Departments of Physics and Chemistry, York University, Toronto, ON, Canada (EDM Collaboration)
| | - Gregory K. Koyanagi
- Departments of Physics and Chemistry, York University, Toronto, ON, Canada (EDM Collaboration)
| | - Anita Ragyanszki
- Departments of Physics and Chemistry, York University, Toronto, ON, Canada (EDM Collaboration)
| | - Marko Horbatsch
- Departments of Physics and Chemistry, York University, Toronto, ON, Canada (EDM Collaboration)
| | - Rene Fournier
- Departments of Physics and Chemistry, York University, Toronto, ON, Canada (EDM Collaboration)
| | - Eric A. Hessels
- Departments of Physics and Chemistry, York University, Toronto, ON, Canada (EDM Collaboration)
| |
Collapse
|
9
|
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
|
10
|
Ströker P, Hellmann R, Meier K. Thermodynamic properties of krypton from Monte Carlo simulations using ab initio potentials. J Chem Phys 2022; 157:114504. [PMID: 36137797 DOI: 10.1063/5.0107851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ten different thermodynamic properties of the noble gas krypton were calculated by Monte Carlo simulations in the isothermal-isobaric ensemble using a highly accurate ab initio pair potential, Feynman-Hibbs corrections for quantum effects, and an extended Axilrod-Teller-Muto potential to account for nonadditive three-body interactions. Fourteen state points at a liquid and a supercritical isotherm were simulated. To obtain results representative for macroscopic systems, simulations with several particle numbers were carried out and extrapolated to the thermodynamic limit. Our results agree well with experimental data from the literature, an accurate equation of state for krypton, and a recent virial equation of state (VEOS) for krypton in the region where the VEOS has converged. These results demonstrate that very good agreement between simulation and experiment can only be achieved if nonadditive three-body interactions and quantum effects are taken into account.
Collapse
Affiliation(s)
- Philipp Ströker
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Robert Hellmann
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Karsten Meier
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| |
Collapse
|
11
|
Ströker P, Hellmann R, Meier K. Thermodynamic properties of argon from Monte Carlo simulations using ab initio potentials. Phys Rev E 2022; 105:064129. [PMID: 35854585 DOI: 10.1103/physreve.105.064129] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2022] [Accepted: 05/16/2022] [Indexed: 06/15/2023]
Abstract
Ten different thermodynamic properties of the noble gas argon in the liquid and supercritical regions were obtained from semiclassical Monte Carlo simulations in the isothermal-isobaric ensemble using ab initio potentials for the two-body and nonadditive three-body interactions. Our results for the density and speed of sound agree with the most accurate experimental data for argon almost within the uncertainty of these data, a level of agreement unprecedented for many-particle simulations. This demonstrates the high predictive but yet unexploited power of ab initio potentials in the field of molecular modeling and simulation for thermodynamic properties of fluids.
Collapse
Affiliation(s)
- Philipp Ströker
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Robert Hellmann
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| | - Karsten Meier
- Institut für Thermodynamik, Helmut-Schmidt-Universität/Universität der Bundeswehr Hamburg, Holstenhofweg 85, 22043 Hamburg, Germany
| |
Collapse
|
12
|
Nitzke I, Pohl S, Thol M, Span R, Vrabec J. How well does the Tang-Toennies potential represent the thermodynamic properties of argon? Mol Phys 2022. [DOI: 10.1080/00268976.2022.2078240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Isabel Nitzke
- Thermodynamics and Process Engineering, Technische Universtität Berlin, Berlin, Germany
| | - Sven Pohl
- Thermodynamics, Ruhr-Universität Bochum, Bochum, Germany
| | - Monika Thol
- Thermodynamics, Ruhr-Universität Bochum, Bochum, Germany
| | - Roland Span
- Thermodynamics, Ruhr-Universität Bochum, Bochum, Germany
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technische Universtität Berlin, Berlin, Germany
| |
Collapse
|
13
|
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
|
14
|
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
|
15
|
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
|
16
|
Walker PJ, Zhao T, Haslam AJ, Jackson G. Ab initio development of generalized Lennard-Jones (Mie) force fields for predictions of thermodynamic properties in advanced molecular-based SAFT equations of state. J Chem Phys 2022; 156:154106. [PMID: 35459299 DOI: 10.1063/5.0087125] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A methodology for obtaining molecular parameters of a modified statistical associating fluid theory for variable-range interactions of Mie form (SAFT-VR Mie) equation of state (EoS) from ab initio calculations is proposed for non-associative species that can be modeled as single spherical segments. The methodology provides a strategy to map interatomic or intermolecular potentials obtained from ab initio quantum-chemistry calculations to the corresponding Mie potentials that can be used within the SAFT-VR Mie EoS. The inclusion of corrections for quantum and many-body effects allows for an excellent, fully predictive description of the vapor-liquid envelope and other bulk thermodynamic properties of noble gases; this description is of similar or superior quality to that obtained using SAFT-VR Mie with parameters regressed in the traditional way using experimental thermodynamic-property data. The methodology is extended to an anisotropic species, methane, where similar levels of accuracy are obtained. The efficacy of using less-accurate quantum-chemistry methods in this methodology is explored, showing that these methods do not provide satisfactory results, although we note that the description is nevertheless substantially better than those obtained using the conductor-like screening model for describing real solvents (COSMO-RS), the only other fully predictive ab initio method currently available. Overall, the reliance on thermophysical data is completely dispensed with, providing the first extensible, wholly predictive SAFT-type EoSs.
Collapse
Affiliation(s)
- Pierre J Walker
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Tianpu Zhao
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Andrew J Haslam
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - George Jackson
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
17
|
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]
|
18
|
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
|
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
|
Sadus RJ. Effect of the range of particle cohesion on the phase behavior and thermodynamic properties of fluids. J Chem Phys 2020; 153:244502. [PMID: 33380081 DOI: 10.1063/5.0031517] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Molecular simulations are performed for the (m + 1, m) potential to systematically investigate the effect of changing the range of particle cohesion on both vapor-liquid equilibria and thermodynamic properties of fluids. The results are reported for m = 4-11, which represent a progressive narrowing of the potential energy well. The conventional Lennard-Jones potential is used as a reference point for normal fluid behavior. Small values of m result in a broadening of the phase envelope compared with the Lennard-Jones potential, whereas a contraction is observed in other cases. The critical properties are reported, and a relationship between the critical temperature and the Boyle temperature is determined. The low values of the critical compressibility factor when m < 6 reflect the behavior observed for real fluids such as n-alkanes. The results for supercritical thermodynamic properties are much more varied. Properties such as pressure, potential energy, isochoric thermal pressure coefficient, and thermal expansion coefficient vary consistently with m, whereas other properties such as the Joule-Thomson coefficient exhibit much more nuanced behavior. Maximum and minimum values are reported for both the isochoric heat capacity and isothermal compressibility. A minimum in the speed of sound is also observed.
Collapse
Affiliation(s)
- Richard J Sadus
- Centre for Computational Innovations, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
21
|
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
|
22
|
Hogan A, Space B. Next-Generation Accurate, Transferable, and Polarizable Potentials for Material Simulations. J Chem Theory Comput 2020; 16:7632-7644. [PMID: 33251798 DOI: 10.1021/acs.jctc.0c00837] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
PHAHST (potentials with high accuracy, high speed, and transferability) intermolecular potential energy functions have been developed from first principles for H2, N2, the noble gases, and a metal-organic material, HKUST-1. The potentials are designed from the outset to be transferable to heterogeneous environments including porous materials, interfaces, and material simulations. This is accomplished by theoretically justified choices for all functional forms, parameters, and mixing rules, including explicit polarization in every environment and fitting to high quality electronic structure calculations using methods that are tractable for real systems. The models have been validated in neat systems by comparison to second virial coefficients and bulk pressure-density isotherms. For inhomogeneous applications, our main target, comparisons are presented to previously published experimental studies on the metal-organic material HKUST-1 including adsorption, isosteric heats of adsorption, binding site locations, and binding site energies. A systematic prescription is provided for developing compatible potentials for additional small molecules and materials. The resulting models are recommended for use in complex heterogeneous simulations where existing potentials may be inadequate.
Collapse
Affiliation(s)
- Adam Hogan
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, Florida 33620-5250, United States
| | - Brian Space
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, Florida 33620-5250, United States
| |
Collapse
|
23
|
Sadus RJ. Combining intermolecular potentials for the prediction of fluid properties: Two-body and three-body interactions. J Chem Phys 2020; 153:214509. [DOI: 10.1063/5.0033466] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [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
|
24
|
Bell IH. Effective hardness of interaction from thermodynamics and viscosity in dilute gases. J Chem Phys 2020; 152:164508. [PMID: 32357769 DOI: 10.1063/5.0007583] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The hardness of the effective inverse power law (IPL) potential, which can be obtained from thermodynamics or collision integrals, can be used to demonstrate similarities between thermodynamic and transport properties. This link is investigated for systems of increasing complexity (i.e., the EXP, square-well, Lennard-Jones, and Stockmayer potentials; ab initio results for small molecules; and rigid linear chains of Lennard-Jones sites). These results show that while the two approaches do not yield precisely the same values of effective IPL exponent, their qualitative behavior is intriguingly similar, offering a new way of understanding the effective interactions between molecules, especially at high temperatures. In both approaches, the effective hardness is obtained from a double-logarithmic temperature derivative of an effective area.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| |
Collapse
|
25
|
Altun A, Neese F, Bistoni G. HFLD: A Nonempirical London Dispersion-Corrected Hartree-Fock Method for the Quantification and Analysis of Noncovalent Interaction Energies of Large Molecular Systems †. J Chem Theory Comput 2019; 15:5894-5907. [PMID: 31538779 DOI: 10.1021/acs.jctc.9b00425] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A nonempirical quantum mechanical method for the efficient and accurate quantification and analysis of intermolecular interactions is presented and tested on existing benchmark sets. The leading idea here is to focus on the intermolecular part of the correlation energy that contains the all-important London dispersion (LD) interaction. To keep the cost of the method low, essentially at the level of a Hartree-Fock (HF) calculation, the intramolecular part of the correlation energy is neglected. We also neglect the nondispersive parts of the intermolecular correlation energy. This scheme that we denote as Hartree-Fock plus London dispersion (HFLD) can be readily realized on the basis of the recently reported multilevel implementation of the domain-based local pair natural orbital coupled-cluster (DLPNO-CC) theory in conjunction with the well-established local energy decomposition (LED) analysis. The accuracy and efficiency of the HFLD method are evaluated on rare gas dimers, on the S66 and L7 benchmark sets of noncovalent interactions, and on an additional set (LP14) consisting of bulky Lewis pairs held together by intermolecular interactions of various strengths, with interaction energies ranging from -8 to -107 kcal/mol. It is first shown that the LD energy calculated with this approach is essentially identical to that obtained from the full DLPNO-CCSD(T)/LED calculation, with a mean absolute error of 0.2 kcal/mol on the S66 benchmark set. Moreover, in terms of the overall interaction energies, the HFLD method shows an efficiency that is comparable to that of the HF method, while retaining an accuracy between that of the DLPNO-CCSD and DLPNO-CCSD(T) schemes. Since the underlying DLPNO-CCSD method is linear scaling with respect to the system size, the HFLD approach also does not lead to new bottlenecks for large systems. As an illustrative example of its efficiency, the HFLD scheme was applied to the interaction between the substrate and the residues in the active site of the cyclohexanone monooxygenase enzyme. The excellent cost/performance ratio indicates that the HFLD method opens new avenues for the accurate calculation and analysis of noncovalent interaction energies in large molecular systems.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1 , D-45470 Mülheim an der Ruhr , Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1 , D-45470 Mülheim an der Ruhr , Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1 , D-45470 Mülheim an der Ruhr , Germany
| |
Collapse
|
26
|
Abstract
Symmetry Adapted Perturbation Theory (SAPT) has become an important tool when predicting and analyzing intermolecular interactions. Unfortunately, Density Functional Theory (DFT)-SAPT, which uses DFT for the underlying monomers, has some arbitrariness concerning the exchange-correlation potential and the exchange-correlation kernel involved. By using ab initio Brueckner Doubles densities and constructing Kohn-Sham orbitals via the Zhao-Morrison-Parr (ZMP) method, we are able to lift the dependence of DFT-SAPT on DFT exchange-correlation potential models in first order. This way, we can compute the monomers at the coupled-cluster level of theory and utilize SAPT for the intermolecular interaction energy. The resulting ZMP-SAPT approach is tested for small dimer systems involving rare gas atoms, cations, and anions and shown to compare well with the Tang-Toennies model and coupled cluster results.
Collapse
Affiliation(s)
- A Daniel Boese
- Institute of Chemistry, Physical and Theoretical Chemistry, University of Graz, Heinrichstrasse 28/IV, 8010 Graz, Austria
| | - Georg Jansen
- Faculty of Chemistry, University of Duisburg-Essen, Universitätsstraße 5, 45117 Essen, Germany
| |
Collapse
|
27
|
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
|
28
|
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
|
29
|
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
|
30
|
Piel H, Chrysos M. A shortcut to quantum-mechanical absolute scattering phase-shift computations in van der Waals systems. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1484951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Henri Piel
- Laboratoire MOLTECH-Anjou, LUNAM Université, Université d'Angers, CNRS UMR 6200, Angers, France
| | - Michael Chrysos
- Laboratoire MOLTECH-Anjou, LUNAM Université, Université d'Angers, CNRS UMR 6200, Angers, France
| |
Collapse
|
31
|
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
|
32
|
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
|
33
|
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
|
34
|
Affiliation(s)
- Frederike Jaeger
- Department of Physics, Imperial College London, London SW7 2AZ, United Kingdom
| | - Omar K. Matar
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Erich A. Müller
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
35
|
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
|
36
|
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
|
37
|
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
|
38
|
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
|
39
|
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
|
40
|
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
|
41
|
Schwerdtfeger P. Toward an accurate description of solid-state properties of superheavy elements. EPJ WEB OF CONFERENCES 2016. [DOI: 10.1051/epjconf/201613107004] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
42
|
Walton JR, Rivera-Rivera LA, Lucchese RR, Bevan JW. Morse, Lennard-Jones, and Kratzer Potentials: A Canonical Perspective with Applications. J Phys Chem A 2016; 120:8347-8359. [DOI: 10.1021/acs.jpca.6b05371] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jay R. Walton
- Department of Mathematics, Texas A&M University, College Station, Texas 77843-3368, United States
| | - Luis A. Rivera-Rivera
- Department of Chemistry, Texas A&M University, College Station, Texas 77843-3255, United States
| | - Robert R. Lucchese
- Department of Chemistry, Texas A&M University, College Station, Texas 77843-3255, United States
| | - John W. Bevan
- Department of Chemistry, Texas A&M University, College Station, Texas 77843-3255, United States
| |
Collapse
|
43
|
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
|
44
|
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]
|
45
|
Sun H, Kang D, Dai J, Ma W, Zhou L, Zeng J. First-principles study on equation of states and electronic structures of shock compressed Ar up to warm dense regime. J Chem Phys 2016; 144:124503. [DOI: 10.1063/1.4943767] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Huayang Sun
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Dongdong Kang
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Jiayu Dai
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Wen Ma
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Liangyuan Zhou
- Academy of Ocean Science and Engineering, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | - Jiaolong Zeng
- Department of Physics, College of Science, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
- IFSA Collaborative Innovation Center, Shanghai Jiao Tong University, Shanghai 200240, People’s Republic of China
| |
Collapse
|
46
|
Jäger B, Hellmann R, Bich E, Vogel E. State-of-the-art ab initio potential energy curve for the krypton atom pair and thermophysical properties of dilute krypton gas. J Chem Phys 2016; 144:114304. [DOI: 10.1063/1.4943959] [Citation(s) in RCA: 72] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Benjamin Jäger
- Institut für Chemie, Universität Rostock, D-18059 Rostock, Germany
| | - Robert Hellmann
- Institut für Chemie, Universität Rostock, D-18059 Rostock, Germany
| | - Eckard Bich
- Institut für Chemie, Universität Rostock, D-18059 Rostock, Germany
| | - Eckhard Vogel
- Institut für Chemie, Universität Rostock, D-18059 Rostock, Germany
| |
Collapse
|
47
|
Řezáč J, Hobza P. Benchmark Calculations of Interaction Energies in Noncovalent Complexes and Their Applications. Chem Rev 2016; 116:5038-71. [DOI: 10.1021/acs.chemrev.5b00526] [Citation(s) in RCA: 281] [Impact Index Per Article: 35.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Řezáč
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
| | - Pavel Hobza
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
48
|
Tian C, Liu F, Cai L, Yuan H, Chen H, Zhong M. Ab initio calculations of many-body interactions for compressed solid argon. J Chem Phys 2015; 143:174506. [PMID: 26547175 DOI: 10.1063/1.4935050] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
An investigation on many-body effects of solid argon at high pressure was conducted based on a many-body expansion of interaction energy. The three- and four-body terms in the expansion were calculated using the coupled-cluster method with single, double, and noniterative triple theory and incremental method, in which the configurations of argon trimers and tetramers were chosen as the same as those in the actual lattice. The four-body interactions in compressed solid argon were estimated for the first time, and the three-body interaction ab initio calculations were extended to a small distance. It shows that the four-body contribution is repulsive at high densities and effectively cancels the three-body lattice energy. The dimer potential plus three-body interaction can well reproduce the measurements of equation of state at pressure approximately lower than ∼60 GPa, when including the four-body effects extends the agreement up to the maximum experimental pressure of 114 GPa.
Collapse
Affiliation(s)
- Chunling Tian
- School of Physical Science and Technology, Southwest University, Chongqing 400715, China
| | - Fusheng Liu
- Institute of High Pressure Physics, Southwest Jiaotong University, Chengdu 610031, China
| | - Lingcang Cai
- Laboratory for Shock Wave and Detonation Physics Research, China Academy of Engineering Physics, Mianyang 621900, China
| | - Hongkuan Yuan
- School of Physical Science and Technology, Southwest University, Chongqing 400715, China
| | - Hong Chen
- School of Physical Science and Technology, Southwest University, Chongqing 400715, China
| | - Mingmin Zhong
- School of Physical Science and Technology, Southwest University, Chongqing 400715, China
| |
Collapse
|
49
|
Senn F, Wiebke J, Schwerdtfeger P, Pahl E. Long-range contributions for the use of truncated pair potentials of molecular systems – application to nitrogen N 2. Mol Phys 2015. [DOI: 10.1080/00268976.2014.994568] [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]
Affiliation(s)
- F. Senn
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study and the Institute of Natural and Mathematical Sciences, Massey University Albany, Auckland, New Zealand
| | - J. Wiebke
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study and the Institute of Natural and Mathematical Sciences, Massey University Albany, Auckland, New Zealand
| | - P. Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study and the Institute of Natural and Mathematical Sciences, Massey University Albany, Auckland, New Zealand
| | - E. Pahl
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study and the Institute of Natural and Mathematical Sciences, Massey University Albany, Auckland, New Zealand
| |
Collapse
|
50
|
Waldrop JM, Song B, Patkowski K, Wang X. Accurate ab initio potential for the krypton dimer and transport properties of the low-density krypton gas. J Chem Phys 2015; 142:204307. [DOI: 10.1063/1.4921623] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jonathan M. Waldrop
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| | - 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, Shaanxi 710049, China
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| | - Xiaopo Wang
- Key Laboratory of Thermo-Fluid Science and Engineering, Ministry of Education, School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi 710049, China
| |
Collapse
|