Damped response theory description of two-photon absorption

Non-linear light-matter interactions from the Bethe-Salpeter equation

A route to assess non-linear light-matter interactions from the increasingly popular GW-Bethe-Salpeter equation (GW-BSE) method is outlined. In the present work, the necessary analytic expressions within the static-screened exchange approximation of the BSE are derived. This enables a straightforward implementation of the computation of the first hyperpolarizability as well as two-photon absorption processes for molecular systems. Benchmark calculations on small molecular systems reveal that the GW-BSE method is intriguingly accurate for predicting both first hyperpolarizabilities and two-photon absorption strengths. Using state-of-the-art Kohn-Sham references as a starting point, the accuracy of the GW-BSE method rivals that of the coupled-cluster singles-and-doubles method, outperforming both second-order coupled-cluster and time-dependent density-functional theory.

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.

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.

In Silico Ultrafast Nonlinear Spectroscopy Meets Experiments: The Case of Perylene Bisimide Dye

Spectroscopy simulations are of paramount importance for the interpretation of experimental electronic spectra, the disentangling of overlapping spectral features, and the tracing of the microscopic origin of the observed signals. Linear and nonlinear simulations are based on the results drawn from electronic structure calculations that provide the necessary parameterization of the molecular systems probed by light. Here, we investigate the applicability of excited-state properties obtained from linear-response time-dependent density functional theory (TDDFT) in the description of nonlinear spectra by employing the pseudowavefunction approach and compare them with benchmarks from highly accurate RASSCF/RASPT2 calculations and with high temporal resolution experimental results. As a test case, we consider the prediction of femtosecond transient absorption and two-dimensional electronic spectroscopy of a perylene bisimide dye in solution. We find that experimental signals are well reproduced by both theoretical approaches, showing that the computationally cheaper TDDFT can be a suitable option for the simulation of nonlinear spectroscopy of molecular systems that are too large to be treated with higher-level RASSCF/RASPT2 methods.

Strong Coupling between Localized Surface Plasmons and Molecules by Coupled Cluster Theory

Plasmonic nanocavities enable the confinement of molecules and electromagnetic fields within nanometric volumes. As a consequence, the molecules experience a remarkably strong interaction with the electromagnetic field to such an extent that the quantum states of the system become hybrids between light and matter: polaritons. Here, we present a nonperturbative method to simulate the emerging properties of such polaritons: it combines a high-level quantum chemical description of the molecule with a quantized description of the localized surface plasmons in the nanocavity. We apply the method to molecules of realistic complexity in a typical plasmonic nanocavity, featuring also a subnanometric asperity (picocavity). Our results disclose the effects of the mutual polarization and correlation of plasmons and molecular excitations, disregarded so far. They also quantify to what extent the molecular charge density can be manipulated by nanocavities and stand as benchmarks to guide the development of methods for molecular polaritonics.

The orbital picture of the first dipole hyperpolarizability from many-body response theory

We present an approach for obtaining a molecular orbital picture of the first dipole hyperpolarizability (β) from correlated many-body electronic structure methods. Ab initio calculations of β rely on quadratic response theory, which recasts the sum-over-all-states expression of β into a closed-form expression by calculating a handful of first- and second-order response states; for resonantly enhanced β, damped response theory is used. These response states are then used to construct second-order response reduced one-particle density matrices (1PDMs), which, upon visualization in terms of natural orbitals (NOs), facilitate a rigorous and black-box mapping of the underlying electronic structure with β. We explain the interpretation of different components of the response 1PDMs and the corresponding NOs within both the undamped and damped response theory framework. We illustrate the utility of this new tool by deconstructing β for cis-difluoroethene, para-nitroaniline, and hemibonded OH^· + H₂O complex, computed within the framework of coupled-cluster singles and doubles response theory, in terms of the underlying response 1PDMs and NOs for a range of frequencies.

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.

How to stay out of trouble in RIXS calculations within equation-of-motion coupled-cluster damped response theory? Safe hitchhiking in the excitation manifold by means of core–valence separation

We present a novel approach with robust convergence of the response equations for computing resonant inelastic X-ray scattering (RIXS) cross sections within the equation-of-motion coupled-cluster (EOM-CC) framework.

One- and two-photon absorption properties of quadrupolar thiophene-based dyes with acceptors of varying strengths

