Quasienergy formulation of damped response theory

Subspace Methods for the Simulation of Molecular Response Properties on a Quantum Computer

We explore Davidson methods for obtaining excitation energies and other linear response properties within the recently developed quantum self-consistent linear response (q-sc-LR) method. Davidson-type methods allow for obtaining only a few selected excitation energies without explicitly constructing the electronic Hessian since they only require the ability to perform Hessian-vector multiplications. We apply the Davidson method to calculate the excitation energies of hydrogen chains (up to H10) and analyze aspects of statistical noise for computing excitation energies on quantum simulators. Additionally, we apply Davidson methods for computing linear response properties such as static polarizabilities for H2, LiH, H2O, OH-, and NH3, and show that unitary coupled cluster outperforms classical projected coupled cluster for molecular systems with strong correlation. Finally, we formulate the Davidson method for damped (complex) linear response, with application to the nitrogen K-edge X-ray absorption of ammonia, and the C6 coefficients of H2, LiH, H2O, OH-, and NH3.

Quantum Simulation of Molecular Response Properties in the NISQ Era

Accurate modeling of the response of molecular systems to an external electromagnetic field is challenging on classical computers, especially in the regime of strong electronic correlation. In this article, we develop a quantum linear response (qLR) theory to calculate molecular response properties on near-term quantum computers. Inspired by the recently developed variants of the quantum counterpart of equation of motion (qEOM) theory, the qLR formalism employs "killer condition" satisfying excitation operator manifolds that offer a number of theoretical advantages along with reduced quantum resource requirements. We also used the qEOM framework in this work to calculate the state-specific response properties. Further, through noiseless quantum simulations, we show that response properties calculated using the qLR approach are more accurate than the ones obtained from the classical coupled-cluster-based linear response models due to the improved quality of the ground-state wave function obtained using the ADAPT-VQE algorithm.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Efficient Protocol for Computing MCD Spectra in a Broad Frequency Range Combining Resonant and Damped CC2 Quadratic Response Theory

Coupled cluster response theory offers a path to high-accuracy calculations of spectroscopic properties, such as magnetic circular dichroism (MCD). However, divergence or slow convergence issues are often encountered for electronic transitions in high-energy regions with a high density of states. This is here addressed for MCD by an implementation of damped quadratic response theory for resolution-of-identity coupled cluster singles-and-approximate-doubles (RI-CC2), along with an implementation of the MCD A term from resonant response theory. Combined, damped and resonant response theory calculations provide an efficient strategy to obtain MCD spectra over a broad frequency range and for systems that include highly symmetric molecules with degenerate excited states. The protocol is illustrated by application to zinc tetrabenzoporphyrin in the energy region of 2-8 eV and comparison to experimental data. Timings are reported for the resonant and damped approaches, showing that a greater part of the calculation time is consumed by the construction of the building blocks for the final MCD ellipticity. A recommendation on how to use the procedure is outlined.

A method to capture the large relativistic and solvent effects on the UV-vis spectra of photo-activated metal complexes

We have recently developed a method based on relativistic time-dependent density functional theory (TD-DFT) that allows the calculation of electronic spectra in solution (Creutzberg, Hedegård, J. Chem. Theory Comput.18, 2022, 3671). This method treats the solvent explicitly with a classical, polarizable embedding (PE) description. Furthermore, it employs the complex polarization propagator (CPP) formalism which allows calculations on complexes with a dense population of electronic states (such complexes are known to be problematic for conventional TD-DFT). Here, we employ this method to investigate both the dynamic and electronic effects of the solvent for the excited electronic states of trans-trans-trans-[Pt(N₃)₂(OH)₂(NH₃)₂] in aqueous solution. This complex decomposes into species harmful to cancer cells under light irradiation. Thus, understanding its photo-physical properties may lead to a more efficient method to battle cancer. We quantify the effect of the underlying structure and dynamics by classical molecular mechanics simulations, refined with a subsequent DFT or semi-empirical optimization on a cluster. Moreover, we quantify the effect of employing different methods to set up the solvated system, e.g., how sensitive the results are to the method used for the refinement, and how large a solvent shell that is required. The electronic solvent effect is always included through a PE potential.

Cherry-Picking Resolvents: Recovering the Valence Contribution in X-ray Two-Photon Absorption within the Core-Valence-Separated Equation-of-Motion Coupled-Cluster Response Theory

