Hydrogen-Bonding-Induced Shifts of the Excitation Energies in Nucleic Acid Bases: An Interplay between Electrostatic and Electron Density Overlap Effects

Calculating Molecular Polarizabilities Using Exact Frozen Density Embedding with External Orthogonality

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.

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.

Assessment of Approximations to the Embedding Potential in Frozen-Density Embedding Theory for the Calculation of Electric Field Gradients

The approximations to the embedding potential in frozen-density embedding theory (FDET) have been assessed for the first time for the calculation of the electric field gradient (EFG) at a nucleus. FDET-based methods using a hierarchy of approximations are applied to evaluate the EFG at the nuclei of an HCl molecule in several noncovalently bound clusters chosen to represent potential liquid or molecular crystal systems. A detailed assessment of such approximations is made for the Hartree-Fock treatment of electron-electron correlation (both in FDET and in the reference calculations for the whole cluster). The emerging choice of the optimal set of approximations is reconfirmed in calculations in which electron-electron calculations are treated at the MP2 level. Our optimized protocol produces average errors in the complexation-induced EFG shift on the order of 25% relative to conventional quantum mechanical calculations for the whole cluster. This protocol is shown to be numerically robust and leads to enormous computational savings compared to a complete quantum mechanical treatment of the embedded species and its environment. For a cluster comprising a Na+ cation and up to 24 water molecules, the computation time is reduced by a factor of 30,000 at the expense of introducing an error in the environment-induced EFG shift of 22%.

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.

Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations

Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.

Symmetrized non-decomposable approximations of the non-additive kinetic energy functional

In subsystem density functional theory (DFT), the bottom-up strategy to approximate the multivariable functional of the non-additive kinetic energy (NAKE) makes it possible to impose exact properties on the corresponding NAKE potential (NAKEP). Such a construction might lead to a non-symmetric and non-homogeneous functional, which excludes the use of such approximations for the evaluation of the total energy. We propose a general formalism to construct a symmetric version based on a perturbation theory approach of the energy expression for the asymmetric part. This strategy is then applied to construct a symmetrized NAKE corresponding to the NAKEP developed recently [Polak et al., J. Chem. Phys. 156, 044103 (2022)], making it possible to evaluate consistently the energy. These functionals were used to evaluate the interaction energy in several model intermolecular complexes using the formal framework of subsystem DFT. The new symmetrized energy expression shows a superior qualitative performance over common decomposable models.

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.

Excitation Energies of Embedded Chromophores from Frozen-Density Embedding Theory Using State-Specific Electron Densities of the Environment

Starting from the Perdew-Levy theorem on extrema of the Hohenberg-Kohn functional, the expression for the vertical excitation energy is derived within the formal framework of Frozen-Density Embedding Theory (FDET) that makes it possible to use state-specific electron densities of the environment (ρ_B) of an embedded species. The derived general expression involves the embedded wave functions for ground and excited states that are orthogonal and is exact up to quadratic terms in the appropriate density expansion. It can be applied in practice using various methods differing in the treatment of the electron-electron correlation for embedded electrons, the method to evaluate different contributions to the excitation energy, the method to generate state-specific ρ_B, and the approximation used for the non-electrostatic component of the FDET embedding potential. The derived expression is applied for 47 local excitations in 10 embedded organic chromophores. The explicit treatment of the differential polarization of ρ_B improves indeed the accuracy of the excitation energy as compared to the implicit treatment in which the same ρ_B is used for all states of embedded chromophore. For 47 local excitations in 10 embedded organic chromophores, the average absolute errors in excitation energies drop from 0.04 to 0.03 eV and their standard deviations from 0.032 to 0.025 eV, respectively. The maximal errors show similar trends.

Efficient Implementation of Density Functional Theory Based Embedding for Molecular and Periodic Systems Using Gaussian Basis Functions

