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

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.

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.

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 (ρ) 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 ρ, and the approximation used for the non-electrostatic component of the FDET embedding potential.

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

The T 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, [Im][BF], to provide a framework for future interpretation of NMR experiments. Our calculations demonstrate that the Sternheimer approximation holds for Na in [Im][BF], and the anti-shielding coefficient is comparable to its value in water.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

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.

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 E [ρ,ρ] component of the FDET energy functional and (ii) the choice of the density ρ(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.

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 (ρ → 0 and ∫ρ = 2), as introduced in the work of Lastra et al. [J. Chem.

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.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods.

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 [ 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 [ 121101]. Local density approximation is used for [ρ,ρ] in both cases.

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 CH upon complexation with HO. 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.

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 [ρ] 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 , 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, [ρ] and its functional derivative (correlation potential) are entirely neglected.

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 E[ρ] performed using the auxiliary functional E_{v_{AB}}^{\rm FDET}[\Psi _A, \rho _B], where Ψ is the embedded N-electron wavefunction and ρ(r) is a non-negative function in real space integrating to a given number of electrons N. 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 ρ(r) reconstructed from X-ray diffraction data for a molecular crystal is demonstrated for the first time.

OpenMolcas: From Source Code to Insight.

In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform.

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.

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.

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

In the original formulation, frozen-density embedding theory [T. A. Wesolowski and A.

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.

BioTIME: A database of biodiversity time series for the Anthropocene.

Motivation: The BioTIME database contains raw data on species identities and abundances in ecological assemblages through time. These data enable users to calculate temporal trends in biodiversity within and amongst assemblages using a broad range of metrics. BioTIME is being developed as a community-led open-source database of biodiversity time series.

Benchmark of Excitation Energy Shifts from Frozen-Density Embedding Theory: Introduction of a Density-Overlap-Based Applicability Threshold.

We present a thorough investigation of the errors in results obtained with the combination of frozen-density embedding theory and the algebraic diagrammatic construction scheme for the polarization propagator of second order (FDE-ADC(2)). The study was carried out on a set of 52 intermolecular complexes with varying interaction strength, each consisting of a chromophore of fundamental interest and a few small molecules in its environment. The errors emerging in frozen-density embedding theory-based methods originate from (a) the solver of the quantum many-body problem used to obtain the embedded wave function (Ψ), (b) the approximation for the explicit density functional for the embedding potential, and (c) the choice of the density representing the environment (ρ( r⃗)).