1
|
Finney JM, McCoy AB. Correlations between the Structures and Spectra of Protonated Water Clusters. J Phys Chem A 2024; 128:868-879. [PMID: 38265889 DOI: 10.1021/acs.jpca.3c07338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2024]
Abstract
Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.
Collapse
Affiliation(s)
- Jacob M Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
2
|
Santis GD, Herman KM, Heindel JP, Xantheas SS. Descriptors of water aggregation. J Chem Phys 2024; 160:054306. [PMID: 38341703 DOI: 10.1063/5.0179815] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2023] [Accepted: 01/05/2024] [Indexed: 02/13/2024] Open
Abstract
We rely on a total of 23 (cluster size, 8 structural, and 14 connectivity) descriptors to investigate structural patterns and connectivity motifs associated with water cluster aggregation. In addition to the cluster size n (number of molecules), the 8 structural descriptors can be further categorized into (i) one-body (intramolecular): covalent OH bond length (rOH) and HOH bond angle (θHOH), (ii) two-body: OO distance (rOO), OHO angle (θOHO), and HOOX dihedral angle (ϕHOOX), where X lies on the bisector of the HOH angle, (iii) three-body: OOO angle (θOOO), and (iv) many-body: modified tetrahedral order parameter (q) to account for two-, three-, four-, five-coordinated molecules (qm, m = 2, 3, 4, 5) and radius of gyration (Rg). The 14 connectivity descriptors are all many-body in nature and consist of the AD, AAD, ADD, AADD, AAAD, AAADD adjacencies [number of hydrogen bonds accepted (A) and donated (D) by each water molecule], Wiener index, Average Shortest Path Length, hydrogen bond saturation (% HB), and number of non-short-circuited three-membered cycles, four-membered cycles, five-membered cycles, six-membered cycles, and seven-membered cycles. We mined a previously reported database of 4 948 959 water cluster minima for (H2O)n, n = 3-25 to analyze the evolution and correlation of these descriptors for the clusters within 5 kcal/mol of the putative minima. It was found that rOH and % HB correlated strongly with cluster size n, which was identified as the strongest predictor of energetic stability. Marked changes in the adjacencies and cycle count were observed, lending insight into changes in the hydrogen bond network upon aggregation. A Principal Component Analysis (PCA) was employed to identify descriptor dependencies and group clusters into specific structural patterns across different cluster sizes. The results of this study inform our understanding of how water clusters evolve in size and what appropriate descriptors of their structural and connectivity patterns are with respect to system size, stability, and similarity. The approach described in this study is general and can be easily extended to other hydrogen-bonded systems.
Collapse
Affiliation(s)
- Garrett D Santis
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MSIN J7-10, Richland, Washington 99352, USA
- Computational and Theoretical Chemistry Institute (CTCI), Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
3
|
LaCour RA, Heindel JP, Head-Gordon T. Predicting the Raman Spectra of Liquid Water with a Monomer-Field Model. J Phys Chem Lett 2023; 14:11742-11749. [PMID: 38116782 DOI: 10.1021/acs.jpclett.3c02873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2023]
Abstract
The Raman spectrum of liquid water is quite complex, reflecting its strong sensitivity to the local environment of the individual waters. The OH-stretch region of the spectrum, which captures the influence of hydrogen bonding, has only just begun to be unraveled. Here we develop a model for predicting the Raman spectra of the OH-stretch region by considering how local electric fields distort the energy surface of each water monomer. We find that our model is capable of reproducing the bimodal nature of the main peak, with the shoulder at 3250 cm-1 resulting almost entirely from Fermi resonance. Furthermore, we capture the temperature and polarization dependence of the shoulder, which has proven to be difficult to obtain with previous methods, and analyze the origin of this dependence. We expect our model to be generally useful for understanding and predicting how Raman spectra change under different conditions and with different probe reporters beyond water.
Collapse
Affiliation(s)
- R Allen LaCour
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Joseph P Heindel
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Departments of Bioengineering and Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, United States
| |
Collapse
|
4
|
Foreman MM, Stanton JF, Weber JM. Relation Between Bond Angle and Carbon-Oxygen Stretching Frequencies in CO 2-Containing Compounds. J Phys Chem A 2023; 127:9717-9722. [PMID: 37944122 DOI: 10.1021/acs.jpca.3c05082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2023]
Abstract
The symmetric (νs) and antisymmetric (νas) O-C-O stretching modes of CO2-containing compounds encode structural information that can be difficult to decipher, due to the sensitivity of these spectral features to small shifts in charge distribution and structure, as well as the anharmonicities of these two vibrational modes. In this work, we discuss the relation between the frequency of these modes and the geometry of the O-C-O group, showing that the splitting between νs and νas (Δνas-s = νas - νs) can be predicted based only on the O-C-O bond angle obtained from quantum chemical calculations with reasonable accuracy (±46 cm-1, R2 = 0.994). The relationship is shown to hold for the infrared spectra of a variety of CO2-containing molecules measured in vacuo. The origins of this model are discussed in the framework of elementary mode analysis.
Collapse
Affiliation(s)
- Madison M Foreman
- JILA and Department of Chemistry, University of Colorado-Boulder, 440 UCB, Boulder, Colorado 80309-0440, United States
| | - John F Stanton
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - J Mathias Weber
- JILA and Department of Chemistry, University of Colorado-Boulder, 440 UCB, Boulder, Colorado 80309-0440, United States
| |
Collapse
|
5
|
Kato T, Fujii A. Experimental confirmation of the Badger-Bauer rule in the protonated methanol clusters: weak hydrogen bond formation as a measure of terminal OH acidity in hydrogen bond networks. Phys Chem Chem Phys 2023; 25:30188-30192. [PMID: 37920966 DOI: 10.1039/d3cp04644c] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
We report a linear correlation between the OH stretch frequency shift of the protonated methanol cluster, H+(MeOH)n, upon the π-hydrogen bond formation with benzene and the enthalpy change in clustering of H+(MeOH)n to H+(MeOH)n+1. This result suggests a new method to explore hydrogen bond strength in hydrogen bond networks.
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
|
6
|
Liberka M, Zychowicz M, Hooper J, Nakabayashi K, Ohkoshi SI, Chorazy S. Synchronous Switching of Dielectric Constant and Photoluminescence in Cyanidonitridorhenate-Based Crystals. Angew Chem Int Ed Engl 2023; 62:e202308284. [PMID: 37615930 DOI: 10.1002/anie.202308284] [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: 06/12/2023] [Revised: 08/05/2023] [Accepted: 08/24/2023] [Indexed: 08/25/2023]
Abstract
Switching of multiple physical properties by external stimuli in dynamic materials enables applications in, e.g., smart sensors, biomedical tools, as well as data-storage devices. Among stimuli-responsive materials, inorganic-organic molecular hybrids exhibiting thermal order-disorder phase transitions were tested as promising molecular switches of electrical characteristics, including dielectric constant. We aimed at broadening the multifunctional potential of such hybrid materials towards the switching of not only electrical but also other physical properties, e.g., light emission. We report two ionic salts based on luminescent tetracyanidonitridorhenate(V) anions bearing two different diamine ligands, 1,2-diaminoethane (1) and 1,3-diaminopropane (2), both crystallizing with polar N-methyl-dabconium cations. They exhibit an order-disorder phase transition related to the heating-induced turning-on of the rotation of polar cations. This leads to a unique synchronous switching of the dielectric constant as well as metal-complex-centered photoluminescence, as demonstrated by changes in, e.g., emission lifetime. The roles of organic cations, non-trivial Re(V) complexes, and their interaction in achieving the coupled thermal switching of electrical and optical properties are discussed utilizing experimental and theoretical approaches.
Collapse
Affiliation(s)
- Michal Liberka
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387, Kraków, Poland
- Doctoral School of Exact and Natural Sciences, Jagiellonian University, Łojasiewicza 11, 30-348, Kraków, Poland
| | - Mikolaj Zychowicz
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387, Kraków, Poland
- Doctoral School of Exact and Natural Sciences, Jagiellonian University, Łojasiewicza 11, 30-348, Kraków, Poland
| | - James Hooper
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387, Kraków, Poland
| | - Koji Nakabayashi
- Department of Chemistry, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan
| | - Shin-Ichi Ohkoshi
- Department of Chemistry, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan
| | - Szymon Chorazy
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387, Kraków, Poland
| |
Collapse
|
7
|
Kameda Y, Amo Y, Usuki T, Ikeda K, Honda T, Otomo T. Direct Determination of the Relationship between the Intramolecular Oxygen-Hydrogen Bond Length and Its Stretching Vibrational Frequency of the Methanol Molecule in the Liquid State. J Phys Chem B 2023; 127:7758-7763. [PMID: 37671846 DOI: 10.1021/acs.jpcb.3c05090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2023]
Abstract
Time-of-flight (TOF) neutron diffraction measurements on pure liquid deuterated methanol and concentrated methanolic LiClO4 and LiTFSA solutions have been carried out to investigate the effect of intermolecular hydrogen bonds on the intramolecular O-D distance (rOD) of the methanol molecule in the liquid state. Intramolecular parameters for the methanol molecule are determined by the least-squares fitting analysis of the neutron total interference term observed in the high-Q region. Attenuated total reflection (ATR) IR spectra have been measured for methanolic solutions of natural abundance to determine the gravitational center frequency (νOH) of the stretching vibrational band of the methanol molecule. The relationship between rOD and νOH is approximated well by a linear function. The value dνOH/drOD = -17000 ± 3000 cm-1 Å-1 has been derived from the slope of the fitted function. It has been revealed that the O-D bond length of the methanol molecule is sensitively affected by the intermolecular hydrogen bonding interaction.
Collapse
Affiliation(s)
- Yasuo Kameda
- Department of Material and Biological Chemistry, Faculty of Science, Yamagata University, 1-4-12, Kojirakawa-machi, Yamagata City, Yamagata 990-8560, Japan
| | - Yuko Amo
- Department of Material and Biological Chemistry, Faculty of Science, Yamagata University, 1-4-12, Kojirakawa-machi, Yamagata City, Yamagata 990-8560, Japan
| | - Takeshi Usuki
- Department of Material and Biological Chemistry, Faculty of Science, Yamagata University, 1-4-12, Kojirakawa-machi, Yamagata City, Yamagata 990-8560, Japan
| | - Kazutaka Ikeda
- Neutron Industrial Application Promotion Center, Comprehensive Research Organization for Science and Society (CROSS), Shirakata, Tokai, Naka, Ibaraki 319-1106, Japan
| | - Takashi Honda
- Institute of Materials Structure Science, KEK, 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan
| | - Toshiya Otomo
- Institute of Materials Structure Science, KEK, 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan
| |
Collapse
|
8
|
Ahirwar MB, Deshmukh MM. Fragments-in-fragments method for efficient and reliable estimates of individual hydrogen bond energies in large molecular clusters. J Comput Chem 2023. [PMID: 37191018 DOI: 10.1002/jcc.27133] [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/04/2023] [Revised: 05/01/2023] [Accepted: 05/03/2023] [Indexed: 05/17/2023]
Abstract
The knowledge of individual hydrogen bond (HB) strength in molecular clusters is indispensable to get insights into the bulk properties of condensed systems. Recently, we have developed the molecular tailoring approach based (MTA-based) method for the estimation of individual HB energy in molecular clusters. However, the direct use of this MTA-based method to large molecular clusters becomes progressively difficult with the increase in the size of a cluster. To overcome this caveat, herein, we propose the use of linear scaling method (such as the original MTA method) for the estimation of single-point (SP) energies of large-sized parent molecular cluster and their respective fragments. Because the fragments of the MTA-based method, for the estimation of HB energy, are further fragmented, this proposed strategy is called as Fragments-in-Fragments (Frags-in-Frags) method. The SP energies of fragments and parent cluster calculated by the Frags-in-Frags approach were utilized to estimate the individual HB energy. The estimated individual HB energies, in various molecular clusters, by Frags-in-Frags method are found to be in excellent linear agreement with their MTA-based counterparts (R2 = 0.9975 of 348 data points). The difference being less than 0.5 kcal/mol in most of the cases. Furthermore, RMSD is 0.43 kcal/mol, MAE is 0.33 kcal/mol, and the standard deviation is 0.44 kcal/mol. Importantly, the Frags-in-Frags method not only enables the reliable estimation of HB energy in large molecular clusters but also requires less computational time and can be possible even with off-the-shelf hardware.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| |
Collapse
|
9
|
Roget SA, Heck TR, Carter-Fenk KA, Fayer MD. Ion/Water Network Structural Dynamics in Highly Concentrated Lithium Chloride and Lithium Bromide Solutions Probed with Ultrafast Infrared Spectroscopy. J Phys Chem B 2023; 127:4532-4543. [PMID: 37172191 DOI: 10.1021/acs.jpcb.2c08792] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
The structural dynamics of highly concentrated LiCl and LiBr aqueous solutions were observed from 1-4 to 1-16 water molecules per ion pair using ultrafast polarization-selective pump-probe (PSPP) experiments on the OD stretch of dilute HOD. At these high salt concentrations, an extended ion/water network exists with complex structural dynamics. Population decays from PSPP experiments highlight two distinct water components. From the frequency-dependent amplitudes of the decays, the spectra of hydroxyls bound to halides and to water oxygens are obtained, which are not observable in the FT-IR spectra. PSPP experiments also measure frequency-dependent water orientational relaxation. At short times, wobbling dynamics within a restricted angular cone occurs. At high concentrations, the cone angles are dependent on frequency (hydrogen bond strength), but at higher water concentrations (>10 waters per ion pair), there is no frequency dependence. The average cone angle increases as the ion concentration decreases. The slow time constant for complete HOD orientational relaxation is independent of concentration but slower in LiCl than in LiBr. Comparison to structural MD simulations of LiCl from the literature indicates that the loss of the cone angle wavelength dependence and the increase in the cone angles as the concentration decreases occur as the prevalence of large ion/water clusters gives way to contact ion pairs.
Collapse
Affiliation(s)
- Sean A Roget
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Tristan R Heck
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | | | - Michael D Fayer
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
10
|
Herman KM, Xantheas SS. An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2-25. Phys Chem Chem Phys 2023; 25:7120-7143. [PMID: 36853239 DOI: 10.1039/d2cp03241d] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2-11, 16-17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (De's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2-25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83-119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (De's) throughout the entire cluster range (n = 2-25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ∼4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6-3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.
Collapse
Affiliation(s)
- Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA.
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA. .,Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, WA, 99352, USA.
| |
Collapse
|
11
|
Kemper M, Drost DA, Engelage E, Merten C. Stereochemistry Controls Dihydrogen Bonding Strengths in Chiral Amine Boranes Adducts. Angew Chem Int Ed Engl 2022; 61:e202213859. [PMID: 36245340 PMCID: PMC10099978 DOI: 10.1002/anie.202213859] [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: 09/20/2022] [Indexed: 11/07/2022]
Abstract
The growing interest in exploiting novel concepts of non-covalent interactions in catalysts and supramolecular chemistry made us revisit a special kind of hydrogen bonding: the dihydrogen bond (DHB), formed between a classical hydrogen bond donor and a hydridic hydrogen as acceptor. Herein, we investigate how the strength of the N-Hδ+ ⋅⋅⋅δ- H-B interaction and hence the DHB-driven self-aggregation of amine-borane adducts is governed by steric effects by comparing the structures and binding enthalpies of various chiral derivatives. For a diastereomeric pair of amine-boranes prepared from a chiral secondary amine, we show that the stereochemistry at the nitrogen has significant influence on the interaction enthalpy. Based on this finding, N-chiral amine boranes can be envisioned to become interesting building blocks in supramolecular chemistry to fine-tune the formation dynamics of assemblies.
Collapse
Affiliation(s)
- Michael Kemper
- Ruhr Universität Bochum, Fakultät für Chemie und Biochemie, Organische Chemie II, Universitätsstraße 150, 44801, Bochum, Germany
| | - Deborah A Drost
- Ruhr Universität Bochum, Fakultät für Chemie und Biochemie, Organische Chemie II, Universitätsstraße 150, 44801, Bochum, Germany
| | - Elric Engelage
- Ruhr Universität Bochum, Fakultät für Chemie und Biochemie, Organische Chemie II, Universitätsstraße 150, 44801, Bochum, Germany
| | - Christian Merten
- Ruhr Universität Bochum, Fakultät für Chemie und Biochemie, Organische Chemie II, Universitätsstraße 150, 44801, Bochum, Germany
| |
Collapse
|
12
|
Nandi A, Martin JML. Thermally-Activated Tunneling in the Two-Water Bridge Catalyzed Tautomerization of Phosphinylidene Compounds. Chemphyschem 2022; 23:e202200396. [PMID: 35867911 PMCID: PMC9804263 DOI: 10.1002/cphc.202200396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Revised: 07/15/2022] [Indexed: 01/05/2023]
Abstract
Phosphinylidenes are an important class of organophosphorus compounds that can exhibit tautomerization between tricoordinated P(III) hydroxide (R1 R2 POH) and a pentacoordinated P(V) oxide (R1 R2 P(O)H) form. Herein we show, using the canonical variational transition state theory combined with multidimensional small-curvature tunneling approximation, the dominance of proton tunneling in the two-water-bridged tautomerizations of phosphinous acid and model phosphinylidenes comprising phosphosphinates, H-phosphonates, H-phosphinates and secondary phosphine oxides. Based on the studied system, the contribution of thermally-activated tunneling is predicted to speed up the semiclassical reaction rate by ca. threefold to as large as two orders of magnitude at 298.15 K in the gas phase. The large KIE and the concavity in the Arrhenius plots are further fingerprints of tunneling. The simulations also predicted that the rapid tunneling rate and short half-life span for the forward reaction, as opposed to the reverse reaction in fluorinated secondary phosphine oxides, would result in P(V) being elusive and only P(III) being isolable, which agrees with previous experiments where only P(III) was detected by IR and NMR spectroscopy. We also explored the role of solvent and predicted tunneling to be substantial.
Collapse
Affiliation(s)
- Ashim Nandi
- Department of Molecular Chemistry and Materials ScienceWeizmann Institute of Science7610001RehovotIsrael
| | - Jan M. L. Martin
- Department of Molecular Chemistry and Materials ScienceWeizmann Institute of Science7610001RehovotIsrael
| |
Collapse
|
13
|
Schröder M, Gatti F, Lauvergnat D, Meyer HD, Vendrell O. The coupling of the hydrated proton to its first solvation shell. Nat Commun 2022; 13:6170. [PMID: 36257946 PMCID: PMC9579203 DOI: 10.1038/s41467-022-33650-w] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2022] [Accepted: 09/28/2022] [Indexed: 11/30/2022] Open
Abstract
The Zundel ([Formula: see text]) and Eigen ([Formula: see text]) cations play an important role as intermediate structures for proton transfer processes in liquid water. In the gas phase they exhibit radically different infrared (IR) spectra. The question arises: is there a least common denominator structure that explains the IR spectra of both, the Zundel and Eigen cations, and hence of the solvated proton? Full dimensional quantum simulations of these protonated cations demonstrate that two dynamical water molecules and an excess proton constitute this fundamental subunit. Embedded in the static environment of the parent Eigen cation, this subunit reproduces the positions and broadenings of its main excess-proton bands. In isolation, its spectrum reverts to the well-known Zundel ion. Hence, the dynamics of this subunit polarized by an environment suffice to explain the spectral signatures and anharmonic couplings of the solvated proton in its first solvation shell.
Collapse
Affiliation(s)
- Markus Schröder
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120, Heidelberg, Germany.
| | - Fabien Gatti
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay UMR 8214, 91405, Orsay, France
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR 8000, 91405, Orsay, France
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120, Heidelberg, Germany
| | - Oriol Vendrell
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120, Heidelberg, Germany.
| |
Collapse
|
14
|
Yang N, Huchmala RM, McCoy AB, Johnson MA. Character of the OH Bend-Stretch Combination Band in the Vibrational Spectra of the "Magic" Number H 3O +(H 2O) 20 and D 3O +(D 2O) 20 Cluster Ions. J Phys Chem Lett 2022; 13:8116-8121. [PMID: 35998327 DOI: 10.1021/acs.jpclett.2c02318] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The fundamental transitions that contribute to the diffuse OH stretching spectrum of water are known to increase in width and intensity with increasing red shift from the free OH frequency. In contrast, the profile of the higher-energy combination band involving the OH stretching and the intramolecular HOH bending modes displays a qualitatively different spectral shape with a much faster falloff on the lower-energy side. We elucidate the molecular origin of this difference by analyzing the shapes of the combination bands in the IR spectra of cryogenically cooled H3O+(H2O)20 and D3O+(D2O)20 clusters. The difference in the shapes of the bands is traced to differences in the dependence of their transition dipole matrix elements on the hydrogen-bonding environment. The fact that individual transitions across the combination band envelope have similar intensities makes it a useful way to determine the participation of various sites in extended H-bonding networks.
Collapse
Affiliation(s)
- Nan Yang
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| | - Rachel M Huchmala
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Mark A Johnson
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
15
|
Palchowdhury S, Mukherjee K, Maroncelli M. Rapid Water Dynamics Structures the OH-Stretching Spectra of Solitary Water in Ionic Liquids and Dipolar Solvents. J Chem Phys 2022; 157:084502. [DOI: 10.1063/5.0107348] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In a recent study [ J. Phys. Chem. B 126, 4584 (2022)] we used infrared spectroscopy to investigate the solvation and dynamics of solitary water in ionic liquids and dipolar solvents. Complex shapes observed for water OH-stretching bands common to all high-polarity solvents were assigned to water in several solvation states. In the present study, classical molecular dynamics simulations of a single water molecule in four ionic liquids and three dipolar solvents were used to test and refine this interpretation. Consistent with past assignments, simulations show solitary water usually donates two hydrogen bonds to distinct solvent molecules. Such symmetrically solvated water produces the primary pair of peaks identified in the OH spectra of water in nearly all solvents. We had further proposed that additional features flanking this main peak are due to asymmetric solvation states, states in which only one OH group makes a hydrogen bond to solvent. Such states were found in significant concentrations in all of the systems simulated. Simulations of the OH stretching spectra using a semiclassical description and the vibrational map developed by Auer and Skinner [ J. Chem. Phys. 128, 224511 (2008)] provided semi-quantitative agreement with experiment. Analysis of species-specific spectra also confirmed assignment of the additional features in the experimental spectra to asymmetrically solvated water. The simulations also showed that rapid water motions cause a marked motional narrowing compared to the inhomogeneous limit, and that this narrowing is largely responsible for making the additional features due to minority solvation states manifest in the spectra.
Collapse
Affiliation(s)
- Souarv Palchowdhury
- The Pennsylvania State University - University Park Campus, United States of America
| | - Kallol Mukherjee
- The Pennsylvania State University - University Park Campus, United States of America
| | - Mark Maroncelli
- Department of Chemsitry, The Pennsylvania State University - University Park Campus, United States of America
| |
Collapse
|
16
|
Weinhold F. High-Density “Windowpane” Coordination Patterns of Water Clusters and Their NBO/NRT Characterization. Molecules 2022; 27:molecules27134218. [PMID: 35807463 PMCID: PMC9268199 DOI: 10.3390/molecules27134218] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Revised: 06/25/2022] [Accepted: 06/27/2022] [Indexed: 02/05/2023] Open
Abstract
Cluster mixture models for liquid water at higher pressures suggest the need for water clusters of higher coordination and density than those commonly based on tetrahedral H-bonding motifs. We show here how proton-ordered water clusters of increased coordination and density can assemble from a starting cyclic tetramer or twisted bicyclic (Möbius-like) heptamer to form extended Aufbau sequences of stable two-, three-, and four-coordinate “windowpane” motifs. Such windowpane clusters exhibit sharply reduced (~90°) bond angles that differ appreciably from the tetrahedral angles of idealized crystalline ice Ih. Computed free energy and natural resonance theory (NRT) bond orders provide quantitative descriptors for the relative stabilities of clusters and strengths of individual coordinative linkages. The unity and consistency of NRT description is demonstrated to extend from familiar supra-integer bonds of the molecular regime to the near-zero bond orders of the weakest linkages in the present H-bond clusters. Our results serve to confirm that H-bonding exemplifies resonance–covalent (fractional) bonding in the sub-integer range and to further discount the dichotomous conceptions of “electrostatics” for intermolecular bonding vs. “covalency” for intramolecular bonding that still pervade much of freshman-level pedagogy and force-field methodology.
Collapse
Affiliation(s)
- Frank Weinhold
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, Madison, WI 53706, USA
| |
Collapse
|
17
|
Mitra S, Khuu T, Choi TH, Huchmala RM, Jordan KD, McCoy AB, Johnson MA. Vibrational Signatures of HNO 3 Acidity When Complexed with Microhydrated Alkali Metal Ions, M +·(HNO 3)(H 2O) n=5 (M = Li, K, Na, Rb, Cs), at 20 K. J Phys Chem A 2022; 126:1640-1647. [PMID: 35249322 DOI: 10.1021/acs.jpca.1c10352] [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
The speciation of strong acids like HNO3 under conditions of restricted hydration is an important factor in the rates of chemical reactions at the air-water interface. Here, we explore the trade-offs at play when HNO3 is attached to alkali ions (Li+-Cs+) with four water molecules in their primary hydration shells. This is achieved by analyzing the vibrational spectra of the M+·(HNO3)(H2O)5 clusters cooled to about 20 K in a cryogenic photofragmentation mass spectrometer. The local acidity of the acidic OH group is estimated by the extent of the red shift in its stretching frequency when attached to a single water molecule. The persistence of this local structural motif (HNO3-H2O) in all of these alkali metal clusters enables us to determine the competition between the effect of the direct complexation of the acid with the cation, which acts to enhance acidity, and the role of the water network in the first hydration shell around the ions, which acts to counter (screen) the intrinsic effect of the ion. Analysis of the vibrational features associated with the acid molecule, as well as those of the water network, reveals how cooperative interactions in the microhydration regime conspire to effectively offset the intrinsic enhancement of HNO3 acidity afforded by attachment to the smaller cations.
Collapse
Affiliation(s)
- Sayoni Mitra
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| | - Thien Khuu
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| | - Tae Hoon Choi
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Rachel M Huchmala
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Kenneth D Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Mark A Johnson
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
18
|
Brünig FN, Geburtig O, Canal AV, Kappler J, Netz RR. Time-Dependent Friction Effects on Vibrational Infrared Frequencies and Line Shapes of Liquid Water. J Phys Chem B 2022; 126:1579-1589. [PMID: 35167754 PMCID: PMC8883462 DOI: 10.1021/acs.jpcb.1c09481] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
From ab initio simulations
of liquid water, the time-dependent
friction functions and time-averaged nonlinear effective bond potentials
for the OH stretch and HOH bend vibrations are extracted. The obtained
friction exhibits not only adiabatic contributions at and below the
vibrational time scales but also much slower nonadiabatic contributions,
reflecting homogeneous and inhomogeneous line broadening mechanisms,
respectively. Intermolecular interactions in liquid water soften both
stretch and bend potentials compared to the gas phase, which by itself
would lead to a red-shift of the corresponding vibrational bands.
In contrast, nonadiabatic friction contributions cause a spectral
blue shift. For the stretch mode, the potential effect dominates,
and thus, a significant red shift when going from gas to the liquid
phase results. For the bend mode, potential and nonadiabatic friction
effects are of comparable magnitude, so that a slight blue shift results,
in agreement with well-known but puzzling experimental findings. The
observed line broadening is shown to be roughly equally caused by
adiabatic and nonadiabatic friction contributions for both the stretch
and bend modes in liquid water. Thus, the quantitative analysis of
the time-dependent friction that acts on vibrational modes in liquids
advances the understanding of infrared vibrational frequencies and
line shapes.
Collapse
|
19
|
Nanayakkara S, Tao Y, Kraka E. Capturing Individual Hydrogen Bond Strengths in Ices via Periodic Local Vibrational Mode Theory: Beyond the Lattice Energy Picture. J Chem Theory Comput 2021; 18:562-579. [PMID: 34928619 DOI: 10.1021/acs.jctc.1c00357] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Local stretching force constants derived from periodic local vibrational modes at the vdW-DF2 density functional level have been employed to quantify the intrinsic hydrogen bond strength of 16 ice polymorphs, ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI, XII, XIII, XIV, XV, XVII, and XIX, that are stable under ambient to elevated pressures. Based on this characterization on 1820 hydrogen bonds, relationships between local stretching force constants and structural parameters such as hydrogen bond length and angle were identified. Moreover, different bond strength distributions, from uniform to inhomogeneous, were observed for the 16 ices and could be explained in relation to different local structural elements within ices, that is, rings, that consist of different hydrogen bond types. In addition, criteria for the classification of hydrogen bonds as strong, intermediate, and weak were introduced. The latter was used to explore a different dimension of the water-ice phase diagram. These findings will provide important guidelines for assessing the credibility of new ice structures.
Collapse
Affiliation(s)
- Sadisha Nanayakkara
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Yunwen Tao
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, United States
| |
Collapse
|
20
|
Investigation of intermolecular interactions in organic solutions by combining two-dimensional correlation Raman spectroscopy and DFT simulation: Example of methanol and Chloralkane. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117147] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
21
|
Kameda Y, Kowaguchi M, Tsutsui K, Amo Y, Usuki T, Ikeda K, Otomo T. Experimental Determination of Relationship between Intramolecular O-D Bond Length and Its Stretching Vibrational Frequency of D 2O Molecule in the Liquid State. J Phys Chem B 2021; 125:11285-11291. [PMID: 34605237 DOI: 10.1021/acs.jpcb.1c07527] [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
Experimental evidence has been obtained for the structure-spectra relationship of hydrogen bonds in aqueous solutions. Intramolecular O-D distance, rOD, has been determined by the least-squares fitting analysis of the neutron interference term in the high-Q region observed for pure D2O and concentrated aqueous solutions. The average O-D stretching frequency, νOD, has been obtained from the position of the center of gravity of the observed ATR-IR O-D stretching band. The linear relationship between rOD and νOD has been confirmed in the liquid state. The slope of dνOD/drOD is evaluated to be -21 000 ± 1000 cm-1 Å-1.
Collapse
Affiliation(s)
- Yasuo Kameda
- Department of Material and Biological Chemistry, Faculty of Science, Yamagata University, 1-4-12, Kojirakawa-machi, Yamagata City, Yamagata 990-8560, Japan
| | - Misaki Kowaguchi
- Department of Material and Biological Chemistry, Faculty of Science, Yamagata University, 1-4-12, Kojirakawa-machi, Yamagata City, Yamagata 990-8560, Japan
| | - Kana Tsutsui
- Department of Material and Biological Chemistry, Faculty of Science, Yamagata University, 1-4-12, Kojirakawa-machi, Yamagata City, Yamagata 990-8560, Japan
| | - Yuko Amo
- Department of Material and Biological Chemistry, Faculty of Science, Yamagata University, 1-4-12, Kojirakawa-machi, Yamagata City, Yamagata 990-8560, Japan
| | - Takeshi Usuki
- Department of Material and Biological Chemistry, Faculty of Science, Yamagata University, 1-4-12, Kojirakawa-machi, Yamagata City, Yamagata 990-8560, Japan
| | - Kazutaka Ikeda
- Institute of Material Structure Science, KEK, Tsukuba, Ibaraki 305-080, Japan
| | - Toshiya Otomo
- Institute of Material Structure Science, KEK, Tsukuba, Ibaraki 305-080, Japan
| |
Collapse
|
22
|
DiRisio RJ, Finney JM, Dzugan LC, Madison LR, McCoy AB. Using Diffusion Monte Carlo Wave Functions to Analyze the Vibrational Spectra of H 7O 3+ and H 9O 4. J Phys Chem A 2021; 125:7185-7197. [PMID: 34433268 DOI: 10.1021/acs.jpca.1c05025] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An approach for evaluating spectra from ground state probability amplitudes (GSPA) obtained from diffusion Monte Carlo (DMC) simulations is extended to improve the description of excited state energies and allow for coupling among vibrational excited states. This approach is applied to studies of the protonated water trimer and tetramer, and their deuterated analogs. These ions provide models for solvated hydronium, and analysis of these spectra provides insights into spectral signatures of proton transfer in aqueous environments. In this approach, we obtain a separable set of internal coordinates from the DMC ground state probability amplitude. A basis is then developed from products of the DMC ground state wave function and low-order polynomials in these internal coordinates. This approach provides a compact basis in which the Hamiltonian and dipole moment matrix are evaluated and used to obtain the spectrum. The resulting spectra are in good agreement with experiment and in many cases provide comparable agreement to the results obtained using much larger basis sets. In addition, the compact basis allows for interpretation of the spectral features and how they evolve with cluster size and deuteration.
Collapse
Affiliation(s)
- Ryan J DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Jacob M Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Laura C Dzugan
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Lindsey R Madison
- Department of Chemistry, Colby College, Waterville, Maine 04901, United States
| | - Anne B McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
23
|
Wang R, Wang J, Zhu Y, Yu F, Yang Y, Wang Z. A Covalent‐Like Feature of Intermolecular Hydrogen Bonding in Energetic Molecules 3,6‐Dihydrazino‐s‐tetrazine (DHT). ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202100179] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Rui Wang
- Institute of Atomic and Molecular Physics Jilin University Changchun 130012 P.R. China
| | - Jia Wang
- College of Information Technology Jilin Normal University Siping 136000 P.R. China
| | - Yu Zhu
- Institute of Atomic and Molecular Physics Jilin University Changchun 130012 P.R. China
| | - Famin Yu
- Institute of Atomic and Molecular Physics Jilin University Changchun 130012 P.R. China
| | - Yanqiang Yang
- Institute of Fluid Physics China Academy of Engineering Physics Mianyang 621900 P.R. China
| | - Zhigang Wang
- Institute of Atomic and Molecular Physics Institute of Theoretical Chemistry Jilin University Changchun 130012 P.R. China
| |
Collapse
|
24
|
Glendening ED, Weinhold F. Pauling's Conceptions of Hybridization and Resonance in Modern Quantum Chemistry. Molecules 2021; 26:4110. [PMID: 34299384 PMCID: PMC8303469 DOI: 10.3390/molecules26144110] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2021] [Revised: 06/30/2021] [Accepted: 07/02/2021] [Indexed: 12/05/2022] Open
Abstract
We employ the tools of natural bond orbital (NBO) and natural resonance theory (NRT) analysis to demonstrate the robustness, consistency, and accuracy with which Linus Pauling's qualitative conceptions of directional hybridization and resonance delocalization are manifested in all known variants of modern computational quantum chemistry methodology.
Collapse
Affiliation(s)
- Eric D. Glendening
- Department of Chemistry and Physics, Indiana State University, Terre Haute, IN 47809, USA;
| | - Frank Weinhold
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, Madison, WI 53706, USA
| |
Collapse
|
25
|
Chakraborty A, Schmahl S, Asmis KR. Isomer-Specific Vibrational Spectroscopy of Microhydrated Lithium Dichloride Anions: Spectral Fingerprint of Solvent-Shared Ion Pairs. Chemphyschem 2021; 22:1036-1041. [PMID: 33783947 PMCID: PMC8252531 DOI: 10.1002/cphc.202100170] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Revised: 03/28/2021] [Indexed: 12/21/2022]
Abstract
The vibrational spectroscopy of lithium dichloride anions microhydrated with one to three water molecules, [LiCl2 (H2 O)1-3 ]- , is studied in the OH stretching region (3800-2800 cm-1 ) using isomer-specific IR/IR double-resonance population labelling experiments. The spectroscopic fingerprints of individual isomers can only be unambiguously assigned after anharmonic effects are considered, but then yield molecular level insight into the onset of salt dissolution in these gas phase model systems. Based on the extent of the observed frequency shifts ΔνOH of the hydrogen-bonded OH stretching oscillators solvent-shared ion pair motifs (<3200 cm-1 ) can be distinguished from intact-core structures (>3200 cm-1 ). The characteristic fingerprint of a water molecule trapped directly in-between two ions of opposite charge provides an alternative route to evaluate the extent of ion pairing in aqueous electrolyte solutions.
Collapse
Affiliation(s)
- Arghya Chakraborty
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig, Linnéstrasse 2, D-04103, Leipzig, Germany
| | - Sonja Schmahl
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig, Linnéstrasse 2, D-04103, Leipzig, Germany
| | - Knut R Asmis
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig, Linnéstrasse 2, D-04103, Leipzig, Germany
| |
Collapse
|
26
|
Zeng HJ, Johnson MA. Demystifying the Diffuse Vibrational Spectrum of Aqueous Protons Through Cold Cluster Spectroscopy. Annu Rev Phys Chem 2021; 72:667-691. [PMID: 33646816 DOI: 10.1146/annurev-physchem-061020-053456] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The ease with which the pH is routinely determined for aqueous solutions masks the fact that the cationic product of Arrhenius acid dissolution, the hydrated proton, or H+(aq), is a remarkably complex species. Here, we review how results obtained over the past 30 years in the study of H+⋅(H2O)n cluster ions isolated in the gas phase shed light on the chemical nature of H+(aq). This effort has also revealed molecular-level aspects of the Grotthuss relay mechanism for positive-charge translocation in water. Recently developed methods involving cryogenic cooling in radiofrequency ion traps and the application of two-color, infrared-infrared (IR-IR) double-resonance spectroscopy have established a clear picture of how local hydrogen-bond topology drives the diverse spectral signatures of the excess proton. This information now enables a new generation of cluster studies designed to unravel the microscopic mechanics underlying the ultrafast relaxation dynamics displayed by H+(aq).
Collapse
Affiliation(s)
- Helen J Zeng
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, USA;
| | - Mark A Johnson
- Sterling Chemistry Laboratory, Yale University, New Haven, Connecticut 06520, USA;
| |
Collapse
|
27
|
Torii H, Ukawa R. Role of Intermolecular Charge Fluxes in the Hydrogen-Bond-Induced Frequency Shifts of the OH Stretching Mode of Water. J Phys Chem B 2021; 125:1468-1475. [PMID: 33506673 DOI: 10.1021/acs.jpcb.0c11461] [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/29/2022]
Abstract
The relation between the vibrational properties and the electrostatic situations of the vibrating functional group is useful to predict vibrational spectroscopic features based on, for example, classical molecular dynamics of liquids or biomolecular systems, but to pursue its generality or the extent of applicability, it is required to understand the mechanisms giving rise to it. Here such an analysis is carried out for the OH stretching mode of water. By examining the correlations among various (structural, vibrational, and electrostatic) properties and by analyzing the spatial characteristics of the behavior of electrons occurring upon the vibration, it is shown that the dependence of the vibrational frequency and the dipole derivative of the OH stretching mode on the electric field is not of purely electrostatic origin, and the delocalized electronic motions occurring with this mode, called intermolecular charge fluxes, related to both the dipole first and second derivatives play important roles.
Collapse
Affiliation(s)
- Hajime Torii
- Department of Applied Chemistry and Biochemical Engineering, Faculty of Engineering, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan.,Department of Optoelectronics and Nanostructure Science, Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
| | - Ryota Ukawa
- Department of Applied Chemistry and Biochemical Engineering, Faculty of Engineering, Shizuoka University, 3-5-1 Johoku, Naka-ku, Hamamatsu 432-8561, Japan
| |
Collapse
|
28
|
Herbert JM, Carter-Fenk K. Electrostatics, Charge Transfer, and the Nature of the Halide-Water Hydrogen Bond. J Phys Chem A 2021; 125:1243-1256. [PMID: 33502859 DOI: 10.1021/acs.jpca.0c11356] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Binary halide-water complexes X-(H2O) are examined by means of symmetry-adapted perturbation theory, using charge-constrained promolecular reference densities to extract a meaningful charge-transfer component from the induction energy. As is known, the X-(H2O) potential energy surface (for X = F, Cl, Br, or I) is characterized by symmetric left and right hydrogen bonds separated by a C2v-symmetric saddle point, with a tunneling barrier height that is <2 kcal/mol except in the case of F-(H2O). Our analysis demonstrates that the charge-transfer energy is correspondingly small (<2 kcal/mol except for X = F), considerably smaller than the electrostatic interaction energy. Nevertheless, charge transfer plays a crucial role determining the conformational preferences of X-(H2O) and provides a driving force for the formation of quasi-linear X··· H-O hydrogen bonds. Charge-transfer energies correlate well with measured O-H vibrational redshifts for the halide-water complexes and also for OH-(H2O) and NO2-(H2O), providing some indication of a general mechanism.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
29
|
Dereka B, Yu Q, Lewis NHC, Carpenter WB, Bowman JM, Tokmakoff A. Crossover from hydrogen to chemical bonding. Science 2021; 371:160-164. [DOI: 10.1126/science.abe1951] [Citation(s) in RCA: 66] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Accepted: 10/20/2020] [Indexed: 11/02/2022]
Affiliation(s)
- Bogdan Dereka
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637, USA
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, GA 30322, USA
| | - Nicholas H. C. Lewis
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637, USA
| | - William B. Carpenter
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637, USA
| | - Joel M. Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, GA 30322, USA
| | - Andrei Tokmakoff
- Department of Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, IL 60637, USA
| |
Collapse
|
30
|
Knorke H, Li H, Warneke J, Liu ZF, Asmis KR. Cryogenic ion trap vibrational spectroscopy of the microhydrated sulfate dianions SO 42-(H 2O) 3-8. Phys Chem Chem Phys 2020; 22:27732-27745. [PMID: 33242322 DOI: 10.1039/d0cp04386a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Infrared photodissociation spectra of the D2-tagged microhydrated sulfate dianions with three to eight water molecules are presented over a broad spectral range that covers the OH stretching and H2O bending modes of the solvent molecules at higher energies, the sulfate stretching modes of the solute at intermediate energies and the intermolecular solute librational modes at the lowest energies. A low ion temperature combined with messenger-tagging ensures well-resolved vibrational spectra that allow for structure assignments based on a comparison to harmonic and anharmonic IR spectra from density functional theory (DFT) calculations. DFT ab initio molecular dynamics simulations are required to disentangle the broad and complex spectral signatures of microhydrated sulfate dianions in the OH stretching region and to identify systematic trends in the correlation of the strength and evolution of the solute-solvent and solvent-solvent interactions with cluster size. The onset for the formation of the second solvation shell is observed for n = 8.
Collapse
Affiliation(s)
- Harald Knorke
- Wilhelm-Ostwald-Institut für Physikalische und Theoretische Chemie, Universität Leipzig, Linnéstr. 2, 04103 Leipzig, Germany.
| | | | | | | | | |
Collapse
|
31
|
Mandal N, Datta A. Molecular designs for expanding the limits of ultralong C-C bonds and ultrashort HH non-bonded contacts. Chem Commun (Camb) 2020; 56:15377-15386. [PMID: 33210669 DOI: 10.1039/d0cc06690g] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recent experiments have reported the formation of very long C-C bonds (dC-C > 1.80 Å) and very short HH non-bonded contacts (dHH < 1.5 Å) in several sets of molecules. Both these rare phenomena arise due to specific donor-acceptor interactions and London dispersion interactions respectively. Favorable negative hyperconjugation, namely H2N(lone-pair) →σ*(C-C), creates an ultralong C-C bond in diamino-o-carborane with dC-C > 1.829 Å and a planar amine reminiscent of a transition-state like structure for ammonia inversion. The small and narrow barrier favours rapid inversion through quantum mechanical tunnelling (QMT) and produces a translationally averaged planar amine as observed in the experiments. On the other hand, designing specific confined molecular cavities or chambers like in,in-bis(hydrosilane) or its germanane analogs furnishes an ultrashort HH distance = 1.47 Å and 1.38 Å respectively. The predisposition of such closely placed HH contacts arises from the rather effective attractive dispersion interactions between them. Controlling the strength of the dispersion interactions provides a rich landscape for realizing such close HH distances. Molecular design ably assisted by computational modeling to further tune these interactions provides new avenues to break the glass-ceilings of ultralong C-C bonds or ultrashort HH contacts. Dispersion-corrected DFT calculations and ab initio molecular dynamics simulations generate a large library of such unique features in a diverse class of molecules. This feature article highlights the design principles to realize hitherto longest C-C bonds/shortest HH contacts.
Collapse
Affiliation(s)
- Nilangshu Mandal
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A and 2B Raja S. C. Mullick Road, Jadavpur - 700032, Kolkata, West Bengal, India.
| | | |
Collapse
|
32
|
Finney JM, DiRisio RJ, McCoy AB. Guided Diffusion Monte Carlo: A Method for Studying Molecules and Ions That Display Large Amplitude Vibrational Motions. J Phys Chem A 2020; 124:9567-9577. [DOI: 10.1021/acs.jpca.0c07181] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jacob M. Finney
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Ryan J. DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B. McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
33
|
Biswas S, Mallik BS. Probing the vibrational dynamics of amide bands of N-methylformamide, N, N-dimethylacetamide, and N-methylacetamide in water. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.113001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
34
|
Carpenter WB, Yu Q, Hack JH, Dereka B, Bowman JM, Tokmakoff A. Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations. J Chem Phys 2020; 153:124506. [PMID: 33003749 DOI: 10.1063/5.0020279] [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/18/2023] Open
Abstract
The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H5O2 + and Eigen H9O4 + ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure-spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm-1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. Finally, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.
Collapse
Affiliation(s)
- William B Carpenter
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - John H Hack
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Bogdan Dereka
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
35
|
Kananenka AA, Skinner JL. Unusually strong hydrogen bond cooperativity in particular (H 2O) 20 clusters. Phys Chem Chem Phys 2020; 22:18124-18131. [PMID: 32761035 DOI: 10.1039/d0cp02343d] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Drawing upon an intuitive charge-transfer-based picture of hydrogen bonding, we demonstrate that cooperativity effects acting in concert can lead to unusually strong hydrogen bonds in neutral water clusters. The structure, vibrational, and NMR properties of a (H2O)20 pentagonal dodecahedron cluster containing such a strong hydrogen bond were studied using second-order perturbation theory and density functional theory. The hydrogen bond length was found to be shorter than 2.50 Å. A large redshift of over 2000 cm-1 with respect to the isolated water molecule was predicted for the OH stretching frequency of the donor water molecule. A large downfield shift to 13.5 ppm of the isotropic part of the 1H magnetic shielding tensor together with an unusually large shielding anisotropy of 49.9 ppm was obtained. The hydrogen bond energy was calculated using symmetry-adapted perturbation theory and was found to be more than three times stronger than a typical hydrogen bond in liquid water.
Collapse
Affiliation(s)
- Alexei A Kananenka
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA.
| | - J L Skinner
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
36
|
Tao Y, Qiu Y, Zou W, Nanayakkara S, Yannacone S, Kraka E. In Situ Assessment of Intrinsic Strength of X-I⋯OA-Type Halogen Bonds in Molecular Crystals with Periodic Local Vibrational Mode Theory. Molecules 2020; 25:molecules25071589. [PMID: 32235623 PMCID: PMC7181175 DOI: 10.3390/molecules25071589] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2020] [Revised: 03/20/2020] [Accepted: 03/24/2020] [Indexed: 12/03/2022] Open
Abstract
Periodic local vibrational modes were calculated with the rev-vdW-DF2 density functional to quantify the intrinsic strength of the X-I⋯OA-type halogen bonding (X = I or Cl; OA: carbonyl, ether and N-oxide groups) in 32 model systems originating from 20 molecular crystals. We found that the halogen bonding between the donor dihalogen X-I and the wide collection of acceptor molecules OA features considerable variations of the local stretching force constants (0.1–0.8 mdyn/Å) for I⋯O halogen bonds, demonstrating its powerful tunability in bond strength. Strong correlations between bond length and local stretching force constant were observed in crystals for both the donor X-I bonds and I⋯O halogen bonds, extending for the first time the generalized Badger’s rule to crystals. It is demonstrated that the halogen atom X controlling the electrostatic attraction between the σ-hole on atom I and the acceptor atom O dominates the intrinsic strength of I⋯O halogen bonds. Different oxygen-containing acceptor molecules OA and even subtle changes induced by substituents can tweak the n→σ∗(X-I) charge transfer character, which is the second important factor determining the I⋯O bond strength. In addition, the presence of the second halogen bond with atom X of the donor X-I bond in crystals can substantially weaken the target I⋯O halogen bond. In summary, this study performing the in situ measurement of halogen bonding strength in crystalline structures demonstrates the vast potential of the periodic local vibrational mode theory for characterizing and understanding non-covalent interactions in materials.
Collapse
Affiliation(s)
- Yunwen Tao
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, USA; (Y.T.); (S.N.); (S.Y.)
| | - Yue Qiu
- Grimwade Centre for Cultural Materials Conservation, School of Historical and Philosophical Studies, Faculty of Arts, University of Melbourne, Parkville, VIC 3052, Australia;
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi’an 710127, China;
| | - Sadisha Nanayakkara
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, USA; (Y.T.); (S.N.); (S.Y.)
| | - Seth Yannacone
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, USA; (Y.T.); (S.N.); (S.Y.)
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, USA; (Y.T.); (S.N.); (S.Y.)
- Correspondence:
| |
Collapse
|
37
|
Morawietz T, Urbina AS, Wise PK, Wu X, Lu W, Ben-Amotz D, Markland TE. Hiding in the Crowd: Spectral Signatures of Overcoordinated Hydrogen-Bond Environments. J Phys Chem Lett 2019; 10:6067-6073. [PMID: 31549833 DOI: 10.1021/acs.jpclett.9b01781] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Molecules with an excess number of hydrogen-bonding partners play a crucial role in fundamental chemical processes, ranging from anomalous diffusion in supercooled water to transport of aqueous proton defects and ordering of water around hydrophobic solutes. Here we show that overcoordinated hydrogen-bond environments can be identified in both the ambient and supercooled regimes of liquid water by combining experimental Raman multivariate curve resolution measurements and machine learning accelerated quantum simulations. In particular, we find that OH groups appearing in spectral regions usually associated with non-hydrogen-bonded species actually correspond to hydrogen bonds formed in overcoordinated environments. We further show that only these species exhibit a turnover in population as a function of temperature, which is robust and persists under both constant pressure and density conditions. This work thus provides a new tool to identify, interpret, and elucidate the spectral signatures of crowded hydrogen-bond networks.
Collapse
Affiliation(s)
- Tobias Morawietz
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | - Andres S Urbina
- Department of Chemistry , Purdue University , West Lafayette , Indiana 47907 , United States
| | - Patrick K Wise
- Department of Chemistry , Purdue University , West Lafayette , Indiana 47907 , United States
| | - Xiangen Wu
- College of Marine Science and Technology , China University of Geosciences , Wuhan 430074 , China
| | - Wanjun Lu
- College of Marine Science and Technology , China University of Geosciences , Wuhan 430074 , China
| | - Dor Ben-Amotz
- Department of Chemistry , Purdue University , West Lafayette , Indiana 47907 , United States
| | - Thomas E Markland
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| |
Collapse
|
38
|
Lacroix MR, Liu Y, Strauss SH. Hydrated Metal Ion Salts of the Weakly Coordinating Fluoroanions PF 6-, TiF 62-, B 12F 122-, Ga(C 2F 5) 4-, B(3,5-C 6H 3(CF 3) 2) 4-, and Al(OC(CF 3) 3) 4-. In Search of the Weakest HOH···F Hydrogen Bonds. Inorg Chem 2019; 58:14900-14911. [PMID: 31617354 DOI: 10.1021/acs.inorgchem.9b02646] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
FTIR spectra of microcrystalline samples of 11 metal ion salt hydrates of a variety of weakly coordinating fluoroanions are reported. The compounds studied were Li(H2O)4(Al(OC(CF3)3)4), Li(H2O)(B(3,5-C6H3(CF3)2)4), Li(H2O)n(Ga(C2F5)4), Li(H2O)(PF6), Li2(H2O)2(TiF6), Li2(H2O)4(B12F12), Na(H2O)(PF6), Na2(H2O)2(B12F12), K2(H2O)2(B12F12), Rb2(H2O)2(B12F12), Cs2(H2O)(B12F12), and their partially or completely deuterated isotopologs and isotopomers. The O-D···F hydrogen bonds in Li(HOD)(H2O)3(Al(OC(CF3)3)4) (ν(OD) = 2706 cm-1), Li(HOD)(B(3,5-C6H3(CF3)2)4) (ν(OD) = 2705 cm-1), and Li(HOD)(H2O)n(Ga(C2F5)4) (ν(OD) = 2697 cm-1) rival HOD absorbed in polyvinylidene difluoride (ν(OD) = 2696 cm-1) and HOD···FCH3 in a frozen Ar matrix (ν(OD) = 2685 cm-1) for the weakest hydrogen bonds between a water molecule and an F atom in any compound. As weak as they are, minor differences in O-H···F hydrogen bonds in the same fluoroanion salt can be distinguished spectroscopically. Uncoupled HOD molecules in asymmetric F···HOD···F' hydrogen bonding environments in Rb+, Cs+, Mg2+, and Co2+ hydrates of B12F122- gave rise to two observable ν(OD) bands even though the two R(O···F) distances differ by only 0.010(4) Å (Mg2+), 0.033(2) Å (Co2+), 0.074(4) Å (Rb+), and 0.106(6) Å (Cs+). A plot of ν(OD) for hydrates with a single uncoupled HOD molecule per metal ion (e.g., Li(HOD)(H2O)3(Al(OC(CF3)3)4)) vs R(O···F) distance from single-crystal X-ray or neutron diffraction structures was prepared. The ν(OD) values range from 2305 to 2706 cm-1 and the R(O···F) distances range from 2.58 to 3.17 Å. The plot consists of 53 {ν(OD), R(O···F)} data points, 23 of which are new and have ν(OD) > 2600 cm-1, in contrast to a 1994 ν(OD) vs R(O···F) plot with 28 data points, none of which had ν(OD) > 2600 cm-1. There is a clear and significant difference between the new ν(OD) vs R(O···F) plot and a literature ν(OD) vs R(O···O) plot for hydrates containing O-D···O hydrogen bonds. For a given ν(OD) stretching frequency, the exponential regression curves show that R(O···F) is typically 0.1-0.2 Å shorter than R(O···O), in harmony with the lower basicity and smaller size of F atoms vs O atoms.
Collapse
Affiliation(s)
- Matthew R Lacroix
- Department of Chemistry , Colorado State University , Fort Collins , Colorado 80523 , United States
| | - Yong Liu
- Department of Chemistry , University of Colorado at Denver , Denver , Colorado 80217 , United States
| | - Steven H Strauss
- Department of Chemistry , Colorado State University , Fort Collins , Colorado 80523 , United States
| |
Collapse
|
39
|
Michalczyk M, Zierkiewicz W, Wysokiński R, Scheiner S. Theoretical Studies of IR and NMR Spectral Changes Induced by Sigma-Hole Hydrogen, Halogen, Chalcogen, Pnicogen, and Tetrel Bonds in a Model Protein Environment. Molecules 2019; 24:E3329. [PMID: 31547416 PMCID: PMC6767630 DOI: 10.3390/molecules24183329] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Revised: 09/09/2019] [Accepted: 09/10/2019] [Indexed: 01/27/2023] Open
Abstract
Various types of σ-hole bond complexes were formed with FX, HFY, H2FZ, and H3FT (X = Cl, Br, I; Y = S, Se, Te; Z = P, As, Sb; T = Si, Ge, Sn) as Lewis acid. In order to examine their interactions with a protein, N-methylacetamide (NMA), a model of the peptide linkage was used as the base. These noncovalent bonds were compared by computational means with H-bonds formed by NMA with XH molecules (X = F, Cl, Br, I). In all cases, the A-F bond, which lies opposite the base and is responsible for the σ-hole on the A atom (A refers to the bridging atom), elongates and its stretching frequency undergoes a shift to the red with a band intensification, much as what occurs for the X-H bond in a H-bond (HB). Unlike the NMR shielding decrease seen in the bridging proton of a H-bond, the shielding of the bridging A atom is increased. The spectroscopic changes within NMA are similar for H-bonds and the other noncovalent bonds. The C=O bond of the amide is lengthened and its stretching frequency red-shifted and intensified. The amide II band shifts to higher frequency and undergoes a small band weakening. The NMR shielding of the O atom directly involved in the bond rises, whereas the C and N atoms both undergo a shielding decrease. The frequency shifts of the amide I and II bands of the base as well as the shielding changes of the three pertinent NMA atoms correlate well with the strength of the noncovalent bond.
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.
| | - Wiktor Zierkiewicz
- Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland.
| | - Rafał Wysokiński
- 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, UT 84322-0300, USA.
| |
Collapse
|
40
|
Yu Q, Carpenter WB, Lewis NHC, Tokmakoff A, Bowman JM. High-Level VSCF/VCI Calculations Decode the Vibrational Spectrum of the Aqueous Proton. J Phys Chem B 2019; 123:7214-7224. [PMID: 31361141 DOI: 10.1021/acs.jpcb.9b05723] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The hydrated excess proton is a common species in aqueous chemistry, which complexes with water in a variety of structures. The infrared spectrum of the aqueous proton is particularly sensitive to this array of structures, which manifests as continuous IR absorption from 1000 to 3000 cm-1 known as the "proton continuum". Because of the extreme breadth of the continuum and strong anharmonicity of the involved vibrational modes, this spectrum has eluded straightforward interpretation and simulation. Using protonated water hexamer clusters from reactive molecular dynamics trajectories, and focusing on their central H+(H2O)2 structures' spectral contribution, we reproduce the linear IR spectrum of the aqueous proton with a high-level local monomer quantum method and highly accurate many-body potential energy surface. The accuracy of this approach is first verified in the vibrational spectra of the two isomers of the protonated water hexamer in the gas phase. We then apply this approach to 800 H+(H2O)6 clusters, also written as [H+(H2O)2](H2O)4, drawn from multistate empirical valence bond simulations of the bulk liquid to calculate the infrared spectrum of the aqueous proton complex. Incorporation of anharmonic effects to the vibrational potential and quantum mechanical treatment of the proton produces a better agreement to the infrared spectrum compared to that of the double-harmonic approximation. We assess the correlation of the proton stretching mode with different atomistic coordinates, finding the best correlation with ⟨ROH⟩, the expectation value of the proton-oxygen distance ROH. We also decompose the IR spectrum based on normal mode vibrations and ⟨ROH⟩ to provide insight on how different frequency regions in the continuum report on different configurations, vibrational modes, and mode couplings.
Collapse
Affiliation(s)
- Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| | - William B Carpenter
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Nicholas H C Lewis
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Andrei Tokmakoff
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| |
Collapse
|
41
|
Lee VGM, McCoy AB. An Efficient Approach for Studies of Water Clusters Using Diffusion Monte Carlo. J Phys Chem A 2019; 123:8063-8070. [DOI: 10.1021/acs.jpca.9b06444] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Victor G. M. Lee
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Anne B. McCoy
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|