1
|
Nanni L. Investigating proton tunneling dynamics in the time-dependent Schrödinger equation. J Comput Chem 2024; 45:1614-1623. [PMID: 38523084 DOI: 10.1002/jcc.27348] [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: 12/29/2023] [Revised: 03/02/2024] [Accepted: 03/07/2024] [Indexed: 03/26/2024]
Abstract
Understanding the temporal evolution of the wave function in the proton tunneling reactions allows us to make theoretical predictions on the possible femtosecond spectroscopy patterns. However, the analytical solution of the time-dependent Schrödinger equation of complex molecules is challenging and requires a high computational cost. In this study, we solve the time-dependent Schrödinger equation using the Fourier grid Hamiltonian method, highlighting its simplicity of calculation, even for multidimensional tunneling reactions. The obtained model is applied to studying malonaldehyde's two-dimensional intramolecular proton tunneling, comparing the results with those obtained using other computational methods.
Collapse
Affiliation(s)
- Luca Nanni
- Department of Physics, University of Ferrara, Ferrara, Italy
| |
Collapse
|
2
|
Asatryan R, Hudzik J, Swihart M. Intramolecular Catalytic Hydrogen Atom Transfer (CHAT). J Phys Chem A 2024; 128:2169-2190. [PMID: 38451855 DOI: 10.1021/acs.jpca.3c06794] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2024]
Abstract
Intramolecular catalysis (IntraCat) is the acceleration of a process at one site of a molecule catalyzed by a functional group in the same molecule; an external agent such as a solvent typically facilitates it. Here, we report a general first-principles-based IntraCat mechanism, which strictly occurs within a single molecule with no coreagent being involved─we call it intramolecular catalytic transfer of hydrogen atoms (CHAT). A reactive part of a molecule (chat catalyst moiety or chat agent, represented by -OOH, -COOH, -SH, -CH2OH, -HPO4, or another bifunctional H-donor/acceptor group) catalyzes an interconversion process, such as keto-enol or amino-imino tautomerization, and cyclization in the same molecule, while being regenerated in the process. It can thus be regarded as an intramolecular version of the intermolecular H atom transfer processes mediated by an external molecular catalyst, e.g., dihydrogen, water, or a carboxylic acid. Earlier, we proposed a general mechanistic systematization of intermolecular processes, illustrated in the simplest case of the H2-mediated reactions classified as dihydrogen catalysis [Asatryan, R.; et al. Catal. Rev.: Sci. Eng., 2014, 56, 403-475]. Following this systematization, the CHAT catalysis belongs to the category of relay transfer of H atoms, albeit in an intramolecular manner. A broader class of intramolecular processes includes all types of H-transfer reactions stimulated by an H-migration, which we call self-catalyzed H atom transfer (SC-HAT). The CHAT mechanism comprises a subset of SC-HAT in which the catalytic moiety is regenerated (i.e., acts as a true catalyst and not a reagent). We provide several characteristic examples of CHAT mechanism based on detailed analysis of the corresponding potential energy surfaces. All such cases showed a dramatically reduced activation barrier relative to the corresponding uncatalyzed H-transfer reactions. For example, we show that CHAT can facilitate long-range H-migration in larger molecules and can occur multiple times in one molecule with multiple interconverting groups. It also facilitates amino-imino tautomerization of unsaturated GABA-analogues and peptides, as well as intramolecular cyclization processes to form heterocycles, e.g., oxygenated rings. CHAT pathways may also explain the pH-dependent increase of mutarotation rate of glucose-6-phosphate demonstrated in pioneering experiments that introduced the classical IntraCat concept. In addition, we identify a ground electronic state CHAT pathway as an alternative to the UV-promoted long-range molecular crane keto-enol conversion with a remarkably low activation energy. To initially assess the possible impact of the new keto-enol conversion pathway on combustion of n-alkanes, we present a detailed kinetic analysis of isomerization and decomposition of pentane-2,4-ketohydroperoxide (2,4-KHP). The results are compared with key alternative reactions, including direct dissociation and Korcek channels (for which a new alkyl group migration channel is also identified), revealing the competitiveness of the CHAT pathway across a range of conditions. Taken together, this work provides insight into a general class of reaction pathways that has not previously being systematically considered and that may occur in a broad range of contexts from combustion to atmospheric chemistry to biochemistry.
Collapse
Affiliation(s)
- Rubik Asatryan
- Department of Chemical and Biological Engineering, and Center for Hybrid Rocket Exascale Simulation Technology (CHREST), University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Jason Hudzik
- Department of Chemical and Biological Engineering, and Center for Hybrid Rocket Exascale Simulation Technology (CHREST), University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Mark Swihart
- Department of Chemical and Biological Engineering, and Center for Hybrid Rocket Exascale Simulation Technology (CHREST), University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
3
|
Saha J, Banerjee S, Malo S, Das AK, Das I. A Torquoselective Thermal 6π-Electrocyclization Approach to 1,4-Cyclohexadienes via Solvent-Aided Proton Transfer: Experimental and Theoretical Studies. Chemistry 2024; 30:e202304009. [PMID: 38179806 DOI: 10.1002/chem.202304009] [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: 12/01/2023] [Revised: 12/26/2023] [Accepted: 01/05/2024] [Indexed: 01/06/2024]
Abstract
The thermal 6π-electrocyclization of hexatriene typically delivers 1,3-cyclohexadiene (1,3-CHD). However, there is only limited success in directly synthesizing 1,4-cyclohexadiene (1,4-CHD) using such an approach, probably due to the difficulty in realizing thermally-forbidden 1,3-hydride shift after electrocyclic ring closure. The present study shows that by heating (2E,4E,6E)-hexatrienes bearing ester or ketone substituents at the C1-position in a mixture of toluene/MeOH or EtOH (2 : 1) solvents at 90-100 °C, 1,4-CHDs can be selectively synthesized. This is achieved through a torquoselective disrotatory 6π-electrocyclic ring closure followed by a proton-transfer process. The success of this method depends on the polar protic solvent-assisted intramolecular proton transfer from 1,3-CHD to 1,4-CHD, which has been confirmed by deuterium-labeling experiments. There are no reports to date for such a solvent-assisted isomerization. Density functional theory (DFT) studies have suggested that forming 1,3-CHD and subsequent isomerization is a thermodynamically feasible process, regardless of the functional groups involved. Two possible successive polar solvent-assisted proton-transfer pathways have been identified for isomerization.
Collapse
Affiliation(s)
- Jayanta Saha
- Organic and Medicinal Chemistry Division, CSIR-Indian Institute of Chemical BiologyJadavpur, Kolkata, 700032, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, 201002, India
| | - Soumadip Banerjee
- School of Mathematical and Computational Sciences, Indian Association for the Cultivation of Science Jadavpur, Kolkata, 700032, India
| | - Sidhartha Malo
- Organic and Medicinal Chemistry Division, CSIR-Indian Institute of Chemical BiologyJadavpur, Kolkata, 700032, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, 201002, India
| | - Abhijit Kumar Das
- School of Mathematical and Computational Sciences, Indian Association for the Cultivation of Science Jadavpur, Kolkata, 700032, India
| | - Indrajit Das
- Organic and Medicinal Chemistry Division, CSIR-Indian Institute of Chemical BiologyJadavpur, Kolkata, 700032, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, 201002, India
| |
Collapse
|
4
|
Mauksch M. Spontaneous emergence of enantioenriched chiral aldol reaction products from Achiral precursors in solution and origin of biological homochirality of sugars: a first-principles study. Phys Chem Chem Phys 2023; 25:1734-1754. [PMID: 36594779 DOI: 10.1039/d2cp04285a] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Experimental reports about observation of spontaneous mirror symmetry breaking and chiral amplification in stereoselective Mannich and aldol reactions, run under fully achiral initial conditions, have drawn a lot of attention, fuelled partly by the role these reactions could have played in chemical evolution as a cause for still puzzling observed homochirality of biomolecules, often considered a prerequisite for the origin of life. We have now revisited this still unresolved problem, using DFT computation of all combinatorially possible transition states and numerical solution of complete set of resulting coupled kinetic rate equations to model the aldol reaction rigorously "from the first principles" and without making any a priori assumptions. Spontaneous mirror symmetry breaking in this autocatalytic, reversible, closed and homogenous system is explained by a supercritical pitchfork bifurcation, occurring in concentrations of enantiomers due to time-delayed kinetic instability of racemic composition of reaction mixture, when reactants are initially provided in non-stoichiometric quantities. Same process, taking place under similar conditions in primordial "soup" of chemicals, might conceivably explain origin of biological homochirality of sugar molecules on early earth billions of years ago. Our results suggest that seemingly innocuous chemical reactions could exhibit unexpected and counter-intuitive emergent behaviour, when initial conditions are appropriately chosen. Chiral amplification in self-catalyzed aldol reaction occurs during approach of thermodynamic equilibrium in accord with principle of microscopic reversibility and second law of thermodynamics.
Collapse
Affiliation(s)
- Michael Mauksch
- Department of Chemistry and Pharmacy, Institute of Theoretical Chemistry, Computer Chemistry Center, Nägelsbachstrasse 25a, 91052 Erlangen, Germany.
| |
Collapse
|
5
|
Constantin T, Górski B, Tilby MJ, Chelli S, Juliá F, Llaveria J, Gillen KJ, Zipse H, Lakhdar S, Leonori D. Halogen-atom and group transfer reactivity enabled by hydrogen tunneling. Science 2022; 377:1323-1328. [DOI: 10.1126/science.abq8663] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
Abstract
The generation of carbon radicals by halogen-atom and group transfer reactions is generally achieved using tin and silicon reagents that maximize the interplay of enthalpic (thermodynamic) and polar (kinetic) effects. In this work, we demonstrate a distinct reactivity mode enabled by quantum mechanical tunneling that uses the cyclohexadiene derivative γ-terpinene as the abstractor under mild photochemical conditions. This protocol activates alkyl and aryl halides as well as several alcohol and thiol derivatives. Experimental and computational studies unveiled a noncanonical pathway whereby a cyclohexadienyl radical undergoes concerted aromatization and halogen-atom or group abstraction through the reactivity of an effective H atom. This activation mechanism is seemingly thermodynamically and kinetically unfavorable but is rendered feasible through quantum tunneling.
Collapse
Affiliation(s)
| | - Bartosz Górski
- Department of Chemistry, University of Manchester, Manchester M13 9PL, UK
| | - Michael J. Tilby
- Department of Chemistry, University of Manchester, Manchester M13 9PL, UK
| | - Saloua Chelli
- CNRS/Université Toulouse III—Paul Sabatier, Laboratoire Hétérochimie Fondamentale et Appliquée, LHFA UMR 5069, 31062 Toulouse Cedex 09, France
| | - Fabio Juliá
- Department of Chemistry, University of Manchester, Manchester M13 9PL, UK
| | - Josep Llaveria
- Global Discovery Chemistry, Therapeutics Discovery, Janssen Research & Development, Janssen-Cilag S.A., 45007 Toledo, Spain
| | - Kevin J. Gillen
- LifeArc, Accelerator Building, Open Innovation Campus, Stevenage SG1 2FX, UK
| | - Hendrik Zipse
- Department Chemie, LMU München, D-81377 München, Germany
| | - Sami Lakhdar
- CNRS/Université Toulouse III—Paul Sabatier, Laboratoire Hétérochimie Fondamentale et Appliquée, LHFA UMR 5069, 31062 Toulouse Cedex 09, France
| | - Daniele Leonori
- Institute of Organic Chemistry, RWTH Aachen University, 52056 Aachen, Germany
| |
Collapse
|
6
|
|
7
|
Jana K, Wakchaure PD, Hussain N, Mukherjee D, Ganguly B. The mechanism of conversion of substituted glycals to chiral acenes via Diels-Alder reaction: a computational study. Org Biomol Chem 2021; 19:6353-6367. [PMID: 34231642 DOI: 10.1039/d1ob00408e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Synthesis of linearly fused aromatic systems using a glycal-based diene with an aryne is a long-standing topic of interest in glycal chemistry. We have examined the mechanistic pathways for the transformation of substituted glycals to chiral fused aromatic cores via Diels-Alder (DA) reaction using the SMDACN-M06-2X/6-31G(d) level of theory. The DA reactions of E (1a) and Z (1a') forms of C-2 alkenyl glycal and an aryl glycal (1b) as a diene were examined with a benzyne intermediate generated as a dienophile. The computational results reveal that 1a and 1b can only be transformed into the fused aromatic cores by the base-catalyzed reaction because a [1,5] sigmatropic hydrogen shift is not feasible. The activation free energy barrier for the base-catalyzed proton abstraction process is 4.2 kcal mol-1 and there is almost no barrier for stereoisomeric 1a DA-complexes. The activation free energy barrier values for stereoisomeric 1b DA-complexes for the base-catalyzed proton abstraction process are 10.8 and 12.4 kcal mol-1. The appropriate orientation of glycal-ring-oxygen and hydrogen at the 5th position of Z (1a') forms of C-2 alkenyl glycal facilitates the [1,5] sigmatropic hydrogen shift; however, the base-catalyzed reaction is energetically more favored than the former case. The rate-determining step for 1a and 1a' is the ring-opening step (18.2 and 19.5 kcal mol-1 for the S-stereoisomer), whereas the DA adduct formation step is the rate-determining step for 1b (16.1 kcal mol-1 for the S-stereoisomer). The structural analysis reveals the formation of the preferred S-stereoisomer over the R-stereoisomer with the respective dienes.
Collapse
Affiliation(s)
- Kalyanashis Jana
- CSIR-Central Salt and Marine Chemicals Research Institute, Bhavnagar 364002, India. and Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh-201 002, India
| | - Padmaja D Wakchaure
- CSIR-Central Salt and Marine Chemicals Research Institute, Bhavnagar 364002, India. and Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh-201 002, India
| | - Nazar Hussain
- Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh-201 002, India and Natural Product Chemistry Division, Indian Institute of Integrative Medicine (IIIM), Jammu 180001, India
| | - Debaraj Mukherjee
- Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh-201 002, India and Natural Product Chemistry Division, Indian Institute of Integrative Medicine (IIIM), Jammu 180001, India
| | - Bishwajit Ganguly
- CSIR-Central Salt and Marine Chemicals Research Institute, Bhavnagar 364002, India. and Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh-201 002, India
| |
Collapse
|