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%.

Embedding Nonrigid Solutes in an Averaged Environment: A Case Study on Rhodopsins.

Many simulation methods concerning solvated molecules are based on the assumption that the solvated species and the solvent can be characterized by some representative structures of the solute and some embedding potential corresponding to this structure. While the averaging of the solvent configurations to obtain an embedding potential has been studied in great detail, this hinges on a single solute structure representation. This assumption is re-examined and generalized for conformationally flexible solutes and tested on 4 nonrigid systems. In this generalized approach, the solute is characterized by a set of representative structures and the corresponding embedding potentials. The representative structures are identified by means of subdividing the statistical ensemble, which in this work is generated by a constant-temperature molecular dynamics simulation. The embedding potential defined in the Frozen-Density Embedding Theory is used to characterize the average effect of the solvent in each subensemble. The numerical examples concern the vertical excitation energies of protonated retinal Schiff bases in protein environments. It is comprehensively shown that subensemble averaging leads to huge computational savings compared with explicit averaging of the excitation energies in the whole ensemble while introducing only minor errors in the case of the systems examined.

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.

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.

Simulations of electric field gradient fluctuations and dynamics around sodium ions in ionic liquids.

The T1 relaxation time measured in nuclear magnetic resonance experiments contains information about electric field gradient (EFG) fluctuations around a nucleus, but computer simulations are typically required to interpret the underlying dynamics. This study uses classical molecular dynamics (MD) simulations and quantum chemical calculations, to investigate EFG fluctuations around a Na+ ion dissolved in the ionic liquid 1-ethyl 3-methylimidazolium tetrafluoroborate, [Im21][BF4], to provide a framework for future interpretation of NMR experiments. Our calculations demonstrate that the Sternheimer approximation holds for Na+ in [Im21][BF4], and the anti-shielding coefficient is comparable to its value in water. EFG correlation functions, CEFG(t), calculated using quantum mechanical methods or from force field charges are roughly equivalent after 200 fs, supporting the use of classical MD for estimating T1 times of monatomic ions in this ionic liquid. The EFG dynamics are strongly bi-modal, with 75%-90% of the de-correlation attributable to inertial solvent motion and the remainder to a highly distributed diffusional processes. Integral relaxation times, ⟨τEFG⟩, were found to deviate from hydrodynamic predictions and were non-linearly coupled to solvent viscosity. Further investigation showed that Na+ is solvated by four tetrahedrally arranged [BF4]- anions and directly coordinated by ∼6 fluorine atoms. Exchange of [BF4]- anions is rare on the 25-50 ns timescale and suggests that motion of solvent-shell [BF4]- is the primary mechanism for the EFG fluctuations. Different couplings of [BF4]- translational and rotational diffusion to viscosity are shown to be the source of the non-hydrodynamic scaling of ⟨τEFG⟩.

Is the non-additive kinetic potential always equal to the difference of effective potentials from inverting the Kohn-Sham equation?

The relation used frequently in the literature according to which the non-additive kinetic potential which is a functional depending on a pair of electron densities is equal (up to a constant) to the difference of two potentials obtained from inverting two Kohn-Sham equations, is examined. The relation is based on a silent assumption that the two densities can be obtained from two independent Kohn-Sham equations, i.e., are vs-representable. It is shown that this assumption does not hold for pairs of densities: ρtot being the Kohn-Sham density in some system and ρB obtained from such partitioning of ρtot that the difference ρtot - ρB vanishes on a Lebesgue measurable volume element. The inversion procedure is still applicable for ρtot - ρB but cannot be interpreted as the inversion of the Kohn-Sham equation. It is rather the inversion of a Kohn-Sham-like equation. The effective potential in the latter equation comprises a "contaminant" that might even not be unique. It is shown that the construction of the non-additive kinetic potential based on the examined relation is not applicable for such pairs.

N-representability of the target density in Frozen-Density Embedding Theory based methods: Numerical significance and its relation to electronic polarization.

The accuracy of any observable derived from multi-scale simulations based on Frozen-Density Embedding Theory (FDET) is affected by two inseparable factors: (i) the approximation for the ExcT nad[ρA,ρB] component of the FDET energy functional and (ii) the choice of the density ρB(r) for which the FDET eigenvalue equation for the embedded wavefunction is solved. A procedure is proposed to estimate the relative significance of these two factors. Numerical examples are given for four weakly bound intermolecular complexes. It is shown that the violation of the non-negativity condition is the principal source of error in the FDET energy if ρB is the density of the isolated environment, i.e., it is generated without taking into account the interactions with the embedded species. Reduction of both the magnitude of the violation of the non-negativity condition and the error in the FDET energy can be pragmatically achieved by means of the explicit treatment of the electronic polarization of the environment.

