1
|
Slipchenko LV. Detangling Solvatochromic Effects by the Effective Fragment Potential Method. J Phys Chem A 2024; 128:656-669. [PMID: 38193780 DOI: 10.1021/acs.jpca.3c06194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
Understanding molecular interactions in complex systems opens avenues for the efficient design of new materials with target properties. Energy decomposition methods provide a means to obtain a detailed picture of intermolecular interactions. This work introduces a molecular modeling approach for decomposing the solvatochromic shifts of the electronic excited states into the contributions of the individual molecular fragments of the environment surrounding the chromophore. The developed approach is implemented for the QM/EFP (quantum mechanics/effective fragment potential) model that provides a rigorous first-principles-based description of the electronic states of the chromophores in complex polarizable environments. On the example of two model systems, water pentamer and hydrated uracil, we show how the decomposition of the solvatochromic shifts into the contributions of individual solvent water molecules provides a detailed picture of the intermolecular interactions in the ground and excited states of these systems. The analysis also demonstrates the nonadditivity of solute-solvent interactions and the significant contribution of solute polarization to the total values of solvatochromic shifts.
Collapse
Affiliation(s)
- Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47906, United States
| |
Collapse
|
2
|
Pan X, Van R, Pu J, Nam K, Mao Y, Shao Y. Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions. J Chem Theory Comput 2023; 19:8234-8244. [PMID: 37943896 PMCID: PMC10835707 DOI: 10.1021/acs.jctc.3c00973] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2023]
Abstract
In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.
Collapse
Affiliation(s)
- Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Richard Van
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
- Laboratory of Computational Biology, National, Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, United States
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182, United States
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| |
Collapse
|
3
|
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
|
4
|
Anstöter CS, Abou-Hatab S, Thodika M, Matsika S. Effective Fragment Potentials for Microsolvated Excited and Anionic States. J Phys Chem A 2022; 126:8508-8518. [DOI: 10.1021/acs.jpca.2c06122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Cate S. Anstöter
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania19122, United States
| | - Salsabil Abou-Hatab
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania19122, United States
| | - Mushir Thodika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania19122, United States
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania19122, United States
| |
Collapse
|
5
|
Liu X, Humeniuk A, Glover WJ. Conical Intersections in Solution with Polarizable Embedding: Integral-Exact Direct Reaction Field. J Chem Theory Comput 2022; 18:6826-6839. [PMID: 36251342 DOI: 10.1021/acs.jctc.2c00662] [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
A common strategy to exploring the properties and reactivity of complex systems is to use quantum mechanics/molecular mechanics (QM/MM) embedding, wherein a QM region is defined and treated with electronic structure theory, and the remainder of the system is treated with a force field. Important to the description of electronic excited states, especially those of charge-transfer character, is the treatment of the coupling between the QM and MM subsystems. The state of the art is to use a polarizable force field for the MM region and mutually couple the QM wavefunction and MM induced dipoles, in addition to the usual electrostatic embedding, yielding a polarizable embedding (QM/MM-Pol) approach. However, we showed previously that current popular QM/MM-Pol approaches exhibit issues of root flipping and/or incorrect descriptions of electronic crossings in multistate calculations [J. Chem. Theory Comput. 14, 2137 (2018)]. Here, we demonstrate a solution to these problems with an integral-exact reformulation of the direct reaction field approach of Thole and Van Duijnen (QM/MM-IEDRF). The resulting embedding potential includes one- and two-electron operators and many-body dipole-induced dipole interactions and thus includes a natural description of the screening of electron-electron interactions by the MM induced dipoles. Pauli repulsion from the environment is mimicked by effective core potentials on the MM atoms. Inherent to the DRF approach is the assumption that MM dipoles respond instantaneously to the positions of the QM electrons; therefore, dispersion interactions are captured approximately. All electronic states are eigenfunctions of the same Hamiltonian, while the polarization induced in the environment and the associated energetic stabilization are unique to each state. This allows for a consistent definition of transition properties and state crossings. We demonstrate QM/MM-IEDRF by exploring the influence of a (polarizable) inert xenon matrix environment on the conical intersection underlying the photoisomerization of ethylene.
Collapse
Affiliation(s)
- Xiao Liu
- NYU Shanghai, 1555 Century Avenue, Shanghai200122, China
| | - Alexander Humeniuk
- NYU Shanghai, 1555 Century Avenue, Shanghai200122, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai200062, China
| | - William J Glover
- NYU Shanghai, 1555 Century Avenue, Shanghai200122, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai200062, China.,Department of Chemistry, New York University, New York, New York10003, United States
| |
Collapse
|
6
|
Barcza B, Szirmai ÁB, Szántó KJ, Tajti A, Szalay PG. Comparison of approximate intermolecular potentials for ab initio fragment calculations on medium sized N-heterocycles. J Comput Chem 2022; 43:1079-1093. [PMID: 35478353 PMCID: PMC9321956 DOI: 10.1002/jcc.26866] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 03/26/2022] [Accepted: 03/29/2022] [Indexed: 01/15/2023]
Abstract
The ground state intermolecular potential of bimolecular complexes of N-heterocycles is analyzed for the impact of individual terms in the interaction energy as provided by various, conceptually different theories. Novel combinations with several formulations of the electrostatic, Pauli repulsion, and dispersion contributions are tested at both short- and long-distance sides of the potential energy surface, for various alignments of the pyrrole dimer as well as the cytosine-uracil complex. The integration of a DFT/CCSD density embedding scheme, with dispersion terms from the effective fragment potential (EFP) method is found to provide good agreement with a reference CCSD(T) potential overall; simultaneously, a quantum mechanics/molecular mechanics approach using CHELPG atomic point charges for the electrostatic interaction, augmented by EFP dispersion and Pauli repulsion, comes also close to the reference result. Both schemes have the advantage of not relying on predefined force fields; rather, the interaction parameters can be determined for the system under study, thus being excellent candidates for ab initio modeling.
Collapse
Affiliation(s)
- Bónis Barcza
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| | - Ádám B. Szirmai
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| | - Katalin J. Szántó
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| | - Attila Tajti
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| | - Péter G. Szalay
- Institute of Chemistry, Laboratory of Theoretical ChemistryELTE Eötvös Loránd UniversityBudapestHungary
| |
Collapse
|
7
|
Liang W, Pei Z, Mao Y, Shao Y. Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges. J Chem Phys 2022; 156:210901. [PMID: 35676148 PMCID: PMC9162785 DOI: 10.1063/5.0088271] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2022] [Accepted: 05/15/2022] [Indexed: 01/04/2023] Open
Abstract
Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin-orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm-Dancoff approximation to TDDFT, spin-flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine-Ag20 complex), we emphasize that caution is needed to properly account for system-environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.
Collapse
Affiliation(s)
- WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Zheng Pei
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| |
Collapse
|
8
|
Queizán M, Sánchez-Lozano M, Mandado M, Hermida-Ramón JM. A Highly Efficient Neutral Anion Receptor in Polar Environments by Synergy of Anion-π Interactions and Hydrogen Bonding. J Chem Inf Model 2021; 61:4455-4461. [PMID: 34396775 DOI: 10.1021/acs.jcim.1c00595] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Herein, it is shown how anion recognition in highly polar solvents by neutral metal-free receptors is feasible when multiple hydrogen bonding and anion-π interactions are suitably combined. A neutral aromatic molecular tweezer functionalized with azo groups is shown to merge these two kinds of interactions in a unique system and its efficiency as an anion catcher in water is evaluated using first-principles quantum methods. Theoretical calculations unequivocally prove the high thermodynamic stability in water of a model anion, bromide, captured within the tweezer's cavity. Thus, static calculations indicate anion-tweezer interaction energies within the range of covalent or ionic bonds and stability constants in water of more than 10 orders of magnitude. First-principles molecular dynamics calculations also corroborate the stability through the time of the anion-tweezer complex in water. It shows that the anion is always found within the tweezer's cavity due to the combination of the tweezer-anion interactions plus a hydrogen bond between the anion and a water molecule that is inside the tweezer's cavity.
Collapse
Affiliation(s)
- Marta Queizán
- Department of Physical Chemistry, Faculty of Chemistry, University of Vigo, Marcosende (As Lagoas) sn, 36310 Vigo, Galicia, Spain
| | - Marta Sánchez-Lozano
- Department of Physical Chemistry, Faculty of Chemistry, University of Vigo, Marcosende (As Lagoas) sn, 36310 Vigo, Galicia, Spain
| | - Marcos Mandado
- Department of Physical Chemistry, Faculty of Chemistry, University of Vigo, Marcosende (As Lagoas) sn, 36310 Vigo, Galicia, Spain
| | - Jose M Hermida-Ramón
- Department of Physical Chemistry, Faculty of Chemistry, University of Vigo, Marcosende (As Lagoas) sn, 36310 Vigo, Galicia, Spain
| |
Collapse
|
9
|
Pan X, Nam K, Epifanovsky E, Simmonett AC, Rosta E, Shao Y. A simplified charge projection scheme for long-range electrostatics in ab initio QM/MM calculations. J Chem Phys 2021; 154:024115. [PMID: 33445891 DOI: 10.1063/5.0038120] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
In a previous work [Pan et al., Molecules 23, 2500 (2018)], a charge projection scheme was reported, where outer molecular mechanical (MM) charges [>10 Å from the quantum mechanical (QM) region] were projected onto the electrostatic potential (ESP) grid of the QM region to accurately and efficiently capture long-range electrostatics in ab initio QM/MM calculations. Here, a further simplification to the model is proposed, where the outer MM charges are projected onto inner MM atom positions (instead of ESP grid positions). This enables a representation of the long-range MM electrostatic potential via augmentary charges (AC) on inner MM atoms. Combined with the long-range electrostatic correction function from Cisneros et al. [J. Chem. Phys. 143, 044103 (2015)] to smoothly switch between inner and outer MM regions, this new QM/MM-AC electrostatic model yields accurate and continuous ab initio QM/MM electrostatic energies with a 10 Å cutoff between inner and outer MM regions. This model enables efficient QM/MM cluster calculations with a large number of MM atoms as well as QM/MM calculations with periodic boundary conditions.
Collapse
Affiliation(s)
- Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Pkwy, Norman, Oklahoma 73019, USA
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, USA
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Andrew C Simmonett
- National Institutes of Health-National Heart, Lung and Blood Institute, Laboratory of Computational Biology, Bethesda, Maryland 20892, USA
| | - Edina Rosta
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Pkwy, Norman, Oklahoma 73019, USA
| |
Collapse
|
10
|
Kim Y, Bui Y, Tazhigulov RN, Bravaya KB, Slipchenko LV. Effective Fragment Potentials for Flexible Molecules: Transferability of Parameters and Amino Acid Database. J Chem Theory Comput 2020; 16:7735-7747. [PMID: 33236635 DOI: 10.1021/acs.jctc.0c00758] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
An accurate but efficient description of noncovalent interactions is a key to predictive modeling of biological and materials systems. The effective fragment potential (EFP) is an ab initio-based force field that provides a physically meaningful decomposition of noncovalent interactions of a molecular system into Coulomb, polarization, dispersion, and exchange-repulsion components. An EFP simulation protocol consists of two steps, preparing parameters for molecular fragments by a series of ab initio calculations on each individual fragment, and calculation of interaction energy and properties of a total molecular system based on the prepared parameters. As the fragment parameters (distributed multipoles, polarizabilities, localized wave function, etc.) depend on a fragment geometry, straightforward application of the EFP method requires recomputing parameters of each fragment if its geometry changes, for example, during thermal fluctuations of a molecular system. Thus, recomputing fragment parameters can easily become both computational and human bottlenecks and lead to a loss of efficiency of a simulation protocol. An alternative approach, in which fragment parameters are adjusted to different fragment geometries, referred to as "flexible EFP", is explored here. The parameter adjustment is based on translations and rotations of local coordinate frames associated with fragment atoms. The protocol is validated on extensive benchmark of amino acid dimers extracted from molecular dynamics snapshots of a cryptochrome protein. A parameter database for standard amino acids is developed to automate flexible EFP simulations in proteins. To demonstrate applicability of flexible EFP in large-scale protein simulations, binding energies and vertical electron ionization and electron attachment energies of a lumiflavin chromophore of the cryptochrome protein are computed. The results obtained with flexible EFP are in a close agreement with the standard EFP procedure but provide a significant reduction in computational cost.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Yen Bui
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ruslan N Tazhigulov
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Ksenia B Bravaya
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
11
|
Viquez Rojas CI, Slipchenko LV. Exchange Repulsion in Quantum Mechanical/Effective Fragment Potential Excitation Energies: Beyond Polarizable Embedding. J Chem Theory Comput 2020; 16:6408-6417. [PMID: 32786899 DOI: 10.1021/acs.jctc.9b01156] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Hybrid quantum mechanical and molecular mechanical (QM/MM) approaches facilitate computational modeling of large biological and materials systems. Typically, in QM/MM, a small region of the system is modeled with an accurate quantum mechanical method and its surroundings with a more efficient alternative, such as a classical force field or the effective fragment potential (EFP). The reliability of QM/MM calculations depends largely on the treatment of interactions between the two subregions, also known as embedding. The polarizable embedding, which allows mutual polarization between solvent and solute, is considered to be essential for describing electronic excitations in polar solvents. In this work, we employ the QM/EFP model and extend the polarizable embedding by incorporating two short-range terms-a charge penetration correction to the electrostatic term and the exchange-repulsion term-both of which are modeled with one-electron contributions to the quantum Hamiltonian. We evaluate the accuracy of these terms by computing excitation energies across 37 molecular clusters consisting of biologically relevant chromophores surrounded by polar solvent molecules. QM/EFP excitation energies are compared to the fully quantum mechanical calculations with the configuration interaction singles (CIS) method. We find that the charge penetration correction diminishes the accuracy of the QM/EFP calculations. On the other hand, while the effect of exchange-repulsion is negligible for most ππ* transitions, the exchange-repulsion significantly improves description of nπ* transitions with blue solvatochromic shifts. As a result, addition of the exchange-repulsion term improves the overall accuracy of QM/EFP. Performances of QM/EFP models remain similar when excitation energies are modeled with cc-pVDZ and aug-cc-pVDZ basis sets.
Collapse
Affiliation(s)
- Claudia I Viquez Rojas
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47906, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47906, United States
| |
Collapse
|
12
|
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
|
13
|
Sattasathuchana T, Xu P, Gordon MS. An Accurate Quantum-Based Approach to Explicit Solvent Effects: Interfacing the General Effective Fragment Potential Method with Ab Initio Electronic Structure Theory. J Phys Chem A 2019; 123:8460-8475. [PMID: 31365250 DOI: 10.1021/acs.jpca.9b05801] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
An interface between ab initio quantum mechanics (QM) methods and the general effective fragment potential (EFP2) method, QM-EFP2, is implemented in which the intermolecular interactions between a QM molecule and EFP fragments consist of Coulomb, polarization, exchange repulsion (exrep), and dispersion components. In order to ensure accuracy in the QM-EFP2 exrep interaction energy, the EFP2-EFP2 spherical Gaussian overlap (SGO) approximation is abandoned and replaced with the exact electron repulsion integrals (ERI) that are evaluated with a direct method to reduce disk usage. A Gaussian damping function for the QM-EFP2 Coulomb component damps both the EFP nuclear and electronic charges. A new overlap damping function has been implemented for the QM-EFP2 dispersion component. The current QM-EFP2 implementation has been benchmarked with the S22 and S66 data sets and demonstrates excellent agreement with symmetry-adapted perturbation theory (SAPT) for component energies and with coupled cluster theory [CCSD(T)] for the total interaction energies. Water clusters of various sizes (up to 256 water molecules) have been tested; it is shown that the QM-EFP2 method has an accuracy that is comparable to that of EFP2-EFP2. It has been shown previously that the accuracy of EFP2-EFP2 intermolecular interactions is comparable to that of second-order perturbation theory (MP2) or better. The implementation of the distributed data interface (DDI) parallelization scheme significantly improves the efficiency of QM-EFP2 calculations. The time to form the QM-EFP2 Fock operator per SCF iteration for water clusters scales linearly with the number EFP basis functions.
Collapse
Affiliation(s)
- Tosaporn Sattasathuchana
- Department of Chemistry , Iowa State University and Ames Laboratory Ames , Iowa 50011 , United States
| | - Peng Xu
- Department of Chemistry , Iowa State University and Ames Laboratory Ames , Iowa 50011 , United States
| | - Mark S Gordon
- Department of Chemistry , Iowa State University and Ames Laboratory Ames , Iowa 50011 , United States
| |
Collapse
|
14
|
Conrad JA, Kim S, Gordon MS. Ionic liquids from a fragmented perspective. Phys Chem Chem Phys 2019; 21:16878-16888. [PMID: 31359024 DOI: 10.1039/c9cp02836f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The efficacy of using fragmentation methods, such as the effective fragment potential, the fragment molecular orbital and the effective fragment molecular orbital methods is discussed. The advantages and current limitations of these methods are considered, potential improvements are suggested, and a prognosis for the future is provided.
Collapse
Affiliation(s)
- Justin A Conrad
- Department of Chemistry, Iowa State University, Ames, IA 50014, USA.
| | - Shinae Kim
- Department of Chemistry, Iowa State University, Ames, IA 50014, USA.
| | - Mark S Gordon
- Department of Chemistry, Iowa State University, Ames, IA 50014, USA.
| |
Collapse
|
15
|
Tazhigulov RN, Gurunathan PK, Kim Y, Slipchenko LV, Bravaya KB. Polarizable embedding for simulating redox potentials of biomolecules. Phys Chem Chem Phys 2019; 21:11642-11650. [PMID: 31116217 PMCID: PMC6611676 DOI: 10.1039/c9cp01533g] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Redox reactions play a key role in various biological processes, including photosynthesis and respiration. Quantitative and predictive computational characterization of redox events is therefore highly desirable for enriching our knowledge on mechanistic features of biological redox-active macromolecules. Here, we present a computational protocol exploiting polarizable embedding hybrid quantum-classical approach and resulting in accurate estimates of redox potentials of biological macromolecules. A special attention is paid to fundamental aspects of the theoretical description such as the effects of environment polarization and of the long-range electrostatic interactions on the computed energetic parameters. Environment (protein and the solvent) polarization is shown to be crucial for accurate estimates of the redox potential: hybrid quantum-classical results with and without account for environment polarization differ by 1.4 V. Long-range electrostatic interactions are shown to contribute significantly to the computed redox potential value even at the distances far beyond the protein outer surface. The approach is tested on simulating reduction potential of cryptochrome 1 protein from Arabidopsis thaliana. The theoretical estimate (0.07 V) of the midpoint reduction potential is in good agreement with available experimental data (-0.15 V).
Collapse
Affiliation(s)
- Ruslan N Tazhigulov
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA.
| | | | | | | | | |
Collapse
|
16
|
Hapka M, Przybytek M, Pernal K. Second-Order Dispersion Energy Based on Multireference Description of Monomers. J Chem Theory Comput 2019; 15:1016-1027. [PMID: 30525591 DOI: 10.1021/acs.jctc.8b01058] [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/22/2022]
Abstract
We propose a method for calculating a second-order dispersion energy for weakly interacting multireference systems in arbitrary electronic states. It is based on response properties obtained from extended random phase approximation equations. The introduced formalism is general and requires only one- and two-particle reduced density matrices of monomers. We combine the new method with either generalized valence bond perfect pairing (GVB) or complete active space (CAS) self-consistent field description of the interacting systems. In addition to a general scheme, three approximations, leading to significant reduction of the computational cost, are developed by exploiting Dyall partitioning of the monomer Hamiltonians. For model multireference systems (H2···H2 and Be···Be) the method is accurate, unlike its single-reference-based counterpart. Neither GVB nor CAS description of single-reference monomers improves the dispersion energy with respect to the Hartree-Fock-based results.
Collapse
Affiliation(s)
- Michał Hapka
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland.,Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - Michał Przybytek
- Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - Katarzyna Pernal
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland
| |
Collapse
|
17
|
Abstract
We review recent work on property decomposition techniques using quantum chemical methods and discuss some topical applications in terms of quantum mechanics-molecular mechanics calculations and the constructing of properties of large molecules and clusters. Starting out from the so-called LoProp decomposition scheme [Gagliardi et al., J. Chem. Phys., 2004, 121, 4994] for extracting atomic and inter-atomic contributions to molecular properties we show how this method can be generalized to localized frequency-dependent polarizabilities, to localized hyperpolarizabilities and to localized dispersion coefficients. Some applications of the generalized decomposition technique are reviewed - calculations of frequency-dependent polarizabilities, Rayleigh scattering of large clusters, and calculations of hyperpolarizabilities of proteins.
Collapse
Affiliation(s)
- Hans Ågren
- KTH Royal Institute of Technology, School of Engineering Sciences in Chemistry, Biotechnology and Health, Department of Theoretical Chemistry and Biology, SE-106 91 Stockholm, Sweden.
| | | | | |
Collapse
|
18
|
Nanda KD, Krylov AI. The effect of polarizable environment on two-photon absorption cross sections characterized by the equation-of-motion coupled-cluster singles and doubles method combined with the effective fragment potential approach. J Chem Phys 2018; 149:164109. [DOI: 10.1063/1.5048627] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Kaushik D. Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
19
|
Viquez Rojas CI, Fine J, Slipchenko LV. Exchange-repulsion energy in QM/EFP. J Chem Phys 2018; 149:094103. [PMID: 30195305 DOI: 10.1063/1.5043107] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The effective fragment potential (EFP) is a quantum mechanics (QM)-based model designed to accurately describe intermolecular interactions. Hybrid QM/EFP calculations combine quantum mechanical methods with an EFP embedding to study complex systems in which many-body effects are relevant. As in EFP-only calculations, non-bonded interactions between the QM region and EFP fragments are computed as a sum of electrostatic, polarization, dispersion, and exchange-repulsion energies. The exchange-repulsion term is a computational bottleneck of the EFP calculations. Here, we present a general procedure for computing the QM/EFP exchange-repulsion interactions based on one-electron contributions to the QM Hamiltonian, by using Gaussian functions to represent localized molecular orbitals of the effective fragments. The accuracy of the exchange-repulsion and total QM/EFP interaction energies is evaluated on a diverse set of dimers, including complexes from the S22 dataset of non-covalent interactions. In most cases, the QM/EFP energies are at least as accurate as corresponding EFP energies. A simple and computationally efficient form of the introduced QM/EFP exchange-repulsion term will facilitate further developments and applications of QM/EFP methods.
Collapse
Affiliation(s)
| | - Jonathan Fine
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | | |
Collapse
|
20
|
Xu P, Guidez EB, Bertoni C, Gordon MS. Perspective:Ab initioforce field methods derived from quantum mechanics. J Chem Phys 2018. [DOI: 10.1063/1.5009551] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Peng Xu
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Emilie B. Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, USA
| | - Colleen Bertoni
- Argonne Leadership Computing Facility, Argonne, Illinois 60439, USA
| | - Mark S. Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
21
|
Curutchet C, Cupellini L, Kongsted J, Corni S, Frediani L, Steindal AH, Guido CA, Scalmani G, Mennucci B. Density-Dependent Formulation of Dispersion–Repulsion Interactions in Hybrid Multiscale Quantum/Molecular Mechanics (QM/MM) Models. J Chem Theory Comput 2018; 14:1671-1681. [DOI: 10.1021/acs.jctc.7b00912] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Carles Curutchet
- Departament de Farmàcia i Tecnologia Farmacèutica i Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia i Ciències de l’Alimentació, Universitat de Barcelona, Av. Joan XXIII s/n, 08028 Barcelona, Spain
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Stefano Corni
- Dipartimento di Scienze Chimiche, Università di Padova, v. F. Marzolo 1, 35131 Padova, Italy
- CNR-NANO Istituto Nanoscienze, 41125 Modena, Italy
| | - Luca Frediani
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT, The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Arnfinn Hykkerud Steindal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT, The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Ciro A. Guido
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, Nantes 44322 Cedex 3, France
| | - Giovanni Scalmani
- Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| |
Collapse
|