The one-photon (1P) and two-photon (2P) absorption properties of three quadrupolar dyes, featuring thiophene as a donor and acceptors of varying strengths, are determined by a combination of experimental and computational methods employing the density functional theory (DFT). The emission shifts in different solvents are well reproduced by time-dependent DFT calculations with the linear response and state specific approaches in the framework of the polarizable continuum model. The calculations show that the energies of both 1P- and 2P-active states decrease with an increase of the strength of the acceptor. The 2P absorption cross-sections predicted by the response theory are accounted for by considering just one intermediate state (S1) in the sum-over-states formulation. For the chromophore featuring the stronger acceptor, the energetic positions of the 1P- and 2P-active states prevent the exploitation of the theoretically predicted very high 2P activity due to the competing 1P absorption into the S1 state.

Excited-state absorption in tetrapyridyl porphyrins: comparing real-time and quadratic-response time-dependent density functional theory

Three meso-substituted tetrapyridyl porphyrins (free base, Ni(ii), and Cu(ii)) were investigated for their optical limiting (OL) capabilities using real-time (RT-), linear-response (LR-), and quadratic-response (QR-) time-dependent density functional theory (TDDFT) methods. These species are experimentally known to display a prominent reverse saturable absorption feature between the Q and B bands of the ground-state absorption (GSA), which has been attributed to increased excited-state absorption (ESA) relative to GSA. A recently developed RT-TDDFT based method for calculating ESA from a LR-TDDFT density was utilized with eight exchange-correlation functionals (BLYP, PBE, B3LYP, CAM-B3LYP, PBE0, M06, BHLYP, and BHandH) and contrasted with calculations of ESA using QR-TDDFT with five exchange-correlation functionals (BLYP, B3LYP, CAM-B3LYP, BHLYP, and BHandH). This allowed for comparison between functionals with varying amounts of exact exchange as well as between the ability of RT-TDDFT and QR-TDDFT to reproduce OL behavior in porphyrin systems. The absorption peak positions and intensities for GSA and ESA are significantly impacted by the choice of DFT functional, with the most critical factor identified as the amount of exact exchange in the functional form. Calculating ESA with QR-TDDFT is found to be significantly more sensitive to the amount of exact exchange than GSA and ESA with RT-TDDFT, as well as GSA with LR-TDDFT. An analogous behavior is also demonstrated for the polycyclic aromatic hydrocarbon coronene. This is problematic when using the same approximate functional for calculation of both GSA and ESA, as the LR- and QR-TDDFT excitation energies will not have similar errors. Overall, the RT-TDDFT method with hybrid functionals reproduces the OL features for the porphyrin systems studied here and is a viable computational approach for efficient screening of molecular complexes for OL properties.

Open-ended formulation of self-consistent field response theory with the polarizable continuum model for solvation

The study of high-order absorption properties of molecules is a field of growing importance. Quantum-chemical studies can help design chromophores with desirable characteristics. Given that most experiments are performed in solution, it is important to devise a cost-effective strategy to include solvation effects in quantum-chemical studies of these properties. We here present an open-ended formulation of self-consistent field (SCF) response theory for a molecular solute coupled to a polarizable continuum model (PCM) description of the solvent. Our formulation relies on the open-ended, density matrix-based quasienergy formulation of SCF response theory of Thorvaldsen, et al., [J. Chem. Phys., 2008, 129, 214108] and the variational formulation of the PCM, as presented by Lipparini et al., [J. Chem. Phys., 2010, 133, 014106]. Within the PCM approach to solvation, the mutual solute-solvent polarization is represented by means of an apparent surface charge (ASC) spread over the molecular cavity defining the solute-solvent boundary. In the variational formulation, the ASC is an independent, variational degree of freedom. This allows us to formulate response theory for molecular solutes in the fixed-cavity approximation up to arbitrary order and with arbitrary perturbation operators. For electric dipole perturbations, pole and residue analyses of the response functions naturally lead to the identification of excitation energies and transition moments. We document the implementation of this approach in the Dalton program package using a recently developed open-ended response code and the PCMSolver libraries and present results for one-, two-, three-, four- and five-photon absorption processes of three small molecules in solution.

Importance of double-resonance effects in two-photon absorption properties of Au25(SR)18. Chem Sci 2017. [PMID: 28626569 PMCID: PMC5471455 DOI: 10.1039/c7sc00968b]

We show that double-resonance effects for Au₂₅(SR)₁₈^– are less pronounced and do not lead to significantly enhanced two-photon absorption cross-sections.

