Multistate QM/QM Extrapolation of UV/Vis Absorption Spectra with Point Charge Embedding

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.

Fast and Accurate Calculation of the UV-Vis Spectrum with the Modified Local Excitation Approximation

The local excitation approximation (LEA), a method for the calculation of electronic excitations localized in a specific region of a molecule, has been modified with new approaches to enhance the accuracy of the original method. The primary concept behind LEA involves isolating the region of interest as a submolecule from the full molecule using a localization method, followed by calculating electronic excitations solely within this submolecule. In this study, we examined approaches that improve the accuracy in describing the region of interest, particularly its molecular orbital energies. Additionally, the localization method was extended with a new projection technique to accelerate calculations. These approaches were studied in time-dependent density functional theory (TDDFT) calculations applied to four testing systems with a chromophore as the region of interest: two basic linear molecules, acrolein surrounded by 24 water molecules, and a model of a green fluorescent protein. For all studied systems, the results of TDDFT calculations combined with LEA exhibited near-zero error when groups of atoms adjacent to the chromophore were explicitly included in the submolecule. This was achieved with at least a quadratic speedup of the calculation time as a function of the submolecule size.

Water Maintains the UV-Vis Spectral Features During the Insertion of Anionic Naproxen and Ibuprofen into Model Cell Membranes

UV-vis spectra of anionic ibuprofen and naproxen in a model lipid bilayer of the cell membrane are investigated using computational techniques in combination with a comparative analysis of drug spectra in purely aqueous environments. The simulations aim at elucidating the intricacies behind the negligible changes in the maximum absorption wavelength in the experimental spectra. A set of configurations of the systems constituted by lipid, water, and drugs or just water and drugs are obtained from classical Molecular Dynamics simulations. UV-vis spectra are computed in the framework of atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) approaches together with Time-Dependent Density Functional Theory (TD-DFT). Our results suggest that the molecular orbitals involved in the electronic transitions are the same, regardless of the chemical environment. A thorough analysis of the contacts between the drug and water molecules reveals that no significant changes in UV-vis spectra are a consequence of ibuprofen and naproxen molecules being permanently microsolvated by water molecules, despite the presence of lipid molecules. Water molecules microsolvate the charged carboxylate group as expected but also microsolvate the aromatic regions of the drugs.

Integrated Multiscale Multilevel Approach to Open Shell Molecular Systems

We present a novel multiscale approach to study the electronic structure of open shell molecular systems embedded in an external environment. The method is based on the coupling of multilevel Hartree-Fock (MLHF) and Density Functional Theory (MLDFT), suitably extended to the unrestricted formalism, to Molecular Mechanics (MM) force fields (FF). Within the ML region, the system is divided into active and inactive parts, thus describing the most relevant interactions (electrostatic, polarization, and Pauli repulsion) at the quantum level. The surrounding MM part, which is formulated in terms of nonpolarizable or polarizable FFs, permits a physically consistent treatment of long-range electrostatics and polarization effects. The approach is extended to the calculation of hyperfine coupling constants and applied to selected nitroxyl radicals in an aqueous solution.

Performance of Multilevel Methods for Excited States

The performance of multilevel quantum chemical approaches, which utilize an atom-based system partitioning scheme to model various electronic excited states, is studied. The considered techniques include the mechanical-embedding (ME) of “our own N-layered integrated molecular orbital and molecular mechanics” (ONIOM) method, the point charge embedding (PCE), the electronic-embedding (EE) of ONIOM, the frozen density-embedding (FDE), the projector-based embedding (PbE), and our local domain-based correlation method. For the investigated multilevel approaches, the second-order algebraic-diagrammatic construction [ADC(2)] approach was utilized as the high-level method, which was embedded in either Hartree–Fock or a density functional environment. The XH-27 test set of Zech et al. [J. Chem. Theory Comput., 2018, 14, 402829906111] was used for the assessment, where organic dyes interact with several solvent molecules. With the selection of the chromophores as active subsystems, we conclude that the most reliable approach is local domain-based ADC(2) [L-ADC(2)], and the least robust schemes are ONIOM-ME and ONIOM-EE. The PbE, FDE, and PCE techniques often approach the accuracy of the L-ADC(2) scheme, but their precision is far behind. The results suggest that a more conservative subsystem selection algorithm or the inclusion of subsystem charge-transfers is required for the atom-based cost-efficient methods to produce high-accuracy excitation energies.

An efficient and robust procedure to calculate absorption spectra of aqueous charged species applied to NO2. Phys Chem Chem Phys 2021;23:14857-14872. [PMID: 34223573 DOI: 10.1039/d1cp00652e]

Accurate calculation of absorption spectra of aqueous NO2- requires rigorously sampling the quantum potential energy surfaces for microsolvation of NO2- with at least five explicit water molecules and embedding the resulting clusters in a continuum solvent accounting for the statistical weighted contributions of individual isomers. This method, which we address as ASCEC + PCM, introduces several desired features when compared against MD simulations derived QM/MM spectra: comparatively fewer explicit solvent molecules to be treated with expensive QM methods, the identification of equilibrium structures in the quantum PES to be used in further vibrational spectroscopy, and the unequivocal identification of cluster orbitals undergoing electronic transitions and charge transfer that originate the spectral bands.

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

We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.

Multilevel Density Functional Theory

Following recent developments in multilevel embedding methods, we introduce a novel density matrix-based multilevel approach within the framework of density functional theory (DFT). In this multilevel DFT, the system is partitioned in an active and an inactive fragment, and all interactions are retained between the two parts. The decomposition of the total system is performed upon the density matrix. The orthogonality between the two parts is maintained by solving the Kohn–Sham equations in the MO basis for the active part only, while keeping the inactive density matrix frozen. This results in the reduction of computational cost. We outline the theory and implementation and discuss the differences and similarities with state-of-the-art DFT embedding methods. We present applications to aqueous solutions of methyloxirane and glycidol.

Energy-Based Molecular Orbital Localization in a Specific Spatial Region

We present a novel energy-based localization procedure able to localize molecular orbitals into predefined spatial regions. The method is defined in a multiscale framework based on the multilevel Hartree-Fock approach. In particular, the system is partitioned into active and inactive fragments. The localized molecular orbitals are obtained maximizing the repulsion between the two fragments. The method is applied to several cases including both conjugated and non-conjugated systems. Our multiscale approach is compared with reference values for both ground-state properties, such as dipole moments, and local excitation energies. The proposed approach is useful to extend the application range of high-level electron correlation methods. In fact, the reduced number of molecular orbitals can lead to a large reduction in the computational cost of correlated calculations.

Combining multilevel Hartree–Fock and multilevel coupled cluster approaches with molecular mechanics: a study of electronic excitations in solutions

We present the coupling of different quantum-embedding approaches with a third molecular-mechanics layer, which can be either polarizable or non-polarizable.