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

A variational formulation of the Harris functional as a correction to approximate Kohn-Sham density functional theory

Accurate descriptions of intermolecular interactions are of great importance in simulations of molecular liquids. We present an electronic structure method that combines the accuracy of the Harris functional approach with the computational efficiency of approximately linear-scaling density functional theory (DFT). The non-variational nature of the Harris functional has been addressed by constructing a Lagrangian energy functional, which restores the variational condition by imposing stationarity with respect to the reference density. The associated linear response equations may be treated with linear-scaling efficiency in an atomic orbital based scheme. Key ingredients to describe the structural and dynamical properties of molecular systems are the forces acting on the atoms and the stress tensor. These first-order derivatives of the Harris Lagrangian have been derived and implemented in consistence with the energy correction. The proposed method allows for simulations with accuracies close to the Kohn-Sham DFT reference. Embedded in the CP2K program package, the method is designed to enable ab initio molecular dynamics simulations of molecular solutions for system sizes of several thousand atoms. Available subsystem DFT methods may be used to provide the reference density required for the energy correction at near linear-scaling efficiency. As an example of production applications, we applied the method to molecular dynamics simulations of the binary mixtures cyclohexane-methanol and toluene-methanol, performed within the isobaric-isothermal ensemble, to investigate the hydrogen bonding network in these non-ideal mixtures.

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 E_xcT^nad[ρ_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 E_xcT^nad[ρ_A,ρ_B] was determined based on deviations from the interaction energies from the conventional MP2 calculations. The local density approximation for E_xcT^nad[ρ_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.

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.

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.