1
|
Hu G, Liu P, Jensen L. Calculating Molecular Polarizabilities Using Exact Frozen Density Embedding with External Orthogonality. J Chem Theory Comput 2024. [PMID: 39105755 DOI: 10.1021/acs.jctc.4c00692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/07/2024]
Abstract
Frozen density embedding (FDE) with freeze-thaw cycles is a formally exact embedding scheme. In practice, this method is limited to systems with small density overlaps when approximate nonadditive kinetic energy functionals are used. It has been shown that the use of approximate nonadditive kinetic energy functionals can be avoided when external orthogonality (EO) is enforced, and FDE can then generate exact results even for strongly overlapping subsystems. In this work, we present an implementation of exact FDEc-EO (coupled FDE TDDFT with EO) for the calculation of polarizabilities in the Amsterdam density functional program package. EO is enforced using the level-shift projection operator method, which ensures that orbitals between fragments are orthogonal. For pure functionals, we show that only the symmetric EO contributions to the induced density matrix are needed. This leads to a simplified implementation for the calculation of polarizability that can exactly reproduce the supermolecular TDDFT results. We further discuss the limitation of exact FDEc-EO in interpreting subsystem polarizabilities due to the nonunique partitioning of the total density. We show that this limitation is due to the fact that subsystem polarizability partitioning is dependent on how the subsystems are initially polarized. As supermolecular virtual orbitals are exactly reproduced, this dependence is attributed to the description of the occupied orbitals. In contrast, for excitations of subsystems that are localized within one subsystem, we show that the excitation energies are stable with respect to the order of polarization. This observation shows that impacts from the nonunique nature of exact FDE on subsystem properties can be minimized by better fragmentation of the supermolecular systems if the property is localized. For global properties like polarizability, this is not the case, and nonuniqueness remains independent of the fragmentation used.
Collapse
Affiliation(s)
- Gaohe Hu
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| | - Pengchong Liu
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| | - Lasse Jensen
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| |
Collapse
|
2
|
Giovannini T, Scavino M, Koch H. Time-Dependent Multilevel Density Functional Theory. J Chem Theory Comput 2024; 20:3601-3612. [PMID: 38648031 DOI: 10.1021/acs.jctc.4c00041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
We present a novel three-layer approach based on multilevel density functional theory (MLDFT) and polarizable molecular mechanics to simulate the electronic excitations of chemical systems embedded in an external environment within the time-dependent DFT formalism. In our method, the electronic structure of a target system, the chromophore, is determined in the field of an embedded inactive layer, which is treated as frozen. Long-range interactions are described by employing the polarizable fluctuating charge (FQ) force field. The resulting MLDFT/FQ thus accurately describes both electrostatics (and polarization) and non-electrostatic target-environment interactions. The robustness and reliability of the approach are demonstrated by comparing our results with experimental data reported for various organic molecules in solution.
Collapse
Affiliation(s)
| | - Marco Scavino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| |
Collapse
|
3
|
Lafiosca P, Rossi F, Egidi F, Giovannini T, Cappelli C. Multiscale Frozen Density Embedding/Molecular Mechanics Approach for Simulating Magnetic Response Properties of Solvated Systems. J Chem Theory Comput 2024; 20:266-279. [PMID: 38109486 PMCID: PMC10782454 DOI: 10.1021/acs.jctc.3c00850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 11/06/2023] [Accepted: 11/10/2023] [Indexed: 12/20/2023]
Abstract
We present a three-layer hybrid quantum mechanical/quantum embedding/molecular mechanics approach for calculating nuclear magnetic resonance (NMR) shieldings and J-couplings of molecular systems in solution. The model is based on the frozen density embedding (FDE) and polarizable fluctuating charges (FQ) and fluctuating dipoles (FQFμ) force fields and permits the accurate ab initio description of short-range nonelectrostatic interactions by means of the FDE shell and cost-effective treatment of long-range electrostatic interactions through the polarizable force field FQ(Fμ). Our approach's accuracy and potential are demonstrated by studying NMR spectra of Brooker's merocyanine in aqueous and nonaqueous solutions.
Collapse
Affiliation(s)
- Piero Lafiosca
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Federico Rossi
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Franco Egidi
- Software
for Chemistry and Materials BV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | | | - Chiara Cappelli
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
4
|
Tran VA, Teucher M, Galazzo L, Sharma B, Pongratz T, Kast SM, Marx D, Bordignon E, Schnegg A, Neese F. Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory. J Phys Chem A 2023; 127:6447-6466. [PMID: 37524058 PMCID: PMC10424240 DOI: 10.1021/acs.jpca.3c02879] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 07/01/2023] [Indexed: 08/02/2023]
Abstract
Nitroxides are common EPR sensors of microenvironmental properties such as polarity, numbers of H-bonds, pH, and so forth. Their solvation in an aqueous environment is facilitated by their high propensity to form H-bonds with the surrounding water molecules. Their g- and A-tensor elements are key parameters to extracting the properties of their microenvironment. In particular, the gxx value of nitroxides is rich in information. It is known to be characterized by discrete values representing nitroxide populations previously assigned to have different H-bonds with the surrounding waters. Additionally, there is a large g-strain, that is, a broadening of g-values associated with it, which is generally correlated with environmental and structural micro-heterogeneities. The g-strain is responsible for the frequency dependence of the apparent line width of the EPR spectra, which becomes evident at high field/frequency. Here, we address the molecular origin of the gxx heterogeneity and of the g-strain of a nitroxide moiety (HMI: 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water. To treat the solvation effect on the g-strain, we combined a multi-frequency experimental approach with ab initio molecular dynamics simulations for structural sampling and quantum chemical EPR property calculations at the highest realistically affordable level, including an explicitly micro-solvated HMI ensemble and the embedded cluster reference interaction site model. We could clearly identify the distinct populations of the H-bonded nitroxides responsible for the gxx heterogeneity experimentally observed, and we dissected the role of the solvation shell, H-bond formation, and structural deformation of the nitroxide in the creation of the g-strain associated with each nitroxide subensemble. Two contributions to the g-strain were identified in this study. The first contribution depends on the number of hydrogen bonds formed between the nitroxide and the solvent because this has a large and well-understood effect on the gxx-shift. This contribution can only be resolved at high resonance frequencies, where it leads to distinct peaks in the gxx region. The second contribution arises from configurational fluctuations of the nitroxide that necessarily lead to g-shift heterogeneity. These contributions cannot be resolved experimentally as distinct resonances but add to the line broadening. They can be quantitatively analyzed by studying the apparent line width as a function of microwave frequency. Interestingly, both theory and experiment confirm that this contribution is independent of the number of H-bonds. Perhaps even more surprisingly, the theoretical analysis suggests that the configurational fluctuation broadening is not induced by the solvent but is inherently present even in the gas phase. Moreover, the calculations predict that this broadening decreases upon solvation of the nitroxide.
Collapse
Affiliation(s)
- Van Anh Tran
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Markus Teucher
- Max-Planck-Institut
für Chemische Energiekonversion, Stiftstraße 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Laura Galazzo
- Department
of Physical Chemistry, University of Geneva, Quai Ernest Ansermet 30, 1211 Geneva, Switzerland
- Faculty
of Chemistry and Biochemistry, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Bikramjit Sharma
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Tim Pongratz
- Fakultät
für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227 Dortmund, Germany
| | - Stefan M. Kast
- Fakultät
für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227 Dortmund, Germany
| | - Dominik Marx
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Enrica Bordignon
- Department
of Physical Chemistry, University of Geneva, Quai Ernest Ansermet 30, 1211 Geneva, Switzerland
- Faculty
of Chemistry and Biochemistry, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Alexander Schnegg
- Max-Planck-Institut
für Chemische Energiekonversion, Stiftstraße 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
5
|
Van den Heuvel W, Reinholdt P, Kongsted J. Embedding Beyond Electrostatics: The Extended Polarizable Density Embedding Model. J Phys Chem B 2023; 127:3248-3256. [PMID: 37002869 DOI: 10.1021/acs.jpcb.2c08721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
Abstract
The polarizable density embedding (PDE) model is a focused QM/QM fragment-based embedding model designed to model solvation effects on molecular properties. We extend the PDE model to include exchange and nonadditive exchange-correlation (for DFT) in the embedding potential in addition to the existing electrostatic, polarization, and nonelectrostatic effects already present. The resulting model, termed PDE-X, yields localized electronic excitation energies that accurately capture the range dependence of the solvent interaction and gives close agreement with full quantum mechanical (QM) results, even when using minimal QM regions. We show that the PDE-X embedding description consistently improves the accuracy of excitation energies for a diverse set of organic chromophores. The improved embedding description leads to systematic solvent effects that do not average out when applying configurational sampling.
Collapse
Affiliation(s)
- Willem Van den Heuvel
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| |
Collapse
|
6
|
Gómez S, Giovannini T, Cappelli C. Multiple Facets of Modeling Electronic Absorption Spectra of Systems in Solution. ACS PHYSICAL CHEMISTRY AU 2022; 3:1-16. [PMID: 36718266 PMCID: PMC9881242 DOI: 10.1021/acsphyschemau.2c00050] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Revised: 11/08/2022] [Accepted: 11/08/2022] [Indexed: 11/24/2022]
Abstract
In this Perspective, we outline the essential physicochemical aspects that need to be considered when building a reliable approach to describe absorption properties of solvated systems. In particular, we focus on how to properly model the complexity of the solvation phenomenon, arising from dynamical aspects and specific, strong solute-solvent interactions. To this end, conformational and configurational sampling techniques, such as Molecular Dynamics, have to be coupled to accurate fully atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) methodologies. By exploiting different illustrative applications, we show that an effective reproduction of experimental spectral signals can be achieved by delicately balancing exhaustive sampling, hydrogen bonding, mutual polarization, and nonelectrostatic effects.
Collapse
|
7
|
Tölle J, Neugebauer J. The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory. J Phys Chem Lett 2022; 13:1003-1018. [PMID: 35061387 DOI: 10.1021/acs.jpclett.1c04023] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
8
|
Ambrosetti M, Skoko S, Giovannini T, Cappelli C. Quantum Mechanics/Fluctuating Charge Protocol to Compute Solvatochromic Shifts. J Chem Theory Comput 2021; 17:7146-7156. [PMID: 34619965 PMCID: PMC8582258 DOI: 10.1021/acs.jctc.1c00763] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
Despite the potentialities
of the quantum mechanics (QM)/fluctuating
charge (FQ) approach to model the spectral properties of solvated
systems, its extensive use has been hampered by the lack of reliable
parametrizations of solvents other than water. In this paper, we substantially
extend the applicability of QM/FQ to solvating environments of different
polarities and hydrogen-bonding capabilities. The reliability and
robustness of the approach are demonstrated by challenging the model
to simulate solvatochromic shifts of four organic chromophores, which
display large shifts when dissolved in apolar, aprotic or polar, protic
solvents.
Collapse
Affiliation(s)
| | - Sulejman Skoko
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|