A practical and effective implementation of density functional theory based embedding is reported, which allows us to treat both periodic and aperiodic systems on an equal footing. Its essence is the expansion of orbitals and electron density of the periodic system using Gaussian basis functions, rather than plane-waves, which provides a unique all-electron direct-space representation, thus avoiding the need for pseudopotentials. This makes the construction of embedding potential for a molecular active subsystem due to a periodic environment quite convenient, as transformation between representations is far from trivial. The three flavors of embedding, molecule-in-molecule, molecule-in-periodic, and periodic-in-periodic embedding, are implemented using embedding potentials based on non-additive kinetic energy density functionals (approximate) and level-shift projection operator (exact). The embedding scheme is coupled with a variety of correlated wave function theory (WFT) methods, thereby providing an efficient way to study the ground and excited state properties of low-dimensional systems using high-level methods for the region of interest. Finally, an implementation of real time-time-dependent density functional embedding theory (RT-TDDFET) is presented that uses a projection operator-based embedding potential and provides accurate results compared to full RT-TDDFT for systems with uncoupled excitations. The embedding potential is calculated efficiently using a combination of density fitting and continuous fast multipole method for the Coulomb term. The applicability of (i) WFT-in-DFT embedding, in predicting the adsorption and excitation energies, and (ii) RT-TDDFET, in predicting the absorption spectra, is explored for various test systems.

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.

Protein Response Effects on Cofactor Excitation Energies from First Principles: Augmenting Subsystem Time-Dependent Density-Functional Theory with Many-Body Expansion Techniques

We investigate the possibility of describing protein response effects on a chromophore excitation by means of subsystem time-dependent density-functional theory (sTDDFT) in combination with a many-body expansion (MBE) approach. While sTDDFT is in principle intrinsically able to include such contributions, addressing cofactor excitations in protein models or entire proteins with full environment-response treatments is currently out of reach. Taking different model structures of the green fluorescent protein (GFP) and bovine rhodopsin as examples, we demonstrate that an embedded-MBE approach based on sTDDFT in its simplest version leads to a good agreement of the predicted protein response effect already at second order. To reproduce reference response effects from nonsubsystem TDDFT calculations quantitatively (error ≤ 5%), however, a third- or even fourth-order MBE may be required. For the latter case, we explore a selective inclusion of fourth-order terms that drastically reduces the computational burden. In addition, we demonstrate how this sTDDFT-MBE treatment can be utilized as an analysis tool to identify residues with dominant response contributions. This, in turn, can be employed to arrive at smaller structural models for light-absorbing proteins, which still feature all of the main characteristics in terms of photoresponse properties.

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.

Subsystem-Based GW/Bethe-Salpeter Equation

Subsystem Density-Functional Theory and its extension to excited states, namely, subsystem Time-Dependent Density-Functional Theory, have been proven to be efficient and accurate fragmentation approaches for ground and excited states. In the present study we extend this approach to the subsystem-based description of total systems by means of GW and the Bethe-Salpeter equation (BSE). For this, we derive the working equations starting from a subsystem-based partitioning of the screened-Coulomb interaction for an arbitrary number of subsystems. Making use of certain approximations, we develop a parameter-free approach in which environmental screening contributions are effectively included for each subsystem. We demonstrate the applicability of these approximations by comparing quasi-particle energies and excitation energies from subsystem-based GW/BSE calculations to the supermolecular reference. Furthermore, we demonstrate the computational efficiency and the usefulness of this method for the description of photoinduced processes in complex chemical environments.

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.

Electronic couplings for photo-induced processes from subsystem time-dependent density-functional theory: The role of the diabatization

