A polarizable three-layer frozen density embedding/molecular mechanics approach

Time-Dependent Multilevel Density Functional Theory

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.

Multiscale Frozen Density Embedding/Molecular Mechanics Approach for Simulating Magnetic Response Properties of Solvated Systems

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.

Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory

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.

Embedding Beyond Electrostatics: The Extended Polarizable Density Embedding Model

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.

Multiple Facets of Modeling Electronic Absorption Spectra of Systems in Solution

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.

The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory

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.

Quantum Mechanics/Fluctuating Charge Protocol to Compute Solvatochromic Shifts

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.