A non-decomposable approximation on the complete density function space for the non-additive kinetic potential.

A new non-decomposable approximation of the non-additive kinetic energy potential is constructed starting from the same exact property in the limit (ρA → 0 and ∫ρB = 2), as introduced in the work of Lastra et al. [J. Chem. Phys. 129, 074107 (2008)]. In order to cover the complete function space for exponentially decaying densities, the kernel of a differential operator Dγ[ρ] is introduced and analyzed in dependence of γ. The conclusive choice of γ = 1 assures that the solution functions span the complete space of molecular electron densities. As a result, the new approximant preserves the desired feature of the older approximation, which is the reciprocal singularity if the electron density decays exponentially, and eliminates artificial shallow wells (holes), which are responsible for an artificial "charge leak." Numerical considerations using the standard validation procedure introduced by Wesolowski and Weber [Chem. Phys. Lett. 248, 71-76 (1996)] demonstrate the numerical performance of the developed approximation, which increases the range of applicability of semilocal functionals.

Quantifying Fluctuations of Average Solvent Environments for Embedding Calculations.

The viability and effectiveness of replacing an ensemble of embedded solute calculations by a single calculation using an average description of the solvent environment are evaluated. This work explores the fluctuations of the average description of the system obtained in two ways: from calculations on an ensemble of geometries and from an average environment constructed from the same ensemble. To this end, classical molecular dynamics simulations of a rigid acetone solute in SPCE water are performed in order to generate an ensemble of solvent environments. From this ensemble of solvent configurations, a number of different approaches for constructing an average solvent environment are employed. We perform a thorough numerical analysis of the fluctuations of the electrostatic potential experienced by the solute, as well as the resulting fluctuations of the solute's electronic density, measured through its dipole moment and fitted atomic point charges. At the same time, we inspect the accuracy of the methods used to construct average environments. Finally, the proposed method for generating the embedding potential from an average environment density is applied to estimate the solvatochromic shift of the first excitation of acetone. In order to account for quantum confinement effects, which may be important in certain cases, the fluctuations in the shift due to the interaction with the solvent are evaluated using frozen-density-embedding theory. Our results demonstrate that, for normally distributed environments, the constructed average environment is a reasonably good representation of a fluctuating molecular solvent environment. We then provide guidance for future comparisons between these theoretical treatments of solute/solvent systems to experimental measurements.

Benchmark of the Extension of Frozen-Density Embedding Theory to Nonvariational Correlated Methods: The Embedded-MP2 Case.

The extension of the frozen-density embedding theory for nonvariational methods [J. Chem. Theory Comput. 2020, 16, 6880] was utilized to evaluate intermolecular interaction energies for complexes in the Zhao-Truhlar basis set. In the applied method (FDET-MP2-FAT-LDA), the same auxiliary system is used to evaluate the correlation energy by means of the second-order Møller-Plesset perturbation theory (MP2), as in our previous work [J. Chem. Phys. 2019, 150, 121101]. Local density approximation is used for ExcTnad[ρA,ρB] in both cases. Additionally, the contribution to the energy due to the neglected correlation potential was evaluated and analyzed. The domain of applicability of the local density approximation for ExcTnad[ρA,ρB] was determined based on deviations from the interaction energies from the conventional MP2 calculations. The local density approximation for ExcTnad[ρA,ρB] performs well for hydrogen- or dipole-bound complexes. The relative errors in the interaction energy lie within 3-30%. While for charge-transfer complexes, this approximation fails consistently, and for other types of complexes, the performance of this approximation is not systematic. The sources of error are discussed in detail.

The Challenge of Accurate Computation of Two-Photon Absorption Properties of Organic Chromophores in the Condensed Phase.