Subsystem time-dependent density-functional theory (sTDDFT) making use of approximate non-additive kinetic energy (NAKE) functionals is known to be capable of describing excitation energy transfer processes in a variety of applications. Here, we show that sTDDFT, especially when combined with projection-based embedding (PbE), can be employed for the entire range of photo-induced electronic couplings essential for modeling photophysical properties of complex chemical and biological systems and therefore represents a complete toolbox for this class of problems. This means that it is capable of capturing the interaction/coupling associated with local- and charge-transfer (CT) excitons. However, this requires the choice of a reasonable diabatic basis. We therefore propose different diabatization strategies of the virtual orbital space in PbE-sTDDFT and show how CT excitations can be included in sTDDFT using NAKE functionals via a phenomenological approach. Finally, these electronic couplings are compared to couplings from a multistate fragment excitation difference (FED)-fragment charge difference (FCD) diabatization procedure. We show that both procedures, multistate FED-FCD and sTDDFT (with the right diabatization procedure chosen), lead to an overall good agreement for the electronic couplings, despite differences in their general diabatization strategy. We conclude that the entire range of photo-induced electronic couplings can be obtained using sTDDFT (with the right diabatization procedure chosen) in a black-box manner.

Environmental Effects with Frozen-Density Embedding in Real-Time Time-Dependent Density Functional Theory Using Localized Basis Functions

Frozen-density embedding (FDE) represents a versatile embedding scheme to describe the environmental effect on electron dynamics in molecular systems. The extension of the general theory of FDE to the real-time time-dependent Kohn-Sham method has previously been presented and implemented in plane waves and periodic boundary conditions [Pavanello, M.; J. Chem. Phys. 2015, 142, 154116]. In the current paper, we extend our recent formulation of the real-time time-dependent Kohn-Sham method based on localized basis set functions and developed within the Psi4NumPy framework to the FDE scheme. The latter has been implemented in its "uncoupled" flavor (in which the time evolution is only carried out for the active subsystem, while the environment subsystems remain at their ground state), using and adapting the FDE implementation already available in the PyEmbed module of the scripting framework PyADF. The implementation was facilitated by the fact that both Psi4NumPy and PyADF, being native Python API, provided an ideal framework of development using the Python advantages in terms of code readability and reusability. We employed this new implementation to investigate the stability of the time-propagation procedure, which is based on an efficient predictor/corrector second-order midpoint Magnus propagator employing an exact diagonalization, in combination with the FDE scheme. We demonstrate that the inclusion of the FDE potential does not introduce any numerical instability in time propagation of the density matrix of the active subsystem, and in the limit of the weak external field, the numerical results for low-lying transition energies are consistent with those obtained using the reference FDE calculations based on the linear-response TDDFT. The method is found to give stable numerical results also in the presence of a strong external field inducing nonlinear effects. Preliminary results are reported for high harmonic generation (HHG) of a water molecule embedded in a small water cluster. The effect of the embedding potential is evident in the HHG spectrum reducing the number of the well-resolved high harmonics at high energy with respect to the free water. This is consistent with a shift toward lower ionization energy passing from an isolated water molecule to a small water cluster. The computational burden for the propagation step increases approximately linearly with the size of the surrounding frozen environment. Furthermore, we have also shown that the updating frequency of the embedding potential may be significantly reduced, much less than one per time step, without jeopardizing the accuracy of the transition energies.

Embedding-theory-based simulations using experimental electron densities for the environment

The basic idea of frozen-density embedding theory (FDET) is the constrained minimization of the Hohenberg-Kohn density functional EHK[ρ] performed using the auxiliary functional E_{v_{AB}}^{\rm FDET}[\Psi _A, \rho _B], where ΨA is the embedded NA-electron wavefunction and ρB(r) is a non-negative function in real space integrating to a given number of electrons NB. This choice of independent variables in the total energy functional E_{v_{AB}}^{\rm FDET}[\Psi _A, \rho _B] makes it possible to treat the corresponding two components of the total density using different methods in multi-level simulations. The application of FDET using ρB(r) reconstructed from X-ray diffraction data for a molecular crystal is demonstrated for the first time. For eight hydrogen-bonded clusters involving a chromophore (represented as ΨA) and the glycylglycine molecule [represented as ρB(r)], FDET is used to derive excitation energies. It is shown that experimental densities are suitable for use as ρB(r) in FDET-based simulations.