Calculations of first-order response wave functions in the X-ray regime often diverge within correlated frameworks such as equation-of-motion coupled-cluster singles and doubles (EOM-CCSD), a consequence of the coupling with the valence ionization continuum. Here, we extend our strategy of introducing a hierarchy of approximations to the EOM-EE-CCSD resolvent (or, inversely, the model Hamiltonian) involved in the response equations for the calculation of X-ray two-photon absorption (X2PA) cross sections. We exploit the frozen-core core-valence separation (fc-CVS) scheme to first decouple the core and valence Fock spaces, followed by a separate approximate treatment of the valence resolvent. We demonstrate the robust convergence of X-ray response calculations within this framework and compare X2PA spectra of small benchmark molecules with the previously reported density functional theory results.

Quantum chemistry study of the multiphoton absorption in enhanced green fluorescent protein at the single amino acid residue level

The chromophore (CRO) of fluorescent proteins (FPs) is embedded in a complex environment that is a source of specific interactions with the CRO. Understanding how these interactions influence FPs spectral properties is important for a directed design of novel markers with desired characteristics. In this work, we apply computational chemistry methods to gain insight into one-, two- and three-photon absorption (1PA, 2PA, 3PA) tuning in enhanced green fluorescent protein (EGFP). To achieve this goal, we built EGFP models differing in: i) number and position of hydrogen-bonds (h-bonds) donors to the CRO and ii) the electric field, as approximated by polarizable force  field, acting on the CRO. We  find that h-bonding to the CRO's phenolate oxygen results in stronger one- and multiphoton absorption. The brighter absorption can be also achieved by creating more positive electric field near the CRO's phenolate moiety. Interestingly, while individual CRO-environment h-bonds usually enhance 1PA and 2PA, it takes a few h-bond donors to enhance 3PA. Clearly, response of the absorption intensity to many-body effects depends on the excitation mechanism.  We further employ symmetry-adapted perturbation theory (SAPT) to reveal excellent (2PA) and good (3PA) correlation of multiphoton intensity with electrostatic and induction interaction energies. This points to importance of accounting for mutual CRO-environment polarization in quantitative calculations of absorption spectra in FPs.

Polarizable Embedding Complex Polarization Propagator in Four- and Two-Component Frameworks

Explicit embedding methods combined with the complex polarization propagator (CPP) enable the modeling of spectroscopy for increasingly complex systems with a high density of states. We present the first derivation and implementation of the CPP in four- and exact-two-component (X2C) polarizable embedding (PE) frameworks. We denote the developed methods PE-4c-CPP and PE-X2C-CPP, respectively. We illustrate the methods by estimating the solvent effect on ultraviolet-visible (UV-vis) and X-ray atomic absorption (XAS) spectra of [Rh(H₂O)₆]³⁺ and [Ir(H₂O)₆]³⁺ immersed in aqueous solution. We moreover estimate solvent effects on UV-vis spectra of a platinum complex that can be photochemically activated (in water) to kill cancer cells. Our results clearly show that the inclusion of the environment is required: UV-vis and (to a lesser degree) XAS spectra can become qualitatively different from vacuum calculations. Comparison of PE-4c-CPP and PE-X2C-CPP methods shows that X2C essentially reproduces the solvent effect obtained with the 4c methods.

Origin invariant electronic circular dichroism in the length dipole gauge without London atomic orbitals

We present a method for obtaining origin-independent electronic circular dichroism (ECD) in the length-gauge representation LG(OI) without the usage of London atomic orbitals. This approach builds upon the work by Caricato [J. Chem. Phys. 153, 151101 (2020)] and is applied to rotatory strengths and ECD spectra from damped response theory. Numerical results are presented for time-dependent Hartree-Fock and density-functional theory, the second-order algebraic diagrammatic construction method, and linear-response coupled-cluster theory with singles and approximate doubles. We can support the finding that the common choice of placing the gauge origin in the center of mass of a molecule in conventional length-gauge calculations involving chiroptical properties might not be optimal and show that LG(OI) is a valuable alternative for the origin-independent calculation of ECD spectra. We show that, for a limited test set, the convergence of the rotatory strengths calculated with the LG(OI) approach toward the basis-set limit tends to be faster than for the established velocity gauge representation. Relationships between the sum-over-states expression of the optical rotation in the LG(OI) framework and its representation in terms of response functions are analyzed.

