1
|
Alessandrini S, Melosso M, Bizzocchi L, Barone V, Puzzarini C. The Semiexperimental Approach at Work: Equilibrium Structure of Radical Species. J Phys Chem A 2024; 128:5833-5855. [PMID: 38991181 DOI: 10.1021/acs.jpca.4c01758] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
The so-called semiexperimental (SE) approach is a powerful technique for obtaining highly accurate equilibrium structures for isolated systems. This Featured Article describes its extension to open-shell species, thus providing the first systematic investigation on radical equilibrium geometries to be used for benchmarking purposes. The small yet significant database obtained demonstrates that there is no reduction in accuracy when moving from closed-shell species to radicals. We also provide an extension of the applicability of the SE approach to medium-/large-sized radicals by exploiting the so-called "Lego-brick" approach, which is based on the assumption that a molecular system can be seen as formed by smaller fragments for which the SE equilibrium structure is available. In this Featured Article we show that this model can be successfully applied also to open-shell species.
Collapse
Affiliation(s)
- Silvia Alessandrini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Mattia Melosso
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Luca Bizzocchi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | | | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| |
Collapse
|
2
|
Jaworski A, Hedin N. Electron correlation and vibrational effects in predictions of paramagnetic NMR shifts. Phys Chem Chem Phys 2022; 24:15230-15244. [PMID: 35703010 DOI: 10.1039/d2cp01206e] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Electronic structure calculations are fundamentally important for the interpretation of nuclear magnetic resonance (NMR) spectra from paramagnetic systems that include organometallic and inorganic compounds, catalysts, or metal-binding sites in proteins. Prediction of induced paramagnetic NMR shifts requires knowledge of electron paramagnetic resonance (EPR) parameters: the electronic g tensor, zero-field splitting D tensor, and hyperfine A tensor. The isotropic part of A, called the hyperfine coupling constant (HFCC), is one of the most troublesome properties for quantum chemistry calculations. Yet, even relatively small errors in calculations of HFCC tend to propagate into large errors in the predicted NMR shifts. The poor quality of A tensors that are currently calculated using density functional theory (DFT) constitutes a bottleneck in improving the reliability of interpretation of the NMR spectra from paramagnetic systems. In this work, electron correlation effects in calculations of HFCCs with a hierarchy of ab initio methods were assessed, and the applicability of different levels of DFT approximations and the coupled cluster singles and doubles (CCSD) method was tested. These assessments were performed for the set of selected test systems comprising an organic radical, and complexes with transition metal and rare-earth ions, for which experimental data are available. Severe deficiencies of DFT were revealed but the CCSD method was able to deliver good agreement with experimental data for all systems considered, however, at substantial computational costs. We proposed a more computationally tractable alternative, where the A was computed with the coupled cluster theory exploiting locality of electron correlation. This alternative is based on the domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) method. In this way the robustness and reliability of the coupled cluster theory were incorporated into the modern formalism for the prediction of induced paramagnetic NMR shifts, and became applicable to systems of chemical interest. This approach was verified for the bis(cyclopentadienyl)vanadium(II) complex (Cp2V; vanadocene), and the metal-binding site of the Zn2+ → Co2+ substituted superoxide dismutase (SOD) metalloprotein. Excellent agreement with experimental NMR shifts was achieved, which represented a substantial improvement over previous theoretical attempts. The effects of vibrational corrections to orbital shielding and hyperfine tensor were evaluated and discussed within the second-order vibrational perturbation theory (VPT2) framework.
Collapse
Affiliation(s)
- Aleksander Jaworski
- Department of Materials and Environmental Chemistry, Stockholm University, SE-106 91 Stockholm, Sweden.
| | - Niklas Hedin
- Department of Materials and Environmental Chemistry, Stockholm University, SE-106 91 Stockholm, Sweden.
| |
Collapse
|
3
|
Mukhopadhyay DP, Gerlach M, Hartweg S, Fischer I, Loison JC. Photoelectron spectroscopy of low valent organophosphorus compounds, P-CH 3, H-PCH 2 and PCH 2. Phys Chem Chem Phys 2022; 24:10993-10999. [PMID: 35467677 DOI: 10.1039/d2cp01082h] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report the mass-selected slow photoelectron spectra of three reactive organophosphorus species, PCH2, and the two isomers, methylenephosphine or phosphaethylene, HPCH2 and methylphosphinidine, P-CH3. All spectra were recorded by double imaging photoelectron-photoion coincidence spectroscopy (i2PEPICO) using synchrotron radiation and all species were generated in a flow reactor by the reaction of trimethyl phosphine with fluorine atoms. Adiabatic ionisation energies of 8.80 ± 0.02 eV (PCH2), 10.07 ± 0.03 eV (H-PCH2) and 8.91 ± 0.04 eV (P-CH3) were determined and the vibronic structure was simulated by calculating Franck-Condon factors from optimised structures based on quantum chemical methods. Observation of biradicalic P-CH3 isomer with its triplet ground state is surprising because it is less stable than H-PCH2.
Collapse
Affiliation(s)
- Deb Pratim Mukhopadhyay
- Institute of Physical and Theoretical Chemistry, University of Würzburg, Am Hubland, D-97074 Würzburg, Germany.
| | - Marius Gerlach
- Institute of Physical and Theoretical Chemistry, University of Würzburg, Am Hubland, D-97074 Würzburg, Germany.
| | - Sebastian Hartweg
- Synchrotron SOLEIL, L'Orme des Merisiers, St Aubin, B. P. 48, F-91192 Gif sur Yvette, France
| | - Ingo Fischer
- Institute of Physical and Theoretical Chemistry, University of Würzburg, Am Hubland, D-97074 Würzburg, Germany.
| | | |
Collapse
|
4
|
Bizzocchi L, Alessandrini S, Melosso M, Rivilla VM, Puzzarini C. Ab Initio Study of Fine and Hyperfine Interactions in Triplet POH. Molecules 2022; 27:302. [PMID: 35011533 PMCID: PMC8746536 DOI: 10.3390/molecules27010302] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Revised: 12/28/2021] [Accepted: 12/30/2021] [Indexed: 11/16/2022] Open
Abstract
Phosphorous-containing molecules have a great relevance in prebiotic chemistry in view of the fact that phosphorous is a fundamental constituent of biomolecules, such as RNA, DNA, and ATP. Its biogenic importance has led astrochemists to investigate the possibility that P-bearing species could have formed in the interstellar medium (ISM) and subsequently been delivered to early Earth by rocky bodies. However, only two P-bearing molecules have been detected so far in the ISM, with the chemistry of interstellar phosphorous remaining poorly understood. Here, in order to shed further light on P-carriers in space, we report a theoretical spectroscopic characterisation of the rotational spectrum of POH in its 3A″ ground electronic state. State-of-the-art coupled-cluster schemes have been employed to derive rotational constants, centrifugal distortion terms, and most of the fine and hyperfine interaction parameters, while the electron spin-spin dipolar coupling has been investigated using the multi-configuration self-consistent-field method. The computed spectroscopic parameters have been used to simulate the appearance of triplet POH rotational and ro-vibrational spectra in different conditions, from cold to warm environments, either in gas-phase experiments or in molecular clouds. Finally, we point out that the predicted hyperfine structures represent a key pattern for the recognition of POH in laboratory and interstellar spectra.
Collapse
Affiliation(s)
- Luca Bizzocchi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy; (L.B.); (S.A.)
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy;
| | - Silvia Alessandrini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy; (L.B.); (S.A.)
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy;
| | - Mattia Melosso
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy;
- Scuola Superiore Meridionale, Università di Napoli Federico II, Largo San Marcellino 10, 80138 Naples, Italy
| | - Víctor M. Rivilla
- Centro de Astrobiología (CSIC-INTA), Ctra. de Ajalvir Km. 4, Torrejón de Ardoz, 28850 Madrid, Spain;
- INAF—Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Florence, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy;
| |
Collapse
|
5
|
Chen G, Qin Z, Li J, Liu L. A global CHIPR potential energy surface of PH 2(X 2B 1) via extrapolation to the complete basis set limit and the dynamics of P( 2D) + H 2(X 1Σ+g) → PH(X 3Σ −) + H( 2S). Phys Chem Chem Phys 2022; 24:19371-19381. [DOI: 10.1039/d2cp02690b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
A global PEC of the ground state PH2 was constructed using the CHIPR method based on accurate MRCI(Q)/CBS(T, Q) energy points. The ICS and k(T) of P(2D) + H2(X1Σ+g) → PH(X3Σ−) + H(2S) were calculated based on the QCT method.
Collapse
Affiliation(s)
- Guangan Chen
- School of Energy and Power Engineering, Shandong University, 250061, Jinan, China
| | - Zhi Qin
- School of Energy and Power Engineering, Shandong University, 250061, Jinan, China
- Optics and Thermal Radiation Research Center, Institute of Frontier and Interdisciplinary Science, Shandong University, 266237, Qingdao, China
| | - Jing Li
- School of Physics and Physical Engineering, Qufu Normal University, 273165, Qufu, China
| | - Linhua Liu
- School of Energy and Power Engineering, Shandong University, 250061, Jinan, China
- Optics and Thermal Radiation Research Center, Institute of Frontier and Interdisciplinary Science, Shandong University, 266237, Qingdao, China
- School of Energy Science and Engineering, Harbin Institute of Technology, 150001, Harbin, China
| |
Collapse
|
6
|
Sharma B, Tran VA, Pongratz T, Galazzo L, Zhurko I, Bordignon E, Kast SM, Neese F, Marx D. A Joint Venture of Ab Initio Molecular Dynamics, Coupled Cluster Electronic Structure Methods, and Liquid-State Theory to Compute Accurate Isotropic Hyperfine Constants of Nitroxide Probes in Water. J Chem Theory Comput 2021; 17:6366-6386. [PMID: 34516119 PMCID: PMC8515807 DOI: 10.1021/acs.jctc.1c00582] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2021] [Indexed: 01/11/2023]
Abstract
The isotropic hyperfine coupling constant (HFCC, Aiso) of a pH-sensitive spin probe in a solution, HMI (2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water, is computed using an ensemble of state-of-the-art computational techniques and is gauged against X-band continuous wave electron paramagnetic resonance (EPR) measurement spectra at room temperature. Fundamentally, the investigation aims to delineate the cutting edge of current first-principles-based calculations of EPR parameters in aqueous solutions based on using rigorous statistical mechanics combined with correlated electronic structure techniques. In particular, the impact of solvation is described by exploiting fully atomistic, RISM integral equation, and implicit solvation approaches as offered by ab initio molecular dynamics (AIMD) of the periodic bulk solution (using the spin-polarized revPBE0-D3 hybrid functional), embedded cluster reference interaction site model integral equation theory (EC-RISM), and polarizable continuum embedding (using CPCM) of microsolvated complexes, respectively. HFCCs are obtained from efficient coupled cluster calculations (using open-shell DLPNO-CCSD theory) as well as from hybrid density functional theory (using revPBE0-D3). Re-solvation of "vertically desolvated" spin probe configuration snapshots by EC-RISM embedding is shown to provide significantly improved results compared to CPCM since only the former captures the inherent structural heterogeneity of the solvent close to the spin probe. The average values of the Aiso parameter obtained based on configurational statistics using explicit water within AIMD and from EC-RISM solvation are found to be satisfactorily close. Using either such explicit or RISM solvation in conjunction with DLPNO-CCSD calculations of the HFCCs provides an average Aiso parameter for HMI in aqueous solution at 300 K and 1 bar that is in good agreement with the experimentally determined one. The developed computational strategy is general in the sense that it can be readily applied to other spin probes of similar molecular complexity, to aqueous solutions beyond ambient conditions, as well as to other solvents in the longer run.
Collapse
Affiliation(s)
- Bikramjit Sharma
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Van Anh Tran
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Tim Pongratz
- Physikalische
Chemie III, Technische Universität
Dortmund, Otto-Hahn-Str. 4a, 44227 Dortmund, Germany
| | - Laura Galazzo
- Faculty
of Chemistry and Biochemistry, Ruhr University
Bochum, 44780 Bochum, Germany
| | - Irina Zhurko
- Laboratory
of Nitrogen Compounds, N.N. Vorozhtsov Novosibirsk Institute of Organic
Chemistry, NIOCH SB RAS, 9 Lavrentiev Avenue, 630090, Novosibirsk, Russia
| | - Enrica Bordignon
- Faculty
of Chemistry and Biochemistry, Ruhr University
Bochum, 44780 Bochum, Germany
| | - Stefan M. Kast
- Physikalische
Chemie III, Technische Universität
Dortmund, Otto-Hahn-Str. 4a, 44227 Dortmund, Germany
| | - Frank Neese
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Dominik Marx
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| |
Collapse
|
7
|
Gromov OI. Performance of the DLPNO-CCSD and recent DFT methods in the calculation of isotropic and dipolar contributions to 14N hyperfine coupling constants of nitroxide radicals. J Mol Model 2021; 27:194. [PMID: 34075533 DOI: 10.1007/s00894-021-04807-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Accepted: 05/26/2021] [Indexed: 10/21/2022]
Abstract
In the present study, the performance of a set of density functionals: BP86, PBE, OLYP, BEEF, PBEpow, TPSS, SCAN, PBEGXPBE, M06L, MN15L, B3LYP, PBE0, mPW1PW, B97, BHandHLYP, mPW1PW, B98, TPSS0, PBE1KCIS, SCAN0, M06, M06-2X, MN15, CAM-B3LYP, ωB97x, B2PLYP, and the B3LYP/N07D and PBE/N07D schemes in the calculation of the 14N anisotropic hyperfine coupling (HFC) constants of a set of 23 nitroxide radicals is evaluated. The results are compared with those obtained with the DLPNO-CCSD method and experimental HFC values. Harmonic contribution to the 14N HFC vibrational correction was calculated at the revPBE0/def2-TZVPP level and included in the evaluation. With the vibrational correction, the DLPNO-CCSD method yielded HFC values in good agreement with the experiment (mean absolute deviation (MAD) = 0.3 G for the dipole-dipole contribution and MAD = 0.8 G for the contact coupling contribution). The best DFT results are obtained using the M06 functional with MAD = 0.2 G for the dipole-dipole contribution and MAD = 0.7 G for the contact coupling contribution. In general, vibrational correction significantly improved most DFT functionals' performance but did not change its overall ranking.
Collapse
Affiliation(s)
- Oleg I Gromov
- Chemistry Department, Lomonosov Moscow State University, Leninskiye Gory 1-3, Moscow, 119991, Russia.
| |
Collapse
|
8
|
Chaudhuri RK, Chattopadhyay S. Description of the Methylene Amidogene Radical and Its Anion with an Economical Treatment of Correlation Effects Using Density Functional Theory Orbitals. J Phys Chem A 2021; 125:543-558. [PMID: 33417452 DOI: 10.1021/acs.jpca.0c08635] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The ground and low-lying excited state electronic structural properties (such as equilibrium geometries, harmonic frequencies, excitation energies, barrier energy, and so on) of the methylene amidogene radical (H2CN) and its anion (H2CN-) have been studied using the CASCI (complete active space configuration interaction) and SSMRPT (state-specific multireference Møller-Plesset perturbation theory) methods with density function theory (DFT) orbitals. Here, the span of the active orbitals have been obtained from Kohn-Sham DFT using B3LYP exchange-correlation functionals in the CASCI (DFT-CASCI) approximation to describe nondynamic correlation associated with electronic degeneracies. The DFT-SSMRPT protocol provides an attractive way to deal with both dynamical and nondynamical correlation effects in strongly correlated systems such as H2CN and H2CN-. The present work clearly indicates that the electronic absorption band near 35,050 cm-1 corresponds to the B̃2A1 ← X̃2B2 transition. DFT-SSMRPT findings are in close agreement with high-level theoretical estimates. It is concluded that the transition at 1725 cm-1 could be due to the CN stretching of the trans-HCNH isomer which is originally assigned to the CN stretch of H2CN in the experiment. The present results confirm most of the previous vibrational assignments. It is not possible to definitively assign a transition to the 35,600 cm-1 band with the present estimations, suggesting further experiment is urgently called for.
Collapse
Affiliation(s)
| | - Sudip Chattopadhyay
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| |
Collapse
|
9
|
Ghassemi Tabrizi S, Arbuznikov AV, Jiménez-Hoyos CA, Kaupp M. Hyperfine-Coupling Tensors from Projected Hartree-Fock Theory. J Chem Theory Comput 2020; 16:6222-6235. [PMID: 32841008 DOI: 10.1021/acs.jctc.0c00617] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
We assess the calculation of hyperfine coupling (HFC) tensors by different variants of Projected Hartree-Fock (PHF) theory. For a set of small main-group S = 1/2 radicals (BO, CO+, CN, AlO, vinyl, methyl, ethynyl), spin-symmetry as well as complex-conjugation and point-group symmetry are first broken in a reference determinant, and then variationally restored, in the frame of the modern formulation of PHF theory. Historically, PHF theory was basically restricted to the restoration of spin symmetry from an unrestricted HF determinant (conserving Sz symmetry). This afforded unsatisfactory HFCs. We obtain far better results for isotropic (and anisotropic) HFCs when the variational energy is further lowered by working with generalized determinants that completely break spin symmetry, and when additional symmetries are used. Specifically, complex-conjugation projection recovers a substantial fraction of the dynamical correlation energy in small molecules, and the detailed equations for combined complex-conjugation, spin- and point-group projection in the density-matrix/diagonalization formulation of PHF theory are here reported for the first time. The compact representation of the PHF wave function allows for a straightforward evaluation of the spin-density matrix and of HFC tensors with little effort. The promising performance of PHF theory may motivate the application of post-PHF methods to the calculation of HFC tensors.
Collapse
Affiliation(s)
- Shadan Ghassemi Tabrizi
- Institut für Chemie, Theoretische Chemie, Technische Universität Berlin, Sekr. C7, Strasse des 17. Juni 135, 10623 Berlin, Germany
| | - Alexei V Arbuznikov
- Institut für Chemie, Theoretische Chemie, Technische Universität Berlin, Sekr. C7, Strasse des 17. Juni 135, 10623 Berlin, Germany
| | - Carlos A Jiménez-Hoyos
- Department of Chemistry, Wesleyan University, Middletown, Connecticut 06459, United States
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie, Technische Universität Berlin, Sekr. C7, Strasse des 17. Juni 135, 10623 Berlin, Germany
| |
Collapse
|
10
|
Sanakis Y, Krzystek J, Maganas D, Grigoropoulos A, Ferentinos E, Kostakis MG, Petroulea V, Pissas M, Thirunavukkuarasu K, Wernsdorfer W, Neese F, Kyritsis P. Magnetic Properties and Electronic Structure of the S = 2 Complex [Mn III{(OPPh 2) 2N} 3] Showing Field-Induced Slow Magnetization Relaxation. Inorg Chem 2020; 59:13281-13294. [PMID: 32897702 DOI: 10.1021/acs.inorgchem.0c01636] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The high-spin S = 2 Mn(III) complex [Mn{(OPPh2)2N}3] (1Mn) exhibits field-induced slow relaxation of magnetization (Inorg. Chem. 2013, 52, 12869). Magnetic susceptibility and dual-mode X-band electron paramagnetic resonance (EPR) studies revealed a negative value of the zero-field-splitting (zfs) parameter D. In order to explore the magnetic and electronic properties of 1Mn in detail, a combination of experimental and computational studies is presented herein. Alternating-current magnetometry on magnetically diluted samples (1Mn/1Ga) of 1Mn in the diamagnetic gallium analogue, [Ga{(OPPh2)2N}3], indicates that the slow relaxation behavior of 1Mn is due to the intrinsic properties of the individual molecules of 1Mn. Investigation of the single-crystal magnetization of both 1Mn and 1Mn/1Ga by a micro-SQUID device reveals hysteresis loops below 1 K. Closed hysteresis loops at a zero direct-current magnetic field are observed and attributed to fast quantum tunneling of magnetization. High-frequency and -field EPR (HFEPR) spectroscopic studies reveal that, apart from the second-order zfs terms (D and E), fourth-order terms (B4m) are required in order to appropriately describe the magnetic properties of 1Mn. These studies provide accurate spin-Hamiltonian (sH) parameters of 1Mn, i.e., zfs parameters |D| = 3.917(5) cm-1, |E| = 0.018(4) cm-1, B04 = B42 = 0, and B44 = (3.6 ± 1.7) × 10-3 cm-1 and g = [1.994(5), 1.996(4), 1.985(4)], and confirm the negative sign of D. Parallel-mode X-band EPR studies on 1Mn/1Ga and CH2Cl2 solutions of 1Mn probe the electronic-nuclear hyperfine interactions in the solid state and solution. The electronic structure of 1Mn is investigated by quantum-chemical calculations by employing recently developed computational protocols that are grounded on ab initio wave function theory. From computational analysis, the contributions of spin-spin and spin-orbit coupling to the magnitude of D are obtained. The calculations provide also computed values of the fourth-order zfs terms B4m, as well as those of the g and hyperfine interaction tensor components. In all cases, a very good agreement between the computed and experimentally determined sH parameters is observed. The magnetization relaxation properties of 1Mn are rationalized on the basis of the composition of the ground-state wave functions in the absence or presence of an external magnetic field.
Collapse
Affiliation(s)
- Yiannis Sanakis
- Institute of Nanoscience and Nanotechnolgy, National Centre of Scientific Research "Demokritos", Aghia Paraskevi 15310, Attiki, Greece
| | - J Krzystek
- National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32310, United States
| | - Dimitrios Maganas
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Alexios Grigoropoulos
- Inorganic Chemistry Laboratory, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece
| | - Eleftherios Ferentinos
- Inorganic Chemistry Laboratory, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece
| | - Marios G Kostakis
- Analytical Chemistry Laboratory, Department of Chemistry, National and Kapodistrian University of Athens,15771 Athens, Greece
| | - Vasiliki Petroulea
- Institute of Nanoscience and Nanotechnolgy, National Centre of Scientific Research "Demokritos", Aghia Paraskevi 15310, Attiki, Greece
| | - Michael Pissas
- Institute of Nanoscience and Nanotechnolgy, National Centre of Scientific Research "Demokritos", Aghia Paraskevi 15310, Attiki, Greece
| | | | - Wolfgang Wernsdorfer
- Physikalisches Institut, Karlsruher Institut für Technologie, Wolfgang-Gaede-Strasse 1, 76131 Karlsruhe, Germany.,Institute of Quantum Materials and Technologies, Karlsruher Institut für Technologie, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Panayotis Kyritsis
- Inorganic Chemistry Laboratory, Department of Chemistry, National and Kapodistrian University of Athens, 15771 Athens, Greece
| |
Collapse
|
11
|
Witwicki M, Walencik PK, Jezierska J. How accurate is density functional theory in predicting spin density? An insight from the prediction of hyperfine coupling constants. J Mol Model 2019; 26:10. [PMID: 31834497 DOI: 10.1007/s00894-019-4268-0] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Accepted: 11/25/2019] [Indexed: 01/30/2023]
Abstract
Electron paramagnetic resonance (EPR) spectroscopy has been proven to be an important technique for studying paramagnetic systems. Probably, the most accessible EPR parameter and the one that provides a significant amount of information about molecular structure and spin density is the hyperfine coupling constant (HFCC). Hence, accurate quantum-chemical modeling of HFCCs is frequently essential to the adequate interpretation of EPR spectra. It requires the precise spin density, which is the difference between the densities of α- and β-electrons, and thus, its quality is expected to reflect the quality of the total electron density. The question of which approximate exchange-correlation density functional yields sufficiently accurate HFCCs, and thus, the spin density remains open. To assess the performance of well-established density functionals for calculating HFCCs, we used a series of 26 small paramagnetic species and compared the obtained results to the CCSD reference values. The performance of DFT was also tested on EPR-studied o-semiquinone radical interacting with water molecules and Mg2+ cation. The HFCCs were additionally calculated by the DLPNO-CCSD method, and this wave function-based technique was found superior to all functionals we tested. Although some functionals were found, on average, to be fairly efficient, we found that the most accurate functional is system-dependent, and therefore, the DLPNO-CCSD method should be preferred for theoretical investigations of the HFCCs and spin density.
Collapse
Affiliation(s)
- Maciej Witwicki
- Faculty of Chemistry, Wrocław University, ul. F. Joliot-Curie 14, 50-383, Wrocław, Poland.
| | - Paulina K Walencik
- Faculty of Chemistry, Wrocław University, ul. F. Joliot-Curie 14, 50-383, Wrocław, Poland
| | - Julia Jezierska
- Faculty of Chemistry, Wrocław University, ul. F. Joliot-Curie 14, 50-383, Wrocław, Poland
| |
Collapse
|
12
|
Breidung J, Thiel W. Equilibrium Structures of the Phosphorus Trihalides PF 3 and PCl 3, and the Phosphoranes PH 3F 2, PF 5, PCl 3F 2, and PCl 5. J Phys Chem A 2019; 123:5600-5612. [PMID: 31181157 PMCID: PMC6750832 DOI: 10.1021/acs.jpca.9b04406] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2019] [Revised: 06/07/2019] [Indexed: 11/30/2022]
Abstract
Among the title species, a reliable and accurate equilibrium geometry ( re structure) is available only for PF3, which has been determined experimentally more than 20 years ago. Here, we report accurate re structures for all title molecules, which were obtained using a composite computational approach based on explicitly correlated coupled-cluster theory (CCSD(T)-F12b) in conjunction with a large correlation-consistent basis set (cc-pCVQZ-F12) to take core-valence electron correlation into account. Additional terms were included to correct for the effects of iterative triple excitations (CCSDT), noniterative quadruple excitations (CCSDT(Q)), and scalar relativistic contributions (DKH2-CCSD(T)). The performance of this computational procedure was established through test calculations on selected small molecules (PH, PF, PCl, PH2, PF2, and PH3). For PF3, PCl3, PH3F2, and PF5 sufficiently accurate experimental ground-state rotational constants from the literature were used to determine semiexperimental re structures, which were found to be in excellent agreement with the corresponding best estimates from the current composite approach. The recommended equilibrium structural parameters are for PCl3, re(PCl) = 203.94 pm and θe(ClPCl) = 100.18°; for PH3F2, re(PHeq) = 138.38 pm and re(PFax) = 164.15 pm; for PF5, re(PFeq) = 153.10 pm and re(PFax) = 157.14 pm; for PCl3F2, re(PCleq) = 200.21 pm and re(PFax) = 159.37 pm; and for PCl5, re(PCleq) = 201.29 pm and re(PClax) = 211.83 pm. The associated uncertainties are estimated to be ±0.10 pm and ±0.10°, respectively.
Collapse
Affiliation(s)
- Jürgen Breidung
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
13
|
Datta D, Gauss J. Accurate Prediction of Hyperfine Coupling Tensors for Main Group Elements Using a Unitary Group Based Rigorously Spin-Adapted Coupled-Cluster Theory. J Chem Theory Comput 2019; 15:1572-1592. [PMID: 30698956 DOI: 10.1021/acs.jctc.8b01048] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the development of a perturbative triples correction scheme for the previously reported unitary group based spin-adapted combinatoric open-shell coupled-cluster (CC) singles and doubles (COS-CCSD) approach and report on the applications of the newly developed method, termed "COS-CCSD(T)", to the calculation of hyperfine coupling (HFC) tensors for radicals consisting of hydrogen, second- and third-row elements. The COS-CCSD(T) method involves a single noniterative step with [Formula: see text] scaling of the computational cost for the calculation of triples corrections to the energy. The key feature of this development is the use of spatial semicanonical orbitals generated from standard restricted open-shell Hartree-Fock (ROHF) orbitals, which allows the unperturbed Hamiltonian operator to be defined in terms of a diagonal spin-free Fock operator. The HFC tensors are computed as a first-order property via implementation of an analytic derivative scheme. The required one-particle spin density matrix is computed by using one- and two-particle spin-free density matrices that are obtained from the analytic derivative implementation, in this way avoiding the use of any spin-dependent operator and maintaining spin adaptation of the CC wavefunction. Benchmark calculations of HFC tensors for a set of 21 radicals indicate reasonably good agreement of the COS-CCSD(T) results with experiment and a consistent improvement over the COS-CCSD method. We demonstrate that the accuracies of the isotropic hyperfine coupling constants obtained in unrestricted HF (UHF) reference based spin-orbital CCSD(T) calculations deteriorate when spin contamination in the UHF wavefunction is large, and the results may even become qualitatively incorrect when spin polarization is the driving mechanism. Within a similar noniterative perturbative treatment of triple excitations, the spin-adapted COS-CCSD(T) approach produces accurate results, thus ensuring cost-effectiveness together with reliability.
Collapse
Affiliation(s)
- Dipayan Datta
- Institut für Physikalische Chemie , Johannes Gutenberg-Universität Mainz , Duesbergweg 10-14 , 55128 Mainz , Germany
| | - Jürgen Gauss
- Institut für Physikalische Chemie , Johannes Gutenberg-Universität Mainz , Duesbergweg 10-14 , 55128 Mainz , Germany
| |
Collapse
|
14
|
Abbott AS, Glick ZL, Schaefer HF. Reinterpretation of the electronic absorption spectrum of the methylene amidogen radical (H 2CN). J Chem Phys 2018; 149:094302. [DOI: 10.1063/1.5042763] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Adam S. Abbott
- Center for Computational Quantum Chemistry, The University of Georgia, Athens, Georgia 30602, USA
| | - Zach L. Glick
- Center for Computational Quantum Chemistry, The University of Georgia, Athens, Georgia 30602, USA
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry, The University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
15
|
Samanta PK, Köhn A. First-order properties from internally contracted multireference coupled-cluster theory with particular focus on hyperfine coupling tensors. J Chem Phys 2018; 149:064101. [DOI: 10.1063/1.5040587] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
| | - Andreas Köhn
- Institut für Theoretische Chemie, Universität Stuttgart, D-70569 Stuttgart, Germany
| |
Collapse
|
16
|
Saitow M, Neese F. Accurate spin-densities based on the domain-based local pair-natural orbital coupled-cluster theory. J Chem Phys 2018; 149:034104. [PMID: 30037259 DOI: 10.1063/1.5027114] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
Exploiting locality in the electron correlation reduces the computational cost for solving the Coupled-Cluster (CC) equations. This is important for making CC theory applicable to routine computational chemistry applications where it promises to deliver results of "gold-standard" quality. Recently, we have proposed a series of CC formulations in the domain-based local pair-natural orbital framework [DLPNO-coupled-cluster with singles and doubles (CCSD) and DLPNO-coupled-cluster singles and doubles with perturbative triples] which are designed to reproduce approximately 99.9% of the canonical correlation energy. In our previous work, the DLPNO-CCSD method has been extended to the high-spin open-shell reference and shown to possess comparable accuracy to the closed-shell counterpart [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)]. The so-called Λ-equations have been formulated in the DLPNO framework for the closed-shell species as an exact derivative of the DLPNO-CCSD Lagrangian with respect to the PNO-based cluster amplitudes [D. Datta et al., J. Chem. Phys. 145, 114101 (2016)]. In this paper, we extend the DLPNO-based Lagrangian scheme to the high-spin open-shell reference cases, thus enabling the accurate computation of the electron- and spin-densities for large open-shell species. We apply this newly developed approach to various first-order electronic and magnetic properties such as isotropic and anisotropic components in the hyperfine coupling interactions and the electric field gradient. We demonstrate that the DLPNO-CCSD results converge toward the respective canonical CC density and also that the DLPNO-CCSD-based properties are more accurate than the conventional density functional theory (DFT) results in real-life applications. The additional computational cost is not more than one energy evaluation in the DLPNO-CCSD framework.
Collapse
Affiliation(s)
- Masaaki Saitow
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
17
|
Vogler S, Savasci G, Ludwig M, Ochsenfeld C. Selected-Nuclei Method for the Computation of Hyperfine Coupling Constants within Second-Order Møller-Plesset Perturbation Theory. J Chem Theory Comput 2018; 14:3014-3024. [PMID: 29762028 DOI: 10.1021/acs.jctc.8b00116] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We introduce a new ansatz to compute hyperfine coupling constants of selected nuclei at the level of second-order Møller-Plesset perturbation (MP2) and double-hybrid density functional theory with reduced computational effort, opening the route to the analyis of hyperfine coupling constants of large molecular structures. Our approach is based on a reformulation of the canonical MP2 term in atomic orbitals, thus exploiting the locality of electron correlation. We show that a perturbation-including integral screening reduces the scaling behavior of the number of significant two-electron integrals to sublinear. This selected-nuclei approach allows for an efficient computation within scaled-opposite spin (SOS) RI-MP2 on massively parallelized architectures such as graphical processor units (GPUs), thus enabling studies on the influence of the environment on hyperfine coupling constants.
Collapse
Affiliation(s)
- Sigurd Vogler
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry , University of Munich (LMU) , Butenandtstrasse 7 , 81377 Munich , Germany
| | - Gökcen Savasci
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry , University of Munich (LMU) , Butenandtstrasse 7 , 81377 Munich , Germany.,Max Planck Institute for Solid State Research, Heisenbergstrasse 1 , 70569 Stuttgart , Germany
| | - Martin Ludwig
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry , University of Munich (LMU) , Butenandtstrasse 7 , 81377 Munich , Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry , University of Munich (LMU) , Butenandtstrasse 7 , 81377 Munich , Germany.,Max Planck Institute for Solid State Research, Heisenbergstrasse 1 , 70569 Stuttgart , Germany
| |
Collapse
|
18
|
Giner E, Tenti L, Angeli C, Ferré N. Computation of the Isotropic Hyperfine Coupling Constant: Efficiency and Insights from a New Approach Based on Wave Function Theory. J Chem Theory Comput 2017; 13:475-487. [DOI: 10.1021/acs.jctc.6b00827] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Emmanuel Giner
- Dipartimento
di Scienze Chimiche e Famaceutiche, Universita di Ferrara, Via Fossato
di Mortara 17, I-44121 Ferrara, Italy
| | - Lorenzo Tenti
- Dipartimento
di Scienze Chimiche e Famaceutiche, Universita di Ferrara, Via Fossato
di Mortara 17, I-44121 Ferrara, Italy
| | - Celestino Angeli
- Dipartimento
di Scienze Chimiche e Famaceutiche, Universita di Ferrara, Via Fossato
di Mortara 17, I-44121 Ferrara, Italy
| | | |
Collapse
|
19
|
Cazzoli G, Lattanzi V, Kirsch T, Gauss J, Tercero B, Cernicharo J, Puzzarini C. Laboratory measurements and astronomical search for the HSO radical. ASTRONOMY AND ASTROPHYSICS 2016; 591:A126. [PMID: 27721513 PMCID: PMC5055096 DOI: 10.1051/0004-6361/201628745] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
CONTEXT Despite the fact that many sulfur-bearing molecules, ranging from simple diatomic species up to astronomical complex molecules, have been detected in the interstellar medium, the sulfur chemistry in space is largely unknown and a depletion in the abundance of S-containing species has been observed in the cold, dense interstellar medium (ISM). The chemical form of the missing sulfur has yet to be identified. AIMS For these reasons, in view of the fact that there is a large abundance of triatomic species harbouring sulfur, oxygen, and hydrogen, we decided to investigate the HSO radical in the laboratory to try its astronomical detection. METHODS High-resolution measurements of the rotational spectrum of the HSO radical were carried out within a frequency range well up into the THz region. Subsequently, a rigorous search for HSO in the two most studied high-mass star-forming regions, Orion KL and Sagittarius (Sgr) B2, and in the cold dark cloud Barnard 1 (B1-b) was performed. RESULTS The frequency coverage and the spectral resolution of our measurements allowed us to improve and extend the existing dataset of spectroscopic parameters, thus enabling accurate frequency predictions up to the THz range. These were used to derive the synthetic spectrum of HSO, by means of the MADEX code, according to the physical parameters of the astronomical source under consideration. For all sources investigated, the lack of HSO lines above the confusion limit of the data is evident. CONCLUSIONS The derived upper limit to the abundance of HSO clearly indicates that this molecule does not achieve significant abundances in either the gas phase or in the ice mantles of dust grains.
Collapse
Affiliation(s)
- Gabriele Cazzoli
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via Selmi 2, I-40126 Bologna, Italy
| | - Valerio Lattanzi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via Selmi 2, I-40126 Bologna, Italy; Max-Planck-Institut für Extraterrestrische Physik, Giessenbachstraße 1, D-85748 Garching, Germany
| | - Till Kirsch
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| | - Belén Tercero
- Grupo de Astrofísica Molecular. Instituto de CC. de Materiales de Madrid (ICMM-CSIC). Sor Juana Inés de la Cruz 3, Cantoblanco, 28049 Madrid, Spain
| | - José Cernicharo
- Grupo de Astrofísica Molecular. Instituto de CC. de Materiales de Madrid (ICMM-CSIC). Sor Juana Inés de la Cruz 3, Cantoblanco, 28049 Madrid, Spain
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via Selmi 2, I-40126 Bologna, Italy
| |
Collapse
|
20
|
Datta D, Gauss J. Communication: Spin densities within a unitary group based spin-adapted open-shell coupled-cluster theory: Analytic evaluation of isotropic hyperfine-coupling constants for the combinatoric open-shell coupled-cluster scheme. J Chem Phys 2015; 143:011101. [DOI: 10.1063/1.4923436] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Dipayan Datta
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| |
Collapse
|
21
|
Lan TN, Kurashige Y, Yanai T. Toward Reliable Prediction of Hyperfine Coupling Constants Using Ab Initio Density Matrix Renormalization Group Method: Diatomic 2Σ and Vinyl Radicals as Test Cases. J Chem Theory Comput 2014; 10:1953-67. [DOI: 10.1021/ct400978j] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Tran Nguyen Lan
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Yuki Kurashige
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
- Department
of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
| | - Takeshi Yanai
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
- Department
of Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
22
|
Fišer J, Polák R. On the theoretical description of nuclear quadrupole coupling in Π states of small molecules. Chem Phys 2013. [DOI: 10.1016/j.chemphys.2013.08.007] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
23
|
Rey-Villaverde R, Cybulski H, Flores JR, Fernández B. A high-accuracy theoretical study of the CHnP Systemsn = 1-3. J Comput Chem 2013; 34:2020-31. [DOI: 10.1002/jcc.23357] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2013] [Revised: 05/23/2013] [Accepted: 05/29/2013] [Indexed: 01/06/2023]
Affiliation(s)
- Ringo Rey-Villaverde
- Departamento de Química Física, Facultad de Química; Universidad de Vigo; Vigo; E-36310; Spain
| | | | - Jesús R. Flores
- Departamento de Química Física, Facultad de Química; Universidad de Vigo; Vigo; E-36310; Spain
| | - Berta Fernández
- Center for Research in Biological Chemistry and Molecular Materials (CIQUS), University of Santiago de Compostela; Santiago; E-15782; Spain
| |
Collapse
|
24
|
Barone V, Biczysko M, Bloino J, Egidi F, Puzzarini C. Accurate structure, thermodynamics, and spectroscopy of medium-sized radicals by hybrid coupled cluster/density functional theory approaches: the case of phenyl radical. J Chem Phys 2013; 138:234303. [PMID: 23802956 PMCID: PMC4606980 DOI: 10.1063/1.4810863] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The coupled-cluster singles doubles model with perturbative treatment of triples (CCSD(T)) coupled with extrapolation to the complete basis-set limit and additive approaches represent the "golden standard" for the structural and spectroscopic characterization of building blocks of biomolecules and nanosystems. However, when open-shell systems are considered, additional problems related to both specific computational difficulties and the need of obtaining spin-dependent properties appear. In this contribution, we present a comprehensive study of the molecular structure and spectroscopic (IR, Raman, EPR) properties of the phenyl radical with the aim of validating an accurate computational protocol able to deal with conjugated open-shell species. We succeeded in obtaining reliable and accurate results, thus confirming and, partly, extending the available experimental data. The main issue to be pointed out is the need of going beyond the CCSD(T) level by including a full treatment of triple excitations in order to fulfil the accuracy requirements. On the other hand, the reliability of density functional theory in properly treating open-shell systems has been further confirmed.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | | | | | | | | |
Collapse
|
25
|
Abstract
Rotational spectroscopy is known to be a technique that is widely used to infer information on molecular structure and dynamics. In the last few decades, its role in the field of atmospheric and astrophysical investigations has rapidly grown. However, several are the challenging aspects in rotational spectroscopy, since the detection and analysis of spectra as well as interpretation of obtained results are not at all straightforward. Quantum chemistry has reached such an accuracy that can be used to disentangle these challenging situations by guiding the experimental investigation, assisting in the determination of the spectroscopic parameters, and extracting information of chemical interest. This perspective provides an overview of the theoretical background and computational requirements needed for the accurate evaluation of the spectroscopic parameters of relevance to rotational spectroscopy. The role of theory in guiding and supporting experiment is detailed through a few examples, and the interplay of experiment and theory is discussed in terms of the information of physical and chemical interest that can be derived.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica Giacomo Ciamician, Università di Bologna, Via Selmi 2, I-40126 Bologna, Italy.
| |
Collapse
|
26
|
Barone V, Biczysko M, Bloino J, Puzzarini C. Characterization of the Elusive Conformers of Glycine from State-of-the-Art Structural, Thermodynamic, and Spectroscopic Computations: Theory Complements Experiment. J Chem Theory Comput 2013; 9:1533-47. [DOI: 10.1021/ct3010672] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| | - Malgorzata Biczysko
- Center for Nanotechnology Innovation
@NEST, Istituto Italiano di Tecnologia, Piazza San Silvestro 12, I-56127
Pisa, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
- Consiglio Nazionale
delle Ricerche,
Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), UOS
di Pisa, Area della Ricerca CNR, Via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica “G.
Ciamician,” Università di Bologna, Via F. Selmi 2, 40126
Bologna, Italy
| |
Collapse
|
27
|
Barone V, Biczysko M, Bloino J, Puzzarini C. Glycine conformers: a never-ending story? Phys Chem Chem Phys 2013; 15:1358-63. [DOI: 10.1039/c2cp43884d] [Citation(s) in RCA: 76] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
28
|
Breiten B, Jordan M, Taura D, Zalibera M, Griesser M, Confortin D, Boudon C, Gisselbrecht JP, Schweizer WB, Gescheidt G, Diederich F. Donor-Substituted Octacyano[4]dendralenes: Investigation of π-Electron Delocalization in Their Radical Ions. J Org Chem 2012; 78:1760-7. [DOI: 10.1021/jo301194y] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Benjamin Breiten
- Laboratory of Organic Chemistry, ETH Zurich, Hönggerberg, HCI, CH-8093 Zurich, Switzerland
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
29
|
Puzzarini C. A theoretical study of the CX2N radicals (X = F, Cl, Br): The effect of halogen substitution on structure, isomerization, and energetics. J Chem Phys 2012; 136:044316. [DOI: 10.1063/1.3678006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
|
30
|
Mok DKW, Lee EPF, Chau FT, Dyke JM. Franck-Condon simulation, including anharmonicity, of the photodetachment spectrum of P2H(-): restricted-spin coupled-cluster single-double plus perturbative triple and unrestricted-spin coupled-cluster single-double plus perturbative triple -F12x potential energy functions of P2H and P2H(-). J Chem Phys 2011; 135:124312. [PMID: 21974527 DOI: 10.1063/1.3640037] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Geometry optimization and harmonic vibrational frequency calculations have been carried out on the X̃(2)A(') state of P(2)H and the X̃(1)A(') state of P(2)H(-) using the restricted-spin coupled-cluster single-double plus perturbative triple excitation [RCCSD(T)] and explicitly correlated unrestricted-spin coupled-cluster single-double plus perturbative triple excitation [UCCSD(T)-F12x] methods. For RCCSD(T) calculations, basis sets of up to the augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality were employed, and contributions from extrapolation to the complete basis set limit and from core correlation of the P 2s(2)2p(6) electrons were also included. For UCCSD(T)-F12x calculations, different atomic orbital basis sets of triple-zeta quality with different associated complementary auxiliary basis sets and different geminal Slater exponents were used. When the P 2s(2)2p(6) core electrons were correlated in these F12x calculations, appropriate core-valence basis sets were employed. In addition, potential energy functions (PEFs) of the X̃(2)A(') state of P(2)H and the X̃(1)A(') state of P(2)H(-) were computed at different RCCSD(T) and UCCSD(T)-F12x levels, and were used in variational calculations of anharmonic vibrational wavefunctions, which were then utilized to calculate Franck-Condon factors (FCFs) between these two states, employing a method which includes allowance for anharmonicity and Duschinsky rotation. The photodetachment spectrum of P(2)H(-) was then simulated using the computed FCFs. Simulated spectra obtained using the RCCSD(T)/aug-cc-pV5Z and UCCSD(T)-F12x(x = a or b)/aug-cc-pCVTZ PEFs are compared and found to be essentially identical. Based on the computed FCFs, a more detailed assignment of the observed vibrational structure than previously reported, which includes "hot bands," has been proposed. Comparison between simulated and available experimental spectra has been made, and the currently most reliable sets of equilibrium geometrical parameters for P(2)H and its anion have been derived. The photodetachment spectrum of P(2)D, yet to be recorded, has also been simulated.
Collapse
Affiliation(s)
- Daniel K W Mok
- Department of Applied Biology and Chemical Technology, The Hong Kong Polytechnic University, Hung Hom, Hong Kong.
| | | | | | | |
Collapse
|