Nature inspired solid-liquid phase amphibious adhesive

Here we report a new class of bio-inspired solid-liquid adhesive, obtained by simple mechanical dispersion of PVDF (polyvinylidene fluoride) (solid spheres) into PDMS (polydimethylsiloxane) (liquid). The adhesive behavior arises from strong solid-liquid interactions. This is a chemical reaction free adhesive (no curing time) that can be repeatedly used and is capable of instantaneously joining a large number of diverse materials (metals, ceramic, and polymer) in air and underwater. The current work is a significant advance in the development of amphibious multifunctional adhesives and presents potential applications in a range of sealing applications, including medical ones.

Approximate versus Exact Embedding for Chiroptical Properties: Reconsidering Failures in Potential and Response

We investigate the suitability of subsystem time-dependent density-functional theory (sTDDFT) for describing chiroptical properties with a focus on optical rotation parameters. Our starting point is a new implementation of the recently proposed projection-based, coupled frozen-density embedding (FDEc) framework. We adapt the generalized, non-Hermitian formulation of TDDFT and derive corresponding expressions for regular and damped response properties from subsystem TDDFT. We verify that our implementation of this "exact" formulation allows to reproduce supermolecular results of electronic circular dichroism (ECD) spectra, of optical rotatory dispersion, and of polarizabilities. We present a systematic test of the main approximations typically introduced in practical frozen-density embedding (FDE) calculations of response properties: (i) the use of approximate nonadditive kinetic-energy (NAKE) functionals, which can be avoided through projection techniques, (ii) the use of monomer (subsystem) basis sets rather than supersystem basis sets, and (iii) the neglect of intersubsystem response coupling within the so-called uncoupled FDE (or FDEu) approximation. While approximation (i) is known to generally lead to large errors for covalently bound subsystems, we present cases in which either the basis set or the coupling step are similarly or even (much) more important. In particular, we explicitly demonstrate by comparison to a fully coupled calculation that missing intersubsystem response couplings are responsible for the failure of FDE reported in a previous study [ J. Chem. Theory Comput. 2015, 11, 5305-5315]. We show that good agreement with reference results can be obtained in this case even with standard NAKE approximations for the FDE potentials and efficient monomer basis sets, making calculations for larger systems well accessible.

Embedding Methods for Quantum Chemistry: Applications from Materials to Life Sciences

Quantum mechanical embedding methods hold the promise to transform not just the way calculations are performed, but to significantly reduce computational costs and improve scaling for macro-molecular systems containing hundreds if not thousands of atoms. The field of embedding has grown increasingly broad with many approaches of different intersecting flavors. In this perspective, we lay out the methods into two streams: QM:MM and QM:QM, showcasing the advantages and disadvantages of both. We provide a review of the literature, the underpinning theories including our contributions, and we highlight current applications with select examples spanning both materials and life sciences. We conclude with prospects and future outlook on embedding, and our view on the use of universal test case scenarios for cross-comparisons of the many available (and future) embedding theories.

Synergistic removal of tylosin/sulfamethoxazole and copper by nano-hydroxyapatite modified biochar

Antibiotics and heavy metals are frequently detected simultaneously in aquatic environment. In this study, we investigated the removal performance of biochar modified with nano-hydroxyapatite (nHAP, nHAP@biochar) on tylosin (TYL) /sulfamethoxazole (SMX) and Cu(II) simultaneously. Six nHAP@biochars were prepared with different feedstock and nHAP and biomass ratios. The influences of feedstock and nHAP and biomass ratios, interaction of TYL/SMX and Cu(II) and thermodynamic study were investigated. The adsorption quantities on nHAP@biochars prepared by wood-processing residues were higher than by Chinese medicine residues. The adsorption amounts of TYL decreased with the addition of Cu(II), while the adsorption quantities of SMX increased. The adsorptions of Cu(II) were promoted by TYL and changed slightly with the increasing of SMX. Specific surface area and pore size were two of the main factors influencing the adsorption capacities of nHAP@biochars. According to density functional theory, nHAP@biochar-TYL-Cu and nHAP@biochar-Cu-SMX were more existed in the systems.

