1
|
Sepali C, Lafiosca P, Gómez S, Giovannini T, Cappelli C. Effective fully polarizable QM/MM approaches to compute Raman and Raman Optical Activity spectra in aqueous solution. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 305:123485. [PMID: 37827000 DOI: 10.1016/j.saa.2023.123485] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Revised: 09/22/2023] [Accepted: 10/01/2023] [Indexed: 10/14/2023]
Abstract
Raman and Raman Optical Activity (ROA) signals are amply affected by solvent effects, especially in the presence of strongly solute-solvent interactions such as Hydrogen Bonding (HB). In this work, we extend the fully atomistic polarizable Quantum Mechanics/Molecular Mechanics approach, based on the Fluctuating Charges and Fluctuating Dipoles force field to the calculation of Raman and ROA spectra. Such an approach is able to accurately describe specific HB interactions, by also accounting for anisotropic contributions due to the inclusion of fluctuating dipoles. To highlight the potentiality of the novel approach, Raman and ROA spectra of L-Serine and L-Cysteine dissolved in aqueous solution are computed and compared both with alternative theoretical approaches and experimental measurements.
Collapse
Affiliation(s)
- Chiara Sepali
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, Pisa, 56126, Italy
| | - Piero Lafiosca
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, Pisa, 56126, Italy
| | - Sara Gómez
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, Pisa, 56126, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, Pisa, 56126, Italy.
| |
Collapse
|
2
|
Nicoli L, Giovannini T, Cappelli C. Assessing the quality of QM/MM approaches to describe vacuo-to-water solvatochromic shifts. J Chem Phys 2022; 157:214101. [PMID: 36511555 DOI: 10.1063/5.0118664] [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/13/2022] Open
Abstract
The performance of different quantum mechanics/molecular mechanics embedding models to compute vacuo-to-water solvatochromic shifts is investigated. In particular, both nonpolarizable and polarizable approaches are analyzed and computed results are compared to reference experimental data. We show that none of the approaches outperform the others and that errors strongly depend on the nature of the molecular transition to be described. Thus, we prove that the best choice of embedding model highly depends on the molecular system and that the use of a specific approach as a black box can lead to significant errors and, sometimes, totally wrong predictions.
Collapse
Affiliation(s)
- Luca Nicoli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
3
|
Gómez S, Giovannini T, Cappelli C. Multiple Facets of Modeling Electronic Absorption Spectra of Systems in Solution. ACS PHYSICAL CHEMISTRY AU 2022; 3:1-16. [PMID: 36718266 PMCID: PMC9881242 DOI: 10.1021/acsphyschemau.2c00050] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Revised: 11/08/2022] [Accepted: 11/08/2022] [Indexed: 11/24/2022]
Abstract
In this Perspective, we outline the essential physicochemical aspects that need to be considered when building a reliable approach to describe absorption properties of solvated systems. In particular, we focus on how to properly model the complexity of the solvation phenomenon, arising from dynamical aspects and specific, strong solute-solvent interactions. To this end, conformational and configurational sampling techniques, such as Molecular Dynamics, have to be coupled to accurate fully atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) methodologies. By exploiting different illustrative applications, we show that an effective reproduction of experimental spectral signals can be achieved by delicately balancing exhaustive sampling, hydrogen bonding, mutual polarization, and nonelectrostatic effects.
Collapse
|
4
|
Santra G, Calinsky R, Martin JML. Benefits of Range-Separated Hybrid and Double-Hybrid Functionals for a Large and Diverse Data Set of Reaction Energies and Barrier Heights. J Phys Chem A 2022; 126:5492-5505. [PMID: 35930677 PMCID: PMC9393870 DOI: 10.1021/acs.jpca.2c03922] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
To better understand the thermochemical kinetics and
mechanism
of a specific chemical reaction, an accurate estimation of barrier
heights (forward and reverse) and reaction energies is vital. Because
of the large size of reactants and transition state structures involved
in real-life mechanistic studies (e.g., enzymatically catalyzed reactions),
density functional theory remains the workhorse for such calculations.
In this paper, we have assessed the performance of 91 density functionals
for modeling the reaction energies and barrier heights on a large
and chemically diverse data set (BH9) composed of 449 organic chemistry
reactions. We have shown that range-separated hybrid functionals perform
better than the global hybrids for BH9 barrier heights and reaction
energies. Except for the PBE-based range-separated nonempirical double
hybrids, range separation of the exchange term helps improve the performance
for barrier heights and reaction energies. The 16-parameter Berkeley
double hybrid, ωB97M(2), performs remarkably well for both properties.
However, our minimally empirical range-separated double hybrid functionals
offer marginally better accuracy than ωB97M(2) for BH9 barrier
heights and reaction energies.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Rivka Calinsky
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
5
|
Barone V, Carnimeo I, Mancini G, Pagliai M. Development, Validation, and Pilot Application of a Generalized Fluctuating Charge Model for Computational Spectroscopy in Solution. ACS OMEGA 2022; 7:13382-13394. [PMID: 35474835 PMCID: PMC9026056 DOI: 10.1021/acsomega.2c01110] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/07/2022] [Accepted: 03/23/2022] [Indexed: 06/14/2023]
Abstract
A general approach enforcing nonperiodic boundary conditions for the computation of spectroscopic properties in solution has been improved including an effective description of charge-transfer contributions and coordination number adjustment for explicit solvent molecules. Both contributions are obtained from a continuous description of intermolecular hydrogen bonds, which has been employed also for an effective clustering of molecular dynamics trajectories. Fine tuning of the model has been performed for several water clusters, and then its efficiency and reliability have been demonstrated by computing the absorption spectra of different creatinine tautomers in aqueous solution.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Ivan Carnimeo
- Scuola
Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy
| | - Giordano Mancini
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Pagliai
- Dipartimento
di Chimica “Ugo Schiff”, Università
degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| |
Collapse
|
6
|
Mancini G, Fusè M, Lipparini F, Nottoli M, Scalmani G, Barone V. Molecular Dynamics Simulations Enforcing Nonperiodic Boundary Conditions: New Developments and Application to the Solvent Shifts of Nitroxide Magnetic Parameters. J Chem Theory Comput 2022; 18:2479-2493. [PMID: 35257572 PMCID: PMC9009096 DOI: 10.1021/acs.jctc.2c00046] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Abstract
![]()
Multiscale methods
combining quantum mechanics and molecular mechanics
(QM/MM) have become the most suitable and effective strategies for
the investigation of the spectroscopic properties of medium-to-large
size chromophores in condensed phases. In this context, we are developing
a novel workflow aimed at improving the generality, reliability, and
ease of use of the available computational tools. In this paper, we
report our latest developments with specific reference to a general
protocol based on atomistic simulations, carried out under nonperiodic
boundary conditions (NPBC). In particular, we add to our in house
MD engine a new efficient treatment of mean field electrostatic contributions
to energy and forces, together with the capability of performing the
simulations either in the canonical (NVT) or in the
isothermal–isobaric (NPT) ensemble. Next,
we provide convincing evidence that the NBPC approach enhanced by
specific tweaks for rigid body propagation, allows for the simulation
of solute–solvent systems with a minimum number of degrees
of freedom and large integration time step. After its validation,
this new approach is applied to the challenging case of solvatochromic
effects on the electron paramagnetic resonance (EPR) spectrum of a
prototypical nitroxide radical. To this end, we propose and validate
also an automated protocol to extract and weight simulation snapshots,
making use of a continuous description of the strength of solute–solvent
hydrogen bridges. While further developments are being worked on,
the effectiveness of our approach, even in its present form, is demonstrated
by the accuracy of the results obtained through an unsupervised approach
characterized by a strongly reduced computational cost as compared
to that of conventional QM/MM models, without any appreciable deterioration
of accuracy.
Collapse
Affiliation(s)
- Giordano Mancini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| | - Marco Fusè
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Universitaá di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Universitaá di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Giovanni Scalmani
- Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| |
Collapse
|
7
|
|
8
|
Prioli S, Kongsted J. Modeling environmental effects in two-photon circular dichroism calculations. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02838-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
9
|
Barone V, Puzzarini C, Mancini G. Integration of theory, simulation, artificial intelligence and virtual reality: a four-pillar approach for reconciling accuracy and interpretability in computational spectroscopy. Phys Chem Chem Phys 2021; 23:17079-17096. [PMID: 34346437 DOI: 10.1039/d1cp02507d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The established pillars of computational spectroscopy are theory and computer based simulations. Recently, artificial intelligence and virtual reality are becoming the third and fourth pillars of an integrated strategy for the investigation of complex phenomena. The main goal of the present contribution is the description of some new perspectives for computational spectroscopy, in the framework of a strategy in which computational methodologies at the state of the art, high-performance computing, artificial intelligence and virtual reality tools are integrated with the aim of improving research throughput and achieving goals otherwise not possible. Some of the key tools (e.g., continuous molecular perception model and virtual multifrequency spectrometer) and theoretical developments (e.g., non-periodic boundaries, joint variational-perturbative models) are shortly sketched and their application illustrated by means of representative case studies taken from recent work by the authors. Some of the results presented are already well beyond the state of the art in the field of computational spectroscopy, thereby also providing a proof of concept for other research fields.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | | | | |
Collapse
|
10
|
Nottoli M, Cupellini L, Lipparini F, Granucci G, Mennucci B. Multiscale Models for Light-Driven Processes. Annu Rev Phys Chem 2021; 72:489-513. [PMID: 33561359 DOI: 10.1146/annurev-physchem-090419-104031] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| |
Collapse
|
11
|
Giovannini T, Egidi F, Cappelli C. Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems. Phys Chem Chem Phys 2020; 22:22864-22879. [PMID: 33043930 DOI: 10.1039/d0cp04027d] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.
Collapse
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | | | | |
Collapse
|
12
|
Nottoli M, Mennucci B, Lipparini F. Excited state Born-Oppenheimer molecular dynamics through coupling between time dependent DFT and AMOEBA. Phys Chem Chem Phys 2020; 22:19532-19541. [PMID: 32844823 DOI: 10.1039/d0cp03688a] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We present the implementation of excited state Born-Oppenheimer molecular dynamics (BOMD) using a polarizable QM/MM approach based on a time-dependent density functional theory (TDDFT) formulation and the AMOEBA force field. The implementation relies on an interface between Tinker and Gaussian software and it uses an algorithm for the calculation of QM/MM energy and forces which scales linearly with the number of MM atoms. The resulting code can perform TDDFT/AMOEBA BOMD simulations on real-life systems with standard computational resources. As a test case, the method is applied to the study of the mechanism of locally-excited to charge-transfer conversion in dimethylaminobenzonitrile in a polar solvent. Our simulations confirm that such a conversion is governed by the twisting of the dimethylamino group which is accompanied by an important reorientation of solvent molecules.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy.
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy.
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy.
| |
Collapse
|
13
|
Zhang K, Ren S, Caricato M. Multistate QM/QM Extrapolation of UV/Vis Absorption Spectra with Point Charge Embedding. J Chem Theory Comput 2020; 16:4361-4372. [DOI: 10.1021/acs.jctc.0c00339] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kaihua Zhang
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Sijin Ren
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Marco Caricato
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
14
|
Puzzarini C, Spada L, Alessandrini S, Barone V. The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:343002. [PMID: 32203942 DOI: 10.1088/1361-648x/ab8253] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/27/2019] [Accepted: 03/23/2020] [Indexed: 06/10/2023]
Abstract
In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.
Collapse
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
| | - Lorenzo Spada
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Silvia Alessandrini
- Dipartimento di Chimica 'Giacomo Ciamician', Via F. Selmi 2, I-40126 Bologna, Italy
- 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
|
15
|
Schwinn K, Ferré N, Huix-Rotllant M. Efficient Analytic Second Derivative of Electrostatic Embedding QM/MM Energy: Normal Mode Analysis of Plant Cryptochrome. J Chem Theory Comput 2020; 16:3816-3824. [PMID: 32320612 DOI: 10.1021/acs.jctc.9b01145] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Analytic second derivatives of electrostatic embedding (EE) quantum mechanics/molecular mechanics (QM/MM) energy are important for performing vibrational analysis and simulating vibrational spectra of quantum systems interacting with an environment represented as a classical electrostatic potential. The main bottleneck of EE-QM/MM second derivatives is the solution of coupled perturbed equations for each MM atom perturbation. Here, we exploit the Q-vector method [J. Chem. Phys., 2019, 151, 041102] to workaround this bottleneck. We derive the full analytic second derivative of the EE-QM/MM energy, which allows us to compute QM, MM, and QM-MM Hessian blocks in an efficient and easy to implement manner. To show the capabilities of our method, we compute the normal modes for the full Arabidopsis thaliana plant cryptochrome. We show that the flavin adenine dinucleotide vibrations (QM subsystem) strongly mix with protein modes. We compute approximate vibronic couplings for the lowest bright transition, from which we extract spectral densities and the homogeneous broadening of FAD absorption spectrum in protein using vibrationally resolved electronic spectrum simulations.
Collapse
|
16
|
Bondanza M, Nottoli M, Cupellini L, Lipparini F, Mennucci B. Polarizable embedding QM/MM: the future gold standard for complex (bio)systems? Phys Chem Chem Phys 2020; 22:14433-14448. [DOI: 10.1039/d0cp02119a] [Citation(s) in RCA: 59] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale
- Università di Pisa
- I-56124 Pisa
- Italy
| |
Collapse
|
17
|
Giovannini T, Riso RR, Ambrosetti M, Puglisi A, Cappelli C. Electronic transitions for a fully polarizable QM/MM approach based on fluctuating charges and fluctuating dipoles: Linear and corrected linear response regimes. J Chem Phys 2019; 151:174104. [PMID: 31703497 DOI: 10.1063/1.5121396] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
The fully polarizable Quantum Mechanics/Molecular Mechanics (QM/MM) approach based on fluctuating charges and fluctuating dipoles, named QM/FQFμ [T. Giovannini et al., J. Chem. Theory Comput. 15, 2233 (2019)], is extended to the calculation of vertical excitation energies of solvated molecular systems. Excitation energies are defined within two different solvation regimes, i.e., linear response (LR), where the response of the MM portion is adjusted to the QM transition density, and corrected-Linear Response (cLR) in which the MM response is adjusted to the relaxed QM density, thus being able to account for charge equilibration in the excited state. The model, which is specified in terms of three physical parameters (electronegativity, chemical hardness, and polarizability) is applied to vacuo-to-water solvatochromic shifts of aqueous solutions of para-nitroaniline, pyridine, and pyrimidine. The results show a good agreement with their experimental counterparts, thus highlighting the potentialities of this approach.
Collapse
Affiliation(s)
- Tommaso Giovannini
- 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
|
18
|
Di Remigio R, Giovannini T, Ambrosetti M, Cappelli C, Frediani L. Fully Polarizable QM/Fluctuating Charge Approach to Two-Photon Absorption of Aqueous Solutions. J Chem Theory Comput 2019; 15:4056-4068. [DOI: 10.1021/acs.jctc.9b00305] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Luca Frediani
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway
| |
Collapse
|
19
|
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
|
20
|
Giovannini T, Puglisi A, Ambrosetti M, Cappelli C. Polarizable QM/MM Approach with Fluctuating Charges and Fluctuating Dipoles: The QM/FQFμ Model. J Chem Theory Comput 2019; 15:2233-2245. [DOI: 10.1021/acs.jctc.8b01149] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
| | | | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
21
|
Dziedzic J, Head-Gordon T, Head-Gordon M, Skylaris CK. Mutually polarizable QM/MM model with in situ optimized localized basis functions. J Chem Phys 2019; 150:074103. [PMID: 30795653 DOI: 10.1063/1.5080384] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
We extend our recently developed quantum-mechanical/molecular mechanics (QM/MM) approach [Dziedzic et al., J. Chem. Phys. 145, 124106 (2016)] to enable in situ optimization of the localized orbitals. The quantum subsystem is described with onetep linear-scaling density functional theory and the classical subsystem - with the AMOEBA polarizable force field. The two subsystems interact via multipolar electrostatics and are fully mutually polarizable. A total energy minimization scheme is employed for the Hamiltonian of the coupled QM/MM system. We demonstrate that, compared to simpler models using fixed basis sets, the additional flexibility offered by in situ optimized basis functions improves the accuracy of the QM/MM interface, but also poses new challenges, making the QM subsystem more prone to overpolarization and unphysical charge transfer due to increased charge penetration. We show how these issues can be efficiently solved by replacing the classical repulsive van der Waals term for QM/MM interactions with an interaction of the electronic density with a fixed, repulsive MM potential that mimics Pauli repulsion, together with a modest increase in the damping of QM/MM polarization. We validate our method, with particular attention paid to the hydrogen bond, in tests on water-ion pairs, the water dimer, first solvation shells of neutral and charged species, and solute-solvent interaction energies. As a proof of principle, we determine suitable repulsive potential parameters for water, K+, and Cl-. The mechanisms we employed to counteract the unphysical overpolarization of the QM subsystem are demonstrated to be adequate, and our approach is robust. We find that the inclusion of explicit polarization in the MM part of QM/MM improves agreement with fully QM calculations. Our model permits the use of minimal size QM regions and, remarkably, yields good energetics across the well-balanced QM/MM interface.
Collapse
Affiliation(s)
- Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Teresa Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
22
|
Egidi F, Giovannini T, Del Frate G, Lemler PM, Vaccaro PH, Cappelli C. A combined experimental and theoretical study of optical rotatory dispersion for (R)-glycidyl methyl ether in aqueous solution. Phys Chem Chem Phys 2019; 21:3644-3655. [PMID: 30383044 DOI: 10.1039/c8cp04445g] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The dispersive optical activity for aqueous solutions of non-rigid (R)-glycidyl methyl ether (R-GME) has been explored synergistically from experimental and theoretical perspectives. Density functional theory analyses performed with the polarizable continuum model for implicit solvation identified nine low-lying stable conformers that are interconverted by rotation about two large-amplitude torsional coordinates. The antagonistic chiroptical signatures predicted for these structural isomers were averaged under a Boltzmann-weighting ansatz to estimate the behavior expected for a thermally equilibrated ensemble. This led to optical rotatory dispersion profiles that reproduced the overall shape of observations but failed to achieve uniform agreement with measured specific-rotation values even when anharmonic vibrational corrections were applied. A mixed QM/FQ paradigm, whereby quantum-mechanical (QM) calculations of optical activity were combined with classical molecular dynamics simulations of explicit solvation that included mutual-polarization effects by means of fluctuating charges (FQ), was enlisted to elucidate the microsolvation environment and gauge its impact upon conformer distributions and response properties. Although quantitative accord with experiments remained elusive, this approach revealed strong variations in the magnitude and sign of rotatory powers for R-GME as the configuration of surrounding water molecules evolved, thereby highlighting the inherently dynamical nature of the solvated chiroptical response, calling into question the validity of "static" descriptions based on the presumption of distinct energy minima, and giving insight into the inherent complexity posed by the modeling of such properties for solvated systems.
Collapse
Affiliation(s)
- Franco Egidi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, Italy
| | | | | | | | | | | |
Collapse
|
23
|
Puglisi A, Giovannini T, Antonov L, Cappelli C. Interplay between conformational and solvent effects in UV-visible absorption spectra: curcumin tautomers as a case study. Phys Chem Chem Phys 2019; 21:15504-15514. [PMID: 31259324 DOI: 10.1039/c9cp00907h] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
We present a combined theoretical and experimental study on the UV-vis spectra of enol-keto (EK) and keto-keto (KK) tautomeric forms of curcumin dissolved in aqueous solution. Solvent effects have been investigated by resorting to the implicit polarizable continuum model (QM/PCM) and non-polarizable and fully polarizable QM/MM approaches, the latter based on the fluctuating charges (FQ) force-field. In particular, all methods are challenged to rationalize the contribution of conformational, electrostatic and polarization effects in the calculation of the vertical excitation spectra of curcumin tautomers. The obtained results highlight that for both tautomers specific solute-solvent hydrogen-bond interactions play a minor role with respect to conformational and electrostatic effects.
Collapse
Affiliation(s)
| | | | - Liudmil Antonov
- Bulgarian Academy of Sciences, Institute of Organic Chemistry with Centre of Phytochemistry, Acad. G. Bonchev str., Bldg. 9, Sofia 1113, Bulgaria
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| |
Collapse
|
24
|
Reinholdt P, Nørby MS, Kongsted J. Modeling of Magnetic Circular Dichroism and UV/Vis Absorption Spectra Using Fluctuating Charges or Polarizable Embedding within a Resonant-Convergent Response Theory Formalism. J Chem Theory Comput 2018; 14:6391-6404. [DOI: 10.1021/acs.jctc.8b00660] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Morten S. Nørby
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| |
Collapse
|
25
|
On the nature of charge-transfer excitations for molecules in aqueous solution: a polarizable QM/MM study. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2259-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
26
|
A polarizable embedding approach to second harmonic generation (SHG) of molecular systems in aqueous solutions. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2247-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
27
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
28
|
Macchiagodena M, Mancini G, Pagliai M, Barone V. Accurate prediction of bulk properties in hydrogen bonded liquids: amides as case studies. Phys Chem Chem Phys 2018; 18:25342-25354. [PMID: 27711662 DOI: 10.1039/c6cp04666e] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
In this contribution we show that it is possible to build accurate force fields for small organic molecules allowing the reliable reproduction of a large panel of bulk properties, which are seldom addressed in the same context. Starting from the results obtained in recent studies, we developed a protocol for charge estimation and virtual site generation for the amide class of molecules. The parametrization of electrostatic properties is based on population analysis and orbital localization of quantum mechanical computations rooted in density functional theory and the polarizable continuum model, without any additional external information. The new protocol, coupled to other recent studies in our group targeted at an accurate fitting of internal degrees of freedom, makes available a method for building force fields from scratch (excluding for the moment intermolecular van der Waals interactions) with focus on reproducing the structure and dynamics of hydrogen bonded liquids, yielding results that are in line or better than those delivered by current general force fields. The approach is tested on the demanding series formed by formamide and its two N-methyl derivatives, N-methylformamide and N,N-dimethylformamide. We show that the atomistic structure of the liquids arising from classical molecular dynamics (MD) simulations employing the new force field is in full agreement with X-ray and neutron diffraction experiments and the corresponding spatial distribution functions are in remarkable agreement with the results of ab initio MD simulations. It is noteworthy that the latter result has never been obtained before without using ad hoc (and system dependent) scale factors and that, in addition, our parameter-free procedure is able to reproduce static dielectric constants over a wide range of values without sacrificing the force field accuracy with respect to other observables. Finally, we are able to explain the trend of static dielectric constants followed by the three amides in terms of properties obtained from the simulations, namely hydrogen bond patterns and reorientational lifetimes.
Collapse
Affiliation(s)
| | - Giordano Mancini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126, Pisa, Italy. and Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| | - Marco Pagliai
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126, Pisa, Italy.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126, Pisa, Italy. and Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| |
Collapse
|
29
|
Licari D, Fusè M, Salvadori A, Tasinato N, Mendolicchio M, Mancini G, Barone V. Towards the SMART workflow system for computational spectroscopy. Phys Chem Chem Phys 2018; 20:26034-26052. [DOI: 10.1039/c8cp03417f] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications?
Collapse
Affiliation(s)
- Daniele Licari
- Scuola Normale Superiore
- 56126 Pisa
- Italy
- Istituto Italiano di Tecnologia
- 16163 Genova
| | | | | | | | | | | | | |
Collapse
|
30
|
Giovannini T, Del Frate G, Lafiosca P, Cappelli C. Effective computational route towards vibrational optical activity spectra of chiral molecules in aqueous solution. Phys Chem Chem Phys 2018; 20:9181-9197. [DOI: 10.1039/c8cp00487k] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
A polarizable QM/MM approach to accurately compute the Vibrational Optical Activity (VOA) spectra of chiral systems is proposed and applied to aqueous solutions of (l)-methyl lactate and (S)-glycidol.
Collapse
|
31
|
Giovannini T, Lafiosca P, Cappelli C. A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches. J Chem Theory Comput 2017; 13:4854-4870. [PMID: 28898079 DOI: 10.1021/acs.jctc.7b00776] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
A methodology to account for nonelectrostatic interactions in Quantum Mechanical (QM)/Molecular Mechanics (MM) approaches is developed. Formulations for Pauli repulsion and dispersion energy, explicitly depending on the QM density, are derived. Such expressions are based on the definition of an auxiliary density on the MM portion and the Tkatchenko-Scheffler (TS) approach, respectively. The developed method is general enough to be applied to any QM/MM method and partition, provided an accurate tuning of a small number of parameters is obtained. The coupling of the method with both nonpolarizable and fully polarizable QM/fluctuating charge (FQ) approaches is reported and applied. A suitable parametrization for the aqueous solution, so that its most representative features are well reproduced, is outlined. Then, the obtained parametrization and method are applied to calculate the nonelectrostatic (repulsion and dispersion) interaction energy of nicotine in aqueous solution.
Collapse
Affiliation(s)
| | - Piero Lafiosca
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
32
|
Menger MFSJ, Caprasecca S, Mennucci B. Excited-State Gradients in Polarizable QM/MM Models: An Induced Dipole Formulation. J Chem Theory Comput 2017; 13:3778-3786. [DOI: 10.1021/acs.jctc.7b00475] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Maximilian F. S. J. Menger
- Dipartimento
di Chimica e Chimica Industriale, University of Pisa, Via G. Moruzzi
13, 56124 Pisa, Italy
- Institut
für Theoretische Chemie, Universität Wien, Währinger
Strasse 17, A-1090 Wien, Austria
| | - Stefano Caprasecca
- Dipartimento
di Chimica e Chimica Industriale, University of Pisa, Via G. Moruzzi
13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento
di Chimica e Chimica Industriale, University of Pisa, Via G. Moruzzi
13, 56124 Pisa, Italy
| |
Collapse
|
33
|
Pagliai M, Mancini G, Carnimeo I, De Mitri N, Barone V. Electronic absorption spectra of pyridine and nicotine in aqueous solution with a combined molecular dynamics and polarizable QM/MM approach. J Comput Chem 2016; 38:319-335. [PMID: 27910109 PMCID: PMC6680224 DOI: 10.1002/jcc.24683] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2016] [Revised: 11/04/2016] [Accepted: 11/08/2016] [Indexed: 01/06/2023]
Abstract
The electronic absorption spectra of pyridine and nicotine in aqueous solution have been computed using a multistep approach. The computational protocol consists in studying the solute solvation with accurate molecular dynamics simulations, characterizing the hydrogen bond interactions, and calculating electronic transitions for a series of configurations extracted from the molecular dynamics trajectories with a polarizable QM/MM scheme based on the fluctuating charge model. Molecular dynamics simulations and electronic transition calculations have been performed on both pyridine and nicotine. Furthermore, the contributions of solute vibrational effect on electronic absorption spectra have been taken into account in the so called vertical gradient approximation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Marco Pagliai
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, I-56126, Italy
| | - Giordano Mancini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, I-56126, Italy
| | - Ivan Carnimeo
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, I-56126, Italy.,Compunet, Istituto Italiano di Tecnologia (IIT), via Morego 30, Genova, I-16163, Italy
| | - Nicola De Mitri
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, I-56126, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, I-56126, Italy
| |
Collapse
|
34
|
Oruganti B, Fang C, Durbeej B. Assessment of a composite CC2/DFT procedure for calculating 0–0 excitation energies of organic molecules. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1235736] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Baswanth Oruganti
- Division of Theoretical Chemistry, IFM, Linköping University, Linköping, Sweden
| | - Changfeng Fang
- Division of Theoretical Chemistry, IFM, Linköping University, Linköping, Sweden
| | - Bo Durbeej
- Division of Theoretical Chemistry, IFM, Linköping University, Linköping, Sweden
| |
Collapse
|
35
|
Dziedzic J, Mao Y, Shao Y, Ponder J, Head-Gordon T, Head-Gordon M, Skylaris CK. TINKTEP: A fully self-consistent, mutually polarizable QM/MM approach based on the AMOEBA force field. J Chem Phys 2016; 145:124106. [DOI: 10.1063/1.4962909] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
- Faculty of Applied Physics and Mathematics, Gdańsk University of Technology, Gdańsk, Poland
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Jay Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
- Department of Bioengineering, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
36
|
Barone V. The Virtual Multifrequency Spectrometer: a new paradigm for spectroscopy. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2016; 6:86-110. [PMID: 29075335 PMCID: PMC5654514 DOI: 10.1002/wcms.1238] [Citation(s) in RCA: 67] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
On going developments of hardware and software are changing computational spectroscopy from a strongly specialized research area to a general tool in the inventory of most researchers. Increased interactions between experimentally-oriented users and theoretically-oriented developers of new methods and models would result in more robust, flexible and reliable tools and studies for the systems of increasing complexity, which are of current scientific and technological interest. This is the philosophy behind this review, which presents the development of a so-called virtual multi-frequency spectrometer (VMS) including state-of-the-art approaches in a user-friendly frame. The current status of the VMS tool will be illustrated by a number of case studies with special reference to infrared and UV-vis regions of the electro-magnetic spectrum including also chiral spectroscopies. Only the basic theoretical background will be provided avoiding explicit equations as far as possible, and pointing out the most recent advancements beyond the standard rigid-rotor harmonic-oscillator model coupled to vertical electronic excitation energies.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
37
|
Carnimeo I, Cappelli C, Barone V. Analytical gradients for MP2, double hybrid functionals, and TD-DFT with polarizable embedding described by fluctuating charges. J Comput Chem 2015; 36:2271-90. [PMID: 26399473 PMCID: PMC5054946 DOI: 10.1002/jcc.24195] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2015] [Revised: 07/31/2015] [Accepted: 08/07/2015] [Indexed: 01/09/2023]
Abstract
A polarizable quantum mechanics (QM)/ molecular mechanics (MM) approach recently developed for Hartree-Fock (HF) and Kohn-Sham (KS) methods has been extended to energies and analytical gradients for MP2, double hybrid functionals, and TD-DFT models, thus allowing the computation of equilibrium structures for excited electronic states together with more accurate results for ground electronic states. After a detailed presentation of the theoretical background and of some implementation details, a number of test cases are analyzed to show that the polarizable embedding model based on fluctuating charges (FQ) is remarkably more accurate than the corresponding electronic embedding based on a fixed charge (FX) description. In particular, a set of electronegativities and hardnesses has been optimized for interactions between QM and FQ regions together with new repulsion-dispersion parameters. After validation of both the numerical implementation and of the new parameters, absorption electronic spectra have been computed for representative model systems including vibronic effects. The results show remarkable agreement with full QM computations and significant improvement with respect to the corresponding FX results. The last part of the article provides some hints about computation of solvatochromic effects on absorption spectra in aqueous solution as a function of the number of FQ water molecules and on the use of FX external shells to improve the convergence of the results. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Ivan Carnimeo
- Istituto Italiano Di Tecnologia (IIT), Compunet, via Morego, 30 I-16163 Genova, Italy and Scuola Normale Superiore, Piazza Dei Cavalieri, 7 I-56126, Pisa, Italy
| | - Chiara Cappelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi, 13, I-56124, Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri, 7 I-56126, Pisa, Italy
| |
Collapse
|