1
|
Szczuka C, Eichel RA, Granwehr J. Gauging the importance of structural parameters for hyperfine coupling constants in organic radicals. RSC Adv 2023; 13:14565-14574. [PMID: 37188254 PMCID: PMC10177955 DOI: 10.1039/d3ra02476h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Accepted: 05/08/2023] [Indexed: 05/17/2023] Open
Abstract
The identification of fundamental relationships between atomic configuration and electronic structure typically requires experimental empiricism or systematic theoretical studies. Here, we provide an alternative statistical approach to gauge the importance of structure parameters, i.e., bond lengths, bond angles, and dihedral angles, for hyperfine coupling constants in organic radicals. Hyperfine coupling constants describe electron-nuclear interactions defined by the electronic structure and are experimentally measurable, for example, by electron paramagnetic resonance spectroscopy. Importance quantifiers are computed with the machine learning algorithm neighborhood components analysis using molecular dynamics trajectory snapshots. Atomic-electronic structure relationships are visualized in matrices correlating structure parameters with coupling constants of all magnetic nuclei. Qualitatively, the results reproduce common hyperfine coupling models. Tools to use the presented procedure for other radicals/paramagnetic species or other atomic structure-dependent parameters are provided.
Collapse
Affiliation(s)
- Conrad Szczuka
- Institute of Energy and Climate Research (IEK-9), Forschungszentrum Jülich GmbH 52425 Jülich Germany
| | - Rüdiger-A Eichel
- Institute of Energy and Climate Research (IEK-9), Forschungszentrum Jülich GmbH 52425 Jülich Germany
- Institute of Physical Chemistry, RWTH Aachen University 52056 Aachen Germany
| | - Josef Granwehr
- Institute of Energy and Climate Research (IEK-9), Forschungszentrum Jülich GmbH 52425 Jülich Germany
- Institute of Technical and Macromolecular Chemistry, RWTH Aachen University 52056 Aachen Germany
| |
Collapse
|
2
|
Giovannini T, Marrazzini G, Scavino M, Koch H, Cappelli C. Integrated Multiscale Multilevel Approach to Open Shell Molecular Systems. J Chem Theory Comput 2023; 19:1446-1456. [PMID: 36780359 PMCID: PMC10018740 DOI: 10.1021/acs.jctc.2c00805] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
Abstract
We present a novel multiscale approach to study the electronic structure of open shell molecular systems embedded in an external environment. The method is based on the coupling of multilevel Hartree-Fock (MLHF) and Density Functional Theory (MLDFT), suitably extended to the unrestricted formalism, to Molecular Mechanics (MM) force fields (FF). Within the ML region, the system is divided into active and inactive parts, thus describing the most relevant interactions (electrostatic, polarization, and Pauli repulsion) at the quantum level. The surrounding MM part, which is formulated in terms of nonpolarizable or polarizable FFs, permits a physically consistent treatment of long-range electrostatics and polarization effects. The approach is extended to the calculation of hyperfine coupling constants and applied to selected nitroxyl radicals in an aqueous solution.
Collapse
Affiliation(s)
| | - Gioia Marrazzini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Scavino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.,Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
3
|
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
|
4
|
Effects of non-local exchange functionals in the density functional theories for the description of molecular vibrations. J CHEM SCI 2022. [DOI: 10.1007/s12039-022-02061-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
5
|
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
|
6
|
Edeleva M, Van Steenberge PH, Sabbe MK, D’hooge DR. Connecting Gas-Phase Computational Chemistry to Condensed Phase Kinetic Modeling: The State-of-the-Art. Polymers (Basel) 2021; 13:3027. [PMID: 34577928 PMCID: PMC8467432 DOI: 10.3390/polym13183027] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Revised: 09/02/2021] [Accepted: 09/03/2021] [Indexed: 02/06/2023] Open
Abstract
In recent decades, quantum chemical calculations (QCC) have increased in accuracy, not only providing the ranking of chemical reactivities and energy barriers (e.g., for optimal selectivities) but also delivering more reliable equilibrium and (intrinsic/chemical) rate coefficients. This increased reliability of kinetic parameters is relevant to support the predictive character of kinetic modeling studies that are addressing actual concentration changes during chemical processes, taking into account competitive reactions and mixing heterogeneities. In the present contribution, guidelines are formulated on how to bridge the fields of computational chemistry and chemical kinetics. It is explained how condensed phase systems can be described based on conventional gas phase computational chemistry calculations. Case studies are included on polymerization kinetics, considering free and controlled radical polymerization, ionic polymerization, and polymer degradation. It is also illustrated how QCC can be directly linked to material properties.
Collapse
Affiliation(s)
- Mariya Edeleva
- Laboratory for Chemical Technology (LCT), Ghent University, Technologiepark 125, 9052 Zwijnaarde, Belgium; (P.H.M.V.S.); (M.K.S.)
| | - Paul H.M. Van Steenberge
- Laboratory for Chemical Technology (LCT), Ghent University, Technologiepark 125, 9052 Zwijnaarde, Belgium; (P.H.M.V.S.); (M.K.S.)
| | - Maarten K. Sabbe
- Laboratory for Chemical Technology (LCT), Ghent University, Technologiepark 125, 9052 Zwijnaarde, Belgium; (P.H.M.V.S.); (M.K.S.)
- Industrial Catalysis and Adsorption Technology (INCAT), Ghent University, Valentin Vaerwyckweg 1, 9000 Ghent, Belgium
| | - Dagmar R. D’hooge
- Laboratory for Chemical Technology (LCT), Ghent University, Technologiepark 125, 9052 Zwijnaarde, Belgium; (P.H.M.V.S.); (M.K.S.)
- Centre for Textile Science and Engineering (CTSE), Ghent University, Technologiepark 70a, 9052 Zwijnaarde, Belgium
| |
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
|
Vogler S, Dietschreit JCB, Peters LDM, Ochsenfeld C. Important components for accurate hyperfine coupling constants: electron correlation, dynamic contributions, and solvation effects. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1772515] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Sigurd Vogler
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | | | - Laurens D. M. Peters
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| |
Collapse
|
9
|
Auer AA, Tran VA, Sharma B, Stoychev GL, Marx D, Neese F. A case study of density functional theory and domain-based local pair natural orbital coupled cluster for vibrational effects on EPR hyperfine coupling constants: vibrational perturbation theory versus ab initio molecular dynamics. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1797916] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
| | - Van Anh Tran
- MPI für Kohlenforschung, Mülheim an der Ruhr, Germany
| | - Bikramjit Sharma
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum, Germany
| | | | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum, Germany
| | - Frank Neese
- MPI für Kohlenforschung, Mülheim an der Ruhr, Germany
| |
Collapse
|
10
|
A never-ending story in the sky: The secrets of chemical evolution. Phys Life Rev 2020; 32:59-94. [DOI: 10.1016/j.plrev.2019.07.001] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2019] [Revised: 06/26/2019] [Accepted: 07/02/2019] [Indexed: 01/13/2023]
|
11
|
Melosso M, Melli A, Spada L, Zheng Y, Chen J, Li M, Lu T, Feng G, Gou Q, Dore L, Barone V, Puzzarini C. Rich Collection of n-Propylamine and Isopropylamine Conformers: Rotational Fingerprints and State-of-the-Art Quantum Chemical Investigation. J Phys Chem A 2020; 124:1372-1381. [PMID: 31985228 PMCID: PMC7997563 DOI: 10.1021/acs.jpca.9b11767] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
![]()
The
conformational isomerism of isopropylamine and n-propylamine
has been investigated by means of an integrated strategy combining
high-level quantum-chemical calculations and high-resolution rotational
spectroscopy. The equilibrium structures (and thus equilibrium rotational
constants) as well as relative energies of all conformers have been
computed using the so-called “cheap” composite scheme,
which combines the coupled-cluster methodology with second-order Møller–Plesset
perturbation theory for extrapolation to the complete basis set. Methods
rooted in the density functional theory have been instead employed
for computing spectroscopic parameters and for accounting for vibrational
effects. Guided by quantum-chemical predictions, the rotational spectra
of isopropylamine and n-propylamine have been investigated between
2 and 400 GHz with Fourier transform microwave and frequency-modulation
millimeter/submillimeter spectrometers. Spectral assignments confirmed
the presence of several conformers with comparable stability and pointed
out possible Coriolis resonance effects between some of them.
Collapse
Affiliation(s)
- Mattia Melosso
- Dipartimento di Chimica "Giacomo Ciamician" , Università di Bologna , Via Selmi 2 , 40126 Bologna , Italy
| | - Alessio Melli
- Dipartimento di Chimica "Giacomo Ciamician" , Università di Bologna , Via Selmi 2 , 40126 Bologna , Italy
| | - Lorenzo Spada
- Dipartimento di Chimica "Giacomo Ciamician" , Università di Bologna , Via Selmi 2 , 40126 Bologna , Italy.,Scuola Normale Superiore , Piazza dei Cavalieri 7 , 56126 Pisa , Italy
| | - Yang Zheng
- Department of Chemistry, School of Chemistry and Chemical Engineering , Chongqing University , Daxuecheng South Rd. 55 , 401331 Chongqing , China
| | - Junhua Chen
- Department of Chemistry, School of Chemistry and Chemical Engineering , Chongqing University , Daxuecheng South Rd. 55 , 401331 Chongqing , China
| | - Meng Li
- Department of Chemistry, School of Chemistry and Chemical Engineering , Chongqing University , Daxuecheng South Rd. 55 , 401331 Chongqing , China
| | - Tao Lu
- Department of Chemistry, School of Chemistry and Chemical Engineering , Chongqing University , Daxuecheng South Rd. 55 , 401331 Chongqing , China
| | - Gang Feng
- Department of Chemistry, School of Chemistry and Chemical Engineering , Chongqing University , Daxuecheng South Rd. 55 , 401331 Chongqing , China
| | - Qian Gou
- Department of Chemistry, School of Chemistry and Chemical Engineering , Chongqing University , Daxuecheng South Rd. 55 , 401331 Chongqing , China
| | - Luca Dore
- Dipartimento di Chimica "Giacomo Ciamician" , Università di Bologna , Via Selmi 2 , 40126 Bologna , Italy
| | - Vincenzo Barone
- Scuola Normale Superiore , Piazza dei Cavalieri 7 , 56126 Pisa , Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician" , Università di Bologna , Via Selmi 2 , 40126 Bologna , Italy
| |
Collapse
|
12
|
Puzzarini C, Barone V. The challenging playground of astrochemistry: an integrated rotational spectroscopy - quantum chemistry strategy. Phys Chem Chem Phys 2020; 22:6507-6523. [PMID: 32163090 DOI: 10.1039/d0cp00561d] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
While it is now well demonstrated that the interstellar medium (ISM) is characterized by a diverse and complex chemistry, a significant number of features in radioastronomical spectra are still unassigned and call for new laboratory efforts, which are increasingly based on integrated experimental and computational strategies. In parallel, the identification of an increasing number of molecules containing more than five atoms and at least one carbon atom (the so-called "interstellar" complex organic molecules), which can play a relevant role in the chemistry of life, raises the additional issue of how these species can be produced in the typical harsh conditions of the ISM. On these grounds, this perspective aims to present an integrated rotational spectroscopy - quantum chemistry approach for supporting radioastronomical observations and a computational strategy for contributing to the elucidation of chemical reactivity in the interstellar space.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", University of Bologna, via F. Selmi 2, I-40126 Bologna, Italy.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, I-56126, Italy
| |
Collapse
|
13
|
Puzzarini C, Bloino J, Tasinato N, Barone V. Accuracy and Interpretability: The Devil and the Holy Grail. New Routes across Old Boundaries in Computational Spectroscopy. Chem Rev 2019; 119:8131-8191. [DOI: 10.1021/acs.chemrev.9b00007] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
14
|
Giovannini T, Lafiosca P, Chandramouli B, Barone V, Cappelli C. Effective yet reliable computation of hyperfine coupling constants in solution by a QM/MM approach: Interplay between electrostatics and non-electrostatic effects. J Chem Phys 2019; 150:124102. [DOI: 10.1063/1.5080810] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Piero Lafiosca
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Balasubramanian Chandramouli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Compunet, Istituto Italiano di Tecnologia (IIT), Via Morego 30, 16163 Genova, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
15
|
Gromov OI, Kuzin SV, Golubeva EN. Performance of DFT methods in the calculation of isotropic and dipolar contributions to 14N hyperfine coupling constants of nitroxide radicals. J Mol Model 2019; 25:93. [PMID: 30859325 DOI: 10.1007/s00894-019-3966-y] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2018] [Accepted: 02/07/2019] [Indexed: 10/27/2022]
Abstract
In the present study, we tested the widely used density functionals BP86, PBE, OLYP, TPSS, M06-L, B3LYP, PBE0, mPW1PW, B97, BHandHLYP, TPSS0, M06, M06-2X, CAM-B3LYP, ωB97x, and B2PLYP with the cc-pCVQZ basis set in calculations on a set of 23 nitroxide radicals with well-resolved 14N anisotropic hyperfine coupling (HFC) constants. The results were compared with those obtained using the B3LYP/N07D and PBE/N07D methods. The convergence of the HFC values to the complete basis set limit is briefly discussed. The best results were obtained using the M06/COSMO method, with a mean absolute deviation (MAD) of 0.4 G for the dipole-dipole contribution and MAD = 0.6 G for the contact coupling contribution (as compared to 1.1 G and 1.0 G, respectively, for the B3LYP/N07D/COSMO method and 1.7 G and 0.5 G, respectively, for the B3LYP/N07D method). The majority of the functionals yielded satisfactory results for the dipole-dipole contribution, but only the M06 functional yielded similar errors for both the dipole-dipole and isotropic contributions. The RIJCOSX and RI approximations introduced errors equal to or smaller than 0.01 G.
Collapse
Affiliation(s)
- Oleg I Gromov
- Chemistry Department, Lomonosov Moscow State University, Leninskiye Gory 1-3, Moscow, Russia, 119991.
| | - Sergei V Kuzin
- Chemistry Department, Lomonosov Moscow State University, Leninskiye Gory 1-3, Moscow, Russia, 119991
| | - Elena N Golubeva
- Chemistry Department, Lomonosov Moscow State University, Leninskiye Gory 1-3, Moscow, Russia, 119991
| |
Collapse
|
16
|
Baiardi A, Bloino J, Barone V. Time-Dependent Formulation of Resonance Raman Optical Activity Spectroscopy. J Chem Theory Comput 2018; 14:6370-6390. [PMID: 30281300 DOI: 10.1021/acs.jctc.8b00488] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
In this work, we extend the theoretical framework recently developed for the simulation of resonance Raman (RR) spectra of medium-to-large sized systems to its chiral counterpart, namely, resonance Raman optical activity (RROA). The theory is based on a time-dependent (TD) formulation, with the transition tensors obtained as half-Fourier transforms of the appropriate cross-correlation functions. The implementation has been kept as general as possible, supporting adiabatic and vertical models for the PES representation, both in Cartesian and internal coordinates, with the possible inclusion of Herzberg-Teller (HT) effects. Thanks to the integration of this TD-RROA procedure within a general-purpose quantum-chemistry program, both solvation and leading anharmonicity effects can be included in an effective way. The implementation is validated on one of the smallest chiral molecule (methyloxirane). Practical applications are illustrated with three medium-size organic molecules (naproxen-OCD3, quinidine and 2-Br-hexahelicene), whose simulated spectra are compared to the corresponding experimental data.
Collapse
Affiliation(s)
- Alberto Baiardi
- Scuola Normale Superiore , piazza dei Cavalieri 7 , I-56126 Pisa , Italy
| | - Julien Bloino
- Scuola Normale Superiore , piazza dei Cavalieri 7 , I-56126 Pisa , Italy
| | - Vincenzo Barone
- Scuola Normale Superiore , piazza dei Cavalieri 7 , I-56126 Pisa , Italy
| |
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
|
Baiardi A, Paoloni L, Barone V, Zakrzewski VG, Ortiz JV. Assessment of Electron Propagator Methods for the Simulation of Vibrationally Resolved Valence and Core Photoionization Spectra. J Chem Theory Comput 2017; 13:3120-3135. [PMID: 28521087 PMCID: PMC5732571 DOI: 10.1021/acs.jctc.6b00958] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The analysis of photoelectron spectra is usually facilitated by quantum mechanical simulations. Because of the recent improvement of experimental techniques, the resolution of experimental spectra is rapidly increasing, and the inclusion of vibrational effects is usually mandatory to obtain a reliable reproduction of the spectra. With the aim of defining a robust computational protocol, a general time-independent formulation to compute different kinds of vibrationally resolved electronic spectra has been generalized to also support photoelectron spectroscopy. The electronic structure data underlying the simulation are computed using different electron propagator approaches. In addition to the more standard approaches, a new and robust implementation of the second-order self-energy approximation of the electron propagator based on a transition operator reference (TOEP2) is presented. To validate our implementation, a series of molecules has been used as test cases. The result of the simulations shows that, for ultraviolet photoionization spectra, the more accurate nondiagonal approaches are needed to obtain a reliable reproduction of vertical ionization energies but that diagonal approaches are sufficient for energy gradients and pole strengths. For X-ray photoelectron spectroscopy, the TOEP2 approach, besides being more efficient, is also the most accurate in the reproduction of both vertical ionization energies and vibrationally resolved bandshapes.
Collapse
Affiliation(s)
- A Baiardi
- Scuola Normale Superiore , Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - L Paoloni
- Scuola Normale Superiore , Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - V Barone
- Scuola Normale Superiore , Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - V G Zakrzewski
- Department of Chemistry and Biochemistry, Auburn University , Auburn, Alabama 36849-5312, United States
| | - J V Ortiz
- Department of Chemistry and Biochemistry, Auburn University , Auburn, Alabama 36849-5312, United States
| |
Collapse
|
19
|
Baiardi A, Bloino J, Barone V. Simulation of Vibronic Spectra of Flexible Systems: Hybrid DVR-Harmonic Approaches. J Chem Theory Comput 2017; 13:2804-2822. [PMID: 28467703 DOI: 10.1021/acs.jctc.7b00236] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Our general framework for the simulation of vibrational signatures in electronic spectra has been extended to treat one large-amplitude motion (LAM) at the anharmonic level, coupled to the other small-amplitude motions (SAM) treated as harmonic. The coupling between LAM and SAM is minimized thanks to the use of delocalized internal coordinates, which are built automatically from the molecular topology. General LAMs can be employed, ranging from intrinsic reaction coordinates to rigid or flexible paths based on the distinguished coordinate approach. The anharmonic model is based on a fully numerical method based on the discrete variable representation (DVR) theory, supporting different types of boundary conditions. The inclusion of this model in a general-purpose electronic structure code makes available to the user a large panel of quantum chemistry models, for both isolated and condensed phases. The flexibility and reliability of the new framework are illustrated by some case studies, covering various types of LAMs, ranging from a small test case, the photoelectron spectrum of ammonia, to larger systems, such as phenylanthracene and cyclobutanone.
Collapse
Affiliation(s)
- Alberto Baiardi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 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
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
20
|
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
|
21
|
Kirste B. DFT calculations of hyperfine coupling constants of organic π radicals and comparison with empirical equations and experiment. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2016; 54:835-841. [PMID: 27334701 DOI: 10.1002/mrc.4467] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2015] [Revised: 05/26/2016] [Accepted: 05/27/2016] [Indexed: 06/06/2023]
Abstract
Density functional theory calculations (UB3LYP/EPR-III) for a series of radicals and radical ions were performed to check the internal consistency of the method and its implications to the theoretical concepts of electron paramagnetic resonance such as π-σ spin polarization, hyperconjugation and phenyl hyperconjugation. In the second part, experimental data for seven radicals (43 hyperfine coupling constants) are compared with calculations, yielding a correlation of r2 = 0.97. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Burkhard Kirste
- Institut für Chemie und Biochemie, Freie Universität Berlin, Takustrasse 3, 14195, Berlin, Germany
| |
Collapse
|
22
|
Baiardi A, Bloino J, Barone V. General formulation of vibronic spectroscopy in internal coordinates. J Chem Phys 2016; 144:084114. [PMID: 26931688 DOI: 10.1063/1.4942165] [Citation(s) in RCA: 51] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Our general platform integrating time-independent and time-dependent evaluations of vibronic effects at the harmonic level for different kinds of absorption and emission one-photon, conventional and chiral spectroscopies has been extended to support various sets of internal coordinates. Thanks to the implementation of analytical first and second derivatives of different internal coordinates with respect to cartesian ones, both vertical and adiabatic models are available, with the inclusion of mode mixing and, possibly, Herzberg-Teller contributions. Furthermore, all supported non-redundant sets of coordinates are built from a fully automatized algorithm using only a primitive redundant set derived from a bond order-based molecular topology. Together with conventional stretching, bending, and torsion coordinates, the availability of additional coordinates (including linear and out-of-plane bendings) allows a proper treatment of specific systems, including, for instance, inter-molecular hydrogen bridges. A number of case studies are analysed, showing that cartesian and internal coordinates are nearly equivalent for semi-rigid systems not experiencing significant geometry distortions between initial and final electronic states. At variance, delocalized (possibly weighted) internal coordinates become much more effective than their cartesian counterparts for flexible systems and/or in the presence of significant geometry distortions accompanying electronic transitions.
Collapse
Affiliation(s)
- Alberto Baiardi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Julien Bloino
- 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
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
23
|
Hodecker M, Biczysko M, Dreuw A, Barone V. Simulation of Vacuum UV Absorption and Electronic Circular Dichroism Spectra of Methyl Oxirane: The Role of Vibrational Effects. J Chem Theory Comput 2016; 12:2820-33. [PMID: 27159495 PMCID: PMC5612404 DOI: 10.1021/acs.jctc.6b00121] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Vibrationally resolved one-photon absorption and electronic circular dichroism spectra of (R)-methyl oxirane were calculated with different electronic and vibronic models selecting, through an analysis of the convergence of the results, the best compromise between reliability and computational cost. Linear-response TD-DFT/CAM-B3LYP/SNST electronic computations in conjunction with the simple vertical gradient vibronic model were chosen and employed for systematic comparison with the available experimental data. Remarkable agreement between simulated and experimental spectra was achieved for both one-photon absorption and circular dichroism concerning peak positions, relative intensities, and general spectral shapes considering the computational efficiency of the chosen theoretical approach. The significant improvement of the results with respect to smearing of vertical electronic transitions by phenomenological Gaussian functions and the possible inclusion of solvent effects by polarizable continuum models at a negligible additional cost paves the route toward the simulation and analysis of spectral shapes of complex molecular systems in their natural environment.
Collapse
Affiliation(s)
- Manuel Hodecker
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 368, 69120 Heidelberg, Germany
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, College of Sciences, Shanghai University, 99 Shangda Road, Shanghai, 200444 China
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 368, 69120 Heidelberg, Germany
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
24
|
Fornaro T, Biczysko M, Bloino J, Barone V. Reliable vibrational wavenumbers for C=O and N-H stretchings of isolated and hydrogen-bonded nucleic acid bases. Phys Chem Chem Phys 2016; 18:8479-90. [PMID: 26940362 PMCID: PMC5612391 DOI: 10.1039/c5cp07386c] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The accurate prediction of vibrational wavenumbers for functional groups involved in hydrogen-bonded bridges remains an important challenge for computational spectroscopy. For the specific case of the C=O and N-H stretching modes of nucleobases and their oligomers, the paucity of experimental reference values needs to be compensated by reliable computational data, which require the use of approaches going beyond the standard harmonic oscillator model. Test computations performed for model systems (formamide, acetamide and their cyclic homodimers) in the framework of the second order vibrational perturbation theory (VPT2) confirmed that anharmonic corrections can be safely computed by global hybrid (GHF) or double hybrid (DHF) functionals, whereas the harmonic part is particularly challenging. As a matter of fact, GHFs perform quite poorly and even DHFs, while fully satisfactory for C=O stretchings, face unexpected difficulties when dealing with N-H stretchings. On these grounds, a linear regression for N-H stretchings has been obtained and validated for the heterodimers formed by 4-aminopyrimidine with 6-methyl-4-pyrimidinone (4APM-M4PMN) and by uracil with water. In view of the good performance of this computational model, we have built a training set of B2PLYP-D3/maug-cc-pVTZ harmonic wavenumbers (including linear regression scaling for N-H) for six-different uracil dimers and a validation set including 4APM-M4PMN, one of the most stable hydrogen-bonded adenine homodimers, as well as the adenine-uracil, adenine-thymine, guanine-cytosine and adenine-4-thiouracil heterodimers. Because of the unfavourable scaling of DHF harmonic wavenumbers with the dimensions of the investigated systems, we have optimized a linear regression of B3LYP-D3/N07D harmonic wavenumbers for the training set, which has been next checked against the validation set. This relatively cheap model, which shows very good agreement with experimental data (average errors of about 10 cm(-1)), paves the route toward a reliable analysis of spectroscopic signatures for larger polynucleotides.
Collapse
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, College of Sciences, Shanghai University, 99 Shangda Road, Shanghai, 200444, China.
| | - Julien Bloino
- 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
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| |
Collapse
|
25
|
Barone V, Cimino P, Stendardo E. Development and Validation of the B3LYP/N07D Computational Model for Structural Parameter and Magnetic Tensors of Large Free Radicals. J Chem Theory Comput 2015; 4:751-64. [PMID: 26621090 DOI: 10.1021/ct800034c] [Citation(s) in RCA: 213] [Impact Index Per Article: 23.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Extensive calculations on a large set of free radicals containing atoms of the second and third row show that the B3LYP/N07D computational model provides remarkably accurate structural parameters and magnetic tensors at reasonable computational costs. The key of this success is the optimization of core-valence s functions for hyperfine coupling constants, while retaining (and even improving) the good performances of the parent 6-31+G(d,p) basis set for valence properties through reoptimization of polarization and diffuse p functions.
Collapse
Affiliation(s)
- Vincenzo Barone
- LSDM and INSTM-Village, Dipartimento di Chimica 'Paolo Corradini', Complesso Universitario Monte S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Paola Cimino
- LSDM and INSTM-Village, Dipartimento di Chimica 'Paolo Corradini', Complesso Universitario Monte S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Emiliano Stendardo
- LSDM and INSTM-Village, Dipartimento di Chimica 'Paolo Corradini', Complesso Universitario Monte S. Angelo, via Cintia, I-80126 Napoli, Italy
| |
Collapse
|
26
|
Baiardi A, Bloino J, Barone V. Accurate Simulation of Resonance-Raman Spectra of Flexible Molecules: An Internal Coordinates Approach. J Chem Theory Comput 2015; 11:3267-80. [DOI: 10.1021/acs.jctc.5b00241] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Alberto Baiardi
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 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
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
27
|
Fornaro T, Burini D, Biczysko M, Barone V. Hydrogen-Bonding Effects on Infrared Spectra from Anharmonic Computations: Uracil–Water Complexes and Uracil Dimers. J Phys Chem A 2015; 119:4224-36. [DOI: 10.1021/acs.jpca.5b01561] [Citation(s) in RCA: 99] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| | - Diletta Burini
- Dipartimento
di Matematica e Informatica, Università di Perugia, INFN Sezione
Perugia Via Vanvitelli, I-106123 Perugia, Italy
| | - Malgorzata Biczysko
- Physics
Department, and International Centre for Quantum and Molecular Structures, Shanghai University, 99 Shangda Road, Shanghai, 200444 China
- 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
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
28
|
Bloino J. A VPT2 Route to Near-Infrared Spectroscopy: The Role of Mechanical and Electrical Anharmonicity. J Phys Chem A 2015; 119:5269-87. [DOI: 10.1021/jp509985u] [Citation(s) in RCA: 76] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Julien Bloino
- Consiglio
Nazionale delle Ricerche, Istituto di Chimica dei Composti OrganoMetallici
(ICCOM-CNR), UOS di Pisa,
Area della Ricerca CNR, Via G. Moruzzi
1, 56124 Pisa, Italy
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
29
|
Belogolova EF, Vakul'skaya TI, Sidorkin VF. Radical anions of hypervalent silicon compounds: 1-substituted silatranes. Phys Chem Chem Phys 2015; 17:12735-46. [DOI: 10.1039/c4cp06046f] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The first representatives of the radical anions of silatranes XSi(OCH2CH2)3N were obtained and studied using the EPR and quantum chemistry methods. Depending on X, the addition of an electron to silatranes is accompanied by the retention of the dative Si–N contact, its significant weakening, or a change in the nature of the coordination center.
Collapse
Affiliation(s)
- Elena F. Belogolova
- A. E. Favorsky Irkutsk Institute of Chemistry
- Siberian Branch of the Russian Academy of Sciences
- Irkutsk 664033
- Russian Federation
| | - Tamara I. Vakul'skaya
- A. E. Favorsky Irkutsk Institute of Chemistry
- Siberian Branch of the Russian Academy of Sciences
- Irkutsk 664033
- Russian Federation
| | - Valery F. Sidorkin
- A. E. Favorsky Irkutsk Institute of Chemistry
- Siberian Branch of the Russian Academy of Sciences
- Irkutsk 664033
- Russian Federation
| |
Collapse
|
30
|
Hermosilla L, Calle P, García de la Vega JM. Modeling EPR parameters of nitrogen containing conjugated radical cations. RSC Adv 2015. [DOI: 10.1039/c5ra08758a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
DFT investigation on conjugated radical cations containing14N nucleus to obtain accurate isotropic hyperfine coupling constants.
Collapse
Affiliation(s)
- L. Hermosilla
- Departamento de Química Física Aplicada
- Facultad de Ciencias
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
| | - P. Calle
- Departamento de Química Física Aplicada
- Facultad de Ciencias
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
| | - J. M. García de la Vega
- Departamento de Química Física Aplicada
- Facultad de Ciencias
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
| |
Collapse
|
31
|
Barnes L, Schindler B, Allouche AR, Simon D, Chambert S, Oomens J, Compagnon I. Anharmonic simulations of the vibrational spectrum of sulfated compounds: application to the glycosaminoglycan fragment glucosamine 6-sulfate. Phys Chem Chem Phys 2015; 17:25705-13. [DOI: 10.1039/c5cp02079d] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Anharmonic behavior of sulfated glucosamine resolved by hybrid GVPT2 approach.
Collapse
Affiliation(s)
- Loïc Barnes
- Université de Lyon
- Lyon
- France
- Université Lyon 1
- Villeurbanne
| | | | | | - Daniel Simon
- Université de Lyon
- Lyon
- France
- Université Lyon 1
- Villeurbanne
| | | | - Jos Oomens
- Radboud University
- Institute for Molecules and Materials
- FELIX Laboratory
- 6525ED Nijmegen
- The Netherlands
| | | |
Collapse
|
32
|
Fornaro T, Carnimeo I, Biczysko M. Toward Feasible and Comprehensive Computational Protocol for Simulation of the Spectroscopic Properties of Large Molecular Systems: The Anharmonic Infrared Spectrum of Uracil in the Solid State by the Reduced Dimensionality/Hybrid VPT2 Approach. J Phys Chem A 2014; 119:5313-26. [DOI: 10.1021/jp510101y] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| | - Ivan Carnimeo
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
- Compunet, Istituto Italiano di Tecnologia (IIT), Via Morego 30, I-16163 Genova, Italy
| | - Malgorzata Biczysko
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), UOS di Pisa, Area della
Ricerca CNR, Consiglio Nazionale delle Ricerche, Via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
33
|
Benassi E, Cappelli C, Carlotti B, Barone V. An integrated computational tool to model the broadening of the absorption bands of flexible dyes in solution: cationic chromophores as test cases. Phys Chem Chem Phys 2014; 16:26963-73. [DOI: 10.1039/c4cp03419h] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
34
|
Kesharwani MK, Brauer B, Martin JML. Frequency and Zero-Point Vibrational Energy Scale Factors for Double-Hybrid Density Functionals (and Other Selected Methods): Can Anharmonic Force Fields Be Avoided? J Phys Chem A 2014; 119:1701-14. [DOI: 10.1021/jp508422u] [Citation(s) in RCA: 349] [Impact Index Per Article: 34.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Manoj K. Kesharwani
- Department
of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Brina Brauer
- Department
of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Jan M. L. Martin
- Department
of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| |
Collapse
|
35
|
Baiardi A, Bloino J, Barone V. A general time-dependent route to resonance-Raman spectroscopy including Franck-Condon, Herzberg-Teller and Duschinsky effects. J Chem Phys 2014; 141:114108. [PMID: 25240346 PMCID: PMC4608049 DOI: 10.1063/1.4895534] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a new formulation of the time-dependent theory of Resonance-Raman spectroscopy (TD-RR). Particular attention has been devoted to the generality of the framework and to the possibility of including different effects (Duschinsky mixing, Herzberg-Teller contributions). Furthermore, the effects of different harmonic models for the intermediate electronic state are also investigated. Thanks to the implementation of the TD-RR procedure within a general-purpose quantum-chemistry program, both solvation and leading anharmonicity effects have been included in an effective way. The reliability and stability of our TD-RR implementation are validated against our previously proposed and well-tested time-independent procedure. Practical applications are illustrated with some closed- and open-shell medium-size molecules (anthracene, phenoxyl radical, benzyl radical) and the simulated spectra are compared to the experimental results. More complex and larger systems, not limited to organic compounds, can be also studied, as shown for the case of Tris(bipyridine)ruthenium(II) chloride.
Collapse
Affiliation(s)
- Alberto Baiardi
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Julien Bloino
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
36
|
Barone V, Baiardi A, Bloino J. New developments of a multifrequency virtual spectrometer: stereo-electronic, dynamical, and environmental effects on chiroptical spectra. Chirality 2014; 26:588-600. [PMID: 24839096 PMCID: PMC4604663 DOI: 10.1002/chir.22325] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2013] [Revised: 01/29/2014] [Accepted: 03/03/2014] [Indexed: 11/09/2022]
Abstract
Computational spectroscopy has recently evolved from a field reserved for specialists toward a general tool allowing interpretations and analyses of experimental results. However, the current practice of providing tables of transitions for rigid geometries, possibly tuned by phenomenological broadening, is by far too naive. In order to improve this situation in the last few years we have been developing a general, robust, and user-friendly virtual spectrometer (VS) able to complement experimental studies for complex systems in condensed phases. The VS is based on flexible graphical pre- and postprocessing tools interfaced with general number-crunching software. This last tool is rooted in several electronic structure methodologies (DFT, TD-DFT, post-Hartree-Fock), powerful discrete/continuum models for describing environmental effects, and general vibrational and vibronic models. These last topics are the main focus of this work, which sketches our latest developments related to effective inclusion of anharmonic contributions, together with time-independent and/or time-dependent descriptions of vibronic transitions including Franck-Condon, Herzberg-Teller, and Duschinsky effects. Some test cases are described in some detail with the aim of showing the role of different effects in ruling vibrational (VCD) and electronic (ECD, CPL) chiral spectroscopies.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Alberto Baiardi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organometallici, UOS di Pisa, Via G. Moruzzi 1, 56124 Pisa, Italy
| |
Collapse
|
37
|
De Mitri N, Prampolini G, Monti S, Barone V. Structural, dynamic and photophysical properties of a fluorescent dye incorporated in an amorphous hydrophobic polymer bundle. Phys Chem Chem Phys 2014; 16:16573-87. [PMID: 24988373 PMCID: PMC4618303 DOI: 10.1039/c4cp01828a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The properties of a low molecular weight organic dye, namely 4-naphthyloxy-1-methoxy-2,2,6,6-tetramethylpiperidine, covalently bound to an apolar polyolefin were investigated by means of a multi-level approach, combining classical molecular dynamics simulations, based on purposely parameterized force fields, and quantum mechanical calculations based on density functional theory (DFT) and its time-dependent extension (TD-DFT). The structure and dynamics of the dye in its embedding medium were analyzed and discussed taking the entangling effect of the surrounding polymer into account, and also by comparing the results to those obtained for a different environment, i.e. toluene solution. Finally, the influence was investigated of long lived cages found in the polymeric embedding on photophysical properties, in terms of the slow and fast dye's internal dynamics, by comparing computed IR and UV spectra with their experimental counterparts.
Collapse
Affiliation(s)
- N De Mitri
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | | | | | | |
Collapse
|
38
|
Barone V, Biczysko M, Bloino J, Carta L, Pedone A. Reprint of “Environmental and dynamical effects on the optical properties of molecular systems by time-independent and time-dependent approaches: Coumarin derivatives as test cases”. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.05.004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
39
|
Fornaro T, Biczysko M, Monti S, Barone V. Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers. Phys Chem Chem Phys 2014; 16:10112-28. [PMID: 24531740 PMCID: PMC4612423 DOI: 10.1039/c3cp54724h] [Citation(s) in RCA: 68] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.
Collapse
Affiliation(s)
- Teresa Fornaro
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | | | | | | |
Collapse
|
40
|
Barone V, Biczysko M, Bloino J, Carta L, Pedone A. Environmental and dynamical effects on the optical properties of molecular systems by time-independent and time-dependent approaches: Coumarin derivatives as test cases. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.03.027] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
41
|
Carta L, Biczysko M, Bloino J, Licari D, Barone V. Environmental and complexation effects on the structures and spectroscopic signatures of organic pigments relevant to cultural heritage: the case of alizarin and alizarin-Mg(II)/Al(III) complexes. Phys Chem Chem Phys 2014; 16:2897-911. [PMID: 24424261 PMCID: PMC4627506 DOI: 10.1039/c3cp50499a] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
An integrated computational approach allowed an unbiased analysis of optical and structural properties of alizarin-based pigments, which can be directly compared with experimental results. Madder lake pigments have been modeled by Mg(II)- and Al(III)-coordinated alizarin taking into account solvation and metal-linkage effects, responsible for colour modifications. Moreover, different environmental conditions have been analyzed for free alizarin, showing in all cases semi-quantitative agreement with experimental spectroscopic data (UV-VIS). Our results point out the ability of in silico approaches to unravel the subtle interplay of stereo-electronic, dynamic, and environmental effects in tuning the physico-chemical properties of pigments relevant to cultural heritage.
Collapse
Affiliation(s)
- Luciano Carta
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | | | | | | | | |
Collapse
|
42
|
Barone V, Biczysko M, Bloino J. Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation. Phys Chem Chem Phys 2014; 16:1759-87. [PMID: 24346191 PMCID: PMC4604664 DOI: 10.1039/c3cp53413h] [Citation(s) in RCA: 295] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | | | | |
Collapse
|
43
|
De Mitri N, Monti S, Prampolini G, Barone V. Absorption and Emission Spectra of a Flexible Dye in Solution: a Computational Time-Dependent Approach. J Chem Theory Comput 2013; 9:4507-4516. [PMID: 26504457 PMCID: PMC4618302 DOI: 10.1021/ct4005799] [Citation(s) in RCA: 67] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The spectroscopic properties of the organic chromophore 4-naphthoyloxy-1-methoxy-2,2,6,6-tetramethylpiperidine (NfO-TEMPO-Me) in toluene solution are explored through an integrated computational strategy combining a classical dynamic sampling with a quantum mechanical description within the framework of the time-dependent density functional theory (TDDFT) approach. The atomistic simulations are based on an accurately parametrized force field, specifically designed to represent the conformational behavior of the molecule in its ground and bright excited states, whereas TDDFT calculations are performed through a selected combination of hybrid functionals and basis sets to obtain optical spectra closely matching the experimental findings. Solvent effects, crucial to obtain good accuracy, are taken into account through explicit molecules and polarizable continuum descriptions. Although, in the case of toluene, specific solvation is not fundamental, the detailed conformational sampling in solution has confirmed the importance of a dynamic description of the molecular geometry for a reliable description of the photophysical properties of the dye. The agreement between theoretical and experimental data is established and a robust protocol for the prediction of the optical behaviour of flexible fluorophores in solution is set.
Collapse
Affiliation(s)
- Nicola De Mitri
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Susanna Monti
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy ; Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Giacomo Prampolini
- Istituto per i Processi Chimico-Fisici (IPCF-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
44
|
Baiardi A, Bloino J, Barone V. General Time Dependent Approach to Vibronic Spectroscopy Including Franck-Condon, Herzberg-Teller, and Duschinsky Effects. J Chem Theory Comput 2013; 9:4097-115. [PMID: 26592403 DOI: 10.1021/ct400450k] [Citation(s) in RCA: 258] [Impact Index Per Article: 23.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
An effective time-dependent (TD) approach to compute vibrationally resolved optical spectra from first principles is presented for the computation of one-photon electronic spectra induced by either electric or magnetic transition dipoles or by their mutual interaction, namely absorption, emission, and circular dichroism. Particular care has been devoted to generality, modularity, and numerical stability including all the contributions that play a role at the harmonic level of approximation, namely Franck-Condon, Herzberg-Teller, and Dushinsky (i.e., mode mixing) effects. The implementation shares the same general framework of our previous time-independent (TI) model, thus allowing an effective integration between both approaches with the consequent enhancement of their respective strengths (e.g., spectrum completeness and straightforward account of temperature effects for the TD route versus band resolution and assignment for the TI route) using a single set of starting data. Implementation of both models in the same general computer program allows comprehensive studies using several levels of electronic structure description together with effective account of environmental effects by atomistic and/or continuum models of different sophistication. A few medium-size molecules (furan, phenyl radical, anthracene, dimethyloxirane, coumarin 339) have been studied in order to fully validate the approach.
Collapse
Affiliation(s)
- Alberto Baiardi
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 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
| | - Vincenzo Barone
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
45
|
Carnimeo I, Puzzarini C, Tasinato N, Stoppa P, Charmet AP, Biczysko M, Cappelli C, Barone V. Anharmonic theoretical simulations of infrared spectra of halogenated organic compounds. J Chem Phys 2013; 139:074310. [PMID: 23968095 PMCID: PMC4604659 DOI: 10.1063/1.4817401] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The recent implementation of the computation of infrared (IR) intensities beyond the double-harmonic approximation [J. Bloino and V. Barone, J. Chem. Phys. 136, 124108 (2012)] paved the route to routine calculations of infrared spectra for a wide set of molecular systems. Halogenated organic compounds represent an interesting class of molecules, from both an atmospheric and computational point of view, due to the peculiar chemical features related to the halogen atoms. In this work, we simulate the IR spectra of eight halogenated molecules (CH2F2, CHBrF2, CH2DBr, CF3Br, CH2CHF, CF2CFCl, cis-CHFCHBr, cis-CHFCHI), using two common hybrid and double-hybrid density functionals in conjunction with both double- and triple-ζ quality basis sets (SNSD and cc-pVTZ) as well as employing the coupled-cluster theory with basis sets of at least triple-ζ quality. Finally, we compare our results with available experimental spectra, with the aim of checking the accuracy and the performances of the computational approaches.
Collapse
Affiliation(s)
- Ivan Carnimeo
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, I-56126 Pisa, Italy.
| | | | | | | | | | | | | | | |
Collapse
|
46
|
Pedone A. Role of Solvent on Charge Transfer in 7-Aminocoumarin Dyes: New Hints from TD-CAM-B3LYP and State Specific PCM Calculations. J Chem Theory Comput 2013; 9:4087-96. [DOI: 10.1021/ct4004349] [Citation(s) in RCA: 119] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Alfonso Pedone
- Department of Chemical and Geological
Sciences, University of Modena and Reggio Emilia, via G. Campi 183, I-41125, Modena, Italy
| |
Collapse
|
47
|
Adhikary A, Kumar A, Palmer BJ, Todd AD, Sevilla MD. Formation of S-Cl phosphorothioate adduct radicals in dsDNA S-oligomers: hole transfer to guanine vs disulfide anion radical formation. J Am Chem Soc 2013; 135:12827-38. [PMID: 23885974 DOI: 10.1021/ja406121x] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In phosphorothioate-containing dsDNA oligomers (S-oligomers), one of the two nonbridging oxygen atoms in the phosphate moiety of the sugar-phosphate backbone is replaced by sulfur. In this work, electron spin resonance (ESR) studies of one-electron oxidation of several S-oligomers by Cl2(•-) at low temperatures are performed. Electrophilic addition of Cl2(•-) to phosphorothioate with elimination of Cl(-) leads to the formation of a two-center three-electron σ(2)σ*(1)-bonded adduct radical (-P-S-̇Cl). In AT S-oligomers with multiple phosphorothioates, i.e., d[ATATAsTsAsT]2, -P-S-̇Cl reacts with a neighboring phosphorothioate to form the σ(2)σ*(1)-bonded disulfide anion radical ([-P-S-̇S-P-](-)). With AT S-oligomers with a single phosphorothioate, i.e., d[ATTTAsAAT]2, reduced levels of conversion of -P-S-̇Cl to [-P-S-̇S-P-](-) are found. For guanine-containing S-oligomers containing one phosphorothioate, -P-S-̇Cl results in one-electron oxidation of guanine base but not of A, C, or T, thereby leading to selective hole transfer to G. The redox potential of -P-S-̇Cl is thus higher than that of G but is lower than those of A, C, and T. Spectral assignments to -P-S-̇Cl and [-P-S-̇S-P-](-) are based on reaction of Cl2(•-) with the model compound diisopropyl phosphorothioate. The results found for d[TGCGsCsGCGCA]2 suggest that [-P-S-̇S-P-](-) undergoes electron transfer to the one-electron-oxidized G, healing the base but producing a cyclic disulfide-bonded backbone with a substantial bond strength (50 kcal/mol). Formation of -P-S-̇Cl and its conversion to [-P-S-̇S-P-](-) are found to be unaffected by O2, and this is supported by the theoretically calculated electron affinities and reduction potentials of [-P-S-S-P-] and O2.
Collapse
Affiliation(s)
- Amitava Adhikary
- Department of Chemistry, Oakland University, Rochester, Michigan 48309, USA
| | | | | | | | | |
Collapse
|
48
|
Hermosilla L, Prampolini G, Calle P, García de la Vega JM, Brancato G, Barone V. Extension of the AMBER Force Field for Nitroxide Radicals and Combined QM/MM/PCM Approach to the Accurate Determination of EPR Parameters of DMPO-H in Solution. J Chem Theory Comput 2013; 9:3626-36. [PMID: 26584116 PMCID: PMC4660035 DOI: 10.1021/ct4003256] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A computational strategy that combines both time-dependent and time-independent approaches is exploited to accurately model molecular dynamics and solvent effects on the isotropic hyperfine coupling constants of the DMPO-H nitroxide. Our recent general force field for nitroxides derived from AMBER ff99SB is further extended to systems involving hydrogen atoms in β-positions with respect to NO-moiety. The resulting force-field has been employed in a series of classical molecular dynamics simulations, comparing the computed EPR parameters from selected molecular configurations to the corresponding experimental data in different solvents. The effect of vibrational averaging on the spectroscopic parameters is also taken into account, by second-order vibrational perturbation theory involving semidiagonal third energy derivatives together first and second property derivatives.
Collapse
Affiliation(s)
- Laura Hermosilla
- Departamento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid (Spain)
| | - Giacomo Prampolini
- Istituto per i Processi Chimico-Fisici, Consiglio Nazionale delle Ricerche, Area della Ricerca, via Moruzzi 1, I-56124, Pisa (Italy)
| | - Paloma Calle
- Departamento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid (Spain)
| | | | - Giuseppe Brancato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa (Italy)
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa (Italy)
| |
Collapse
|
49
|
Charmet AP, Quartarone G, Ronchin L, Tortato C, Vavasori A. Quantum Chemical Investigation on Indole: Vibrational Force Field and Theoretical Determination of Its Aqueous pKa Value. J Phys Chem A 2013; 117:6846-58. [DOI: 10.1021/jp4049692] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Andrea Pietropolli Charmet
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| | - Giuseppe Quartarone
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| | - Lucio Ronchin
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| | - Claudio Tortato
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| | - Andrea Vavasori
- Dipartimento
di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Calle
Larga Santa Marta 2137, I-30123, Venezia, Italy
| |
Collapse
|
50
|
Neshchadin D, Gescheidt G. CIDNP as a tool for rapid rearrangements: New insights into cyclobutanoid 4-center/3-electron radical cations. J PHYS ORG CHEM 2013. [DOI: 10.1002/poc.3153] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Dmytro Neshchadin
- Institute of Physical and Theoretical Chemistry; Graz University of Technology; Stremayrgasse 9 8010 Graz Austria
| | - Georg Gescheidt
- Institute of Physical and Theoretical Chemistry; Graz University of Technology; Stremayrgasse 9 8010 Graz Austria
| |
Collapse
|