Quantum dynamics of vibration-assisted excitation energy transfer in phycobiliprotein light-harvesting complex

Phycobiliprotein is a light-harvesting complex containing linear tetrapyrrole bilin pigments that are responsible for absorption and funneling the sun's energy in cryptophytes algae. In particular, the protein structure determines relative positions and orientations of the pigments and thus controls energy transfer pathways. The present research reveals the impact of molecular vibrations (in the 850-2700 cm^-1 region) on excitation energy transfer in phycobiliprotein. The analysis of the excitation energy transfer pathways indicates a possibility of the coherent mechanism of energy transfer (delocalization) in central dihydrobiliverdin pigments and incoherent vibration-assisted energy transfer to peripheral phycocyanobilin pigments at a sub-picosecond time scale. A computational approach that enables modeling the dynamics of the excitation energy transfer with the quantum master equation formalism employing Huang-Rhys factors to describe electronic-vibrational coupling has been developed. The computational methodology has been implemented in PyFREC software.

Analytical Nuclear Excited-State Gradients for the Second-Order Approximate Coupled-Cluster Singles and Doubles (CC2) Method Employing Uncoupled Frozen-Density Embedding

We report the derivation and implementation of analytical orbital-relaxed properties and nuclear gradients for excited states using the second-order approximate coupled-cluster singles and doubles (CC2) model combined with uncoupled frozen-density embedding (FDEu). An implementation of the algebraic diagrammatic construction through second-order ADC(2), which arises from simplification of RICC2 FDEu, is also presented. In order to ensure a RICC2 FDEu Lagrange functional that is linear in the Lagrange multipliers, the Hartree-Fock density is employed for the target subsystem in the embedding contributions. The accuracy of the new scheme is assessed using the carbon monoxide molecule, 4-aminophthalimide, and a benzonitrile dimer, revealing that the obtained errors are below the method error of RICC2. Using density functional theory for the environment, the efficiency of the new method is illustrated by computing the perturbed excited-state dipole moment of a chromophore in a biological environment. For this system, comprising 32 molecules consisting of 366 atoms in total, the computation requires only a couple of days on a standard compute node. RICC2 FDEu thus enables large-scale calculations of ab initio wave functions for molecules in complex environments as routine applications.

Explicit vs. implicit electronic polarisation of environment of an embedded chromophore in frozen-density embedding theory

A comparison of strategies to account for environment polarisation in Frozen Density Embedding Theory (FDET).

Analytical nuclear excited-state gradients for the Tamm-Dancoff approximation using uncoupled frozen-density embedding

We report the derivation and implementation of analytical nuclear gradients for excited states using time-dependent density functional theory using the Tamm-Dancoff approximation combined with uncoupled frozen-density embedding using density fitting. Explicit equations are presented and discussed. The implementation is able to treat singlet as well as triplet states and functionals using the local density approximation, the generalized gradient approximation, combinations with Hartree-Fock exchange (hybrids), and range-separated functionals such as CAM-B3LYP. The new method is benchmarked against supermolecule calculations in two case studies: The solvatochromic shift of the (vertical) fluorescence energy of 4-aminophthalimide on solvation, and the first local excitation of the benzonitrile dimer. Whereas for the 4-aminophthalimide-water complex deviations of about 0.2 eV are obtained to supermolecular calculations, for the benzonitrile dimer the maximum error for adiabatic excitation energies is below 0.01 eV due to a weak coupling of the subsystems. © 2017 Wiley Periodicals, Inc.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

