1
|
Tu Quyen LT, Tung BN, Thach PN, Tri NN, Trung NT. Characteristics of nonconventional hydrogen bonds and stability of dimers of chalcogenoaldehyde derivatives: a noticeable role of oxygen compared to other chalcogens. RSC Adv 2024; 14:14114-14125. [PMID: 38686288 PMCID: PMC11057360 DOI: 10.1039/d4ra01837k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2024] [Accepted: 04/17/2024] [Indexed: 05/02/2024] Open
Abstract
In this work, twenty-four stable dimers of RCHZ with R = H, F, Cl, Br, CH3 or NH2 and Z = O, S, Se or Te were determined. It was found that the stability of most dimers is primarily contributed by the electrostatic force, except for the dominant role of the induction term in those involving a Te atom, which has been rarely observed. Both electron-donating and -withdrawing groups in substituted formaldehyde cause an increase in the strength of nonconventional Csp2-H⋯Z hydrogen bonds, as well as the dimers, in which the electron donating effect plays a more crucial role. The strength of nonconventional hydrogen bonds decreases in the following order: Csp2-H⋯O ≫ Csp2-H⋯S > Csp2-H⋯Se > Csp2-H⋯Te. Remarkably, a highly significant role of the O atom compared to S, Se and Te in increasing the Csp2-H stretching frequency and strength of the nonconventional hydrogen bonds and dimers is found. A Csp2-H stretching frequency red-shift is observed in Csp2-H⋯S/Se/Te, while a blue-shift is obtained in Csp2-H⋯O. When Z changes from O to S to Se and to Te, the Csp2-H blue-shift tends to decrease and eventually turns to a red-shift, in agreement with the increasing order of the proton affinity at Z in the isolated monomer. The magnitude of the Csp2-H stretching frequency red-shift is larger for Csp2-H⋯Te than Csp2-H⋯S/Se, consistent with the rising trend of proton affinity at the Z site and the polarity of the Csp2-H bond in the substituted chalcogenoaldehydes. The Csp2-H blue-shifting of the Csp2-H⋯O hydrogen bonds is observed in all dimers regardless of the electron effect of the substituents. Following complexation, the electron-donating derivatives exhibit a stronger Csp2-H blue-shift compared to the electron-withdrawing ones. Notably, the stronger Csp2-H blue-shift turns out to involve a less polarized Csp2-H bond and a decrease in the occupation at the σ*(Csp2-H) antibonding orbital in the isolated monomer.
Collapse
Affiliation(s)
- Le Thi Tu Quyen
- Laboratory of Computational Chemistry and Modelling (LCCM), Quy Nhon University 170 An Duong Vuong Street Quy Nhon City 590000 Vietnam
| | - Bui Nhat Tung
- Laboratory of Computational Chemistry and Modelling (LCCM), Quy Nhon University 170 An Duong Vuong Street Quy Nhon City 590000 Vietnam
| | - Pham Ngoc Thach
- Laboratory of Computational Chemistry and Modelling (LCCM), Quy Nhon University 170 An Duong Vuong Street Quy Nhon City 590000 Vietnam
- Faculty of Natural Sciences, Quy Nhon University 170 An Duong Vuong Street Quy Nhon City 590000 Vietnam
| | - Nguyen Ngoc Tri
- Laboratory of Computational Chemistry and Modelling (LCCM), Quy Nhon University 170 An Duong Vuong Street Quy Nhon City 590000 Vietnam
- Faculty of Natural Sciences, Quy Nhon University 170 An Duong Vuong Street Quy Nhon City 590000 Vietnam
| | - Nguyen Tien Trung
- Laboratory of Computational Chemistry and Modelling (LCCM), Quy Nhon University 170 An Duong Vuong Street Quy Nhon City 590000 Vietnam
- Faculty of Natural Sciences, Quy Nhon University 170 An Duong Vuong Street Quy Nhon City 590000 Vietnam
| |
Collapse
|
2
|
de Azevedo Santos L, van der Voort S, Burema SR, Fonseca Guerra C, Bickelhaupt FM. Blueshift in Trifurcated Hydrogen Bonds: A Tradeoff between Tetrel Bonding and Steric Repulsion. Chemphyschem 2024; 25:e202300480. [PMID: 37864778 DOI: 10.1002/cphc.202300480] [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: 08/04/2023] [Revised: 10/17/2023] [Accepted: 10/20/2023] [Indexed: 10/23/2023]
Abstract
We have quantum chemically investigated the origin of the atypical blueshift of the H-C bond stretching frequency in the hydrogen-bonded complex X- •••H3 C-Y (X, Y=F, Cl, Br, I), as compared to the corresponding redshift occurring in Cl- •••H3 N and Cl- •••H3 C-H, using relativistic density functional theory (DFT) at ZORA-BLYP-D3(BJ)/QZ4P. Previously, this blueshift was attributed, among others, to the contraction of the H-C bonds as the H3 C moiety becomes less pyramidal. Herein, we provide quantitative evidence that, instead, the blueshift arises from a direct and strong X- •••C interaction of the HOMO of A- with the backside lobe on carbon of the low-lying C-Y antibonding σ* LUMO of the H3 C-Y fragment. This X- •••C bond, in essence a tetrel bond, pushes the H atoms towards a shorter H-C distance and makes the H3 C moiety more planar. The blueshift may, therefore, serve as a diagnostic for tetrel bonding.
Collapse
Affiliation(s)
- Lucas de Azevedo Santos
- Department of Chemistry and Pharmaceutical Sciences, AIMMS, Vrije Universiteit, Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Storm van der Voort
- Department of Chemistry and Pharmaceutical Sciences, AIMMS, Vrije Universiteit, Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Shiri R Burema
- Department of Chemistry and Pharmaceutical Sciences, AIMMS, Vrije Universiteit, Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Célia Fonseca Guerra
- Department of Chemistry and Pharmaceutical Sciences, AIMMS, Vrije Universiteit, Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - F Matthias Bickelhaupt
- Department of Chemistry and Pharmaceutical Sciences, AIMMS, Vrije Universiteit, Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
- Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ, Nijmegen, The Netherlands
- Department of Chemical Sciences, University of Johannesburg Auckland Park, Johannesburg, 2006, South Africa
| |
Collapse
|
3
|
Czernek J, Brus J, Czerneková V, Kobera L. Quantifying the Intrinsic Strength of C-H⋯O Intermolecular Interactions. Molecules 2023; 28:molecules28114478. [PMID: 37298953 DOI: 10.3390/molecules28114478] [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: 04/26/2023] [Revised: 05/27/2023] [Accepted: 05/29/2023] [Indexed: 06/12/2023] Open
Abstract
It has been recognized that the C-H⋯O structural motif can be present in destabilizing as well as highly stabilizing intermolecular environments. Thus, it should be of interest to describe the strength of the C-H⋯O hydrogen bond for constant structural factors so that this intrinsic strength can be quantified and compared to other types of interactions. This description is provided here for C2h-symmetric dimers of acrylic acid by means of the calculations that employ the coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] together with an extrapolation to the complete basis set (CBS) limit. Dimers featuring the C-H⋯O and O-H⋯O hydrogens bonds are carefully investigated in a wide range of intermolecular separations by the CCSD(T)/CBS approach, and also by the symmetry-adapted perturbation theory (SAPT) method, which is based on the density-functional theory (DFT) treatment of monomers. While the nature of these two types of hydrogen bonding is very similar according to the SAPT-DFT/CBS calculations and on the basis of a comparison of the intermolecular potential curves, the intrinsic strength of the C-H⋯O interaction is found to be about a quarter of its O-H⋯O counterpart that is less than one might anticipate.
Collapse
Affiliation(s)
- Jiří Czernek
- Institute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovsky Square 2, 162 00 Prague, Czech Republic
| | - Jiří Brus
- Institute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovsky Square 2, 162 00 Prague, Czech Republic
| | - Vladimíra Czerneková
- Institute of Physics, Czech Academy of Science, Na Slovance 2, 182 21 Prague, Czech Republic
| | - Libor Kobera
- Institute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovsky Square 2, 162 00 Prague, Czech Republic
| |
Collapse
|
4
|
Aldossary A, Gimferrer M, Mao Y, Hao H, Das AK, Salvador P, Head-Gordon T, Head-Gordon M. Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions. J Phys Chem A 2023; 127:1760-1774. [PMID: 36753558 DOI: 10.1021/acs.jpca.2c08061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.
Collapse
Affiliation(s)
- Abdulrahman Aldossary
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| | - Martí Gimferrer
- Institut de Química Computacional i Catàlsi and Departament de Química, Universitat de Girona, 17003 Girona, Catalonia Spain
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182, United States
| | - Hongxia Hao
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| | - Akshaya K Das
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| | - Pedro Salvador
- Institut de Química Computacional i Catàlsi and Departament de Química, Universitat de Girona, 17003 Girona, Catalonia Spain
| | - Teresa Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| | - Martin Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley California 94720, United States
| |
Collapse
|
5
|
Michalczyk M, Kizior B, Zierkiewicz W, Scheiner S. Factors contributing to halogen bond strength and stretch or contraction of internal covalent bond. Phys Chem Chem Phys 2023; 25:2907-2915. [PMID: 36636920 DOI: 10.1039/d2cp05598h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The halogen bond formed by a series of Lewis acids TF3X (T = C, Si, Ge, Sn, Pb; X = Cl, Br, I) with NH3 is studied by quantum chemical calculations. The interaction energy is closely mimicked by the depth of the σ-hole on the X atom as well as the full electrostatic energy. There is a first trend by which the hole is deepened if the T atom to which X is attached becomes more electron-withdrawing: C > Si > Ge > Sn > Pb. On the other hand, larger more polarizable T atoms are better able to transmit the electron-withdrawing power of the F substituents. The combination of these two opposing factors leaves PbF3X forming the strongest XBs, followed by CF3X, with SiF3X engaging in the weakest bonds. The charge transfer from the NH3 lone pair into the σ*(TX) antibonding orbital tends to elongate the covalent TX bond, and this force is largest for the heavier X and T atoms. On the other hand, the contraction of this bond deepens the σ-hole at the X atom, which would enhance both the electrostatic component and the full interaction energy. This bond-shortening effect is greatest for the lighter X atoms. The combination of these two opposing forces leaves the T-X bond contracting for X = Cl and Br, but lengthening for I.
Collapse
Affiliation(s)
- Mariusz Michalczyk
- Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370, Wrocław, Poland.
| | - Beata Kizior
- Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370, Wrocław, Poland.
| | - Wiktor Zierkiewicz
- Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370, Wrocław, Poland.
| | - Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University Logan, Utah, 84322-0300, USA.
| |
Collapse
|
6
|
Liu R, Wang R, Li D, Zhu Y, Yang X, Wang Z. An ab initio study on boundaries for characterizing cooperative effect of hydrogen bonds by intermolecular compression. CHINESE CHEM LETT 2022. [DOI: 10.1016/j.cclet.2022.107857] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
7
|
Robinson VN, Ghosh R, Egan CK, Riera M, Knight C, Paesani F, Hassanali A. The behavior of methane-water mixtures under elevated pressures from simulations using many-body potentials. J Chem Phys 2022; 156:194504. [PMID: 35597630 DOI: 10.1063/5.0089773] [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/15/2022] Open
Abstract
Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane-water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane-water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane-water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. We find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.
Collapse
Affiliation(s)
- Victor Naden Robinson
- The 'Abdus Salam' International Centre for Theoretical Physics, I-34151 Trieste, Italy
| | - Raja Ghosh
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Colin K Egan
- The 'Abdus Salam' International Centre for Theoretical Physics, I-34151 Trieste, Italy
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Computational Science Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Ali Hassanali
- The 'Abdus Salam' International Centre for Theoretical Physics, I-34151 Trieste, Italy
| |
Collapse
|
8
|
Dutta J, Sahu AK, Bhadauria AS, Biswal HS. Carbon-Centered Hydrogen Bonds in Proteins. J Chem Inf Model 2022; 62:1998-2008. [PMID: 35293733 DOI: 10.1021/acs.jcim.2c00015] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Hydrogen bonding (H-bonding) without lone pair(s) of electrons and π-electrons is a concept developed 2-3 years ago. H-bonds involving less electronegative tetrahedral carbon are beyond the classical concept of H-bonds. Herein, we present the first report on H-bonds with tetravalent carbons in proteins. A special bonding arrangement is needed to increase the negative charge density around the sp3-hybridized carbon atom. Therefore, less electronegative elements such as As and Mg, when bonded to sp3-C, enable the C-atoms as H-bond acceptors. Careful protein structure analysis aided by several quantum chemical calculations suggests that these H-bonds are weak to moderate in strength. We developed an empirical equation to estimate the C-H···C H-bond energy in proteins from the distances between the C- and H-atoms. In proteins, the binding energies range from -5.4 to -14.0 kJ/mol. The C-H···C H-bonds assist the substrate binding in proteins. We also explored the potential role of these carbon-centered H-bonds in C-H bond activation through σ-bond metathesis. To our surprise, contribution from these H-bonds is almost of similar magnitude as that from C-H···π H-bonds for C-H bond activation.
Collapse
Affiliation(s)
- Juhi Dutta
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), Bhimpur-Padanpur, Via-Jatni, Khurda, Bhubaneswar 752050, India.,Training School Complex, Homi Bhabha National Institute, Anushakti Nagar, Mumbai 400094, India
| | - Akshay Kumar Sahu
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), Bhimpur-Padanpur, Via-Jatni, Khurda, Bhubaneswar 752050, India.,Training School Complex, Homi Bhabha National Institute, Anushakti Nagar, Mumbai 400094, India
| | - Abhijeet S Bhadauria
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), Bhimpur-Padanpur, Via-Jatni, Khurda, Bhubaneswar 752050, India.,Training School Complex, Homi Bhabha National Institute, Anushakti Nagar, Mumbai 400094, India
| | - Himansu S Biswal
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), Bhimpur-Padanpur, Via-Jatni, Khurda, Bhubaneswar 752050, India.,Training School Complex, Homi Bhabha National Institute, Anushakti Nagar, Mumbai 400094, India
| |
Collapse
|
9
|
Cuc NTT, An NT, Ngan VT, Chandra AK, Trung NT. Importance of water and intramolecular interaction governs substantial blue shift of Csp2–H stretching frequency in complexes between chalcogenoaldehydes and water. RSC Adv 2022; 12:1998-2008. [PMID: 35425273 PMCID: PMC8979115 DOI: 10.1039/d1ra07444j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2021] [Accepted: 01/04/2022] [Indexed: 11/25/2022] Open
Abstract
Geometrical structure, stability and cooperativity, and contribution of hydrogen bonds to the stability of complexes between chalcogenoaldehydes and water were thoroughly investigated using quantum chemical methods. The stability of the complexes increases significantly when one or more H2O molecules are added to the binary system, whereas it decreases sharply going from O to S, Se, or Te substitution. The O–H⋯O H-bond is twice as stable as Csp2–H⋯O and O–H⋯S/Se/Te H-bonds. It is found that a considerable blue-shift of Csp2–H stretching frequency in the Csp2–H⋯O H-bond is mainly determined by an addition of water into the complexes along with the low polarity of the Csp2–H covalent bond in formaldehyde and acetaldehyde. The Csp2–H stretching frequency shift as a function of net second hyperconjugative energy for the σ*(Csp2–H) antibonding orbital is observed. Remarkably, a considerable Csp2–H blue shift of 109 cm−1 has been reported for the first time. Upon the addition of H2O into the binary systems, halogenated complexes witness a decreasing magnitude of the Csp2–H stretching frequency blue-shift in the Csp2–H⋯O H-bond, whereas CH3-substituted complexes experience the opposite trend. The considerable blue shift of Csp2–H stretching frequency.![]()
Collapse
Affiliation(s)
- Nguyen Thi Thanh Cuc
- Laboratory of Computational Chemistry and Modelling (LCCM), Department of Chemistry, Faculty of Natural Sciences, Quy Nhon University, Vietnam
| | - Nguyen Truong An
- Laboratory of Computational Chemistry and Modelling (LCCM), Department of Chemistry, Faculty of Natural Sciences, Quy Nhon University, Vietnam
| | - Vu Thi Ngan
- Laboratory of Computational Chemistry and Modelling (LCCM), Department of Chemistry, Faculty of Natural Sciences, Quy Nhon University, Vietnam
| | - Asit. K. Chandra
- Department of Chemistry, North-Eastern Hill University, Shillong 793022, Meghalaya, India
| | - Nguyen Tien Trung
- Laboratory of Computational Chemistry and Modelling (LCCM), Department of Chemistry, Faculty of Natural Sciences, Quy Nhon University, Vietnam
| |
Collapse
|
10
|
An NT, Duong NT, Tri NN, Trung NT. Role of O–H⋯O/S conventional hydrogen bonds in considerable C sp2–H blue-shift in the binary systems of acetaldehyde and thioacetaldehyde with substituted carboxylic and thiocarboxylic acids. RSC Adv 2022; 12:35309-35319. [DOI: 10.1039/d2ra05391h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2022] [Accepted: 11/27/2022] [Indexed: 12/13/2022] Open
Abstract
The presence of O–H⋯O/S conventional hydrogen bonds in the complex governs a significant blue shift of Csp2–H bonds.
Collapse
Affiliation(s)
- Nguyen Truong An
- Faculty of Natural Sciences, Quy Nhon University, Quy Nhon, Vietnam
| | - Nguyen Thi Duong
- Faculty of Natural Sciences, Quy Nhon University, Quy Nhon, Vietnam
| | - Nguyen Ngoc Tri
- Faculty of Natural Sciences, Quy Nhon University, Quy Nhon, Vietnam
- Laboratory of Computational Chemistry and Modelling (LCCM), Quy Nhon University, Quy Nhon, Vietnam
| | - Nguyen Tien Trung
- Faculty of Natural Sciences, Quy Nhon University, Quy Nhon, Vietnam
- Laboratory of Computational Chemistry and Modelling (LCCM), Quy Nhon University, Quy Nhon, Vietnam
| |
Collapse
|
11
|
Pullanchery S, Kulik S, Rehl B, Hassanali A, Roke S. Charge transfer across C-H⋅⋅⋅O hydrogen bonds stabilizes oil droplets in water. Science 2021; 374:1366-1370. [PMID: 34882471 DOI: 10.1126/science.abj3007] [Citation(s) in RCA: 66] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
[Figure: see text].
Collapse
Affiliation(s)
- Saranya Pullanchery
- Laboratory for Fundamental BioPhotonics, Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Sergey Kulik
- Laboratory for Fundamental BioPhotonics, Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Benjamin Rehl
- Laboratory for Fundamental BioPhotonics, Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Ali Hassanali
- International Centre for Theoretical Physics, 34100 Trieste, Italy
| | - Sylvie Roke
- Laboratory for Fundamental BioPhotonics, Institute of Bioengineering (IBI), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland.,Institute of Materials Science and Engineering (IMX), School of Engineering (STI), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland.,Lausanne Centre for Ultrafast Science, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
12
|
Cuc NTT, Phan CTD, Nhung NTA, Nguyen MT, Trung NT, Ngan VT. Theoretical Aspects of Nonconventional Hydrogen Bonds in the Complexes of Aldehydes and Hydrogen Chalcogenides. J Phys Chem A 2021; 125:10291-10302. [PMID: 34818019 DOI: 10.1021/acs.jpca.1c06708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Hydrogen bonds (H-bonds) in the complexes between aldehydes and hydrogen chalcogenides, XCHO...nH2Z with X = H, F, Cl, Br, and CH3, Z = O, S, Se, and Te, and n = 1,2, were investigated using high-level ab initio calculations. The Csp2-H...O H-bonds are found to be about twice as strong as the Csp2-H...S/Se/Te counterparts. Remarkably, the S/Se/Te-H...S/Se/Te H-bonds are 4.5 times as weak as the O-H...O ones. The addition of the second H2Z molecule into binary systems induces stronger complexes and causes a positive cooperative effect in ternary complexes. The blue shift of Csp2-H stretching frequency involving the Csp2-H...Z H-bond sharply increases when replacing one H atom in HCHO by a CH3 group. In contrast, when one H atom in HCHO is substituted with a halogen, the magnitude of blue-shifting of the Csp2-H...Z H-bond becomes smaller. The largest blue shift up to 92 cm-1 of Csp2-H stretching frequency in Csp2-H...O H-bond in CH3CHO...2H2O has rarely been observed and is much greater than that in the cases of the Csp2-H...S/Se/Te ones. The Csp2-H blue shift of Csp2-H...Z bonds in the halogenated aldehydes is converted into a red shift when H2O is replaced by a heavier analogue, such as H2S, H2Se, or H2Te. The stability and classification of nonconventional H-bonds including Csp2-H...Se/Te, Te-H...Te, and Se/Te-H...O have been established for the first time.
Collapse
Affiliation(s)
- Nguyen Thi Thanh Cuc
- Laboratory of Computational Chemistry and Modelling (LCCM), Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 55100, Vietnam
| | - Cam-Tu Dang Phan
- Laboratory of Computational Chemistry and Modelling (LCCM), Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 55100, Vietnam
| | - Nguyen Thi Ai Nhung
- Department of Chemistry, University of Sciences, Hue University, Hue 49000, Vietnam
| | | | - Nguyen Tien Trung
- Laboratory of Computational Chemistry and Modelling (LCCM), Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 55100, Vietnam
| | - Vu Thi Ngan
- Laboratory of Computational Chemistry and Modelling (LCCM), Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 55100, Vietnam
| |
Collapse
|
13
|
Zech A, Head-Gordon M. Dissociation of HCl in water nanoclusters: an energy decomposition analysis perspective. Phys Chem Chem Phys 2021; 23:26737-26749. [PMID: 34846396 DOI: 10.1039/d1cp04587c] [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
As known, small HCl-water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl-water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl-water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push-pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.
Collapse
Affiliation(s)
- Alexander Zech
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA.
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
14
|
Cárdenas G, Pérez-Barcia Á, Mandado M, Nogueira JJ. Characterization of cisplatin/membrane interactions by QM/MM energy decomposition analysis. Phys Chem Chem Phys 2021; 23:20533-20540. [PMID: 34505588 DOI: 10.1039/d1cp03382d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We extend for the first time a quantum mechanical energy decomposition analysis scheme based on deformation electron densities to a hybrid electrostatic embedding quantum mechanics/molecular mechanics framework. The implemented approach is applied to characterize the interactions between cisplatin and a dioleyl-phosphatidylcholine membrane, which play a key role in the permeation mechanism of the drug inside the cells. The interaction energy decomposition into electrostatic, induction, dispersion and Pauli repulsion contributions is performed for ensembles of geometries to account for conformational sampling. It is evidenced that the electrostatic and repulsive components are predominant in both polar and non-polar regions of the bilayer.
Collapse
Affiliation(s)
- Gustavo Cárdenas
- Department of Chemistry, Universidad Autónoma de Madrid, Calle Francisco Tomás y Valiente, 7, 28049, Madrid, Spain.
| | - Álvaro Pérez-Barcia
- Department of Physical Chemistry, University of Vigo, Lagoas-Marcosende s/n, ES-36310-Vigo, Galicia, Spain.
| | - Marcos Mandado
- Department of Physical Chemistry, University of Vigo, Lagoas-Marcosende s/n, ES-36310-Vigo, Galicia, Spain.
| | - Juan J Nogueira
- Department of Chemistry, Universidad Autónoma de Madrid, Calle Francisco Tomás y Valiente, 7, 28049, Madrid, Spain. .,IADCHEM, Institute for Advanced Research in Chemistry, Universidad Autónoma de Madrid, Calle Francisco Tomás y Valiente, 7, 28049, Madrid, Spain
| |
Collapse
|
15
|
Enudi OC, Louis H, Edim MM, Agwupuye JA, Ekpen FO, Bisong EA, Utsu PM. Understanding the aqueous chemistry of quinoline and the diazanaphthalenes: insight from DFT study. Heliyon 2021; 7:e07531. [PMID: 34296019 PMCID: PMC8282981 DOI: 10.1016/j.heliyon.2021.e07531] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2021] [Revised: 05/18/2021] [Accepted: 07/06/2021] [Indexed: 11/16/2022] Open
Abstract
The inter-fragment interactions at various binding sites and the overall cluster stability of quinolone (QNOL), cinnoline (CNOL), quinazoline (QNAZ), and quinoxaline (QNOX) complexes with H2O were studied using the density functional theory (DFT) approach. The adsorption and H-bond binding energies, and the energy decomposition mechanism was considered to determine the relative stabilization status of the studied clusters. Scanning tunneling microscopy (STM), natural bonding orbitals (NBO) and charge decomposition were studied to expose the electronic distribution and interaction between fragments. The feasibility of formations of the various complexes were also studied by considering their thermodynamic properties. Results from adsorption studies confirmed the actual adsorption of H2O molecules on the various binding sites studied, with QNOX clusters exhibiting the best adsorptions. Charge decomposition analysis (CDA) revealed significant charge transfer from substrate to H2O fragment in most complexes, except in QNOL, CNOL and QNAZ clusters with H2O at binding position 4, where much charges are back-donated to substrate. The O---H inter-fragment bonds was discovered to be stronger than counterpart N---H bonds in the complexes, whilst polarity indices confirmed N---H as more polar covalent than O---H bonds. Thermodynamic considerations revealed that the formation process of all studied complexes are endothermic (+ve ΔH f ) and non-spontaneous (+ve ΔG f ).
Collapse
Affiliation(s)
- Obieze C. Enudi
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
| | - Hitler Louis
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
| | - Moses M. Edim
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
| | - John A. Agwupuye
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
| | - Francis O. Ekpen
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
| | - Emmanuel A. Bisong
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
| | - Patrick M. Utsu
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
| |
Collapse
|
16
|
Levina EO, Khrenova MG, Tsirelson VG. The explicit role of electron exchange in the hydrogen bonded molecular complexes. J Comput Chem 2021; 42:870-882. [PMID: 33675552 DOI: 10.1002/jcc.26507] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2020] [Revised: 02/16/2021] [Accepted: 02/18/2021] [Indexed: 01/22/2023]
Abstract
We applied a set of advanced bonding descriptors to establish the hidden electron density features and binding energy characteristics of intermolecular DH∙∙∙A hydrogen bonds (OH∙∙∙O, NH∙∙∙O and SH∙∙∙O) in 150 isolated and solvated molecular complexes. The exchange-correlation and Pauli potentials as well as corresponding local one-electron forces allowed us to explicitly ascertain how electron exchange defines the bonding picture in the proximity of the H-bond critical point. The electron density features of DH∙∙∙A interaction are governed by alterations in the electron localization in the H-bond region displaying itself in the exchange hole. At that, they do not depend on the variations in the exchange hole mobility. The electrostatic interaction mainly defines the energy of H-bonds of different types, whereas the strengthening/weakening of H-bonds in complexes with varying substituents depends on the barrier height of the exchange potential near the bond critical point. Energy variations between H-bonds in isolated and solvated systems are also caused the electron exchange peculiarities as follows from the corresponding potential and the interacting quantum atom analyses complemented by electron delocalization index calculations. Our approach is based on the bonding descriptors associated with the characteristics of the observable electron density and can be recommended for in-depth studies of non-covalent bonding.
Collapse
Affiliation(s)
- Elena O Levina
- Federal Research Centre "Fundamentals of Biotechnology" of the Russian Academy of Sciences, Moscow, Russia.,Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | - Maria G Khrenova
- Federal Research Centre "Fundamentals of Biotechnology" of the Russian Academy of Sciences, Moscow, Russia.,Moscow Institute of Physics and Technology, Dolgoprudny, Russia.,Lomonosov Moscow State University, Moscow, Russia
| | | |
Collapse
|
17
|
Mao Y, Loipersberger M, Horn PR, Das A, Demerdash O, Levine DS, Prasad Veccham S, Head-Gordon T, Head-Gordon M. From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis. Annu Rev Phys Chem 2021; 72:641-666. [PMID: 33636998 DOI: 10.1146/annurev-physchem-090419-115149] [Citation(s) in RCA: 47] [Impact Index Per Article: 15.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N2, and BF bound to a [Ru(II)(NH3)5]2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.
Collapse
Affiliation(s)
- Yuezhi Mao
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Matthias Loipersberger
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Paul R Horn
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Akshaya Das
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA; .,Department of Bioengineering and Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
| | - Omar Demerdash
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA; .,Department of Bioengineering and Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
| | - Daniel S Levine
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Srimukh Prasad Veccham
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| | - Teresa Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA; .,Department of Bioengineering and Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA;
| |
Collapse
|
18
|
Veccham SP, Lee J, Mao Y, Horn PR, Head-Gordon M. A non-perturbative pairwise-additive analysis of charge transfer contributions to intermolecular interaction energies. Phys Chem Chem Phys 2021; 23:928-943. [DOI: 10.1039/d0cp05852a] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A non-perturbative scheme for complete decomposition of energy and charge associated with charge transfer interaction into pairwise additive components.
Collapse
Affiliation(s)
| | - Joonho Lee
- Department of Chemistry
- University of California
- Berkeley
- USA
- Chemical Sciences Division
| | - Yuezhi Mao
- Department of Chemistry
- University of California
- Berkeley
- USA
- Chemical Sciences Division
| | - Paul R. Horn
- Department of Chemistry
- University of California
- Berkeley
- USA
- Chemical Sciences Division
| | - Martin Head-Gordon
- Department of Chemistry
- University of California
- Berkeley
- USA
- Chemical Sciences Division
| |
Collapse
|
19
|
Mao Y, Loipersberger M, Kron KJ, Derrick JS, Chang CJ, Sharada SM, Head-Gordon M. Consistent inclusion of continuum solvation in energy decomposition analysis: theory and application to molecular CO 2 reduction catalysts. Chem Sci 2020; 12:1398-1414. [PMID: 34163903 PMCID: PMC8179122 DOI: 10.1039/d0sc05327a] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
To facilitate computational investigation of intermolecular interactions in the solution phase, we report the development of ALMO-EDA(solv), a scheme that allows the application of continuum solvent models within the framework of energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs). In this scheme, all the quantum mechanical states involved in the variational EDA procedure are computed with the presence of solvent environment so that solvation effects are incorporated in the evaluation of all its energy components. After validation on several model complexes, we employ ALMO-EDA(solv) to investigate substituent effects on two classes of complexes that are related to molecular CO2 reduction catalysis. For [FeTPP(CO2-κC)]2- (TPP = tetraphenylporphyrin), we reveal that two ortho substituents which yield most favorable CO2 binding, -N(CH3)3 + (TMA) and -OH, stabilize the complex via through-structure and through-space mechanisms, respectively. The coulombic interaction between the positively charged TMA group and activated CO2 is found to be largely attenuated by the polar solvent. Furthermore, we also provide computational support for the design strategy of utilizing bulky, flexible ligands to stabilize activated CO2 via long-range Coulomb interactions, which creates biomimetic solvent-inaccessible "pockets" in that electrostatics is unscreened. For the reactant and product complexes associated with the electron transfer from the p-terphenyl radical anion to CO2, we demonstrate that the double terminal substitution of p-terphenyl by electron-withdrawing groups considerably strengthens the binding in the product state while moderately weakens that in the reactant state, which are both dominated by the substituent tuning of the electrostatics component. These applications illustrate that this new extension of ALMO-EDA provides a valuable means to unravel the nature of intermolecular interactions and quantify their impacts on chemical reactivity in solution.
Collapse
Affiliation(s)
- Yuezhi Mao
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
| | | | - Kareesa J Kron
- Mork Family Department of Chemical Engineering and Material Science, University of Southern California Los Angeles CA 90089 USA
| | - Jeffrey S Derrick
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA
| | - Christopher J Chang
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA .,Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA.,Department of Molecular and Cell Biology, University of California Berkeley Berkeley CA 94720 USA
| | - Shaama Mallikarjun Sharada
- Mork Family Department of Chemical Engineering and Material Science, University of Southern California Los Angeles CA 90089 USA.,Department of Chemistry, University of Southern California Los Angeles CA 90089 USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California at Berkeley Berkeley CA 94720 USA .,Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
| |
Collapse
|
20
|
Kazim M, Guan L, Chopra A, Sun R, Siegler MA, Lectka T. Switching a HO···π Interaction to a Nonconventional OH···π Hydrogen Bond: A Completed Crystallographic Puzzle. J Org Chem 2020; 85:9801-9807. [PMID: 32633510 DOI: 10.1021/acs.joc.0c01121] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this article, we present crystallographic and spectroscopic evidence of a tunable system wherein a HO···π interaction switches incrementally to a nonconventional OH···π hydrogen bonding (HB) interaction. More specifically, we report the synthesis of substituted forms of model system 1 to study the effects of aryl ring electronic density on the qualitative characteristics of OH···π hydrogen bonds therein. The OH stretch in experimental infrared data, in agreement with density-functional theory (DFT) calculations, shows continuous red-shifts as the adjacent ring becomes more electron rich. For example, the OH stretch of an amino-substituted analogue is red-shifted by roughly 50 cm-1 compared to the same stretch in the CF3 analogue, indicating a significantly stronger HB interaction in the former. Moreover, DFT calculations (ωB97XD/6-311+G**) predict that increasing electronic density on the adjacent top ring reduces the aryl π-OH σ* energy gap with a concomitant enhancement of the OH n-π* energy gap. Consequently, a dominant π-σ* interaction in the amino substituted analogue locks the system in the in-form while a favorable n-π* interaction reverses the orientation of the oxygen-bound hydrogen in its protonated form. Additionally, the 1H NMR data of various analogues reveal that stronger OH···π interactions in systems with electron-rich aromatic rings slow exchange of the alcoholic proton, thereby revealing coupling with the geminal proton. Finally, X-ray crystallographic analyses of a spectrum of analogues clearly visualize the three distinct stages of "switch"-starting with exclusive HO···π, to partitioned HO···π/OH···π, and finally to achieving exclusive OH···π forms. Furthermore, the crystal structure of the amino analogue reveals an interesting feature in which an extended HB network, involving two conventional (NH···O) and two nonconventional (OH···π) HBs, dimerizes and anchors the molecule in the unit cell.
Collapse
Affiliation(s)
- Muhammad Kazim
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Liangyu Guan
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States.,Calibr-A Division of Scripps Research, 11119 N Torrey Pines Rd, La Jolla, San Diego, California 92037, United States
| | - Anant Chopra
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Roy Sun
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Maxime A Siegler
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Thomas Lectka
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| |
Collapse
|
21
|
Mao Y, Levine DS, Loipersberger M, Horn PR, Head-Gordon M. Probing radical-molecule interactions with a second generation energy decomposition analysis of DFT calculations using absolutely localized molecular orbitals. Phys Chem Chem Phys 2020; 22:12867-12885. [PMID: 32510096 DOI: 10.1039/d0cp01933j] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Intermolecular interactions between radicals and closed-shell molecules are ubiquitous in chemical processes, ranging from the benchtop to the atmosphere and extraterrestrial space. While energy decomposition analysis (EDA) schemes for closed-shell molecules can be generalized for studying radical-molecule interactions, they face challenges arising from the unique characteristics of the electronic structure of open-shell species. In this work, we introduce additional steps that are necessary for the proper treatment of radical-molecule interactions to our previously developed unrestricted Absolutely Localized Molecular Orbital (uALMO)-EDA based on density functional theory calculations. A "polarize-then-depolarize" (PtD) scheme is used to remove arbitrariness in the definition of the frozen wavefunction, rendering the ALMO-EDA results independent of the orientation of the unpaired electron obtained from isolated fragment calculations. The contribution of radical rehybridization to polarization energies is evaluated. It is also valuable to monitor the wavefunction stability of each intermediate state, as well as their associated spin density profiles, to ensure the EDA results correspond to a desired electronic state. These radical extensions are incorporated into the "vertical" and "adiabatic" variants of uALMO-EDA for studies of energy changes and property shifts upon complexation. The EDA is validated on two model complexes, H2O˙F and FH˙OH. It is then applied to several chemically interesting radical-molecule complexes, including the sandwiched and T-shaped benzene dimer radical cation, complexes of pyridine with benzene and naphthalene radical cations, binary and ternary complexes of the hydroxyl radical with water (˙OH(H2O) and ˙OH(H2O)2), and the pre-reactive complexes and transition states in the ˙OH + HCHO and ˙OH + CH3CHO reactions. These examples suggest that this second generation uALMO-EDA is a useful tool for furthering one's understanding of both energetic and property changes associated with radical-molecule interactions.
Collapse
Affiliation(s)
- Yuezhi Mao
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, CA 94720, USA.
| | | | | | | | | |
Collapse
|
22
|
Chen X, Gao J. Fragment Exchange Potential for Realizing Pauli Deformation of Interfragment Interactions. J Phys Chem Lett 2020; 11:4008-4016. [PMID: 32308000 DOI: 10.1021/acs.jpclett.0c00933] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
In fragment-based methods, the lack of explicit short-range exchange interactions between monomers can result in unphysical deformation in charge density. In this study, we describe a fragment exchange potential (XFP) to explicitly account for interfragmental Pauli deformation. In our implementation, a Kohn-Sham exchange potential is adopted along with the Yukawa potential. The method has been validated by comparison of the computed exchange energies using the XFP potential with results obtained from antisymmetrized fragmental orbitals on the S66×8 data set containing 528 bimolecular interactions of equilibrium and arbitrary geometries. It was also found that it is only necessary to deploy numerical grids on atoms within their van der Waals contacts, significantly reducing the small, albeit extra, computational cost. We anticipate that the XFP presented here may be applied to molecular dynamics simulations of macromolecules using a fragment-based quantum mechanical potential with improved SCF convergence and computational accuracy.
Collapse
Affiliation(s)
- Xin Chen
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong, China
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
23
|
Kazim M, Siegler MA, Lectka T. Close Amide NH···F Hydrogen Bonding Interactions in 1,8-Disubstituted Naphthalenes. J Org Chem 2020; 85:6195-6200. [PMID: 32227992 DOI: 10.1021/acs.joc.0c00553] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
In this note, we present a series of N-(8-fluoronaphthalen-1-yl)benzamide derivatives designed to maximize amide-NH···F hydrogen bond interactions therein. A combination of IR and NMR spectroscopy indicates a linear correlation between the high energy shift in NH stretching frequency and the electron withdrawing nature of the substituent, consistent with the trend predicted by DFT calculations. Additionally, a limiting case of hydrogen bonding is observed when the benzamide derivatives are replaced with trifluoroacetamide, causing an additional red shift of 44 cm-1 in the NH stretching frequency. Most importantly, 1H-19F coupling constants in this series are among the largest measured for amide-NH···F interactions. X-ray crystallography reveals face-to-face alignment of naphthalene rings in these derivatives resulting in part from the NH···F hydrogen bonds. This motif also dictates the formation of sheets composed of stacked naphthalene rings in the crystal structure as opposed to unfluorinated analogues wherein NH···OC hydrogen-bonding interactions force benzamide and naphthalene rings to engage in T-shaped π-π interactions instead. Additionally, the NH proton in the trifluoroacetamide derivative engages in extended H-bond interactions in its crystal structure.
Collapse
Affiliation(s)
- Muhammad Kazim
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Maxime A Siegler
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| | - Thomas Lectka
- Department of Chemistry, Johns Hopkins University, 3400 North Charles Street, Baltimore, Maryland 21218, United States
| |
Collapse
|
24
|
Loipersberger M, Mao Y, Head-Gordon M. Variational Forward–Backward Charge Transfer Analysis Based on Absolutely Localized Molecular Orbitals: Energetics and Molecular Properties. J Chem Theory Comput 2020; 16:1073-1089. [DOI: 10.1021/acs.jctc.9b01168] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Matthias Loipersberger
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Yuezhi Mao
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
25
|
Rossomme E, Lininger CN, Bell AT, Head-Gordon T, Head-Gordon M. Electronic structure calculations permit identification of the driving forces behind frequency shifts in transition metal monocarbonyls. Phys Chem Chem Phys 2020; 22:781-798. [DOI: 10.1039/c9cp04643g] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Our direct DFT decomposition of CO frequency shifts updates the paradigm for metal carbonyl binding.
Collapse
Affiliation(s)
- Elliot Rossomme
- Kenneth S. Pitzer Center for Theoretical Chemistry
- Department of Chemistry
- University of California
- Berkeley
- USA
| | - Christianna N. Lininger
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Department of Chemical and Biomolecular Engineering
| | - Alexis T. Bell
- Chemical Sciences Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- Department of Chemical and Biomolecular Engineering
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry
- Department of Chemistry
- University of California
- Berkeley
- USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry
- Department of Chemistry
- University of California
- Berkeley
- USA
| |
Collapse
|