Probing Molecular Chirality of Ground and Electronically Excited States in the UV-vis and X-ray Regimes: An EOM-CCSD Study

We present several strategies for computing electronic circular dichroism (CD) spectra across different frequency ranges at the equation-of-motion coupled-cluster singles and doubles level of theory. CD spectra of both ground and electronically excited states are discussed. For selected cases, the approach is compared with coupled-cluster linear response results as well as time-dependent density functional theory. The extension of the theory to include the effect of spin-orbit coupling is presented and illustrated by calculations of X-ray CD spectra at the L-edge.

The role of hydrogen bonds and electrostatic interactions in enhancing two-photon absorption in green and yellow fluorescent proteins

The spectral properties of fluorescent proteins (FPs) depend on the protein environment of the chromophore (CRO). A deeper understanding of the CRO - environment interactions in terms of FPs spectral characteristics will allow for a rational design of novel markers with desired properties. Here, we are taking a step towards achieving this important goal. With the time-dependent density functional theory (TDDFT), we calculate one- and two-photon absorption (OPA and TPA) spectra for 5 green FPs (GFPs) and 3 yellow FPs (YFPs) differing in amino acid sequence. The goal is to reveal a role of: (i) electrostatic interactions, (ii) hydrogen-bonds (h-bonds), and (iii) h-bonds together with distant electrostatic field in absorption spectra tuning. Our results point to design hypothesis towards FPs optimised for TPA-based applications. Both h-bonds and electrostatic interactions co-operate in enhancing TPA cross-section (σ TPA ) for the S 0 ->S 1   transition in GFPs. Furthermore, it seems that details of h-bonds network in the CRO's vicinity influences σ TPA   response to CRO - environment electrostatic interactions in YFPs. We postulate that engineering FPs with more hydrophilic CRO's environment can lead to greater σ TPA . We also find that removing h-bonds formed with the CRO's phenolate leads to TPA enhancement for transition to higher excited states than S 1 . Particularly Y145 and T203 residues are important in this regard.

Self-assembling, structure and nonlinear optical properties of fluorescent organic nanoparticles in water

Owing to their intense emission, low toxicity and solubility in aqueous medium, fluorescent organic nanoparticles (FONs) have emerged as promising alternatives to inorganic ones for the realization of exogenous probes for bioimaging applications. However, the intimate structure of FONs in solution, as well as the role played by intermolecular interactions on their optical properties, remains challenging to study. Following a recent Second-Harmonic Scattering (SHS) investigation led by two of us [Daniel et al., ACS Photonics, 2015, 2, 1209], we report herein a computational study of the structural organization and second-order nonlinear optical (NLO) properties of FONs based on dipolar chromophores incorporating a hydrophobic triphenylamine electron-donating unit and a slightly hydrophilic aldehyde electron-withdrawing unit at their extremities. Molecular dynamics simulations of the FON formation in water are associated with quantum chemical calculations, to provide insight into the molecular aggregation process, the molecular orientation of the dipolar dyes within the nanoparticles, and the dynamical behavior of their NLO properties. Moreover, the impact of intermolecular interactions on the NLO responses of the FONs is investigated by employing the tight-binding version of the recently developed simplified time-dependent density functional theory (sTD-DFT) approach, allowing the all-atom quantum mechanics treatment of nanoparticles.

Removing artifacts in polarizable embedding calculations of one- and two-photon absorption spectra of fluorescent proteins

The multiscale calculations involving excited states may suffer from the electron spill-out (ESO) problem. This seems to be especially the case when the environment of the core region, described with the electronic structure method, is approximated by a polarizable force field. The ESO effect often leads to incorrect physical character of electronic excitations, spreading outside the quantum region, which, in turn, results in erroneous absorption spectra. In this work, we investigate means to remove the artifacts in one-photon absorption (OPA) and two-photon absorption (TPA) spectra of green and yellow fluorescent protein representatives. This includes (i) using different basis sets, (ii) extending the core subsystem beyond the chromophore, (iii) modification of polarization interaction between the core region and its environment, and (iv) including the Pauli repulsion through effective core potentials (ECPs). Our results clearly show that ESO is observed when diffuse functions are used to assemble the multielectron wave function regardless of the exchange-correlation functional used. Furthermore, extending the core region, thus accounting for exchange interactions between the chromophore and its environment, leads to even more spurious excited states. Also, damping the interactions between the core subsystem and the polarizable force field is hardly helpful. In contrast, placing ECPs in the position of sites creating the embedding potential leads to the removal of artificious excited states that presumably should not be observed in the OPA and TPA spectra. We prove that it is a reliable and cost-effective approach for systems where the covalent bond(s) between the core region and its environment must be cut.

