1
|
Joshi P, Voora VK. Generalized perturbative singles corrections to the random phase approximation method: Impact on noncovalent interaction energies of closed- and open-shell dimers. J Chem Phys 2024; 160:044104. [PMID: 38258929 DOI: 10.1063/5.0180526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Accepted: 12/29/2023] [Indexed: 01/24/2024] Open
Abstract
The post-Kohn-Sham (KS) random phase approximation (RPA) method may provide a poor description of interaction energies of weakly bonded molecules due to inherent density errors in approximate KS functionals. To overcome these errors, we develop a generalized formalism to incorporate perturbative singles (pS) corrections to the RPA method using orbital rotations as a perturbation parameter. The pS schemes differ in the choice of orbital-rotation gradient and Hessian. We propose a pS scheme termed RPA singles (RPAS)[Hartree-Fock (HF)] that uses the RPA orbital-rotation gradient and time-dependent HF Hessian. This correction reduces the errors in noncovalent interaction energies of closed- and open-shell dimers. For the open-shell dimers, the RPAS(HF) method leads to a consistent error reduction by 50% or more compared to the RPA method for the cases of hydrogen-bonding, metal-solvent, carbene-solvent, and dispersion interactions. We also find that the pS corrections are more important in error reduction compared to higher-order exchange corrections to the RPA method. Overall, for open shells, the RPAS(HF)-corrected RPA method provides chemical accuracy for noncovalent interactions and is more reliable than other perturbative schemes and dispersion-corrected density functional approximations, highlighting its importance as a reliable beyond-RPA correction.
Collapse
Affiliation(s)
- Pulkit Joshi
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| |
Collapse
|
2
|
Belleflamme F, Hutter J. Radicals in aqueous solution: assessment of density-corrected SCAN functional. Phys Chem Chem Phys 2023; 25:20817-20836. [PMID: 37497572 DOI: 10.1039/d3cp02517a] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
We study self-interaction effects in solvated and strongly-correlated cationic molecular clusters, with a focus on the solvated hydroxyl radical. To address the self-interaction issue, we apply the DC-r2SCAN method, with the auxiliary density matrix approach. Validating our method through simulations of bulk liquid water, we demonstrate that DC-r2SCAN maintains the structural accuracy of r2SCAN while effectively addressing spin density localization issues. Extending our analysis to solvated cationic molecular clusters, we find that the hemibonded motif in the [CH3S∴CH3SH]+ cluster is disrupted in the DC-r2SCAN simulation, in contrast to r2SCAN that preserves the (three-electron-two-center)-bonded motif. Similarly, for the [SH∴SH2]+ cluster, r2SCAN restores the hemibonded motif through spin leakage, while DC-r2SCAN predicts a weaker hemibond formation influenced by solvent-solute interactions. Our findings demonstrate the potential of DC-r2SCAN combined with the auxiliary density matrix method to improve electronic structure calculations, providing insights into the properties of solvated cationic molecular clusters. This work contributes to the advancement of self-interaction corrected electronic structure theory and offers a computational framework for modeling condensed phase systems with intricate correlation effects.
Collapse
Affiliation(s)
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Zurich, Switzerland.
| |
Collapse
|
3
|
Kato T, Fujii A. Infrared Spectroscopy of (Benzene-H 2S-X n) +, X = H 2O ( n = 1 and 2) and CH 3OH ( n = 1), Radical Cation Clusters: Microsolvation Effects on the S-π Hemibond. J Phys Chem A 2023; 127:742-750. [PMID: 36636986 DOI: 10.1021/acs.jpca.2c08324] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
An unconventional covalent bond in which three electrons are shared by two centers is called hemibond. Hemibond formation frequently competes with proton transfer (or ionic hydrogen bond formation), but there have been a few experimental reports on such competition. In the present study, we focus on the (benzene-H2S)+ radical cation cluster, which is a model system of the S-π hemibond. The stability of the S-π hemibond to the microsolvation by water and methanol is explored with infrared spectroscopy of (benzene-H2S-Xn)+, X = H2O (n = 1 and 2) and CH3OH (n = 1), clusters. We also perform energy-optimization and vibrational simulations of (benzene-H2S-Xn)+. By comparison among the observed and simulated spectra, we determine the intermolecular binding motifs in (benzene-H2S-Xn)+. While the S-π hemibonded isomer is exclusively populated in (benzene-H2S-H2O)+, both the hemibonded and proton-transferred isomers coexist in [benzene-H2S-(H2O)2]+ and (benzene-H2S-CH3OH)+. Breaking of the S-π hemibond by the microsolvation is observed, and its solvent and cluster size dependence is interpreted by the proton affinity and the coordination property of the solvent moiety.
Collapse
Affiliation(s)
- Takeru Kato
- Department of Chemistry, Graduate School of Science, Tohoku University Sendai 980-8578, Japan
| | - Asuka Fujii
- Department of Chemistry, Graduate School of Science, Tohoku University Sendai 980-8578, Japan
| |
Collapse
|
4
|
Masaya TW, Goulay F. A Molecular Dynamic Study of the Effects of Surface Partitioning on the OH Radical Interactions with Solutes in Multicomponent Aqueous Aerosols. J Phys Chem A 2023; 127:751-764. [PMID: 36639126 DOI: 10.1021/acs.jpca.2c07419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
The surface-bulk partitioning of small saccharide and amide molecules in aqueous droplets was investigated using molecular dynamics. The air-particle interface was modeled using a 80 Å cubic water box containing a series of organic molecules and surrounded by gaseous OH radicals. The properties of the organic solutes within the interface and the water bulk were examined at a molecular level using density profiles and radial pair distribution functions. Molecules containing only polar functional groups such as urea and glucose are found predominantly in the water bulk, forming an exclusion layer near the water surface. Substitution of a single polar group by an alkyl group in sugars and amides leads to the migration of the molecule toward the interface. Within the first 2 nm from the water surface, surface-active solutes lose their rotational freedom and adopt a preferred orientation with the alkyl group pointing toward the surface. The different packing within the interface leads to different solvation shell structures and enhanced interaction between the organic molecules and absorbed OH radicals. The simulations provide quantitative information about the dimension, composition, and organization of the air-water interface as well as about the nonreactive interaction of the OH radicals with the organic solutes. It suggests that increased concentrations, preferred orientations, and decreased solvation near the air-water surface may lead to differences in reactivities between surface-active and surface-inactive molecules. The results are important to explain how heterogeneous oxidation mechanisms and kinetics within interfaces may differ from those of the bulk.
Collapse
Affiliation(s)
- Tadini Wenyika Masaya
- C. Eugene Bennett Department of Chemistry, West Virginia University, Morgantown, West Virginia26506, United States
| | - Fabien Goulay
- C. Eugene Bennett Department of Chemistry, West Virginia University, Morgantown, West Virginia26506, United States
| |
Collapse
|
5
|
Heindel JP, Hao H, LaCour RA, Head-Gordon T. Spontaneous Formation of Hydrogen Peroxide in Water Microdroplets. J Phys Chem Lett 2022; 13:10035-10041. [PMID: 36264238 DOI: 10.1021/acs.jpclett.2c01721] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
There is accumulating evidence that many chemical reactions are accelerated by several orders of magnitude in micrometer-sized aqueous or organic liquid droplets compared to their corresponding bulk liquid phase. However, the molecular origin of the enhanced rates remains unclear as in the case of spontaneous appearance of 1 μM hydrogen peroxide in water microdroplets. In this Letter, we consider the range of ionization energies and whether interfacial electric fields of a microdroplet can feasibly overcome the high energy step from hydroxide ions (OH-) to hydroxyl radicals (OH•) in a primary H2O2 mechanism. We find that the vertical ionization energies (VIEs) of partially solvated OH- ions are greatly lowered relative to the average VIE in the bulk liquid, unlike the case of the Cl- anion which shows no reduction in the VIEs regardless of solvation environment. Overall reduced hydrogen-bonding and undercoordination of OH- are structural features that are more readily present at the air-water interface, where the energy scale for ionization can be matched by statistically probable electric field values.
Collapse
Affiliation(s)
- Joseph P Heindel
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Hongxia Hao
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California94720, United States
| | - R Allen LaCour
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
- Departments of Bioengineering and Chemical and Biomolecular EngineeringUniversity of California, Berkeley, California94720, United States
| |
Collapse
|
6
|
Sun X, Xie M, Qiu W, Wei C, Chen X, Hu Y. Spectroscopic evidence of S∴N and S∴O hemibonds in heterodimer cations. Phys Chem Chem Phys 2022; 24:19354-19361. [PMID: 35686608 DOI: 10.1039/d2cp00904h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Computational and condensed phase experimental evidence for the existence of S∴N and S∴O hemibonded structures has been reported previously, but no gas phase experimental evidence has been reported. To experimentally explore the existence of the S∴N and S∴O hemibonds in the gas phase, we recorded the infrared photodissociation action spectra of four cationic clusters: [CH3SH-NH3]+, [CH3SCH3-NH3]+, [CH3SCH3-H2O]+, and [CH3OCH3-H2O]+. Combined with the calculation results, it is found that the S∴N hemibonded structure is competitive with the S⋯HN H-bonded structure, though only the latter structure is actually observed in [CH3SH-NH3]+. The spectral and theoretical results show that hemibonds can form between the second- (oxygen or nitrogen) and the third-period elements (sulfur) in the heterodimer clusters of [CH3SCH3-NH3]+ and [CH3SCH3-H2O]+. However, the S∴N and S∴O hemibonded structures are found competitive with the C⋯HN and CH⋯O H-bonded structures, respectively, and both the structures coexist. On the other hand, the O∴O hemibonded structure is much less stable than other hydrogen bonded (H-bonded) structures in [CH3OCH3-H2O]+, and it shows no clear contribution to the observed spectrum. This study provides direct spectroscopic evidence for the existence of S∴N and S∴O hemibonds in the gas phase and their competition with the H-bonds, which may be also fundamentally important in biological processes.
Collapse
Affiliation(s)
- Xiaonan Sun
- MOE Key Laboratory of Laser Life Science & Guangdong Provincial Key Laboratory of Laser Life Science, Guangzhou Key Laboratory of Spectral Analysis and Functional Probes, College of Biophotonics, South China Normal University, Guangzhou, 510631, China.
| | - Min Xie
- MOE Key Laboratory of Laser Life Science & Guangdong Provincial Key Laboratory of Laser Life Science, Guangzhou Key Laboratory of Spectral Analysis and Functional Probes, College of Biophotonics, South China Normal University, Guangzhou, 510631, China.
| | - Wei Qiu
- MOE Key Laboratory of Laser Life Science & Guangdong Provincial Key Laboratory of Laser Life Science, Guangzhou Key Laboratory of Spectral Analysis and Functional Probes, College of Biophotonics, South China Normal University, Guangzhou, 510631, China.
| | - Chengcheng Wei
- MOE Key Laboratory of Laser Life Science & Guangdong Provincial Key Laboratory of Laser Life Science, Guangzhou Key Laboratory of Spectral Analysis and Functional Probes, College of Biophotonics, South China Normal University, Guangzhou, 510631, China.
| | - Xujian Chen
- MOE Key Laboratory of Laser Life Science & Guangdong Provincial Key Laboratory of Laser Life Science, Guangzhou Key Laboratory of Spectral Analysis and Functional Probes, College of Biophotonics, South China Normal University, Guangzhou, 510631, China.
| | - Yongjun Hu
- MOE Key Laboratory of Laser Life Science & Guangdong Provincial Key Laboratory of Laser Life Science, Guangzhou Key Laboratory of Spectral Analysis and Functional Probes, College of Biophotonics, South China Normal University, Guangzhou, 510631, China.
| |
Collapse
|
7
|
Rana B, Coons MP, Herbert JM. Detection and Correction of Delocalization Errors for Electron and Hole Polarons Using Density-Corrected DFT. J Phys Chem Lett 2022; 13:5275-5284. [PMID: 35674719 DOI: 10.1021/acs.jpclett.2c01187] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Modeling polaron defects is an important aspect of computational materials science, but the description of unpaired spins in density functional theory (DFT) often suffers from delocalization error. To diagnose and correct the overdelocalization of spin defects, we report an implementation of density-corrected (DC-)DFT and its analytic energy gradient. In DC-DFT, an exchange-correlation functional is evaluated using a Hartree-Fock density, thus incorporating electron correlation while avoiding self-interaction error. Results for an electron polaron in models of titania and a hole polaron in Al-doped silica demonstrate that geometry optimization with semilocal functionals drives significant structural distortion, including the elongation of several bonds, such that subsequent single-point calculations with hybrid functionals fail to afford a localized defect even in cases where geometry optimization with the hybrid functional does localize the polaron. This has significant implications for traditional workflows in computational materials science, where semilocal functionals are often used for structure relaxation. DC-DFT calculations provide a mechanism to detect situations where delocalization error is likely to affect the results.
Collapse
Affiliation(s)
- Bhaskar Rana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Marc P Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|