The two-photon absorption (TPA) cross-sections of small thiolate-protected gold clusters have been shown to be much larger than typical small organic molecules. In comparison with larger nanoparticles, their TPA cross-sections per gold atom are also found to be larger. Theoretical simulations have suggested that the large enhancement of these TPA cross-sections comes from a one-photon double-resonance mechanism. However, it remains difficult to simulate TPA cross-sections of thiolate-protected gold clusters due to their large system size and a high density of states. In this work, we report a time-dependent density functional theory (TDDFT) study of the TPA spectra of the Au₂₅(SR)₁₈^– cluster based on a damped response theory formalism. Damped response theory enables a consistent treatment of on- and off-resonance molecular properties even for molecules with a high density of states, and thus is well-suited for studying the TPA properties of gold clusters. Our results indicate that the one- and two-photon double-resonance effect is much smaller than previously found, and thus is unlikely to be the main cause of the large TPA cross-sections found experimentally. The effect of symmetry breaking of the Au₂₅(SR)₁₈^– cluster due to the ligands on the TPA cross-sections has been studied and was found to only slightly increase the cross-section. Furthermore, by comparing with larger nanoparticles we find that the TPA cross-section per gold atom scales linearly with the diameter of the particles, and that the Kerr non-linear response of the Au₂₅(SR)₁₈^– cluster is on the same order as that of bulk gold films.

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.

Benchmarking two-photon absorption cross sections: performance of CC2 and CAM-B3LYP

We investigate the performance of CC2 and TDDFT/CAM-B3LYP for the calculation of two-photon absorption (TPA) strengths and cross sections and contrast our results to a recent coupled cluster equation-of-motion (EOM-EE-CCSD) benchmark study [K. D. Nanda and A. I. Krylov, J. Chem. Phys., 2015, 142, 064118]. In particular, we investigate whether CC2 TPA strengths are significantly overestimated compared to higher-level coupled-cluster calculations for fluorescent protein chromophores. Our conclusion is that CC2 TPA strengths are only slightly overestimated compared to the reference EOM-EE-CCSD results and that previously published overestimated cross sections are a result of inconsistencies in the conversion of the TPA strengths to macroscopic units. TDDFT/CAM-B3LYP TPA strengths, on the other hand, are found to be 1.5 to 3 times smaller than the coupled-cluster reference for the molecular systems considered. The unsatisfactory performance of TDDFT/CAM-B3LYP can be linked to an underestimation of excited-state dipole moments predicted by TDDFT/CAM-B3LYP.

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.

Two-photon absorption and two-photon circular dichroism of hexahelicene derivatives: a study of the effect of the nature of intramolecular charge transfer

Effect of the strength and nature of the extension of the π-electronic delocalization on TPA and TPCD of helicenes derivatives.

Chemical control of channel interference in two-photon absorption processes

The two-photon absorption (TPA) process is the simplest and hence the most studied nonlinear optical phenomenon, and various aspects of this process have been explored in the past few decades, experimentally as well as theoretically. Previous investigations have shown that the two-photon (TP) activity of a molecular system can be tuned, and at present, performance-tailored TP active materials are easy to develop by monitoring factors such as length of conjugation, dimensionality of charge-transfer network, strength of donor-acceptor groups, polarity of solvents, self-aggregation, H-bonding, and micellar encapsulation to mention but a few. One of the most intriguing phenomena affecting the TP activity of a molecule is channel interference. The phrase "channel interference" implies that if the TP transition from one electronic state to another involves more than one optical pathway or channel, characterized by the corresponding transition dipole moment (TDM) vectors, the channels may interfere with each other depending upon the angles between the TDM vectors and hence can either increase (constructive interference) or decrease (destructive interference) the overall TP activity of a system to a significant extent. This phenomenon was first pointed out by Cronstrand, Luo, and Ågren [Chem. Phys. Lett. 2002, 352, 262-269] in two-dimensional systems (i.e., only involving two components of the transition moment vectors). For three-dimensional molecules, an extended version of this idea was required. In order to fill this gap, we developed a generalized model for describing and exploring channel interference, valid for systems of any dimensionality. We have in particular applied it to through-bond (TB) and through-space (TS) charge-transfer systems both in gas phase and in solvents with different polarities. In this Account, we will, in addition to briefly describing the concept of channel interference, discuss two key findings of our recent work: (1) how to control the channel interference by chemical means, and (2) the role of channel interference in the anomalous solvent dependence of certain TP chromophores. For example, we will show that simple structurally induced changes in certain dihedral angles of the well-known betaine dye (TB type) will help fine-tune the constructive channel interference and hence increase the overall TP activity of molecules with this general TP channel structure. Another intriguing result we will discuss is observed for a tweezer-trinitrofluorinone complex (TS type) where, on moving from polar to essentially nonpolar solvents, the nature of the channel interference switches from destructive to constructive, leading to a net abnormal solvent dependence of the TP activity of the system. The present Account highlights the usefulness of the channel interference effect and establishes it as a new and unique way of controlling the TP transition probability in different types of three-dimensional molecules.

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.