Damped (linear) response theory within the resolution-of-identity coupled cluster singles and approximate doubles (RI-CC2) method

An implementation of a complex solver for the solution of the linear equations required to compute the complex response functions of damped response theory is presented for the resolution-of-identity (RI) coupled cluster singles and approximate doubles (CC2) method. The implementation uses a partitioned formulation that avoids the storage of double excitation amplitudes to make it applicable to large molecules. The solver is the keystone element for the development of the damped coupled cluster response formalism for linear and nonlinear effects in resonant frequency regions at the RI-CC2 level of theory. Illustrative results are reported for the one-photon absorption cross section of C₆₀, the electronic circular dichroism of n-helicenes (n = 5, 6, 7), and the C₆ dispersion coefficients of a set of selected organic molecules and fullerenes.

Efficient implementation of isotropic cubic response functions for two-photon absorption cross sections within the self-consistent field approximation

Within the self-consistent field approximation, computationally tractable expressions for the isotropic second-order hyperpolarizability have been derived and implemented for the calculation of two-photon absorption cross sections. The novel tensor average formulation presented in this work allows for the evaluation of isotropic damped cubic response functions using only ∼3.3% (one-photon off-resonance regions) and ∼10% (one-photon resonance regions) of the number of auxiliary Fock matrices required when explicitly calculating all the needed individual tensor components. Numerical examples of the two-photon absorption cross section in the one-photon off-resonance and resonance regions are provided for alanine-tryptophan and 2,5-dibromo-1,4-bis(2-(4-diphenylaminophenyl)vinyl)-benzene. Furthermore, a benchmark set of 22 additional small- and medium-sized organic molecules is considered. In all these calculations, a quantitative assessment is made of the reduced and approximate forms of the cubic response function in the one-photon off-resonance regions and results demonstrate a relative error of less than ∼5% when using the reduced expression as compared to the full form of the isotropic cubic response function.

Core–valence-separated coupled-cluster-singles-and-doubles complex-polarization-propagator approach to X-ray spectroscopies

The iterative subspace algorithm to solve the CCSD complex linear response equations has been modified to include a core–valence separation projection step to overcome convergence problems. Illustrative results are reported for XAS, XCD, XES and RIXS.

CC2 oscillator strengths within the local framework for calculating excitation energies (LoFEx)

In a recent work [P. Baudin and K. Kristensen, J. Chem. Phys. 144, 224106 (2016)], we introduced a local framework for calculating excitation energies (LoFEx), based on second-order approximated coupled cluster (CC2) linear-response theory. LoFEx is a black-box method in which a reduced excitation orbital space (XOS) is optimized to provide coupled cluster (CC) excitation energies at a reduced computational cost. In this article, we present an extension of the LoFEx algorithm to the calculation of CC2 oscillator strengths. Two different strategies are suggested, in which the size of the XOS is determined based on the excitation energy or the oscillator strength of the targeted transitions. The two strategies are applied to a set of medium-sized organic molecules in order to assess both the accuracy and the computational cost of the methods. The results show that CC2 excitation energies and oscillator strengths can be calculated at a reduced computational cost, provided that the targeted transitions are local compared to the size of the molecule. To illustrate the potential of LoFEx for large molecules, both strategies have been successfully applied to the lowest transition of the bivalirudin molecule (4255 basis functions) and compared with time-dependent density functional theory.

A Polarization Propagator for Nonlinear X-ray Spectroscopies

A complex polarization propagator approach has been developed to third order and implemented in density functional theory (DFT), allowing for the direct calculation of nonlinear molecular properties in the X-ray wavelength regime without explicitly addressing the excited-state manifold. We demonstrate the utility of this propagator method for the modeling of coherent near-edge X-ray two-photon absorption using, as an example, DFT as the underlying electronic structure model. Results are compared with the corresponding near-edge X-ray absorption fine structure spectra, illuminating the differences in the role of symmetry, localization, and correlation between the two spectroscopies. The ramifications of this new technique for nonlinear X-ray research are briefly discussed.