Revealing electronic open quantum systems with subsystem TDDFT

Open quantum systems (OQSs) are perhaps the most realistic systems one can approach through simulations. In recent years, describing OQSs with Density Functional Theory (DFT) has been a prominent avenue of research with most approaches based on a density matrix partitioning in conjunction with an ad-hoc description of system-bath interactions. We propose a different theoretical approach to OQSs based on partitioning of the electron density. Employing the machinery of subsystem DFT (and its time-dependent extension), we provide a novel way of isolating and analyzing the various terms contributing to the coupling between the system and the surrounding bath. To illustrate the theory, we provide numerical simulations on a toy system (a molecular dimer) and on a condensed phase system (solvated excimer). The simulations show that non-Markovian dynamics in the electronic system-bath interactions are important in chemical applications. For instance, we show that the superexchange mechanism of transport in donor-bridge-acceptor systems is a non-Markovian interaction between the donor-acceptor (OQS) with the bridge (bath) which is fully characterized by real-time subsystem time-dependent DFT.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.

Analytical gradients for excitation energies from frozen-density embedding

Analytical excitation-energy gradients from frozen-density embedding–time-dependent density functional theory are derived and implemented, which are important for photochemistry in complex systems.

Frozen-density embedding theory with average solvent charge densities from explicit atomistic simulations

Besides molecular electron densities obtained within the Born–Oppenheimer approximation (ρ_B(r)) to represent the environment, the ensemble averaged density (〈ρ_B〉(r)) is also admissible in frozen-density embedding theory (FDET) [Wesolowski, Phys. Rev. A, 2008, 77, 11444].

External orthogonality in subsystem time-dependent density functional theory

We present a subsystem TDDFT method which maintains inter-subsystem orbital orthogonality.

Equilibrium Geometries of Noncovalently Bound Intermolecular Complexes Derived from Subsystem Formulation of Density Functional Theory

The subsystem formulation of density functional theory is used to obtain equilibrium geometries and interaction energies for a representative set of noncovalently bound intermolecular complexes. The results are compared with literature benchmark data. The range of applicability of two considered approximations to the exchange-correlation- and nonadditive kinetic energy components of the total energy is determined. Local density approximation, which does not involve any empirical parameters, leads to excellent intermolecular equilibrium distances for hydrogen-bonded complexes (maximal error 0.13 Å for NH3-NH3). It is a method of choice for a wide class of weak intermolecular complexes including also dipole-bound and the ones formed by rare gas atoms or saturated hydrocarbons. The range of applicability of the chosen generalized gradient approximation, which was shown in our previous works to lead to good interaction energies in such complexes, where π-electrons are involved in the interaction, remains limited to this group because it improves neither binding energies nor equilibrium geometries in the wide class of complexes for which local density approximation is adequate. An efficient energy minimization procedure, in which optimization of the geometry and the electron density of each subsystem is made simultaneously, is proposed and tested.

Embedding potentials for excited states of embedded species

Frozen-Density-Embedding Theory (FDET) is a formalism to obtain the upper bound of the ground-state energy of the total system and the corresponding embedded wavefunction by means of Euler-Lagrange equations [T. A. Wesolowski, Phys. Rev. A 77(1), 012504 (2008)]. FDET provides the expression for the embedding potential as a functional of the electron density of the embedded species, electron density of the environment, and the field generated by other charges in the environment. Under certain conditions, FDET leads to the exact ground-state energy and density of the whole system. Following Perdew-Levy theorem on stationary states of the ground-state energy functional, the other-than-ground-state stationary states of the FDET energy functional correspond to excited states. In the present work, we analyze such use of other-than-ground-state embedded wavefunctions obtained in practical calculations, i.e., when the FDET embedding potential is approximated. Three computational approaches based on FDET, that assure self-consistent excitation energy and embedded wavefunction dealing with the issue of orthogonality of embedded wavefunctions for different states in a different manner, are proposed and discussed.