1
|
Rusakov YY, Rusakova IL. Getaway from the Geometry Factor Error in the Molecular Property Calculations: Efficient pecG- n ( n = 1, 2) Basis Sets for the Geometry Optimization of Molecules Containing Light p Elements. J Chem Theory Comput 2024. [PMID: 39075034 DOI: 10.1021/acs.jctc.4c00772] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/31/2024]
Abstract
The basis of all molecular property quantum chemical calculations is the correct equilibrium geometry. In this paper, new efficient pecG-n (n = 1, 2) basis sets for the geometry optimization of molecules containing hydrogen and p elements of 2-3 periods are proposed. These basis sets were optimized via the property-energy consistent (PEC) algorithm directed to the minimization of the molecular energy gradient relative to the bond lengths. New basis sets are compact and give equilibrium geometries of very high quality, which is comparable to that provided by considerably larger energy-optimized basis sets. The equilibrium geometries obtained with the pecG-n (n = 1, 2) basis sets and the other basis sets of diverse quality were tested in the CCSD calculations of different second-order molecular properties, including NMR shielding constants, static polarizabilities, and static magnetizabilities. As a result, new basis sets have demonstrated far superior performance as compared to the other energy-optimized basis sets of the same or close sizes commonly used at the geometry optimization stage.
Collapse
Affiliation(s)
- Yuriy Yu Rusakov
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Irina L Rusakova
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| |
Collapse
|
2
|
Garberoglio G, Lissoni C, Spagnoli L, Harvey AH. Comprehensive quantum calculation of the first dielectric virial coefficient of water. J Chem Phys 2024; 160:024309. [PMID: 38214389 DOI: 10.1063/5.0187774] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2023] [Accepted: 12/20/2023] [Indexed: 01/13/2024] Open
Abstract
We present a complete calculation, fully accounting for quantum effects and for molecular flexibility, of the first dielectric virial coefficient of water and its isotopologues. The contribution of the electronic polarizability is computed from a state-of-the-art intramolecular potential and polarizability surface from the literature, and its small temperature dependence is quantified. The dipolar polarizability is calculated in a similar manner with an accurate literature dipole-moment surface; it differs from the classical result both due to the different molecular geometries sampled at different temperatures and due to the quantization of rotation. We calculate the dipolar contribution independently from spectroscopic information in the HITRAN2020 database and find that the two methods yield consistent results. The resulting first dielectric virial coefficient provides a complete description of the dielectric constant at low density that can be used in humidity metrology and as a boundary condition for new formulations for the static dielectric constant of water and heavy water.
Collapse
Affiliation(s)
- Giovanni Garberoglio
- European Centre for Theoretical Studies in Nuclear Physics and Related Areas (ECT*), Fondazione Bruno Kessler, Trento I-38123, Italy
| | - Chiara Lissoni
- Physics Department, University of Trento, Trento I-38123, Italy
| | - Luca Spagnoli
- Physics Department, University of Trento, Trento I-38123, Italy
| | - Allan H Harvey
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| |
Collapse
|
3
|
Begušić T, Blake GA. Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality. Nat Commun 2023; 14:1950. [PMID: 37029146 PMCID: PMC10082090 DOI: 10.1038/s41467-023-37667-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 03/22/2023] [Indexed: 04/09/2023] Open
Abstract
Two-dimensional spectroscopic techniques combining terahertz (THz), infrared (IR), and visible pulses offer a wealth of information about coupling among vibrational modes in molecular liquids, thus providing a promising probe of their local structure. However, the capabilities of these spectroscopies are still largely unexplored due to experimental limitations and inherently weak nonlinear signals. Here, through a combination of equilibrium-nonequilibrium molecular dynamics (MD) and a tailored spectrum decomposition scheme, we identify a relationship between the tetrahedral order of liquid water and its two-dimensional IR-IR-Raman (IIR) spectrum. The structure-spectrum relationship can explain the temperature dependence of the spectral features corresponding to the anharmonic coupling between low-frequency intermolecular and high-frequency intramolecular vibrational modes of water. In light of these results, we propose new experiments and discuss the implications for the study of tetrahedrality of liquid water.
Collapse
Affiliation(s)
- Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
| | - Geoffrey A Blake
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
- Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, 91125, USA.
| |
Collapse
|
4
|
Yu CC, Seki T, Chiang KY, Tang F, Sun S, Bonn M, Nagata Y. Polarization-Dependent Heterodyne-Detected Sum-Frequency Generation Spectroscopy as a Tool to Explore Surface Molecular Orientation and Ångström-Scale Depth Profiling. J Phys Chem B 2022; 126:6113-6124. [PMID: 35849538 PMCID: PMC9421650 DOI: 10.1021/acs.jpcb.2c02178] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Revised: 07/01/2022] [Indexed: 11/28/2022]
Abstract
Sum-frequency generation (SFG) spectroscopy provides a unique optical probe for interfacial molecules with interface-specificity and molecular specificity. SFG measurements can be further carried out at different polarization combinations, but the target of the polarization-dependent SFG is conventionally limited to investigating the molecular orientation. Here, we explore the possibility of polarization-dependent SFG (PD-SFG) measurements with heterodyne detection (HD-PD-SFG). We stress that HD-PD-SFG enables accurate determination of the peak amplitude, a key factor of the PD-SFG data. Subsequently, we outline that HD-PD-SFG can be used not only for estimating the molecular orientation but also for investigating the interfacial dielectric profile and studying the depth profile of molecules. We further illustrate the variety of combined simulation and PD-SFG studies.
Collapse
Affiliation(s)
- Chun-Chieh Yu
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Takakazu Seki
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kuo-Yang Chiang
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Fujie Tang
- Department
of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Shumei Sun
- Department
of Physics and Applied Optics Beijing Area Major Laboratory, Beijing Normal University, Beijing 100875, China
| | - Mischa Bonn
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Yuki Nagata
- Max
Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
5
|
Omodemi O, Kaledin M, Kaledin AL. Permutationally invariant polynomial representation of polarizability tensor surfaces for linear regression analysis. J Comput Chem 2022; 43:1495-1503. [PMID: 35737590 DOI: 10.1002/jcc.26952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2022] [Accepted: 06/09/2022] [Indexed: 11/07/2022]
Abstract
A linearly parameterized functional form for a Cartesian representation of molecular dipole polarizability tensor surfaces (PTS) is described. The proposed expression for the PTS is a linearization of the recently reported power series ansatz of the original Applequist model, which by construction is non-linear in parameter space. This new approach possesses (i) a unique solution to the least-squares fitting problem; (ii) a low level of the computational complexity of the resulting linear regression procedure, comparable to those of the potential energy and dipole moment surfaces; and (iii) a competitive level of accuracy compared to the non-linear PTS model. Calculations of CH4 PTS, with polarizabilities fitted to 9000 training set points with the energies up to 14,000 cm-1 show an impressive level of accuracy of the linear PTS model obtained with ~1600 parameters: ~1% versus 0.3% RMSE for the non-linear vs. linear model on a test set of 1000 configurations.
Collapse
Affiliation(s)
- Oluwaseun Omodemi
- Department of Chemistry & Biochemistry, Kennesaw State University, Kennesaw, Georgia, USA
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, Kennesaw, Georgia, USA
| | - Alexey L Kaledin
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia, USA
| |
Collapse
|
6
|
Carvalho JR, Vidal LN. Calculation of absolute Raman scattering cross-sections using vibrational self-consistent field/vibrational configuration interaction wave functions. J Comput Chem 2022; 43:1484-1494. [PMID: 35731622 DOI: 10.1002/jcc.26951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2022] [Revised: 05/21/2022] [Accepted: 06/02/2022] [Indexed: 11/12/2022]
Abstract
In the present study, the differential scattering cross-sections, depolarization ratios and Raman shifts of small molecular systems are obtained from configuration iteration wave functions of vibrational self-consistent field (VSCF) states. The transition polarizabilities were modeled using the Placzek approximation, neglecting those contributions not arising from the electric dipole mechanism. This theoretical approach is considered a good approximation for samples that absorb in the UV range if the excitation radiation falls in the visible region, as is the case of the molecules selected for the present study, namely: water, methane, and acetylene. Potential energy and electronic polarizability surfaces are calculated by the CCSD(T) and CC3 methods with aug-cc-p(C)V(T,Q,5)Z basis sets. The vibrational Hamiltonian includes the vibrational angular momentum contribution of the Watson kinetic energy operator. As expected, due to the variational nature of the VSCF and vibrational configuration interaction (VCI) methods, the Raman transition wavenumbers are substantially improved over the harmonic predictions. Surprisingly, the scattering cross-sections obtained using the harmonic approximation or the VSCF method better agrees with the experimental values than those cross-sections predicted using VCI wave functions. The more significant deviations of the VCI results from the experimental reference may be related to the significant uncertainties of the measured cross-sections. Still, it may also indicate that the VCI Raman transition moments may require a more accurate description of the electronic polarizability surface. Finally, the depolarization ratios calculated for H2 O and C2 D2 using harmonic and VCI wave functions have similar accuracy, whereas, for C2 H2 and C2 HD, the VCI results are more accurate.
Collapse
Affiliation(s)
- Jhonatas R Carvalho
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas, USA
| | - Luciano N Vidal
- Departamento Acadêmico de Química e Biologia, Universidade Tecnológica Federal do Paraná, Curitiba, Paraná, Brazil
| |
Collapse
|
7
|
Omodemi O, Sprouse S, Herbert D, Kaledin M, Kaledin AL. On the Cartesian Representation of the Molecular Polarizability Tensor Surface by Polynomial Fitting to Ab Initio Data. J Chem Theory Comput 2021; 18:37-45. [PMID: 34958587 DOI: 10.1021/acs.jctc.1c01015] [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/29/2022]
Abstract
We describe an approach to constructing an analytic Cartesian representation of the molecular dipole polarizability tensor surface in terms of polynomials in interatomic distances with a training set of ab initio data points obtained from a molecular dynamics (MD) simulation or by any other available means. The proposed formulation is based on a perturbation treatment of the unmodified point dipole polarizability model of Applequist [ J. Am. Chem. Soc. 1972, 94, 2952] and is shown here to be, by construction (i) free of short-range or other singularities or discontinuities, (ii) symmetric and translationally invariant, and (iii) nonreliant on a body-fixed coordinate system. Permutational invariance of like nuclei is demonstrated to be readily applicable, making this approach useful for highly fluxional and reactive systems. Derivation of the method is described in detail, adding brief didactic numerical examples of H2 and H2O and concluding with an MD simulation of the Raman spectrum of H5O2+ at 300 K with the polarizability tensor fitted to CCSD(T)/aug-cc-pVTZ data obtained using the HBB-4B potential [ J. Chem. Phys. 2005, 122, 044308].
Collapse
Affiliation(s)
- Oluwaseun Omodemi
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Sarah Sprouse
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Destyni Herbert
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Alexey L Kaledin
- Cherry L. Emerson Center for Scientific Computation, Emory University, 1515 Dickey Drive, Atlanta, Georgia 30322, United States
| |
Collapse
|
8
|
|
9
|
Lao KU, Jia J, Maitra R, DiStasio RA. On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions. J Chem Phys 2018; 149:204303. [PMID: 30501247 DOI: 10.1063/1.5051458] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we investigate how geometric changes influence the static dipole polarizability (α) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-ζ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this ab initio surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order C 6 dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H2O. Using first-order perturbation theory, we compute the zero-point vibrational correction to α at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic (α iso) and nearly 50% to the anisotropic (α aniso) polarizability values. In doing so, we find that α iso = 9.8307 bohr3, which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.
Collapse
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Junteng Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Rahul Maitra
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
10
|
Dey D, Henriksen NE. Non-resonant vibrational excitation of HOD and selective bond breaking. J Chem Phys 2018; 148:234307. [PMID: 29935499 DOI: 10.1063/1.5029548] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper reports a time-dependent quantum mechanical wave packet study for bond-selective excitation and dissociation of HOD into the H + OD and D + OH channels in the first absorption band. Prior to excitation, the HOD molecule is randomly oriented with respect to a linearly polarized laser field and accurate static dipole moment and polarizability surfaces are included in the interaction potential. Vibrational excitation is obtained with intense, non-resonant 800 nm few-cycle excitation using dynamic Stark effect/impulsive Raman scattering. Dissociation is accomplished by another ultrashort vacuum ultraviolet-laser excitation. A laser control scheme is designed with a train of simple, non-resonant laser pulses in order to enhance the selectivity between the fragmentation channels. The effect of the carrier-envelope-phase of the ultrashort laser pulses is also investigated.
Collapse
Affiliation(s)
- Diptesh Dey
- Department of Chemistry, Technical University of Denmark, Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Niels E Henriksen
- Department of Chemistry, Technical University of Denmark, Building 207, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
11
|
Loboda O, Millot C. Geometry-dependent atomic multipole models for the water molecule. J Chem Phys 2017; 147:161718. [DOI: 10.1063/1.4995569] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- O. Loboda
- Karl-Franzens Universität, Institut für Chemie, Heinrichstraße 28/IV, Graz A-8010, Austria
| | - C. Millot
- Université de Lorraine, CNRS, SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes BP 70239, Vandoeuvre-lès-Nancy Cedex F-54506, France
| |
Collapse
|
12
|
Loboda O, Ingrosso F, Ruiz-López MF, Reis H, Millot C. Dipole and quadrupole polarizabilities of the water molecule as a function of geometry. J Comput Chem 2016; 37:2125-32. [DOI: 10.1002/jcc.24431] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2016] [Revised: 05/24/2016] [Accepted: 05/27/2016] [Indexed: 11/11/2022]
Affiliation(s)
- Oleksandr Loboda
- SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239; Université de Lorraine, CNRS; Vandoeuvre-lès-Nancy F-54506 France
| | - Francesca Ingrosso
- SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239; Université de Lorraine, CNRS; Vandoeuvre-lès-Nancy F-54506 France
| | - Manuel F. Ruiz-López
- SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239; Université de Lorraine, CNRS; Vandoeuvre-lès-Nancy F-54506 France
| | - Heribert Reis
- Institute of Biology, Medicinal Chemistry and Biotechnology, National Hellenic Research Foundation; 48 Vas. Constantinou Avenue Athens 11635 Greece
| | - Claude Millot
- SRSMC UMR 7565, Faculté des Sciences et Technologies, Boulevard des Aiguillettes, BP 70239; Université de Lorraine, CNRS; Vandoeuvre-lès-Nancy F-54506 France
| |
Collapse
|
13
|
Loboda O, Ingrosso F, Ruiz-López MF, Szalewicz K, Millot C. Geometry-dependent distributed polarizability models for the water molecule. J Chem Phys 2016; 144:034304. [PMID: 26801031 DOI: 10.1063/1.4939519] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Geometry-dependent distributed polarizability models have been constructed by fits to ab initio calculations at the coupled cluster level of theory with up to noniterative triple excitations in an augmented triple-zeta quality basis set for the water molecule in the field of a point charge. The investigated models include (i) charge-flow polarizabilities between chemically bonded atoms, (ii) isotropic or anisotropic dipolar polarizabilities on oxygen atom or on all atoms, and (iii) combinations of models (i) and (ii). For each model, the polarizability parameters have been optimized to reproduce the induction energy of a water molecule polarized by a point charge successively occupying a grid of points surrounding the molecule. The quality of the models is ascertained by examining their ability to reproduce these induction energies as well as the molecular dipolar and quadrupolar polarizabilities. The geometry dependence of the distributed polarizability models has been explored by changing bond lengths and HOH angle to generate 125 molecular structures (reduced to 75 symmetry-unique ones). For each considered model, the distributed polarizability components have been fitted as a function of the geometry by a Taylor expansion in monomer coordinate displacements up to the sum of powers equal to 4.
Collapse
Affiliation(s)
- Oleksandr Loboda
- Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-les-Nancy F-54506, France
| | - Francesca Ingrosso
- Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-les-Nancy F-54506, France
| | - Manuel F Ruiz-López
- Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-les-Nancy F-54506, France
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Claude Millot
- Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-les-Nancy F-54506, France
| |
Collapse
|
14
|
Medders GR, Paesani F. Infrared and Raman Spectroscopy of Liquid Water through "First-Principles" Many-Body Molecular Dynamics. J Chem Theory Comput 2015; 11:1145-54. [PMID: 26579763 DOI: 10.1021/ct501131j] [Citation(s) in RCA: 163] [Impact Index Per Article: 18.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Vibrational spectroscopy is a powerful technique to probe the structure and dynamics of water. However, deriving an unambiguous molecular-level interpretation of the experimental spectral features remains a challenge due to the complexity of the underlying hydrogen-bonding network. In this contribution, we present an integrated theoretical and computational framework (named many-body molecular dynamics or MB-MD) that, by systematically removing uncertainties associated with existing approaches, enables a rigorous modeling of vibrational spectra of water from quantum dynamical simulations. Specifically, we extend approaches used to model the many-body expansion of interaction energies to develop many-body representations of the dipole moment and polarizability of water. The combination of these "first-principles" representations with centroid molecular dynamics simulations enables the simulation of infrared and Raman spectra of liquid water under ambient conditions that, without relying on any ad hoc parameters, are in good agreement with the corresponding experimental results. Importantly, since the many-body energy, dipole, and polarizability surfaces employed in the simulations are derived independently from accurate fits to correlated electronic structure data, MB-MD allows for a systematic analysis of the calculated spectra in terms of both electronic and dynamical contributions. The present analysis suggests that, while MB-MD correctly reproduces both the shifts and the shapes of the main spectroscopic features, an improved description of quantum dynamical effects possibly combined with a dissociable water potential may be necessary for a quantitative representation of the OH stretch band.
Collapse
Affiliation(s)
- Gregory R Medders
- Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92037, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego , La Jolla, California 92037, United States
| |
Collapse
|
15
|
|
16
|
Sanchéz-Lozano M, Mandado M, Pérez-Juste I, Hermida-Ramón JM. Theoretical vibrational Raman and surface-enhanced Raman scattering spectra of water interacting with silver clusters. Chemphyschem 2014; 15:4067-76. [PMID: 25263101 DOI: 10.1002/cphc.201402454] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2014] [Indexed: 11/11/2022]
Abstract
In this study, we analyzed the Raman spectrum of a water molecule adsorbed on a cluster of 20 silver atoms, and the plasmonic electromagnetic effect of the silver surface was also considered to give a theoretical prediction of the surface-enhanced Raman scattering spectrum. The calculations were performed at the density functional theory (DFT) level by using both frozen and unfrozen silver clusters. Two different models were used to consider the plasmonic enhancement; one of them was a modified classical (dipole) model and the other was the coupled perturbed Hartree-Fock method with excitation frequencies obtained from time-dependent DFT calculations and with proper detuning of these frequencies. The importance of small geometrical distortions of the silver surface in the orientation of the adsorbed water was shown. Moreover, it was shown how the symmetry of the transition dipole moment and the symmetry of the vibrational modes influence the Raman intensities of the SERS spectrum.
Collapse
Affiliation(s)
- Marta Sanchéz-Lozano
- Departamento de Química Física, Facultade de Química, Universidade de Vigo, Campus Lagoas Marcosende s/n, 36310 Vigo (Spain)
| | | | | | | |
Collapse
|
17
|
Babin V, Leforestier C, Paesani F. Development of a “First Principles” Water Potential with Flexible Monomers: Dimer Potential Energy Surface, VRT Spectrum, and Second Virial Coefficient. J Chem Theory Comput 2013; 9:5395-403. [DOI: 10.1021/ct400863t] [Citation(s) in RCA: 318] [Impact Index Per Article: 28.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Volodymyr Babin
- Department
of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Claude Leforestier
- ICG-CTMM
CC 15.01, Universite Montpellier
II , 34095 Montpellier, Cedex 05, France
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
18
|
Monten R, Hajgató B, Deleuze MS. Many-body calculations of molecular electric polarizabilities in asymptotically complete basis sets. Mol Phys 2011. [DOI: 10.1080/00268976.2011.579580] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
19
|
Elking DM, Perera L, Duke R, Darden T, Pedersen LG. A finite field method for calculating molecular polarizability tensors for arbitrary multipole rank. J Comput Chem 2011; 32:3283-95. [PMID: 21915883 DOI: 10.1002/jcc.21914] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2011] [Revised: 07/06/2011] [Accepted: 07/25/2011] [Indexed: 11/07/2022]
Abstract
A finite field method for calculating spherical tensor molecular polarizability tensors α(lm;l'm') = ∂Δ(lm)/∂ϕ(l'm')* by numerical derivatives of induced molecular multipole Δ(lm) with respect to gradients of electrostatic potential ϕ(l'm')* is described for arbitrary multipole ranks l and l'. Interconversion formulae for transforming multipole moments and polarizability tensors between spherical and traceless Cartesian tensor conventions are derived. As an example, molecular polarizability tensors up to the hexadecapole-hexadecapole level are calculated for water using the following ab initio methods: Hartree-Fock (HF), Becke three-parameter Lee-Yang-Parr exchange-correlation functional (B3LYP), Møller-Plesset perturbation theory up to second order (MP2), and Coupled Cluster theory with single and double excitations (CCSD). In addition, intermolecular electrostatic and polarization energies calculated by molecular multipoles and polarizability tensors are compared with ab initio reference values calculated by the Reduced Variation Space method for several randomly oriented small molecule dimers separated by a large distance. It is discussed how higher order molecular polarizability tensors can be used as a tool for testing and developing new polarization models for future force fields.
Collapse
Affiliation(s)
- Dennis M Elking
- University of North Carolina, Department of Chemistry, Chapel Hill, North Carolina 27599, USA
| | | | | | | | | |
Collapse
|
20
|
Hammond JR, Govind N, Kowalski K, Autschbach J, Xantheas SS. Accurate dipole polarizabilities for water clusters n=2–12 at the coupled-cluster level of theory and benchmarking of various density functionals. J Chem Phys 2009; 131:214103. [DOI: 10.1063/1.3263604] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
|
21
|
Páleníková J, Kraus M, Neogrády P, Kellö V, Urban M. Theoretical study of molecular properties of low-lying electronic excited states of H2O and H2S. Mol Phys 2008. [DOI: 10.1080/00268970802454786] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Jana Páleníková
- a Faculty of Natural Sciences, Department of Physical and Theoretical Chemistry , Comenius University , Bratislava, Slovakia
| | - Michal Kraus
- a Faculty of Natural Sciences, Department of Physical and Theoretical Chemistry , Comenius University , Bratislava, Slovakia
| | - Pavel Neogrády
- a Faculty of Natural Sciences, Department of Physical and Theoretical Chemistry , Comenius University , Bratislava, Slovakia
| | - Vladimir Kellö
- a Faculty of Natural Sciences, Department of Physical and Theoretical Chemistry , Comenius University , Bratislava, Slovakia
| | - Miroslav Urban
- a Faculty of Natural Sciences, Department of Physical and Theoretical Chemistry , Comenius University , Bratislava, Slovakia
| |
Collapse
|
22
|
Burnham CJ, Anick DJ, Mankoo PK, Reiter GF. The vibrational proton potential in bulk liquid water and ice. J Chem Phys 2008; 128:154519. [DOI: 10.1063/1.2895750] [Citation(s) in RCA: 108] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Bukowski R, Szalewicz K, Groenenboom GC, van der Avoird A. Polarizable interaction potential for water from coupled cluster calculations. I. Analysis of dimer potential energy surface. J Chem Phys 2008; 128:094313. [DOI: 10.1063/1.2832746] [Citation(s) in RCA: 82] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Wu DY, Duan S, Liu XM, Xu YC, Jiang YX, Ren B, Xu X, Lin SH, Tian ZQ. Theoretical study of binding interactions and vibrational Raman spectra of water in hydrogen-bonded anionic complexes: (H2O)n- (n = 2 and 3), H2O...X- (X = F, Cl, Br, and I), and H2O...M- (M = Cu, Ag, and Au). J Phys Chem A 2008; 112:1313-21. [PMID: 18215023 DOI: 10.1021/jp0722105] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Binding interactions and Raman spectra of water in hydrogen-bonded anionic complexes have been studied by using the hybrid density functional theory method (B3LYP) and ab initio (MP2) method. In order to explore the influence of hydrogen bond interactions and the anionic effect on the Raman intensities of water, model complexes, such as the negatively charged water clusters ((H2O)n-, n = 2 and 3), the water...halide anions (H2O...X-, X = F, Cl, Br, and I), and the water-metal atom anionic complexes (H2O...M-, M = Cu, Ag, and Au), have been employed in the present calculations. These model complexes contained different types of hydrogen bonds, such as O-H...X-, O-H...M-, O-H...O, and O-H...e-. In particular, the last one is a dipole-bound electron involved in the anionic water clusters. Our results showed that there exists a large enhancement in the off-resonance Raman intensities of both the H-O-H bending mode and the hydrogen-bonded O-H stretching mode, and the enhancement factor is more significant for the former than for the latter. The reasons for these spectral properties can be attributed to the strong polarization effect of the proton acceptors (X-, M-, O, and e-) in these hydrogen-bonded complexes. We proposed that the strong Raman signal of the H-O-H bending mode may be used as a fingerprint to address the local microstructures of water molecules in the chemical and biological systems.
Collapse
Affiliation(s)
- De-Yin Wu
- Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, Fujian, People Republic of China.
| | | | | | | | | | | | | | | | | |
Collapse
|
25
|
|
26
|
Elola MD, Ladanyi BM. Intermolecular polarizability dynamics of aqueous formamide liquid mixtures studied by molecular dynamics simulations. J Chem Phys 2007; 126:084504. [PMID: 17343455 DOI: 10.1063/1.2446782] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
A molecular dynamics simulation study is presented for the relaxation of the polarizability anisotropy in liquid mixtures of formamide and water, using a dipolar induction scheme that involves the intrinsic polarizability and first hyperpolarizability tensors of the molecules, and the dipole-quadrupole polarizability of water species. The long time diffusive decay of the collective polarizability anisotropy correlations exhibits a substantial slowing down as the formamide mole fraction increases in the mixture. The diffusive times for the polarizability relaxation obtained from the authors' simulations are in good agreement with optical Kerr effect experimental data, and they are found to correlate nearly linearly with the estimated mean lifetimes of the hydrogen bonds within the mixture, suggesting that the relaxation of the hydrogen bond network is responsible to some extent for the collective relaxation of the polarizability anisotropy of the mixture. The short time behavior of the polarizability anisotropy relaxation was investigated by computing the nuclear response function, R(t), which is very rapidly dominated by the formamide contribution as it is added to water, due to the much larger polarizability anisotropy of formamide molecules compared to that of water. Several contributions to the Raman spectrum were also analyzed as a function of composition, and the dynamical origin of the different bands was determined.
Collapse
Affiliation(s)
- M Dolores Elola
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, USA.
| | | |
Collapse
|