1
|
Shirai S, Iwakiri H, Kanno K, Horiba T, Omiya K, Hirai H, Koh S. Computational Analysis of Chemical Reactions Using a Variational Quantum Eigensolver Algorithm without Specifying Spin Multiplicity. ACS OMEGA 2023; 8:19917-19925. [PMID: 37305284 PMCID: PMC10249088 DOI: 10.1021/acsomega.3c01875] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Accepted: 05/10/2023] [Indexed: 06/13/2023]
Abstract
The analysis of a chemical reaction along the ground-state potential energy surface in conjunction with an unknown spin state is challenging because electronic states must be separately computed several times using different spin multiplicities to find the lowest energy state. However, in principle, the ground state could be obtained with just a single calculation using a quantum computer without specifying the spin multiplicity in advance. In the present work, ground-state potential energy curves for PtCO were calculated as a proof-of-concept using a variational quantum eigensolver (VQE) algorithm. This system exhibits a singlet-triplet crossover as a consequence of the interaction between Pt and CO. VQE calculations using a statevector simulator were found to converge to a singlet state in the bonding region, while a triplet state was obtained at the dissociation limit. Calculations performed using an actual quantum device provided potential energies within ±2 kcal/mol of the simulated energies after error mitigation techniques were adopted. The spin multiplicities in the bonding and dissociation regions could be clearly distinguished even in the case of a small number of shots. The results of this study suggest that quantum computing can be a powerful tool for the analysis of the chemical reactions of systems for which the spin multiplicity of the ground state and variations in this parameter are not known in advance.
Collapse
Affiliation(s)
- Soichi Shirai
- Toyota
Central Research and Development Laboratories, Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | - Hokuto Iwakiri
- QunaSys
Inc., Aqua Hakusan Building
9F, 1-13-7 Hakusan, Bunkyo, Tokyo 113-0001, Japan
| | - Keita Kanno
- QunaSys
Inc., Aqua Hakusan Building
9F, 1-13-7 Hakusan, Bunkyo, Tokyo 113-0001, Japan
| | - Takahiro Horiba
- Toyota
Central Research and Development Laboratories, Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | - Keita Omiya
- QunaSys
Inc., Aqua Hakusan Building
9F, 1-13-7 Hakusan, Bunkyo, Tokyo 113-0001, Japan
| | - Hirotoshi Hirai
- Toyota
Central Research and Development Laboratories, Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | - Sho Koh
- QunaSys
Inc., Aqua Hakusan Building
9F, 1-13-7 Hakusan, Bunkyo, Tokyo 113-0001, Japan
| |
Collapse
|
2
|
Heller ER, Richardson JO. Heavy‐Atom Quantum Tunnelling in Spin Crossovers of Nitrenes**. Angew Chem Int Ed Engl 2022; 61:e202206314. [PMID: 35698730 PMCID: PMC9540336 DOI: 10.1002/anie.202206314] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Indexed: 01/01/2023]
Abstract
We simulate two recent matrix‐isolation experiments at cryogenic temperatures, in which a nitrene undergoes spin crossover from its triplet state to a singlet state via quantum tunnelling. We detail the failure of the commonly applied weak‐coupling method (based on a linear approximation of the potentials) in describing these deep‐tunnelling reactions. The more rigorous approach of semiclassical golden‐rule instanton theory in conjunction with double‐hybrid density‐functional theory and multireference perturbation theory does, however, provide rate constants and kinetic isotope effects in good agreement with experiment. In addition, these calculations locate the optimal tunnelling pathways, which provide a molecular picture of the reaction mechanism. The reactions involve substantial heavy‐atom quantum tunnelling of carbon, nitrogen and oxygen atoms, which unexpectedly even continues to play a role at room temperature.
Collapse
Affiliation(s)
- Eric R. Heller
- Laboratory of Physical Chemistry ETH Zürich 8093 Zürich Switzerland
| | | |
Collapse
|
3
|
Heller ER, Richardson JO. Heavy‐Atom Quantum Tunnelling in Spin Crossovers of Nitrenes. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202206314] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Eric R Heller
- Eidgenossische Technische Hochschule Zurich Lab. Physical Chemistry SWITZERLAND
| | - Jeremy O Richardson
- Eidgenössische Technische Hochschule Zürich Lab. Physical Chemistry Vladimir-Prelog-Weg 2 8093 Zurich SWITZERLAND
| |
Collapse
|
4
|
Ovalle S, Malardier-Jugroot C. Choice of Functional for Iron Porphirin Density Functional Theory Studies: Geometry, Spin-State, and Binding Energy Analysis. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113726] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
5
|
Vitillo JG, Cramer CJ, Gagliardi L. Multireference Methods are Realistic and Useful Tools for Modeling Catalysis. Isr J Chem 2022. [DOI: 10.1002/ijch.202100136] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Jenny G. Vitillo
- Department of Science and High Technology and INSTM Università degli Studi dell'Insubria Via Valleggio 9 I-22100 Como Italy
| | - Christopher J. Cramer
- Underwriters Laboratories Inc. 333 Pfingsten Road Northbrook Illinois 60602 United States
| | - Laura Gagliardi
- Department of Chemistry Pritzker School of Molecular Engineering James Franck Institute University of Chicago Chicago Illinois 60637 United States
| |
Collapse
|
6
|
Sigmund LM, Ehlert C, Enders M, Graf J, Gryn'ova G, Greb L. Disauerstoffaktivierung und Pyrrol‐α‐Spaltung mit Calix[4]pyrrolatoaluminaten: Enzymmodell durch strukturellen Zwang. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202104916] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Lukas Maximilian Sigmund
- Anorganisch-Chemisches Institut Ruprecht-Karls-Universität Heidelberg Im Neuenheimer Feld 270 69120 Heidelberg Deutschland
| | - Christopher Ehlert
- Heidelberger Institut für Theoretische Studien (HITS gGmbH) Schloss-Wolfsbrunnenweg 35 69118 Heidelberg Deutschland
- Interdisziplinäres Zentrum für wissenschaftliches Rechnen (IWR) Ruprecht-Karls-Universität Heidelberg Im Neuenheimer Feld 205 69120 Heidelberg Deutschland
| | - Markus Enders
- Anorganisch-Chemisches Institut Ruprecht-Karls-Universität Heidelberg Im Neuenheimer Feld 270 69120 Heidelberg Deutschland
| | - Jürgen Graf
- Organisch-Chemisches Institut Ruprecht-Karls-Universität Heidelberg Im Neuenheimer Feld 270 69120 Heidelberg Deutschland
| | - Ganna Gryn'ova
- Heidelberger Institut für Theoretische Studien (HITS gGmbH) Schloss-Wolfsbrunnenweg 35 69118 Heidelberg Deutschland
- Interdisziplinäres Zentrum für wissenschaftliches Rechnen (IWR) Ruprecht-Karls-Universität Heidelberg Im Neuenheimer Feld 205 69120 Heidelberg Deutschland
| | - Lutz Greb
- Anorganisch-Chemisches Institut Ruprecht-Karls-Universität Heidelberg Im Neuenheimer Feld 270 69120 Heidelberg Deutschland
| |
Collapse
|
7
|
Sigmund LM, Ehlert C, Enders M, Graf J, Gryn'ova G, Greb L. Dioxygen Activation and Pyrrole α-Cleavage with Calix[4]pyrrolato Aluminates: Enzyme Model by Structural Constraint. Angew Chem Int Ed Engl 2021; 60:15632-15640. [PMID: 33955154 PMCID: PMC8362023 DOI: 10.1002/anie.202104916] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2021] [Indexed: 01/30/2023]
Abstract
The present work describes the reaction of triplet dioxygen with the porphyrinogenic calix[4]pyrrolato aluminates to alkylperoxido aluminates in high selectivity. Multiconfigurational quantum chemical computations disclose the mechanism for this spin‐forbidden process. Despite a negligible spin–orbit coupling constant, the intersystem crossing (ISC) is facilitated by singlet and triplet state degeneracy and spin–vibronic coupling. The formed peroxides are stable toward external substrates but undergo an unprecedented oxidative pyrrole α‐cleavage by ligand aromatization/dearomatization‐initiated O−O σ‐bond scission. A detailed comparison of the calix[4]pyrrolato aluminates with dioxygen‐related enzymology provides insights into the ISC of metal‐ or cofactor‐free enzymes. It substantiates the importance of structural constraint and element–ligand cooperativity for the functions of aerobic life.
Collapse
Affiliation(s)
- Lukas Maximilian Sigmund
- Anorganisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Christopher Ehlert
- Heidelberg Institute for Theoretical Studies (HITS gGmbH), Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing (IWR), Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 205, 69120, Heidelberg, Germany
| | - Markus Enders
- Anorganisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Jürgen Graf
- Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| | - Ganna Gryn'ova
- Heidelberg Institute for Theoretical Studies (HITS gGmbH), Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing (IWR), Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 205, 69120, Heidelberg, Germany
| | - Lutz Greb
- Anorganisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Im Neuenheimer Feld 270, 69120, Heidelberg, Germany
| |
Collapse
|
8
|
Zhao L, Watanabe KJ, Nakatani N, Nakayama A, Xu X, Hasegawa JY. Extending nudged elastic band method to reaction pathways involving multiple spin states. J Chem Phys 2020; 153:134114. [PMID: 33032404 DOI: 10.1063/5.0021923] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
There are diverse reactions including spin-state crossing, especially the reactions catalyzed by transition metal compounds. To figure out the mechanisms of such reactions, the discussion of minimum energy intersystem crossing (MEISC) points cannot be avoided. These points may be the bottleneck of the reaction or inversely accelerate the reactions by providing a better pathway. It is of great importance to reveal their role in the reactions by computationally locating the position of the MEISC points together with the reaction pathway. However, providing a proper initial guess for the structure of the MEISC point is not as easy as that of the transition state. In this work, we extended the nudged elastic band (NEB) method for multiple spin systems, which is named the multiple spin-state NEB method, and it is successfully applied to find the MEISC points while optimizing the reaction pathway. For more precisely locating the MEISC point, a revised approach was adopted. Meanwhile, our examples also suggest that special attention should be paid to the criterion to define an image optimized as the MEISC point.
Collapse
Affiliation(s)
- Liming Zhao
- Institute for Catalysis, Hokkaido University, N21 W10 Kita-ku, Sapporo 001-0021, Hokkaido, Japan
| | - K-Jiro Watanabe
- Institute for Catalysis, Hokkaido University, N21 W10 Kita-ku, Sapporo 001-0021, Hokkaido, Japan
| | - Naoki Nakatani
- Graduate School of Science and Engineering, Tokyo Metropolitan University, Minami-Osawa 1-1, Tokyo 192-0397, Japan
| | - Akira Nakayama
- Department of Chemical System Engineering, The University of Tokyo, Tokyo 113-8656, Japan
| | - Xin Xu
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Jun-Ya Hasegawa
- Institute for Catalysis, Hokkaido University, N21 W10 Kita-ku, Sapporo 001-0021, Hokkaido, Japan
| |
Collapse
|
9
|
Saito K, Watabe Y, Miyazaki T, Takayanagi T, Hasegawa JY. Spin-inversion mechanisms in O 2 binding to a model heme compound: A perspective from nonadiabatic wave packet calculations. J Comput Chem 2020; 41:2527-2537. [PMID: 32841410 DOI: 10.1002/jcc.26409] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Revised: 07/18/2020] [Accepted: 08/08/2020] [Indexed: 12/13/2022]
Abstract
Spin-inversion dynamics in O2 binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe-O distance, Fe-O-O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin-orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin-orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O2 binding dynamics. In addition, the one-dimensional Landau-Zener model underestimated the nonadiabatic transition probability.
Collapse
Affiliation(s)
- Kohei Saito
- Department of Chemistry, Saitama University, Saitama City, Saitama, Japan
| | - Yuya Watabe
- Department of Chemistry, Saitama University, Saitama City, Saitama, Japan
| | - Takaaki Miyazaki
- Department of Chemistry, Saitama University, Saitama City, Saitama, Japan
| | | | - Jun-Ya Hasegawa
- Institute for Catalysis, Hokkaido University, Sapporo, Hokkaido, Japan
| |
Collapse
|
10
|
Saito K, Watabe Y, Fujihara T, Takayanagi T, Hasegawa JY. Spin-inversion mechanisms in O 2 binding to a model heme complex revisited by density function theory calculations. J Comput Chem 2020; 41:1130-1138. [PMID: 32020659 DOI: 10.1002/jcc.26159] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Revised: 01/12/2020] [Accepted: 01/16/2020] [Indexed: 01/10/2023]
Abstract
Spin-inversion mechanisms in O2 binding to a model heme complex, consisting of Fe(II)-porphyrin and imidazole, were investigated using density-functional theory calculations. First, we applied the recently proposed mixed-spin Hamiltonian method to locate spin-inversion structures between different total spin multiplicities. Nine spin-inversion structures were successfully optimized for the singlet-triplet, singlet-quintet, triplet-quintet, and quintet-septet spin-inversion processes. We found that the singlet-triplet spin-inversion points are located around the potential energy surface region at short Fe-O distances, whereas the singlet-quintet and quintet-septet spin-inversion points are located at longer Fe-O distances. This suggests that both narrow and broad crossing models play roles in O2 binding to the Fe-porphyrin complex. To further understand spin-inversion mechanisms, we performed on-the-fly Born-Oppenheimer molecular dynamics calculations. The reaction coordinates, which are correlated to the spin-inversion dynamics between different spin multiplicities, are also discussed.
Collapse
Affiliation(s)
- Kohei Saito
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan
| | - Yuya Watabe
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan
| | - Takashi Fujihara
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama, 338-8570, Japan
| | - Jun-Ya Hasegawa
- Instituteof Catalysis, Hokkaido University, Kita 21, Nishi 10, Kita-ku, Sapporo, Hokkaido, 001-0021, Japan
| |
Collapse
|
11
|
A detailed analysis of the spin-crossover reaction of H 2S binding to heme and the six-coordinated FeP(Im)-HS - porphyrin complex. J Inorg Biochem 2020; 206:111049. [PMID: 32171934 DOI: 10.1016/j.jinorgbio.2020.111049] [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: 10/28/2019] [Revised: 02/21/2020] [Accepted: 02/22/2020] [Indexed: 11/23/2022]
Abstract
The potential energy surfaces of the H2S binding to iron-porphyrin (FeP) with the imidazole (Im) ligand via intersystem crossings are investigated by using density functional theory. The minimum energy intersystem crossing point (MEISCP) between the quintet and triplet states (MEISCPTQ) for the Fe(II)P(Im)-H2S complex is located at a Fe-S distance of 3.39 Å with only 1.1 kcal/mol above the quintet state minimum. The second spin-crossover point, where a change from the triplet to the singlet state occurs, comes at a much shorter Fe-S distance of 2.79 Å, and the MEISCPST is located at 3.7 kcal/mol above the triplet state minimum. The nature of the chemical bonding along the Fe-S reaction coordinate from the ground state singlet to the quintet state along the path to the separated species is analyzed. An inspection of the vibrational modes reveals that the largest contribution to the triplet-quintet transition around the quintet and triplet state minimum comes from the symmetric shrinking of the pyrrole units of the porphyrin ring, indicating that the related reaction coordinate plays a main role in the intersystem crossing. The fully optimized structures of the Fe(II)P(Im)-HS- complex corresponding to three different spin multiplicities (M = 1, 3, 5) are characterized by a bent Fe-H-S conformation. The binding of the hydrosulfide anion to Fe(II)P(Im) in the quintet state induces a 0.2 Å displacement of the Fe atom out of the nitrogen porphyrin (Npyr) plane. The fully optimized structure of the ground state of Fe(II)P(Im)-HS- agrees well with experimental data for the corresponding heme models.
Collapse
|
12
|
Lykhin AO, Varganov SA. Intersystem crossing in tunneling regime: T 1 → S 0 relaxation in thiophosgene. Phys Chem Chem Phys 2020; 22:5500-5508. [PMID: 32101195 DOI: 10.1039/c9cp06956a] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The T1 excited state relaxation in thiophosgene has attracted much attention as a relatively simple model for the intersystem crossing (ISC) transitions in polyatomic molecules. The very short (20-40 ps) T1 lifetime predicted in several theoretical studies strongly disagrees with the experimental values (∼20 ns) indicating that the kinetics of T1 → S0 ISC is not well understood. We use the nonadiabatic transition state theory (NA-TST) with the Zhu-Nakamura transition probability and the multireference perturbation theory (CASPT2) to show that the T1 → S0 ISC occurs in the quantum tunneling regime. We also introduce a new zero-point vibrational energy correction scheme that improves the accuracy of the predicted ISC rate constants at low internal energies. The predicted lifetimes of the T1 vibrational states are between one and two orders of magnitude larger than the experimental values. This overestimation is attributed to the multidimensional nature of quantum tunneling that facilitates ISC transitions along the non-minimum energy path and is not accounted for in the one-dimensional NA-TST.
Collapse
Affiliation(s)
- Aleksandr O Lykhin
- Department of Chemistry, University of Nevada, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA.
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA.
| |
Collapse
|
13
|
Ma Z, Ukaji K, Nakatani N, Fujii H, Hada M. Substitution effects on olefin epoxidation catalyzed by Oxoiron(IV) porphyrin π-cation radical complexes: A dft study. J Comput Chem 2019; 40:1780-1788. [PMID: 30938845 DOI: 10.1002/jcc.25831] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2019] [Revised: 03/12/2019] [Accepted: 03/15/2019] [Indexed: 01/26/2023]
Abstract
The effects of peripheral fluorine atoms on epoxidation reactions of ethylene by oxoiron(IV) porphyrin cation radical complex in the quartet and sextet spin multiplicities are systematically investigated using the DFT method. The overall reaction routes are determined using a model system of ethylene and Fe(IV)OCl-porphyrin with substituted fluorine atoms. By obtaining the energy diagrams and electron- and spin-density difference contour maps of the transition states and intermediate compounds, we confirm that the electron-withdrawing by peripheral fluorine atoms enhances the reactivity as the number of fluorine atoms increases, as is observed experimentally. The intersystem crossing between the quartet and sextet spin multiplicities is discussed by means of the intrinsic reaction coordinate method. We conclude that the rate-determining step is located at the first transition state (TS1) for the activation of CC and FeO bonds, and the ground electronic state changes from quartet to sextet around the TS1. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Zhifeng Ma
- Department of Chemistry, Graduate School of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo 192-0397, Japan
| | - Kasumi Ukaji
- Department of Chemistry, Graduate School of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo 192-0397, Japan
| | - Naoki Nakatani
- Department of Chemistry, Graduate School of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo 192-0397, Japan
| | - Hiroshi Fujii
- Department of Chemistry, Graduate School of Humanities and Science, Nara Women's University, Kitauoyanishi, Nara 630-8506, Japan
| | - Masahiko Hada
- Department of Chemistry, Graduate School of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo 192-0397, Japan
| |
Collapse
|
14
|
Aarabi M, Omidyan R, Soorkia S, Grégoire G, Broquier M, Crestoni ME, de la Lande A, Soep B, Shafizadeh N. The dramatic effect of N-methylimidazole on trans axial ligand binding to ferric heme: experiment and theory. Phys Chem Chem Phys 2019; 21:1750-1760. [PMID: 30623949 DOI: 10.1039/c8cp06210b] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The binding energy of CO, O2 and NO to isolated ferric heme, [FeIIIP]+, was studied in the presence and absence of a σ donor (N-methylimidazole and histidine) as the trans axial ligand. This study combines the experimental determination of binding enthalpies by equilibrium measurements in a low temperature ion trap using the van't Hoff equation and high level DFT calculations. It was found that the presence of N-methylimidazole as the axial ligand on the [FeIIIP]+ porphyrin dramatically weakens the [FeIIIP-ligand]+ bond with an up to sevenfold decrease in binding energy owing to the σ donation by N-methylimidazole to the FeIII(3d) orbitals. This trans σ donor effect is characteristic of ligation to iron in hemes in both ferrous and ferric redox forms; however, to date, this has not been observed for ferric heme.
Collapse
Affiliation(s)
- Mohammad Aarabi
- Department of Chemistry, University of Isfahan, 81746-73441 Isfahan, Iran.
| | | | | | | | | | | | | | | | | |
Collapse
|
15
|
Harada I, Nakayama A, Hasegawa JY. Constraint structure optimization to a specific minimum using ionization energy. J Comput Chem 2019; 40:507-514. [PMID: 30414201 DOI: 10.1002/jcc.25738] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2018] [Revised: 10/02/2018] [Accepted: 10/02/2018] [Indexed: 11/10/2022]
Abstract
A structure optimization method with ionization energy constraint is developed to explore structures with specific ionization energy. The Levine-Coe-Martínez penalty function (J Phys Chem B 2008, 112, 405) was adopted, and the penalty function includes a predefined core-ionization energy and inner-shell ionization energy. For an SN 2 reaction, isomerization of a platinum complex, a proton transfer reaction, and carbon monoxide adsorption on a palladium cluster, the present method was tested, and the targeted energy minima were obtained as designated by the input ionization energy. The shape of the objective function, the parameters in the penalty function, and structural changes during the optimization process were discussed. An automated parameter setting and possible problems are discussed for future direction. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Iori Harada
- Institute for Catalysis, Hokkaido University, Kita-ku, Sapporo, 001-0021, Japan
| | - Akira Nakayama
- Institute for Catalysis, Hokkaido University, Kita-ku, Sapporo, 001-0021, Japan
| | - Jun-Ya Hasegawa
- Institute for Catalysis, Hokkaido University, Kita-ku, Sapporo, 001-0021, Japan
| |
Collapse
|
16
|
Lischka H, Nachtigallová D, Aquino AJA, Szalay PG, Plasser F, Machado FBC, Barbatti M. Multireference Approaches for Excited States of Molecules. Chem Rev 2018; 118:7293-7361. [DOI: 10.1021/acs.chemrev.8b00244] [Citation(s) in RCA: 197] [Impact Index Per Article: 32.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Hans Lischka
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Dana Nachtigallová
- Institute of Organic Chemistry and Biochemistry v.v.i., The Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Palacký University, 78371 Olomouc, Czech Republic
| | - Adélia J. A. Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute for Soil Research, University of Natural Resources and Life Sciences Vienna, Peter-Jordan-Strasse 82, A-1190 Vienna, Austria
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Laboratory of Theoretical Chemistry, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary
| | - Felix Plasser
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Leicestershire LE11 3TU, United Kingdom
| | - Francisco B. C. Machado
- Departamento de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | | |
Collapse
|
17
|
Ostojić B, Schwerdtfeger P, Đorđević D. Modeling the hydrogen sulfide binding to heme. J Inorg Biochem 2018; 184:108-114. [DOI: 10.1016/j.jinorgbio.2018.04.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2017] [Revised: 04/14/2018] [Accepted: 04/17/2018] [Indexed: 10/17/2022]
|
18
|
Ben Amor N, Soupart A, Heitz MC. Methodological CASPT2 study of the valence excited states of an iron-porphyrin complex. J Mol Model 2017; 23:53. [DOI: 10.1007/s00894-017-3226-y] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2016] [Accepted: 01/12/2017] [Indexed: 10/20/2022]
|
19
|
Watanabe KJ, Nakatani N, Nakayama A, Higashi M, Hasegawa JY. Spin-Blocking Effect in CO and H2 Binding Reactions to Molybdenocene and Tungstenocene: A Theoretical Study on the Reaction Mechanism via the Minimum Energy Intersystem Crossing Point. Inorg Chem 2016; 55:8082-90. [DOI: 10.1021/acs.inorgchem.6b01187] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- K-jiro Watanabe
- Institute for Catalysis, Hokkaido University, N21 W10 Kita-ku, Sapporo, Hokkaido 001-0021, Japan
- Department of Chemistry, Biology and Marine Science,
Faculty of Science, University of the Ryukyus, 1 Senbaru, Nishihara, Okinawa 903-0213, Japan
| | - Naoki Nakatani
- Institute for Catalysis, Hokkaido University, N21 W10 Kita-ku, Sapporo, Hokkaido 001-0021, Japan
| | - Akira Nakayama
- Institute for Catalysis, Hokkaido University, N21 W10 Kita-ku, Sapporo, Hokkaido 001-0021, Japan
| | - Masahiro Higashi
- Department of Chemistry, Biology and Marine Science,
Faculty of Science, University of the Ryukyus, 1 Senbaru, Nishihara, Okinawa 903-0213, Japan
| | - Jun-ya Hasegawa
- Institute for Catalysis, Hokkaido University, N21 W10 Kita-ku, Sapporo, Hokkaido 001-0021, Japan
| |
Collapse
|