We report vibrational spectra of the H 2 -tagged, cryogenically cooled X − · HOCl (X = Cl, Br, and I) ion–molecule complexes and analyze the resulting band patterns with electronic structure calculations and an anharmonic theoretical treatment of nuclear motions on extended potential energy surfaces. The complexes are formed by “ligand exchange” reactions of X − · (H 2 O) n clusters with HOCl molecules at low pressure (∼10 −2 mbar) in a radio frequency ion guide. The spectra generally feature many bands in addition to the fundamentals expected at the double harmonic level. These “extra bands” appear in patterns that are similar to those displayed by the X − · HOD analogs, where they are assigned to excitations of nominally IR forbidden overtones and combination bands. The interactions driving these features include mechanical and electronic anharmonicities. Particularly intense bands are observed for the v = 0 → 2 transitions of the out-of-plane bending soft modes of the HOCl molecule relative to the ions. These involve displacements that act to break the strong H-bond to the ion, which give rise to large quadratic dependences of the electric dipoles (electronic anharmonicities) that drive the transition moments for the overtone bands. On the other hand, overtone bands arising from the intramolecular OH bending modes of HOCl are traced to mechanical anharmonic coupling with the v = 1 level of the OH stretch (Fermi resonances). These interactions are similar in strength to those reported earlier for the X − · HOD complexes.
more »
« less
Probing radical–molecule interactions with a second generation energy decomposition analysis of DFT calculations using absolutely localized molecular orbitals
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, H 2 O⋯˙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(H 2 O) and ˙OH(H 2 O) 2 ), and the pre-reactive complexes and transition states in the ˙OH + HCHO and ˙OH + CH 3 CHO 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.
more »
« less
- NSF-PAR ID:
- 10198773
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 22
- Issue:
- 23
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 12867 to 12885
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The discovery of singular organic radical ligands is a formidable challenge due to high reactivity arising from the unpaired electron. Matching radical ligands with metal ions to engender magnetic coupling is crucial for eliciting preeminent physical properties such as conductivity and magnetism that are crucial for future technologies. The metal-radical approach is especially important for the lanthanide ions exhibiting deeply buried 4f-orbitals. The radicals must possess a high spin density on the donor atoms to promote strong coupling. Combining diamagnetic 89 Y ( I = 1/2) with organic radicals allows for invaluable insight into the electronic structure and spin-density distribution. This approach is hitherto underutilized, possibly owing to the challenging synthesis and purification of such molecules. Herein, evidence of an unprecedented bisbenzimidazole radical anion (Bbim 3− ˙) along with its metalation in the form of an yttrium complex, [K(crypt-222)][(Cp* 2 Y) 2 (μ-Bbim˙)] is provided. Access of Bbim 3− ˙ was feasible through double-coordination to the Lewis acidic metal ion and subsequent one-electron reduction, which is remarkable as Bbim 2− was explicitly stated to be redox-inactive in closed-shell complexes. Two molecules containing Bbim 2− (1) and Bbim 3− ˙ (2), respectively, were thoroughly investigated by X-ray crystallography, NMR and UV/Vis spectroscopy. Electrochemical studies unfolded a quasi-reversible feature and emphasize the role of the metal centre for the Bbim redox-activity as neither the free ligand nor the Bbim 2− complex led to analogous CV results. Excitingly, a strong delocalization of the electron density through the Bbim 3− ˙ ligand was revealed via temperature-dependent EPR spectroscopy and confirmed through DFT calculations and magnetometry, rendering Bbim 3− ˙ an ideal candidate for single-molecule magnet design.more » « less
-
null (Ed.)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, N 2 , and BF bound to a [Ru(II)(NH 3 ) 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.more » « less
-
The reactions of thioformaldehyde (H 2 CS) with OH radicals and assisted by a single water molecule have been investigated using high level ab initio quantum chemistry calculations. The H 2 CS + ˙OH reaction can in principle proceed through: (1) abstraction, and (2) addition pathways. The barrier height for the addition reaction in the absence of a catalyst was found to be −0.8 kcal mol −1 , relative to the separated reactants, which has a ∼1.0 kcal mol −1 lower barrier than the abstraction channel. The H 2 CS + ˙OH reaction assisted by a single water molecule reduces the barrier heights significantly for both the addition and abstraction channels, to −5.5 and −6.7 kcal mol −1 respectively, compared to the un-catalyzed H 2 CS + ˙OH reaction. These values suggest that water lowers the barriers by ∼6.0 kcal mol −1 for both reaction paths. The rate constants for the H 2 CS⋯H 2 O + ˙OH and OH⋯H 2 O + H 2 CS bimolecular reaction channels were calculated using Canonical Variational Transition state theory (CVT) in conjunction with the Small Curvature Tunneling (SCT) method over the atmospherically relevant temperatures between 200 and 400 K. Rate constants for the H 2 CS + ˙OH reaction paths for comparison with the H 2 CS + ˙OH + H 2 O reaction in the same temperature range were also computed. The results suggest that the rate of the H 2 CS + ˙OH + H 2 O reaction is slower than that of the H 2 CS + ˙OH reaction by ∼1–4 orders of magnitude in the temperatures between 200 and 400 K. For example, at 300 K, the rates of the H 2 CS + ˙OH + H 2 O and H 2 CS + ˙OH reactions were found to be 2.2 × 10 −8 s −1 and 6.4 × 10 −6 s −1 , respectively, calculated using [OH] = 1.0 × 10 6 molecules cm −3 , and [H 2 O] = 8.2 × 10 17 molecules cm −3 (300 K, RH 100%) atmospheric conditions. Electronic structure calculations on the H 2 C(OH)S˙ product in the presence of 3 O 2 were also performed. The results show that H 2 CS is removed from the atmosphere primarily by reacting with ˙OH and O 2 to form thioformic acid, HO 2 , formaldehyde, and SO 2 as the main end products.more » « less
-
null (Ed.)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 CO 2 reduction catalysis. For [FeTPP(CO 2 -κC)] 2− (TPP = tetraphenylporphyrin), we reveal that two ortho substituents which yield most favorable CO 2 binding, –N(CH 3 ) 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 CO 2 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 CO 2 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 CO 2 , 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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/d0cp01933j