Two strategies are applied to evaluate the effect of the environment on the two-photon absorption (TPA) cross sections for two characteristic excited states of C2H4 upon complexation with H2O. The supermolecular strategy provides the reference complexation-induced shifts and uses either the EOM-CCSD or ADC(2) method. The embedding strategy is based on frozen-density-embedding theory (FDET) and uses only fundamental constants. The TPA cross sections from high-level supermolecular calculations are extremely basis-set-sensitive. Literature data and the present study indicate that accuracy of the absolute TPA cross sections below 100 atomic units and their shifts below 10 atomic units remains a challenge. The obtained FDET results show a similar basis-set behavior. For the largest basis set (d-aug-cc-pVQZ), TPA cross sections obtained from these two strategies are in excellent agreement. The complexation-induced shifts have the correct sign of the effect and a small (12-33%) relative error in magnitude. The deviations of the FDET-derived shifts from the reference are of similar magnitude as the reliability threshold of the reference shifts.

Frontiers in Multiscale Modeling of Photoreceptor Proteins.

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

On the Correlation Potential in Frozen-Density Embedding Theory.

The correlation functional Ec[ρ] known in Levy's constrained search formulation of density functional theory is also one of the components of the energy functional in frozen-density embedding theory (FDET) [Wesolowski Phys. Rev. A, 77, 2008, 012504] if the embedded wave function has the single-determinant form. The relation between the FDET energy and quantities available from an auxiliary system is derived. In the auxiliary system, Ec[ρ] and its functional derivative (correlation potential) are entirely neglected. The relation is exact up to the quadratic term in density changes caused by electron-electron correlation. It is discussed in view of its practical applications in modeling an electronic structure of embedded species.

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.

Origin of the Solvatochromism in Organic Fluorophores with Flexible Side Chains: A Case Study of Flugi-2.

The emission band for Flugi-2 solvated in dimethyl sulfoxide (DMSO) is obtained from the combined quantum-classical simulations in which the quantum mechanics/molecular mechanics excitation energies are evaluated at the equilibrated segment of the classical molecular dynamics trajectory on the lowest-excited-state potential energy surface. The classical force-field parameters were obtained and validated specifically for the purpose of the present work. The calculated gas-phase to DMSO solvatochromic shift amounts to -0.21 eV, which is in line with the experimentally determined difference between the maxima of the emission bands for Flugi-2 in decane and in DMSO (-0.26 eV). The used model describes rather well the effect of DMSO on the broadening of the emission band. The solvatochromic shift in DMSO originates from two competing effects. The structural deformation of Flugi-2 due to the interaction with DMSO, which results in a positive contribution, and the negative contribution of a larger magnitude due to favorable specific interactions with the solvent. The latter is dominated by a single hydrogen bond between the oxygen atom of a DMSO molecule and the N3 hydrogen atom of the Flugi-2 molecule in which the proton of N3 acts as the donor.

Extension of frozen-density embedding theory for non-variational embedded wavefunctions.

In the original formulation, frozen-density embedding theory [T. A. Wesolowski and A. Warshel, J. Phys. Chem. 97, 8050-8053 (1993); T. A. Wesolowski, Phys. Rev. A 77, 012504 (2008)] concerns multi-level simulation methods in which variational methods are used to obtain the embedded NA-electron wavefunction. In this work, an implicit density functional for the total energy is constructed and used to derive a general expression for the total energy in methods in which the embedded NA electrons are treated non-variationally. The formula is exact within linear expansion in density perturbations. Illustrative numerical examples are provided.

The deconvolution analysis of ATR-FTIR spectra of diacetylene during UV exposure.

We performed a detailed deconvolution analysis of ATR-FTIR peaks of a common diacetylene, 10,12-tricosadiynoic acid (TRCDA) during the polymerization and the blue-to-red transition. Based on the analysis and the solvent dependence on the IR signals, we found that the triple peak from CC stretching mode that has been previously suspected as a consequence of Fermi resonance is rather associated with the macromolecular assembly of TRCDA. Besides these CC triple peaks, we found that the background in the region increased during the UV exposure due to the CC signals from polymers. In addition, the anisotropic compression during polymerization was also detected, which supports the proposed interpretation of X-ray data reported previously. These results are the benefits from the deconvolution analysis.

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

Frozen-Density Embedding Theory (FDET) provides a system-independent formal framework for multi-level computational methods. Despite apparent similarity, the interaction energy components commonly used in QM/MM methods do not have their corresponding counterparts in FDET. We show how the effect of the polarisation on the electron distribution in the environment can be (or is) accounted for either explicitly or implicitly within the FDET framework. Numerical examples are provided for vertical excitation energies in four representative cases of embedded chromophores.