1
|
Cerezo J, Gierschner J, Santoro F, Prampolini G. Explicit Modelling of Spectral Bandshapes by a Mixed Quantum-Classical Approach: Solvent Order and Temperature Effects in the Optical Spectra of Distyrylbenzene. Chemphyschem 2024:e202400307. [PMID: 38728539 DOI: 10.1002/cphc.202400307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2024] [Revised: 04/29/2024] [Accepted: 05/08/2024] [Indexed: 05/12/2024]
Abstract
The absorption and emission spectral shapes of a flexible organic probe, the distyrylbenzene (DSB) dye, are simulated accounting for the effect of different environments of increasing complexity, ranging from a homogeneous, low-molecular- weight solvent, to a long-chain alkane, and, eventually, a channel-forming organic matrix. Each embedding is treated explicitly, adopting a mixed quantum-classical approach, the Adiabatic Molecular Dynamics - generalized vertical Hessian (Ad-MD|gVH) model, which allows a direct simulation of the environment-induced constraining effects on the vibronic spectral shapes. In such a theoretical framework, the stiff modes of the dye are described at a quantum level within the harmonic approximation, including Duschinsky mixing effects, while flexible degrees of freedom of the solute (e. g. torsions) and those of the solvent are treated classically by means of molecular dynamics sampling. Such a setup is shown to reproduce the distinct effects exerted by the different environments in varied thermodynamic conditions. Besides allowing for a first-principles rationale on the supramolecular mechanism leading to the experimental spectral features, this result represents the first successful application of the Ad-MD|gVH method to complex embeddings and supports its potential application to other heterogeneous environments, such as for instance, pigment-protein complexes or organic dyes adsorbed into metal-organic frameworks.
Collapse
Affiliation(s)
- Javier Cerezo
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049, Madrid, Spain
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM), Consiglio Nazionale delle Ricerche (CNR), 50019, Sesto Fiorentino, Italien
| | - Johannes Gierschner
- Madrid Institute for Advanced Studies, IMDEA Nanoscience, C/Faraday 9, Ciudad Universitaria de Cantoblanco, 28049, Madrid, Spain
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM), Consiglio Nazionale delle Ricerche (CNR), 50019, Sesto Fiorentino, Italien
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM), Consiglio Nazionale delle Ricerche (CNR), 50019, Sesto Fiorentino, Italien
| |
Collapse
|
2
|
Dellai A, Naim C, Cerezo J, Prampolini G, Castet F. Dynamic effects on the nonlinear optical properties of donor acceptor stenhouse adducts: insights from combined MD + QM simulations. Phys Chem Chem Phys 2024; 26:13639-13654. [PMID: 38511505 DOI: 10.1039/d4cp00310a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/22/2024]
Abstract
The second-order nonlinear optical (NLO) responses of a donor-acceptor stenhouse adduct (DASA) are investigated by using a computational approach combining molecular dynamics simulations and density functional theory (DFT) calculations. Specific force fields for the open and closed photoswitching forms are first parameterized and validated according to the Joyce protocol, in order to finely reproduce the geometrical features and potential energy surfaces of both isomers in chloroform solution. Then, DFT calculations are performed on structural snapshots extracted at regular time steps of the MD trajectories to address the influence of the thermalized conformational dynamics on the NLO responses related to hyper-Rayleigh scattering (HRS) experiments. We show that accounting for the structural dynamics largely enhances the HRS hyperpolarizability (βHRS) compared to DFT calculations considering solely equilibrium geometries, and greatly improves the agreement with experimental measurements. Furthermore, we show that the NLO responses of the NLO-active open form are correlated with the bond order alternation along the triene bridge connecting the donor and acceptor moieties, which is rationalized using simple essential state models.
Collapse
Affiliation(s)
- Angela Dellai
- Univ. Bordeaux, CNRS, Bordeaux INP, Institut des Sciences Moléculaires, UMR 5255, F-33400 Talence, France.
| | - Carmelo Naim
- Univ. Bordeaux, CNRS, Bordeaux INP, Institut des Sciences Moléculaires, UMR 5255, F-33400 Talence, France.
- Donostia International Physics Center (DIPC), Manuel Lardizabal Ibilbidea 4, 20018 Donostia, Euskadi, Spain
| | - Javier Cerezo
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | | | - Frédéric Castet
- Univ. Bordeaux, CNRS, Bordeaux INP, Institut des Sciences Moléculaires, UMR 5255, F-33400 Talence, France.
| |
Collapse
|
3
|
Oliveira MP, Hünenberger PH. Influence of the Lennard-Jones Combination Rules on the Simulated Properties of Organic Liquids at Optimal Force-Field Parametrization. J Chem Theory Comput 2023; 19:2048-2063. [PMID: 36920838 PMCID: PMC10100539 DOI: 10.1021/acs.jctc.2c01170] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
We recently introduced the CombiFF scheme [Oliveira et al., J. Chem. Theory Comput. 2020, 16, 7525], an approach for the automated refinement of force-field parameters against experimental condensed-phase data for large compound families. Using this scheme, once the time-consuming task of target-data selection and curation has been performed, the force-field optimization itself is both straightforward and fast. As a result, CombiFF provides an ideal framework for evaluating the influence of functional-form decisions on the accuracy of a force field at an optimal level of parametrization. We already used this approach to assess the effect of using an all-atom representation compared to united-atom representations in the force field [Oliveira et al., J. Chem. Theory Comput. 2022, 18, 6757]. Here, CombiFF is applied to assess the effect of three Lennard-Jones combination rules, geometric-mean (GM), Lorentz-Berthelot (LB), or Waldman-Hagler (WH), on the simulated properties of organic liquids. The comparison is performed in terms of the experimental liquid density ρliq, vaporization enthalpy ΔHvap, surface-tension coefficient γ, static relative dielectric permittivity ϵ, and self-diffusion coefficient D. The calibrations of the three force-field variants are carried out independently against 2044 experimental values for ρliq, and ΔHvap concerning 1516 compounds. The resulting root-mean-square deviations from experiment are 30.0, 26.9, and 36.7 kg m-3 for ρliq and 2.8, 2.8, and 2.9 kJ mol-1 for ΔHvap, when applying the GM, LB, and WH combination rules, respectively. In terms of these (and the other) properties, the three combination rules perform comparatively well, with the GM and LB results being more similar to each other and slightly more accurate compared to experiment. In contrast, the use of distinct combination rules for the parameter calibration and property calculation leads to much larger errors.
Collapse
Affiliation(s)
- Marina P Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, CH-8093 Zürich, Switzerland
| | | |
Collapse
|
4
|
Janicki TD, Van Vleet MJ, Schmidt JR. Development and Implementation of Atomically Anisotropic First-Principles Force Fields: A Benzene Case Study. J Phys Chem A 2023; 127:1736-1749. [PMID: 36780209 DOI: 10.1021/acs.jpca.2c07244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
Abstract
π-interactions are an important motif in chemical and biochemical systems. However, due to their anisotropic electron densities and complex balance of intermolecular interactions, aromatic molecules represent an ongoing challenge for accurate and transferable force field development. Historically, ab initio force fields for aromatics have not exhibited good accuracy with respect to bulk properties or have only been used to study gas-phase dimers. Using benzene as a proof of concept, herein we show how our own ab initio MASTIFF force field incorporates an atomically anisotropic description of intermolecular interactions to yield an accurate and robust model for aromatic interactions irrespective of phase. Compared to existing models, the MASTIFF benzene force field not only is accurate for liquid phase properties but also offers transferability to the gas and solid phases. Additionally, we introduce a computationally efficient OpenMM plugin which enables customizable anisotropic intermolecular functional forms and which can be generically used in any MD simulation where a model for nonspherical atomic features is required. Overall, our results demonstrate the importance of atomic-level anisotropy in enabling next-generation ab initio force field development.
Collapse
Affiliation(s)
- Tesia D Janicki
- Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Mary J Van Vleet
- Department of Chemistry and Biochemistry, Spelman College, 350 Spelman Ln SW, Atlanta, Georgia 30314, United States
| | - J R Schmidt
- Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
5
|
Csizi K, Reiher M. Universal
QM
/
MM
approaches for general nanoscale applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| |
Collapse
|
6
|
Greff da Silveira L, Livotto PR, Padula D, Vilhena JG, Prampolini G. Accurate Quantum-Mechanically Derived Force-Fields through a Fragment-Based Approach: Balancing Specificity and Transferability in the Prediction of Self-Assembly in Soft Matter. J Chem Theory Comput 2022; 18:6905-6919. [PMID: 36260420 DOI: 10.1021/acs.jctc.2c00747] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The wide range of time/length scales covered by self-assembly in soft matter makes molecular dynamics (MD) the ideal candidate for simulating such a supramolecular phenomenon at an atomistic level. However, the reliability of MD outcomes heavily relies on the accuracy of the adopted force-field (FF). The spontaneous re-ordering in liquid crystalline materials stands as a clear example of such collective self-assembling processes, driven by a subtle and delicate balance between supramolecular interactions and single-molecule flexibility. General-purpose transferable FFs often dramatically fail to reproduce such complex phenomena, for example, the error on the transition temperatures being larger than 100 K. Conversely, quantum-mechanically derived force-fields (QMD-FFs), specifically tailored for the target system, were recently shown (J. Phys. Chem. Lett.2022,13, 243) to allow for the required accuracy as they not only well reproduced transition temperatures but also yielded a quantitative agreement with the experiment on a wealth of structural, dynamic, and thermodynamic properties. The main drawback of this strategy stands in the computational burden connected to the numerous quantum mechanical (QM) calculations usually required for a target-specific parameterization, which has undoubtedly hampered the routine application of QMD-FFs. In this work, we propose a fragment-based strategy to extend the applicability of QMD-FFs, in which the amount of QM calculations is significantly reduced, being a single-molecule-optimized geometry and its Hessian matrix the only QM information required. To validate this route, a new FF is assembled for a large mesogen, exploiting the parameters obtained for two smaller liquid crystalline molecules, in this and previous work. Lengthy MD simulations are carried out with the new transferred QMD-FF, observing again a spontaneous re-orientation in the correct range of temperatures, with good agreement with the available experimental measures. The present results strongly suggest that a partial transfer of QMD-FF parameters can be invoked without a significant loss of accuracy, thus paving the way to exploit the method's intrinsic predictive capabilities in the simulation of novel soft materials.
Collapse
Affiliation(s)
- Leandro Greff da Silveira
- Instituto de Química (Universidade Federal do Rio Grande do Sul), Avenida Bento Gonçalves 9500, CEP 91501-970Porto Alegre, Brazil
| | - Paolo Roberto Livotto
- Instituto de Química (Universidade Federal do Rio Grande do Sul), Avenida Bento Gonçalves 9500, CEP 91501-970Porto Alegre, Brazil
| | - Daniele Padula
- Dipartimento di Biotecnologie, Chimica e Farmacia (Università di Siena), via Aldo Moro 2, 53100Siena, SI, Italy
| | - J G Vilhena
- Departamento de Física Teórica de la Materia Condensada (Universidad Autónoma de Madrid), E-28049Madrid, Spain.,Condensed Matter Physics Center (IFIMAC) (Universidad Autónoma de Madrid), E-28049Madrid, Spain
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124Pisa, Italy
| |
Collapse
|
7
|
Segalina A, Aranda D, Green JA, Cristino V, Caramori S, Prampolini G, Pastore M, Santoro F. How the Interplay among Conformational Disorder, Solvation, Local, and Charge-Transfer Excitations Affects the Absorption Spectrum and Photoinduced Dynamics of Perylene Diimide Dimers: A Molecular Dynamics/Quantum Vibronic Approach. J Chem Theory Comput 2022; 18:3718-3736. [PMID: 35377648 PMCID: PMC9202308 DOI: 10.1021/acs.jctc.2c00063] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
![]()
In this contribution
we present a mixed quantum-classical dynamical
approach for the computation of vibronic absorption spectra of molecular
aggregates and their nonadiabatic dynamics, taking into account the
coupling between local excitations (LE) and charge-transfer (CT) states.
The approach is based on an adiabatic (Ad) separation between the
soft degrees of freedom (DoFs) of the system and the stiff vibrations,
which are described by the quantum dynamics (QD) of wave packets (WPs)
moving on the coupled potential energy surfaces (PESs) of the LE and
CT states. These PESs are described with a linear vibronic coupling
(LVC) Hamiltonian, parameterized by an overlap-based diabatization
on the grounds of time-dependent density functional theory computations.
The WPs time evolution is computed with the multiconfiguration time-dependent
Hartree method, using effective modes defined through a hierarchical
representation of the LVC Hamiltonian. The soft DoFs are sampled with
classical molecular dynamics (MD), and the coupling between the slow
and fast DoFs is included by recomputing the key parameters of the
LVC Hamiltonians, specifically for each MD configuration. This method,
named Ad-MD|gLVC, is applied to a perylene diimide (PDI) dimer in
acetonitrile and water solutions, and it is shown to accurately reproduce
the change in the vibronic features of the absorption spectrum upon
aggregation. Moreover, the microscopic insight offered by the MD trajectories
allows for a detailed understanding of the role played by the fluctuation
of the aggregate structure on the shape of the vibronic spectrum and
on the population of LE and CT states. The nonadiabatic QD predicts
an extremely fast (∼50 fs) energy transfer between the two
LEs. CT states have only a moderate effect on the absorption spectrum,
despite the fact that after photoexcitation they are shown to acquire
a fast and non-negligible population, highlighting their relevance
in dictating the charge separation and transport in PDI-based optical
devices.
Collapse
Affiliation(s)
- Alekos Segalina
- Université de Lorraine and CNRS, LPCT, UMR 7019, F-54000 Nancy, France
| | - Daniel Aranda
- Instituto de Ciencia Molecular (ICMol), Universidad de Valencia, Catedrático J. Beltrán 2, 46980 Paterna, Valencia, Spain
| | - James A Green
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - Vito Cristino
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie, Via Fossato di Mortara 17, 44121 Ferrara, Italy
| | - Stefano Caramori
- Dipartimento di Scienze Chimiche, Farmaceutiche ed Agrarie, Via Fossato di Mortara 17, 44121 Ferrara, Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | | | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
8
|
Prampolini G, Greff da Silveira L, Vilhena JG, Livotto PR. Predicting Spontaneous Orientational Self-Assembly: In Silico Design of Materials with Quantum Mechanically Derived Force Fields. J Phys Chem Lett 2022; 13:243-250. [PMID: 34968058 DOI: 10.1021/acs.jpclett.1c03517] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
De novo design of self-assembled materials hinges upon our ability to relate macroscopic properties to individual building blocks, thus characterizing in such supramolecular architectures a wide range of observables at varied time/length scales. This work demonstrates that quantum mechanical derived force fields (QMD-FFs) do satisfy this requisite and, most importantly, do so in a predictive manner. To this end, a specific FF, built solely based on the knowledge of the target molecular structure, is employed to reproduce the spontaneous transition to an ordered liquid crystal phase. The simulations deliver a multiscale portrait of such self-assembly processes, where conformational changes within the individual building blocks are intertwined with a 200 ns ensemble reorganization. The extensive characterization provided not only is in quantitative agreement with the experiment but also connects the time/length scales at which it was performed. Realizing QMD-FF predictive power and unmatched accuracy stands as an important leap forward for the bottom-up design of advanced supramolecular materials.
Collapse
Affiliation(s)
- Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Leandro Greff da Silveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91 501-970 Porto Alegre, Brazil
| | - J G Vilhena
- Department of Physics, University of Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| | - Paolo Roberto Livotto
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91 501-970 Porto Alegre, Brazil
| |
Collapse
|
9
|
Yang L, Horton JT, Payne MC, Penfold TJ, Cole DJ. Modeling Molecular Emitters in Organic Light-Emitting Diodes with the Quantum Mechanical Bespoke Force Field. J Chem Theory Comput 2021; 17:5021-5033. [PMID: 34264669 DOI: 10.1021/acs.jctc.1c00135] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Combined molecular dynamics (MD) and quantum mechanics (QM) simulation procedures have gained popularity in modeling the spectral properties of functional organic molecules. However, the potential energy surfaces used to propagate long-time scale dynamics in these simulations are typically described using general, transferable force fields designed for organic molecules in their electronic ground states. These force fields do not typically include spectroscopic data in their training, and importantly, there is no general protocol for including changes in geometry or intermolecular interactions with the environment that may occur upon electronic excitation. In this work, we show that parameters tailored for thermally activated delayed fluorescence (TADF) emitters used in organic light-emitting diodes (OLEDs), in both their ground and electronically excited states, can be readily derived from a small number of QM calculations using the QUBEKit (QUantum mechanical BEspoke toolKit) software and improve the overall accuracy of these simulations.
Collapse
Affiliation(s)
- Lupeng Yang
- TCM Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Joshua T Horton
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Michael C Payne
- TCM Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Thomas J Penfold
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Daniel J Cole
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| |
Collapse
|
10
|
Sami S, Menger MFSJ, Faraji S, Broer R, Havenith RWA. Q-Force: Quantum Mechanically Augmented Molecular Force Fields. J Chem Theory Comput 2021; 17:4946-4960. [PMID: 34251194 PMCID: PMC8359013 DOI: 10.1021/acs.jctc.1c00195] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The quality of molecular
dynamics simulations strongly depends
on the accuracy of the underlying force fields (FFs) that determine
all intra- and intermolecular interactions of the system. Commonly,
transferable FF parameters are determined based on a representative
set of small molecules. However, such an approach sacrifices accuracy
in favor of generality. In this work, an open-source and automated
toolkit named Q-Force is presented, which augments these transferable
FFs with molecule-specific bonded parameters and atomic charges that
are derived from quantum mechanical (QM) calculations. The molecular
fragmentation procedure allows treatment of large molecules (>200
atoms) with a low computational cost. The generated Q-Force FFs can
be used at the same computational cost as transferable FFs, but with
improved accuracy: We demonstrate this for the vibrational properties
on a set of small molecules and for the potential energy surface on
a complex molecule (186 atoms) with photovoltaic applications. Overall,
the accuracy, user-friendliness, and minimal computational overhead
of the Q-Force protocol make it widely applicable for atomistic molecular
dynamics simulations.
Collapse
Affiliation(s)
- Selim Sami
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Maximilian F S J Menger
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Shirin Faraji
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Ria Broer
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Remco W A Havenith
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281-(S3), B-9000 Ghent, Belgium
| |
Collapse
|
11
|
Vilhena JG, Greff da Silveira L, Livotto PR, Cacelli I, Prampolini G. Automated Parameterization of Quantum Mechanically Derived Force Fields for Soft Materials and Complex Fluids: Development and Validation. J Chem Theory Comput 2021; 17:4449-4464. [PMID: 34185536 DOI: 10.1021/acs.jctc.1c00213] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The reliability of molecular dynamics (MD) simulations in predicting macroscopic properties of complex fluids and soft materials, such as liquid crystals, colloidal suspensions, or polymers, relies on the accuracy of the adopted force field (FF). We present an automated protocol to derive specific and accurate FFs, fully based on ab initio quantum mechanical (QM) data. The integration of the Joyce and Picky procedures, recently proposed by our group to provide an accurate description of simple liquids, is here extended to larger molecules, capable of exhibiting more complex fluid phases. While the standard Joyce protocol is employed to parameterize the intramolecular FF term, a new automated procedure is here proposed to handle the computational cost of the QM calculations required for the parameterization of the intermolecular FF term. The latter is thus obtained by integrating the old Picky procedure with a fragmentation reconstruction method (FRM) that allows for a reliable, yet computationally feasible sampling of the intermolecular energy surface at the QM level. The whole FF parameterization protocol is tested on a benchmark liquid crystal, and the performances of the resulting quantum mechanically derived (QMD) FF were compared with those delivered by a general-purpose, transferable one, and by the third, "hybrid" FF, where only the bonded terms were refined against QM data. Lengthy atomistic MD simulations are carried out with each FF on extended 5CB systems in both isotropic and nematic phases, eventually validating the proposed protocol by comparing the resulting macroscopic properties with other computational models and with experiments. The QMD-FF yields the best performances, reproducing both phases in the correct range of temperatures and well describing their structure, dynamics, and thermodynamic properties, thus providing a clear protocol that may be explored to predict such properties on other complex fluids or soft materials.
Collapse
Affiliation(s)
- J G Vilhena
- Department of Physics, University of Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| | - Leandro Greff da Silveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Paolo Roberto Livotto
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici, ICCOM-CNR, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
12
|
Schriber JB, Sirianni DA, Smith DGA, Burns LA, Sitkoff D, Cheney DL, Sherrill CD. Optimized damping parameters for empirical dispersion corrections to symmetry-adapted perturbation theory. J Chem Phys 2021; 154:234107. [PMID: 34241276 DOI: 10.1063/5.0049745] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Symmetry-adapted perturbation theory (SAPT) has become an invaluable tool for studying the fundamental nature of non-covalent interactions by directly computing the electrostatics, exchange (steric) repulsion, induction (polarization), and London dispersion contributions to the interaction energy using quantum mechanics. Further application of SAPT is primarily limited by its computational expense, where even its most affordable variant (SAPT0) scales as the fifth power of system size [O(N5)] due to the dispersion terms. The algorithmic scaling of SAPT0 is reduced from O(N5)→O(N4) by replacing these terms with the empirical D3 dispersion correction of Grimme and co-workers, forming a method that may be termed SAPT0-D3. Here, we optimize the damping parameters for the -D3 terms in SAPT0-D3 using a much larger training set than has previously been considered, namely, 8299 interaction energies computed at the complete-basis-set limit of coupled cluster through perturbative triples [CCSD(T)/CBS]. Perhaps surprisingly, with only three fitted parameters, SAPT0-D3 improves on the accuracy of SAPT0, reducing mean absolute errors from 0.61 to 0.49 kcal mol-1 over the full set of complexes. Additionally, SAPT0-D3 exhibits a nearly 2.5× speedup over conventional SAPT0 for systems with ∼300 atoms and is applied here to systems with up to 459 atoms. Finally, we have also implemented a functional group partitioning of the approach (F-SAPT0-D3) and applied it to determine important contacts in the binding of salbutamol to G-protein coupled β1-adrenergic receptor in both active and inactive forms. SAPT0-D3 capabilities have been added to the open-source Psi4 software.
Collapse
Affiliation(s)
- Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Dominic A Sirianni
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Daniel G A Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Doree Sitkoff
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
13
|
Schriber JB, Nascimento DR, Koutsoukas A, Spronk SA, Cheney DL, Sherrill CD. CLIFF: A component-based, machine-learned, intermolecular force field. J Chem Phys 2021; 154:184110. [PMID: 34241025 DOI: 10.1063/5.0042989] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Computation of intermolecular interactions is a challenge in drug discovery because accurate ab initio techniques are too computationally expensive to be routinely applied to drug-protein models. Classical force fields are more computationally feasible, and force fields designed to match symmetry adapted perturbation theory (SAPT) interaction energies can remain accurate in this context. Unfortunately, the application of such force fields is complicated by the laborious parameterization required for computations on new molecules. Here, we introduce the component-based machine-learned intermolecular force field (CLIFF), which combines accurate, physics-based equations for intermolecular interaction energies with machine-learning models to enable automatic parameterization. The CLIFF uses functional forms corresponding to electrostatic, exchange-repulsion, induction/polarization, and London dispersion components in SAPT. Molecule-independent parameters are fit with respect to SAPT2+(3)δMP2/aug-cc-pVTZ, and molecule-dependent atomic parameters (atomic widths, atomic multipoles, and Hirshfeld ratios) are obtained from machine learning models developed for C, N, O, H, S, F, Cl, and Br. The CLIFF achieves mean absolute errors (MAEs) no worse than 0.70 kcal mol-1 in both total and component energies across a diverse dimer test set. For the side chain-side chain interaction database derived from protein fragments, the CLIFF produces total interaction energies with an MAE of 0.27 kcal mol-1 with respect to reference data, outperforming similar and even more expensive methods. In applications to a set of model drug-protein interactions, the CLIFF is able to accurately rank-order ligand binding strengths and achieves less than 10% error with respect to SAPT reference values for most complexes.
Collapse
Affiliation(s)
- Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Daniel R Nascimento
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Steven A Spronk
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| |
Collapse
|
14
|
Fortino M, Cozza C, Bonomi M, Pietropaolo A. Multi-replica biased sampling for photoswitchable π-conjugated polymers. J Chem Phys 2021; 154:174108. [PMID: 34241080 DOI: 10.1063/5.0045944] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In recent years, π-conjugated polymers are attracting considerable interest in view of their light-dependent torsional reorganization around the π-conjugated backbone, which determines peculiar light-emitting properties. Motivated by the interest in designing conjugated polymers with tunable photoswitchable pathways, we devised a computational framework to enhance the sampling of the torsional conformational space and, at the same time, estimate ground- to excited-state free-energy differences. This scheme is based on a combination of Hamiltonian Replica Exchange Method (REM), parallel bias metadynamics, and free-energy perturbation theory. In our scheme, each REM samples an intermediate unphysical state between the ground and the first two excited states, which are characterized by time-dependent density functional theory simulations at the B3LYP/6-31G* level of theory. We applied the method to a 5-mer of 9,9-dioctylfluorene and found that upon irradiation, this system can undergo a dihedral inversion from -155° to 155°, crossing a barrier that decreases from 0.1 eV in the ground state (S0) to 0.05 eV and 0.04 eV in the first (S1) and second (S2) excited states. Furthermore, S1 and even more S2 were predicted to stabilize coplanar dihedrals, with a local free-energy minimum located at ±44°. The presence of a free-energy barrier of 0.08 eV for the S1 state and 0.12 eV for the S2 state can trap this conformation in a basin far from the global free-energy minimum located at 155°. The simulation results were compared with the experimental emission spectrum, showing a quantitative agreement with the predictions provided by our framework.
Collapse
Affiliation(s)
- Mariagrazia Fortino
- Dipartimento di Scienze Della Salute, Università di Catanzaro, Viale Europa, 88100 Catanzaro, Italy
| | - Concetta Cozza
- Dipartimento di Scienze Della Salute, Università di Catanzaro, Viale Europa, 88100 Catanzaro, Italy
| | - Massimiliano Bonomi
- Structural Bioinformatics Unit, Department of Structural Biology and Chemistry; CNRS UMR 3528; C3BI, CNRS USR 3756; Institut Pasteur, Paris, France
| | - Adriana Pietropaolo
- Dipartimento di Scienze Della Salute, Università di Catanzaro, Viale Europa, 88100 Catanzaro, Italy
| |
Collapse
|
15
|
Loeffler JR, Fernández-Quintero ML, Waibl F, Quoika PK, Hofer F, Schauperl M, Liedl KR. Conformational Shifts of Stacked Heteroaromatics: Vacuum vs. Water Studied by Machine Learning. Front Chem 2021; 9:641610. [PMID: 33842433 PMCID: PMC8032969 DOI: 10.3389/fchem.2021.641610] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Accepted: 03/08/2021] [Indexed: 11/13/2022] Open
Abstract
Stacking interactions play a crucial role in drug design, as we can find aromatic cores or scaffolds in almost any available small molecule drug. To predict optimal binding geometries and enhance stacking interactions, usually high-level quantum mechanical calculations are performed. These calculations have two major drawbacks: they are very time consuming, and solvation can only be considered using implicit solvation. Therefore, most calculations are performed in vacuum. However, recent studies have revealed a direct correlation between the desolvation penalty, vacuum stacking interactions and binding affinity, making predictions even more difficult. To overcome the drawbacks of quantum mechanical calculations, in this study we use neural networks to perform fast geometry optimizations and molecular dynamics simulations of heteroaromatics stacked with toluene in vacuum and in explicit solvation. We show that the resulting energies in vacuum are in good agreement with high-level quantum mechanical calculations. Furthermore, we show that using explicit solvation substantially influences the favored orientations of heteroaromatic rings thereby emphasizing the necessity to include solvation properties starting from the earliest phases of drug design.
Collapse
Affiliation(s)
- Johannes R Loeffler
- Center of Molecular Biosciences Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Monica L Fernández-Quintero
- Center of Molecular Biosciences Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Franz Waibl
- Center of Molecular Biosciences Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Patrick K Quoika
- Center of Molecular Biosciences Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Florian Hofer
- Center of Molecular Biosciences Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Michael Schauperl
- Center of Molecular Biosciences Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Klaus R Liedl
- Center of Molecular Biosciences Innsbruck, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
16
|
Li W, Quesada-Moreno MM, Pinacho P, Schnell M. Unlocking the Water Trimer Loop: Isotopic Study of Benzophenone-(H 2 O) 1-3 Clusters with Rotational Spectroscopy. Angew Chem Int Ed Engl 2021; 60:5323-5330. [PMID: 33289239 PMCID: PMC7986920 DOI: 10.1002/anie.202013899] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2020] [Indexed: 12/16/2022]
Abstract
Examined here are the structures of complexes of benzophenone microsolvated with up to three water molecules by using broadband rotational spectroscopy and the cold conditions of a molecular jet. The analysis shows that the water molecules dock sideways on benzophenone for the water monomer and dimer moieties, and they move above one of the aromatic rings when the water cluster grows to the trimer. The rotational spectra shows that the water trimer moiety in the complex adopts an open‐loop arrangement. Ab initio calculations face a dilemma of identifying the global minimum between the open loop and the closed loop, which is only solved when zero‐point vibrational energy correction is applied. An OH⋅⋅⋅π bond and a Bürgi‐Dunitz interaction between benzophenone and the water trimer are present in the cluster. This work shows the subtle balance between water–water and water–solute interactions when the solute molecule offers several different anchor sites for water molecules.
Collapse
Affiliation(s)
- Weixing Li
- Deutsches Elektronen-Synchrotron, Notkestrasse 85, 22607, Hamburg, Germany
| | | | - Pablo Pinacho
- Deutsches Elektronen-Synchrotron, Notkestrasse 85, 22607, Hamburg, Germany
| | - Melanie Schnell
- Deutsches Elektronen-Synchrotron, Notkestrasse 85, 22607, Hamburg, Germany.,Christian-Albrechts-Universität zu Kiel, Institute of Physical Chemistry, Max-Eyth-Str. 1, 24118, Kiel, Germany
| |
Collapse
|
17
|
Li W, Quesada‐Moreno MM, Pinacho P, Schnell M. Unlocking the Water Trimer Loop: Isotopic Study of Benzophenone‐(H
2
O)
1–3
Clusters with Rotational Spectroscopy. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202013899] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Weixing Li
- Deutsches Elektronen-Synchrotron Notkestrasse 85 22607 Hamburg Germany
| | | | - Pablo Pinacho
- Deutsches Elektronen-Synchrotron Notkestrasse 85 22607 Hamburg Germany
| | - Melanie Schnell
- Deutsches Elektronen-Synchrotron Notkestrasse 85 22607 Hamburg Germany
- Christian-Albrechts-Universität zu Kiel Institute of Physical Chemistry Max-Eyth-Str. 1 24118 Kiel Germany
| |
Collapse
|
18
|
Oliveira MP, Andrey M, Rieder SR, Kern L, Hahn DF, Riniker S, Horta BAC, Hünenberger PH. Systematic Optimization of a Fragment-Based Force Field against Experimental Pure-Liquid Properties Considering Large Compound Families: Application to Saturated Haloalkanes. J Chem Theory Comput 2020; 16:7525-7555. [DOI: 10.1021/acs.jctc.0c00683] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Maurice Andrey
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R. Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Leyla Kern
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - David F. Hahn
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
19
|
Segalina A, Cerezo J, Prampolini G, Santoro F, Pastore M. Accounting for Vibronic Features through a Mixed Quantum-Classical Scheme: Structure, Dynamics, and Absorption Spectra of a Perylene Diimide Dye in Solution. J Chem Theory Comput 2020; 16:7061-7077. [DOI: 10.1021/acs.jctc.0c00919] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Alekos Segalina
- Université de Lorraine & CNRS, LPCT, UMR 7019, F-54000 Nancy, France
| | - Javier Cerezo
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR) SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR) SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | | |
Collapse
|
20
|
Campetella M, De Mitri N, Prampolini G. Automated parameterization of quantum-mechanically derived force-fields including explicit sigma holes: A pathway to energetic and structural features of halogen bonds in gas and condensed phase. J Chem Phys 2020; 153:044106. [PMID: 32752684 DOI: 10.1063/5.0014280] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
In classical molecular dynamics, general purpose atomistic force-fields (FFs) often deliver inaccurate results when dealing with halogen bonds (XBs), notwithstanding their crucial role in many fields of science, ranging from material design to drug development. Given the large dimensions of the systems of interest, it would be therefore desirable to increase the FF accuracy maintaining the simplicity of the standard Lennard-Jones (LJ) plus point charge description to avoid an excessive computational cost. A simple yet effective strategy consists in introducing a number of virtual sites able to mimic the so-called "explicit σ-hole." In this work, we present an automated FF parameterization strategy based on a global optimization of both LJ and charge parameters with respect to accurate quantum mechanical data, purposely computed for the system under investigation. As a test case, we report on two homologue series, characterized either by weak or strong XBs, namely, the di-halogenated methanes and the mono-, di-, and tri-substituted acetonitriles, taking into consideration Cl, Br, and I substituents. The resulting quantum mechanically derived FFs are validated for each compound in the gas and in the condensed phase by comparing them to general purpose and specific FFs without virtual sites and to highly accurate reference quantum mechanical data. The results strongly support the adoption of the specific FFs with virtual sites, which overcome the other investigated models in representing both gas phase energetics and the structural patterns of the liquid phase structure related to the presence of XBs.
Collapse
Affiliation(s)
- Marco Campetella
- Institut des Nanosciences de Paris, Sorbonne Université, CNRS, UMR7588, F-75252 Paris, France
| | - Nicola De Mitri
- Enthought Ltd., Broers Building, 21 JJ Thomson Avenue, Cambridge CB3 0FA, United Kingdom
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti Organo Metallici (ICCOM), CNR Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
21
|
Diez-Cabanes V, Prampolini G, Francés-Monerris A, Monari A, Pastore M. Iron's Wake: The Performance of Quantum Mechanical-Derived Versus General-Purpose Force Fields Tested on a Luminescent Iron Complex. Molecules 2020; 25:molecules25133084. [PMID: 32640764 PMCID: PMC7411876 DOI: 10.3390/molecules25133084] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Revised: 06/30/2020] [Accepted: 07/01/2020] [Indexed: 12/29/2022] Open
Abstract
Recently synthetized iron complexes have achieved long-lived excited states and stabilities which are comparable, or even superior, to their ruthenium analogues, thus representing an eco-friendly and cheaper alternative to those materials based on rare metals. Most of computational tools which could help unravel the origin of this large efficiency rely on ab-initio methods which are not able, however, to capture the nanosecond time scale underlying these photophysical processes and the influence of their realistic environment. Therefore, it exists an urgent need of developing new low-cost, but still accurate enough, computational methodologies capable to deal with the steady-state and transient spectroscopy of transition metal complexes in solution. Following this idea, here we focus on the comparison between general-purpose transferable force-fields (FFs), directly available from existing databases, and specific quantum mechanical derived FFs (QMD-FFs), obtained in this work through the Joyce procedure. We have chosen a recently reported FeIII complex with nanosecond excited-state lifetime as a representative case. Our molecular dynamics (MD) simulations demonstrated that the QMD-FF nicely reproduces the structure and the dynamics of the complex and its chemical environment within the same precision as higher cost QM methods, whereas general-purpose FFs failed in this purpose. Although in this particular case the chemical environment plays a minor role on the photo physics of this system, these results highlight the potential of QMD-FFs to rationalize photophysical phenomena provided an accurate QM method to derive its parameters is chosen.
Collapse
Affiliation(s)
- Valentin Diez-Cabanes
- Université de Lorraine & CNRS, LPCT UMR 7019, F-54000 Nancy, France;
- Correspondence: (V.D.-C.); (G.P.); (A.M.); (M.P.)
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
- Correspondence: (V.D.-C.); (G.P.); (A.M.); (M.P.)
| | - Antonio Francés-Monerris
- Université de Lorraine & CNRS, LPCT UMR 7019, F-54000 Nancy, France;
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Antonio Monari
- Université de Lorraine & CNRS, LPCT UMR 7019, F-54000 Nancy, France;
- Correspondence: (V.D.-C.); (G.P.); (A.M.); (M.P.)
| | - Mariachiara Pastore
- Université de Lorraine & CNRS, LPCT UMR 7019, F-54000 Nancy, France;
- Correspondence: (V.D.-C.); (G.P.); (A.M.); (M.P.)
| |
Collapse
|
22
|
Ferretti A, d’Ischia M, Prampolini G. Benchmarking Cation−π Interactions: Assessment of Density Functional Theory and Möller–Plesset Second-Order Perturbation Theory Calculations with Optimized Basis Sets (mp2mod) for Complexes of Benzene, Phenol, and Catechol with Na+, K+, Rb+, and Cs+. J Phys Chem A 2020; 124:3445-3459. [DOI: 10.1021/acs.jpca.0c02090] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Alessandro Ferretti
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Marco d’Ischia
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, I-80126 Napoli, Italy
| | | |
Collapse
|
23
|
Bao Y, Luo Z, Cui S. Environment-dependent single-chain mechanics of synthetic polymers and biomacromolecules by atomic force microscopy-based single-molecule force spectroscopy and the implications for advanced polymer materials. Chem Soc Rev 2020; 49:2799-2827. [PMID: 32236171 DOI: 10.1039/c9cs00855a] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
"The Tao begets the One. One begets all things of the world." This quote from Tao Te Ching is still inspiring for scientists in chemistry and materials science: The "One" can refer to a single molecule. A macroscopic material is composed of numerous molecules. Although the relationship between the properties of the single molecule and macroscopic material is not well understood yet, it is expected that a deeper understanding of the single-chain mechanics of macromolecules will certainly facilitate the development of materials science. Atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) has been exploited extensively as a powerful tool to study the single-chain behaviors of macromolecules. In this review, we summarize the recent advances in the emerging field of environment-dependent single-chain mechanics of synthetic polymers and biomacromolecules by means of AFM-SMFS. First, the single-chain inherent elasticities of several typical linear macromolecules are introduced, which are also confirmed by one of three polymer models with theoretical elasticities of the corresponding macromolecules obtained from quantum mechanical (QM) calculations. Then, the effects of the external environments on the single-chain mechanics of synthetic polymers and biomacromolecules are reviewed. Finally, the impacts of single-chain mechanics of macromolecules on the development of polymer science especially polymer materials are illustrated.
Collapse
Affiliation(s)
- Yu Bao
- Key Laboratory of Advanced Technologies of Materials, Ministry of Education, Southwest Jiaotong University, Chengdu 610031, China.
| | | | | |
Collapse
|
24
|
Development of force fields for binary systems: Application to a dimethylsulfoxide (DMSO) – Oxygen mixture. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.136778] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
25
|
Allen AA, Robertson MJ, Payne MC, Cole DJ. Development and Validation of the Quantum Mechanical Bespoke Protein Force Field. ACS OMEGA 2019; 4:14537-14550. [PMID: 31528808 PMCID: PMC6740169 DOI: 10.1021/acsomega.9b01769] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Accepted: 07/30/2019] [Indexed: 06/10/2023]
Abstract
Molecular mechanics force field parameters for macromolecules, such as proteins, are traditionally fit to reproduce experimental properties of small molecules, and thus, they neglect system-specific polarization. In this paper, we introduce a complete protein force field that is designed to be compatible with the quantum mechanical bespoke (QUBE) force field by deriving nonbonded parameters directly from the electron density of the specific protein under study. The main backbone and sidechain protein torsional parameters are rederived in this work by fitting to quantum mechanical dihedral scans for compatibility with QUBE nonbonded parameters. Software is provided for the preparation of QUBE input files. The accuracy of the new force field, and the derived torsional parameters, is tested by comparing the conformational preferences of a range of peptides and proteins with experimental measurements. Accurate backbone and sidechain conformations are obtained in molecular dynamics simulations of dipeptides, with NMR J coupling errors comparable to the widely used OPLS force field. In simulations of five folded proteins, the secondary structure is generally retained, and the NMR J coupling errors are similar to standard transferable force fields, although some loss of the experimental structure is observed in certain regions of the proteins. With several avenues for further development, the use of system-specific nonbonded force field parameters is a promising approach for next-generation simulations of biological molecules.
Collapse
Affiliation(s)
- Alice
E. A. Allen
- TCM
Group, Cavendish Laboratory, 19 JJ Thomson Ave, Cambridge CB3 0HE, United Kingdom
| | - Michael J. Robertson
- Department of Molecular and Cellular Physiology and Department of Structural Biology Stanford University School of Medicine, 279 Campus Drive, Stanford, California 94305, United States
| | - Michael C. Payne
- TCM
Group, Cavendish Laboratory, 19 JJ Thomson Ave, Cambridge CB3 0HE, United Kingdom
| | - Daniel J. Cole
- School
of Natural and Environmental Sciences, Newcastle
University, Newcastle
upon Tyne NE1 7RU, United
Kingdom
| |
Collapse
|
26
|
Density functional theory study of π-aromatic interaction of benzene, phenol, catechol, dopamine isolated dimers and adsorbed on graphene surface. J Mol Model 2019; 25:302. [PMID: 31486895 DOI: 10.1007/s00894-019-4185-2] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2019] [Accepted: 08/27/2019] [Indexed: 10/26/2022]
Abstract
We analyze the influence of different groups on the intermolecular energy of aromatic homodimers and on the interaction between a single aromatic molecule and a graphene surface. The analysis is performed for benzene, phenol, catechol, and dopamine. For calculating the energies, we employ density functional theory within the local density approximation (LDA-DFT). Our results show that the lowest intermolecular energies between the aromatic molecules are related to the T-shaped configurations. This lower energy results from the quadrupole interaction. In the case of the interaction between the graphene sheet and the aromatic molecules, the lowest energy configuration is the face to face. The adsorption energy of a molecule on a graphene surface involves π - π interactions that explain the face to face arrangement. These results provide insight into the manner by which substituents can be utilized in crystal engineering, supramolecular chemistry, bioinspired materials, formation of various molecular clusters, parameterization of force fields suitable for classical simulations, and design of novel sensing, drug delivery, and filters based on graphene.
Collapse
|
27
|
Cacelli I, Lipparini F, Greff da Silveira L, Jacobs M, Livotto PR, Prampolini G. Accurate interaction energies by spin component scaled Möller-Plesset second order perturbation theory calculations with optimized basis sets (SCS-MP2mod): Development and application to aromatic heterocycles. J Chem Phys 2019; 150:234113. [DOI: 10.1063/1.5094288] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Leandro Greff da Silveira
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Matheus Jacobs
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, 12489 Berlin, Germany
- IRIS Adelrshof, Humboldt-Universität zu Berlin, Zum Großen Windkanal 6, 12489 Berlin, Germany
| | - Paolo Roberto Livotto
- Instituto de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
28
|
Calabrese C, Li W, Prampolini G, Evangelisti L, Uriarte I, Cacelli I, Melandri S, Cocinero EJ. A General Treatment to Study Molecular Complexes Stabilized by Hydrogen‐, Halogen‐, and Carbon‐Bond Networks: Experiment and Theory of (CH
2
F
2
)
n
⋅⋅⋅(H
2
O)
m. Angew Chem Int Ed Engl 2019; 58:8437-8442. [DOI: 10.1002/anie.201902753] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2019] [Revised: 04/05/2019] [Indexed: 11/08/2022]
Affiliation(s)
- Camilla Calabrese
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| | - Weixing Li
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR)Area della Ricerca via G. Moruzzi 1 I-56124 Pisa Italy
| | - Luca Evangelisti
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
- Department of ChemistryUniversity of Virginia McCormick Road VA 22903 Charlottesville USA
| | - Iciar Uriarte
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| | - Ivo Cacelli
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR)Area della Ricerca via G. Moruzzi 1 I-56124 Pisa Italy
- Dipartimento di Chimica e Chimica IndustrialeUniversità di Pisa via Risorgimento 35 I-56126 Pisa Italy
| | - Sonia Melandri
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
| | - Emilio J. Cocinero
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| |
Collapse
|
29
|
Cai W, Xu D, Qian L, Wei J, Xiao C, Qian L, Lu ZY, Cui S. Force-Induced Transition of π–π Stacking in a Single Polystyrene Chain. J Am Chem Soc 2019; 141:9500-9503. [DOI: 10.1021/jacs.9b03490] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Wanhao Cai
- Key Laboratory
of Advanced Technologies of Materials (Ministry of Education), Southwest Jiaotong University, Chengdu 610031, China
| | - Duo Xu
- State Key Laboratory
of Supramolecular Structure and Materials, Institute of Theoretical
Chemistry, Jilin University, Changchun 130023, China
| | - Lu Qian
- Key Laboratory
of Advanced Technologies of Materials (Ministry of Education), Southwest Jiaotong University, Chengdu 610031, China
| | - Junhao Wei
- Key Laboratory
of Advanced Technologies of Materials (Ministry of Education), Southwest Jiaotong University, Chengdu 610031, China
| | - Chen Xiao
- Key Laboratory
of Advanced Technologies of Materials (Ministry of Education), Southwest Jiaotong University, Chengdu 610031, China
| | - Linmao Qian
- Key Laboratory
of Advanced Technologies of Materials (Ministry of Education), Southwest Jiaotong University, Chengdu 610031, China
| | - Zhong-yuan Lu
- State Key Laboratory
of Supramolecular Structure and Materials, Institute of Theoretical
Chemistry, Jilin University, Changchun 130023, China
| | - Shuxun Cui
- Key Laboratory
of Advanced Technologies of Materials (Ministry of Education), Southwest Jiaotong University, Chengdu 610031, China
| |
Collapse
|
30
|
Calabrese C, Li W, Prampolini G, Evangelisti L, Uriarte I, Cacelli I, Melandri S, Cocinero EJ. A General Treatment to Study Molecular Complexes Stabilized by Hydrogen‐, Halogen‐, and Carbon‐Bond Networks: Experiment and Theory of (CH
2
F
2
)
n
⋅⋅⋅(H
2
O)
m
. Angew Chem Int Ed Engl 2019. [DOI: 10.1002/ange.201902753] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Camilla Calabrese
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| | - Weixing Li
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR)Area della Ricerca via G. Moruzzi 1 I-56124 Pisa Italy
| | - Luca Evangelisti
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
- Department of ChemistryUniversity of Virginia McCormick Road VA 22903 Charlottesville USA
| | - Iciar Uriarte
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| | - Ivo Cacelli
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR)Area della Ricerca via G. Moruzzi 1 I-56124 Pisa Italy
- Dipartimento di Chimica e Chimica IndustrialeUniversità di Pisa via Risorgimento 35 I-56126 Pisa Italy
| | - Sonia Melandri
- Dipartimento di Chimica “Giacomo Ciamician”Università degli Studi di Bologna via Selmi 2 I-40126 Bologna Italy
| | - Emilio J. Cocinero
- Departament of Physical ChemistryUniversity of the Basque Country (UPV/EHU) Barrio Sarriena, S/N 48940 Leioa Spain
- Biofisika Institute, (CSIC, UPV/EHU) 48080 Bilbao Spain
| |
Collapse
|
31
|
Horton JT, Allen AEA, Dodda LS, Cole DJ. QUBEKit: Automating the Derivation of Force Field Parameters from Quantum Mechanics. J Chem Inf Model 2019; 59:1366-1381. [DOI: 10.1021/acs.jcim.8b00767] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Joshua T. Horton
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| | - Alice E. A. Allen
- TCM Group, Cavendish Laboratory, 19 JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Leela S. Dodda
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Daniel J. Cole
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU, United Kingdom
| |
Collapse
|
32
|
Mason TO, Buell AK. The Kinetics, Thermodynamics and Mechanisms of Short Aromatic Peptide Self-Assembly. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2019; 1174:61-112. [PMID: 31713197 DOI: 10.1007/978-981-13-9791-2_3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The self-assembly of short aromatic peptides and peptide derivatives into a variety of different nano- and microstructures (fibrillar gels, crystals, spheres, plates) is a promising route toward the creation of bio-compatible materials with often unexpected and useful properties. Furthermore, such simple self-assembling systems have been proposed as model systems for the self-assembly of longer peptides, a process that can be linked to biological function and malfunction. Much effort has been made in the last 15 years to explore the space of peptide sequences, chemical modifications and solvent conditions in order to maximise the diversity of assembly morphologies and properties. However, quantitative studies of the corresponding mechanisms of, and driving forces for, peptide self-assembly have remained relatively scarce until recently. In this chapter we review the current state of understanding of the thermodynamic driving forces and self-assembly mechanisms of short aromatic peptides into supramolecular structures. We will focus on experimental studies of the assembly process and our perspective will be centered around diphenylalanine (FF), a key motif of the amyloid β sequence and a paradigmatic self-assembly building block. Our main focus is the basic physical chemistry and key structural aspects of such systems, and we will also compare the mechanism of dipeptide aggregation with that of longer peptide sequences into amyloid fibrils, with discussion on how these mechanisms may be revealed through detailed analysis of growth kinetics, thermodynamics and other fundamental properties of the aggregation process.
Collapse
Affiliation(s)
- Thomas O Mason
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovot, Israel
| | - Alexander K Buell
- Department of Biotechnology and Biomedicine, Technical University of Denmark, DTU, Lyngby, Denmark.
| |
Collapse
|
33
|
Prampolini G, Ingrosso F, Segalina A, Caramori S, Foggi P, Pastore M. Dynamical and Environmental Effects on the Optical Properties of an Heteroleptic Ru(II)–Polypyridine Complex: A Multilevel Approach Combining Accurate Ground and Excited State QM-Derived Force Fields, MD and TD-DFT. J Chem Theory Comput 2018; 15:529-545. [DOI: 10.1021/acs.jctc.8b01031] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Francesca Ingrosso
- Université de Lorraine, CNRS, Laboratoire de Physique et Chimie Théoriques, F-54000 Nancy, France
| | - Alekos Segalina
- Université de Lorraine, CNRS, Laboratoire de Physique et Chimie Théoriques, F-54000 Nancy, France
| | - Stefano Caramori
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Ferrara, Via Luigi Borsari 46, I-44100, Ferrara, Italy
| | - Paolo Foggi
- European Laboratory for Non Linear Spectroscopy (LENS), Università di Firenze, Via Nello Carrara 1, I-50019 Sesto Fiorentino Florence, Italy
- INO−CNR, Istituto Nazionale di Ottica, Consiglio Nazionale delle Ricerche, Largo Fermi 6, I-50125 Florence, Italy
- Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, Via Elce di Sotto 8, I-06123 Perugia, Italy
| | - Mariachiara Pastore
- Université de Lorraine, CNRS, Laboratoire de Physique et Chimie Théoriques, F-54000 Nancy, France
| |
Collapse
|
34
|
Greff da Silveira L, Jacobs M, Prampolini G, Livotto PR, Cacelli I. Development and Validation of Quantum Mechanically Derived Force-Fields: Thermodynamic, Structural, and Vibrational Properties of Aromatic Heterocycles. J Chem Theory Comput 2018; 14:4884-4900. [PMID: 30040902 DOI: 10.1021/acs.jctc.8b00218] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
A selection of several aromatic molecules, representative of the important class of heterocyclic compounds, has been considered for testing and validating an automated Force Field (FF) parametrization protocol, based only on Quantum Mechanical data. The parametrization is carried out separately for the intra- and intermolecular contributions, employing respectively the Joyce and Picky software packages, previously implemented and refined in our research group. The whole approach is here automated and integrated with a computationally effective yet accurate method, devised very recently ( J. Chem. THEORY Comput., 2018, 14, 543-556) to evaluate a large number of dimer interaction energies. The resulting quantum mechanically derived FFs are then used in extensive molecular dynamics simulations, in order to evaluate a number of thermodynamic, structural, and dynamic properties of the heterocycle's gas and liquid phases. The comparison with the available experimental data is good and furnishes a validation of the presented approach, which can be confidently exploited for the design of novel and more complex materials.
Collapse
Affiliation(s)
- Leandro Greff da Silveira
- Instituto de Química , Universidade Federal do Rio Grande do Sul , Avenida Bento Gonçalves 9500 , CEP 91501-970 Porto , Alegre , Brazil
| | - Matheus Jacobs
- Instituto de Química , Universidade Federal do Rio Grande do Sul , Avenida Bento Gonçalves 9500 , CEP 91501-970 Porto , Alegre , Brazil.,Institut für Physik , Humboldt-Universität zu Berlin , Newtonstrasse 15 , 12489 , Berlin , Germany.,IRIS Adelrshof , Humboldt-Universität zu Berlin , Zum Großen Windkanal 6 , 12489 , Berlin , Germany
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR) , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy
| | - Paolo Roberto Livotto
- Instituto de Química , Universidade Federal do Rio Grande do Sul , Avenida Bento Gonçalves 9500 , CEP 91501-970 Porto , Alegre , Brazil
| | - Ivo Cacelli
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR) , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy.,Dipartimento di Chimica e Chimica Industriale , Università di Pisa , Via G. Moruzzi 13 , I-56124 Pisa , Italy
| |
Collapse
|
35
|
Cerezo J, Prampolini G, Cacelli I. Developing accurate intramolecular force fields for conjugated systems through explicit coupling terms. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2254-8] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
36
|
Jacobs M, Greff Da Silveira L, Prampolini G, Livotto PR, Cacelli I. Interaction Energy Landscapes of Aromatic Heterocycles through a Reliable yet Affordable Computational Approach. J Chem Theory Comput 2018; 14:543-556. [DOI: 10.1021/acs.jctc.7b00602] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Matheus Jacobs
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Leandro Greff Da Silveira
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
- Departamento
de Ciências Exatas e da Terra, Universidade Regional Integrada do Alto Uruguay e da Missões (URI), Avenida Assis Brasil 709, CEP 98400-00 Frederico Westphalen, Brazil
| | - Giacomo Prampolini
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Paolo Roberto Livotto
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Ivo Cacelli
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi
3, I-56124 Pisa, Italy
| |
Collapse
|
37
|
Grimme S, Bannwarth C, Caldeweyher E, Pisarek J, Hansen A. A general intermolecular force field based on tight-binding quantum chemical calculations. J Chem Phys 2017; 147:161708. [DOI: 10.1063/1.4991798] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Christoph Bannwarth
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Eike Caldeweyher
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Jana Pisarek
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| |
Collapse
|
38
|
Assessment of DFT functionals for calculating intermolecular interaction of nitrogen-containing heterocyclic complexes. Theor Chem Acc 2017. [DOI: 10.1007/s00214-017-2154-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
39
|
|
40
|
Jungwirth J, Šebestík J, Šafařík M, Kapitán J, Bouř P. Quantitative Determination of Ala-Ala Conformer Ratios in Solution by Decomposition of Raman Optical Activity Spectra. J Phys Chem B 2017; 121:8956-8964. [DOI: 10.1021/acs.jpcb.7b07154] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jakub Jungwirth
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 16610 Prague, Czech Republic
- Faculty
of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague, Czech Republic
| | - Jaroslav Šebestík
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 16610 Prague, Czech Republic
| | - Martin Šafařík
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 16610 Prague, Czech Republic
| | - Josef Kapitán
- Department
of Optics, Palacký University, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Petr Bouř
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 16610 Prague, Czech Republic
| |
Collapse
|
41
|
Barone V, Cacelli I, Ferretti A, Prampolini G. Noncovalent Interactions in the Catechol Dimer. Biomimetics (Basel) 2017; 2:E18. [PMID: 31105180 PMCID: PMC6352673 DOI: 10.3390/biomimetics2030018] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Revised: 09/04/2017] [Accepted: 09/05/2017] [Indexed: 12/02/2022] Open
Abstract
Noncovalent interactions play a significant role in a wide variety of biological processes and bio-inspired species. It is, therefore, important to have at hand suitable computational methods for their investigation. In this paper, we report on the contribution of dispersion and hydrogen bonds in both stacked and T-shaped catechol dimers, with the aim of delineating the respective role of these classes of interactions in determining the most stable structure. By using second-order Møller⁻Plesset (MP2) calculations with a small basis set, specifically optimized for these species, we have explored a number of significant sections of the interaction potential energy surface and found the most stable structures for the dimer, in good agreement with the highly accurate, but computationally more expensive coupled cluster single and double excitation and the perturbative triples (CCSD(T))/CBS) method.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri, I-56126 Pisa, Italy.
| | - Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy.
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy.
| | - Alessandro Ferretti
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy.
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy.
| |
Collapse
|
42
|
Rubeš M, Trachta M, Koudelková E, Bulánek R, Kasneryk V, Bludský O. Methane adsorption in ADOR zeolites: a combined experimental and DFT/CC study. Phys Chem Chem Phys 2017; 19:16533-16540. [PMID: 28612872 DOI: 10.1039/c7cp02315d] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Physical adsorption of methane in purely siliceous molecular sieves prepared by a recently discovered synthetic pathway using 2D zeolites as nanoscale building blocks has been investigated by means of combined experimental and theoretical approaches. The DFT/CC-based method has been tested on ADOR zeolites of the UTL family and a few experimentally well-characterized siliceous zeolites. Excellent agreement between theoretical and experimental heats of adsorption has been found for OKO, PCR, MFI, CHA and AEI zeolites. The observed discrepancy for the UTL germanosilicate (2 kJ mol-1) has been plausibly explained using a simple model of D4R defects. The proposed methodology can be used as a reliable characterization tool for newly synthesized silica nanomaterials.
Collapse
Affiliation(s)
- M Rubeš
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nám. 2, 162 10 Prague, Czech Republic.
| | | | | | | | | | | |
Collapse
|
43
|
Fraccarollo A, Canti L, Marchese L, Cossi M. Accurate Evaluation of the Dispersion Energy in the Simulation of Gas Adsorption into Porous Zeolites. J Chem Theory Comput 2017; 13:1756-1768. [DOI: 10.1021/acs.jctc.6b01021] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alberto Fraccarollo
- Dipartimento di Scienze e
Innovazione Tecnologica (DISIT), Università del Piemonte Orientale, via T. Michel 11, 15100 Alessandria, Italy
| | - Lorenzo Canti
- Dipartimento di Scienze e
Innovazione Tecnologica (DISIT), Università del Piemonte Orientale, via T. Michel 11, 15100 Alessandria, Italy
| | - Leonardo Marchese
- Dipartimento di Scienze e
Innovazione Tecnologica (DISIT), Università del Piemonte Orientale, via T. Michel 11, 15100 Alessandria, Italy
| | - Maurizio Cossi
- Dipartimento di Scienze e
Innovazione Tecnologica (DISIT), Università del Piemonte Orientale, via T. Michel 11, 15100 Alessandria, Italy
| |
Collapse
|
44
|
Prampolini G, Campetella M, De Mitri N, Livotto PR, Cacelli I. Systematic and Automated Development of Quantum Mechanically Derived Force Fields: The Challenging Case of Halogenated Hydrocarbons. J Chem Theory Comput 2016; 12:5525-5540. [DOI: 10.1021/acs.jctc.6b00705] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Marco Campetella
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi
13, I-56124 Pisa, Italy
| | - Nicola De Mitri
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Paolo Roberto Livotto
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Ivo Cacelli
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, Via G. Moruzzi 1, I-56124 Pisa, Italy
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi
13, I-56124 Pisa, Italy
| |
Collapse
|
45
|
|
46
|
Cerezo J, Santoro F, Prampolini G. Comparing classical approaches with empirical or quantum-mechanically derived force fields for the simulation electronic lineshapes: application to coumarin dyes. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1888-7] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
47
|
Madarász Á, Berta D, Paton RS. Development of a True Transition State Force Field from Quantum Mechanical Calculations. J Chem Theory Comput 2016; 12:1833-44. [PMID: 26925858 DOI: 10.1021/acs.jctc.5b01237] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Transition state force fields (TSFF) treated the TS structure as an artificial minimum on the potential energy surface in the past decades. The necessary parameters were developed either manually or by the Quantum-to-molecular mechanics method (Q2MM). In contrast with these approaches, here we propose to model the TS structures as genuine saddle points at the molecular mechanics level. Different methods were tested on small model systems of general chemical reactions such as protonation, nucleophilic attack, and substitution, and the new procedure led to more accurate models than the Q2MM-type parametrization. To demonstrate the practicality of our approach, transferrable parameters have been developed for Mo-catalyzed olefin metathesis using quantum mechanical properties as reference data. Based on the proposed strategy, any force field can be extended with true transition state force field (TTSFF) parameters, and they can be readily applied in several molecular mechanics programs as well.
Collapse
Affiliation(s)
- Ádám Madarász
- Research Center for Natural Sciences, Hungarian Academy of Sciences , Magyar Tudosok Korutja 2, H-1117 Budapest, Hungary
| | - Dénes Berta
- Research Center for Natural Sciences, Hungarian Academy of Sciences , Magyar Tudosok Korutja 2, H-1117 Budapest, Hungary
| | - Robert S Paton
- Chemistry Research Laboratory, University of Oxford , Mansfield Road, Oxford OX1 3TA, U.K.,Physical and Theoretical Chemistry Laboratory, University of Oxford , South Parks Road, Oxford OX1 3QZ, U.K
| |
Collapse
|
48
|
Reparameterization of 12-6 Lennard-Jones potentials based on quantum mechanism results for novel tetrahedral N4 (Td) explosives. Theor Chem Acc 2016. [DOI: 10.1007/s00214-015-1800-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
49
|
Tafipolsky M, Ansorg K. Toward a Physically Motivated Force Field: Hydrogen Bond Directionality from a Symmetry-Adapted Perturbation Theory Perspective. J Chem Theory Comput 2016; 12:1267-79. [DOI: 10.1021/acs.jctc.5b01057] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maxim Tafipolsky
- Institut
für Physikalische
und Theoretische Chemie, Universität Würzburg, Campus
Hubland Nord, Emil-Fischer-Strasse 42, D-97074 Würzburg, Germany
| | - Kay Ansorg
- Institut
für Physikalische
und Theoretische Chemie, Universität Würzburg, Campus
Hubland Nord, Emil-Fischer-Strasse 42, D-97074 Würzburg, Germany
| |
Collapse
|