Ab initio study of the enantio-selective magnetic-field-induced second harmonic generation in chiral molecules

We present a systematic ab initio study of enantio-selective magnetic-field-induced second harmonic generation (MFISHG) on a set of chiral systems ((l)-alanine, (l)-arginine and (l)-cysteine; 3,4-dehydro-(l)-proline; (S)-α-phellandrene; (R,S)- and (S,S)-cystine disulphide; N-(4-nitrophenyl)-(S)-prolinol, N-(4-(2-nitrovinyl)-phenyl)-(S)-prolinol, N-(4-tricyanovinyl-phenyl)-(S)-prolinol, (R)-BINOL, (S)-BINAM and 6-(M)-helicene). The needed electronic frequency dependent cubic response calculations are performed within a density functional theory (DFT) approach. A study of the dependence of the property on the choice of electron correlation, on one-electron basis set extension and on the choice of magnetic gauge origin is carried out on a prototype system (twisted oxygen peroxide). The magnetic gauge dependence analysis is extended also to the molecules of the set. An attempt to analyze the structure-property relationships is also made, based on the results obtained for biphenyl (in a frozen twisted conformation), for prolinol and for some of their derivatives. The strength of the effect is discussed, in order to establish its measurability with a proposed experimental setup.

4-Component relativistic calculations of L3ionization and excitations for the isoelectronic species UO22+, OUN+and UN2

4-Component relativistic calculations explore uranium 2p_3/2ionization and excitation in the isoelectronic series UO₂²⁺, OUN⁺and UN₂.

A complex-polarization-propagator protocol for magneto-chiral axial dichroism and birefringence dispersion

A schematic representation of magneto-chiral effects.

Excited states in large molecular systems through polarizable embedding

Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.

Two-photon absorption cross sections within equation-of-motion coupled-cluster formalism using resolution-of-the-identity and Cholesky decomposition representations: Theory, implementation, and benchmarks

The equation-of-motion coupled-cluster (EOM-CC) methods provide a robust description of electronically excited states and their properties. Here, we present a formalism for two-photon absorption (2PA) cross sections for the equation-of-motion for excitation energies CC with single and double substitutions (EOM-CC for electronically excited states with single and double substitutions) wave functions. Rather than the response theory formulation, we employ the expectation-value approach which is commonly used within EOM-CC, configuration interaction, and algebraic diagrammatic construction frameworks. In addition to canonical implementation, we also exploit resolution-of-the-identity (RI) and Cholesky decomposition (CD) for the electron-repulsion integrals to reduce memory requirements and to increase parallel efficiency. The new methods are benchmarked against the CCSD and CC3 response theories for several small molecules. We found that the expectation-value 2PA cross sections are within 5% from the quadratic response CCSD values. The RI and CD approximations lead to small errors relative to the canonical implementation (less than 4%) while affording computational savings. RI/CD successfully address the well-known issue of large basis set requirements for 2PA cross sections calculations. The capabilities of the new code are illustrated by calculations of the 2PA cross sections for model chromophores of the photoactive yellow and green fluorescent proteins.

Theoretical study of two-photon circular dichroism on molecular structures simulating aromatic amino acid residues in proteins with secondary structures

Calculation and comparative analysis of the theoretical two-photon circular dichroism (TPCD) spectra of l-His, l-Phe, and l-Tyr simulating residues in proteins with secondary structures (α-helix, β-strand and random coil), down to the far-UV region (FUV).

On the origin of the very strong two-photon activity of squaraine dyes – a standard/damped response theory study

In the present work, we report the mechanism of a very large increase in the two-photon (TP) activity of squaraine based molecules upon changing the substituents.

The Dalton quantum chemistry program system

Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

Comparison of standard and damped response formulations of magnetic circular dichroism

We apply damped response theory to the phenomenon of magnetic circular dichroism (MCD), and we investigate how the numerical instability associated with the simulation of the MCD spectrum from individually calculated A and B terms for close lying states can be remedied by the use of damped response theory. We also present a method for calculating the Faraday A term, formulated as a double residue of the quadratic response function.