Damped response theory description of two-photon absorption

Damped response theory is applied to the calculation of two-photon absorption (TPA) spectra, which are determined directly, at each frequency, from a modified damped cubic response function. The TPA spectrum may therefore be evaluated for selected frequency ranges, making the damped TPA approach attractive for calculations on large molecules with a high density of states, where the calculation of TPA using standard theory is more problematic. Damped response theory can also be applied to the case of intermediate state resonances, where the standard TPA expression is divergent. Both exact damped response theory and its application within density functional theory are discussed. The latter is implemented using an atomic-orbital based density matrix formulation, which makes the approach especially suitable for studies on large systems. A test preliminary study is presented for the TPA spectrum of R-(+)-1,1'-bi(2-naphtol).     

Calculation of two-photon absorption spectra of donor-pi-acceptor compounds in solution using quadratic response time-dependent density functional theory

Linear and quadratic response time-dependent density functional theories have been applied to calculate the photophysical properties of donor-pi-acceptor molecules which are known to have large nonlinear absorption. The linear absorption and two-photon absorption spectra predicted using hybrid functionals, including the Coulomb-attenuated model, with continuum solvation models are reported and compared to experiment and to previous theoretical predictions. While the quadratic response with these functionals overestimated the TPA cross sections relative to experiment when a Gaussian linewidth function was used, a fairly good agreement was obtained when a Lorentzian linewidth function was applied. In addition, the comparison of the TPA cross sections calculated by the sum over states with those calculated by the two-state approximation indicates the importance of the higher energy states in TPA, particularly in nondegenerate experiments.     

Calculations of two-photon charge-transfer excitations using Coulomb-attenuated density-functional theory

In this work, we show that an implementation of Coulomb-attenuated density-functional theory leads to considerably better prospects than hitherto for modeling two-photon absorption cross sections for charge-transfer species. This functional, which corrects for the effect of poor asymptotic dependence of commonly used functionals, essentially brings down the widely different results for larger charge-transfer species between Hartree-Fock and density-functional theory (DFT)-B3LYP into a closer range. The Coulomb-attenuated functional, which retains the best aspects of the Hartree-Fock and DFT-B3LYP methods, proves to be very promising for further modeling design of multiphoton materials with technical applications.     

Time-dependent density-functional theory calculations of triplet-triplet absorption

We present density-functional theory calculations of triplet-triplet absorption by three different approaches based on time-dependent density-functional theory (DFT): unrestricted DFT linear response, open-shell restricted DFT linear response applied to the triplet state, and quadratic response with triplet excitations applied to the ground state. Comparison is also made with corresponding results obtained by Hartree-Fock and multiconfiguration self-consistent-field response theory. Two main conclusions concerning triplet-triplet transitions are drawn in this study: First, the very good agreement between unrestricted and restricted DFT results indicates that spin contamination of the triplet state is not a serious problem when computing triplet-triplet spectra of common organic molecules. Second, DFT response calculations of triplet-triplet transitions can be affected by triplet instability problems, especially for the combination of DFT quadratic response with functionals containing fractional exact Hartree-Fock exchange.     

Real-time linear response for time-dependent density-functional theory

We present a linear-response approach for time-dependent density-functional theories using time-adiabatic functionals. The resulting theory can be performed both in the time and in the frequency domain. The derivation considers an impulsive perturbation after which the Kohn-Sham orbitals develop in time autonomously. The equation describing the evolution is not strictly linear in the wave function representation. Only after going into a symplectic real-spinor representation does the linearity make itself explicit. For performing the numerical integration of the resulting equations, yielding the linear response in time, we develop a modified Chebyshev expansion approach. The frequency domain is easily accessible as well by changing the coefficients of the Chebyshev polynomial, yielding the expansion of a formal symplectic Green's operator.     

Harmonic theory of thermal two-photon absorption in benzene

A correlation function formalism is applied to compute the two-photon absorption spectrum of benzene. Using harmonic Hamiltonians for the ground and excited electronic states, we find that the theory agrees qualitatively with the experimentally observed sparsity of the thermal two-photon absorption spectrum as compared with the single-photon absorption spectrum. An expression for the average vibrational energy in the excited state is derived. We find that cooling of the nascent vibrational energy in the electronically excited state is not as extensive in the two-photon absorption process as compared to the single-photon case.     

On the evaluation of quadratic response functions at the four-component Hartree-Fock level: nonlinear polarization and two-photon absorption in bromo- and iodobenzene

The nonlinear polarization and two-photon absorption parameters have been determined for dibromo- and di-iodobenzene in their meta- and ortho-conformations and with relativistic effects accounted for to a varying degree. By exclusion of small component integrals in the calculations of the first-order hyperpolarizability, results within 1% of fully relativistic four-component Hartree-Fock values are obtained at a cost of 8.7 times the corresponding nonrelativistic calculations. It is shown that the nonlinear absorption in bromobenzene (and even more so in iodobenzene) is broad banded due to spin-orbit interactions among the excited states, and nonrelativistic and scalar relativistic calculations are not to be used in this case.     

Calculation of two-photon absorption strengths with the approximate coupled cluster singles and doubles model CC2 using the resolution-of-identity approximation

An implementation of two-photon absorption matrix elements using the approximate second-order coupled-cluster singles and doubles model CC2 is presented. In this implementation we use the resolution-of-the-identity approximation for the two-electron repulsion integrals to reduce the computational cost. To avoid storage of large arrays we introduce in addition a numerical Laplace transformation of orbital energy denominators for the response of the doubles amplitudes. The error due to the numerical Laplace transformation is found to be negligible. Using this new implementation, we performed a series of benchmark calculations on substituted benzene and azobenzene derivatives to get reference values for TD-DFT results. We show that results obtained with the Coulomb-attenuated B3LYP functional are in reasonable agreement with the coupled-cluster results, whereas other density functionals which do not have a long-range correction give considerably less accurate results. Applications to the AF240 dye molecule and a weakly bound molecular tweezer complex demonstrate that this new RI-CC2 implementation allows for the first time to compute two-photon absorption cross sections with a correlated wave function method for molecules with more than 70 atoms and to apply this method for benchmarking TD-DFT calculations on molecules which are of particular relevance for experimental studies of two-photon absorption.     

Structure determination of small vanadium clusters by density-functional theory in comparison with experimental far-infrared spectra

The far-infrared vibrational spectra for charged vanadium clusters with sizes of 3-15 atoms have been measured using infrared multiple photon dissociation of Vn+Ar-->Vn(+)+Ar. Using density-functional theory calculations, we calculated the ground state energy and vibrational spectra for a large number of stable and metastable geometries of such clusters. Comparison of the calculated vibrational spectra with those obtained in the experiment allows us to deduce the cluster size specific atomic structures. In several cases, a unique atomic structure can be identified, while in other cases our calculations suggest the presence of multiple isomers.     

Two-photon absorption in solution by means of time-dependent density-functional theory and the polarizable continuum model

We present the first study of two-photon absorption (TPA) of solvated molecules based on direct evaluation of TPA cross sections from the quadratic response of time-dependent perturbations. A set of prototypical two-photon (TP) chromophores has been selected and analyzed: a pure pi system (t-stilbene) and its substituted homologs obtained employing a donor (D) and an acceptor (A) group to probe the solvent effects along the series pi, D-pi-D, A-pi-D, and A-pi-A. For the selected systems we have calculated the TPA cross sections in different environments by means of the polarizable continuum model. The data have been analyzed to evaluate how the structural and environmental parameters contribute to the final two-photon absorption cross section. These include molecular structure, geometry relaxation in solution, polarity, and refractive index of the solvent. The performances of the three common functionals SVWN, BLYP, and B3LYP have been compared. The results show a significant solvent dependence of the TPA cross section and an unusual trend when passing from cyclohexane to water. The data have also been rationalized in terms of the main orbital excitations leading to the transitions. Finally, trends along the series have been described and comparison with experiments and previous calculations has been drawn.     

Turbo charging time-dependent density-functional theory with Lanczos chains

We introduce a new implementation of time-dependent density-functional theory which allows the entire spectrum of a molecule or extended system to be computed with a numerical effort comparable to that of a single standard ground-state calculation. This method is particularly well suited for large systems and/or large basis sets, such as plane waves or real-space grids. By using a superoperator formulation of linearized time-dependent density-functional theory, we first represent the dynamical polarizability of an interacting-electron system as an off-diagonal matrix element of the resolvent of the Liouvillian superoperator. One-electron operators and density matrices are treated using a representation borrowed from time-independent density-functional perturbation theory, which permits us to avoid the calculation of unoccupied Kohn-Sham orbitals. The resolvent of the Liouvillian is evaluated through a newly developed algorithm based on the nonsymmetric Lanczos method. Each step of the Lanczos recursion essentially requires twice as many operations as a single step of the iterative diagonalization of the unperturbed Kohn-Sham Hamiltonian. Suitable extrapolation of the Lanczos coefficients allows for a dramatic reduction of the number of Lanczos steps necessary to obtain well converged spectra, bringing such number down to hundreds (or a few thousands, at worst) in typical plane-wave pseudopotential applications. The resulting numerical workload is only a few times larger than that needed by a ground-state Kohn-Sham calculation for a same system. Our method is demonstrated with the calculation of the spectra of benzene, C(60) fullerene, and of chlorophyll a.     

A global investigation of excited state surfaces within time-dependent density-functional response theory

This work investigates the capability of time-dependent density functional response theory to describe excited state potential energy surfaces of conjugated organic molecules. Applications to linear polyenes, aromatic systems, and the protonated Schiff base of retinal demonstrate the scope of currently used exchange-correlation functionals as local, adiabatic approximations to time-dependent Kohn-Sham theory. The results are compared to experimental and ab initio data of various kinds to attain a critical analysis of common problems concerning charge transfer and long range (nondynamic) correlation effects. This analysis goes beyond a local investigation of electronic properties and incorporates a global view of the excited state potential energy surfaces.     

Atomic volumes and polarizabilities in density-functional theory

Becke and Johnson introduced an ad hoc definition of atomic volume [J. Chem. Phys. 124, 014204 (2006)] in order to obtain atom-in-molecule polarizabilities from free-atom polarizabilities in their nonempirical exchange-hole dipole moment model of dispersion interactions. Here we explore the dependence of Becke-Johnson atomic volumes on basis sets and density-functional approximations and provide reference data for all atoms H-Lr. A persuasive theoretical foundation for the Becke-Johnson definition is also provided.     

A superfluorescent fluorenyl probe with efficient two-photon absorption

The linear photophysical, excited state absorption (ESA), superfluorescence, and two-photon absorption (2PA) properties of 4,4'-(1E,1'E)-2,2'-(7,7'(1E,1'E)2,2'(4,4'-sulfonylbis(4,1-phenylene))bis(ethane-2,1-diyl)bis(9,9-didecy-9H-fluorene7,2-diyl))bis(ethane-2,1-diyl)bis(N,N-diphenylaniline) (1) were investigated in organic and aqueous media with respect to its potential application in biological imaging. The analysis of linear photophysical properties revealed a rather complex nature of the main one-photon absorption band, strong solvatochromic effects in the steady-state fluorescence spectra, single-exponential fluorescence decay, and high fluorescence quantum yields in organic solvents (≈1.0). The ESA spectra of 1 suggested potential for light amplification in nonpolar media while efficient superfluorescence in cyclohexane was demonstrated. The degenerate 2PA spectra of 1 were obtained over a broad spectral range (640-900 nm), using a standard two-photon induced fluorescence method under 1 kHz femtosecond excitation. Two well defined 2PA bands with maximum 2PA cross sections up to 1700 GM in the higher energy, short wavelength band and ≈1200 GM in the lower energy, long wavelength band of 1 were shown. The potential use of 1 in bioimaging was demonstrated via one- and two-photon in vitro fluorescence imaging of HCT 116 cells.     

Extended squaraine dyes with large two-photon absorption cross-sections

Extended bis(donor)-substituted squaraine chromophores exhibit very high two-photon cross-sections (as high as 33 000 GM) in the near-IR; these can be attributed to the combination of large transition dipoles with small detuning energies. The modulus of the third-order nonlinear optical susceptibility at 1.3 mum has been found to be 7.0 x 10-11 esu for one of these chromophores.     

Analysis of environment response effects on excitation energies within subsystem-based time-dependent density-functional theory

We analyze the ability of subsystem time-dependent density-functional theory (sTDDFT) to describe environmental response effects. To this end, we utilize the recently proposed “exact” version of sTDDFT relying on projection-based embedding (PbE), which so far was applied only for the special case of two subsystems. We confirm that PbE-sTDDFT in combination with supersystem bases yields results equivalent to those of supermolecular TDDFT calculations for systems solvated by many solvent molecules, using the previously studied system of methylene-cyclopropene⋯(H₂O)₁₇ as an example. By means of this exact reference embedding framework, we are able to disentangle solvent effects introduced in terms of the embedding potential from those caused by solvent response couplings, both for the PbE variant and for sTDDFT with approximate non-additive kinetic energy functionals. Furthermore, we show that the use of a monomer basis introduces significant errors for the environmental response contribution. Employing a virtual-orbital localization strategy on top of PbE-sTDDFT, we can also directly assess the impact of inter-subsystem charge-transfer excitations on the entire solvent effect, which turn out to play a significant role for the environmental response. Finally, we analyze the response effects introduced by the individual solvent molecules and their interdependence, and show that a simple, pair-wise additive correction for solvent response yields excellent results in the present example.     

Calculation of nonadiabatic couplings with restricted open-shell Kohn-Sham density-functional theory

This paper generalizes the recently proposed approaches for calculating the derivative couplings between adiabatic states in density-functional theory (DFT) based on a Slater transition-state density to transitions such as singlet-singlet excitations, where a single-determinant ansatz is insufficient. The proposed approach is based on restricted open-shell Frank et al. [J. Chem. Phys. 108, 4060 (1998)] theory used to describe a spin-adapted Slater transition state. To treat the dependence of electron-electron interactions on the nuclear positions, variational linear-response density-functional perturbation theory is generalized to reference states with an orbital-dependent Kohn-Sham Hamiltonian and nontrivial occupation patterns. The methods proposed in this paper are not limited to the calculation of derivative coupling vectors, but can also be used for the calculation of other transition matrix elements. Moreover, they can be used to calculate the linear response of open-shell systems to arbitrary external perturbations in DFT.     

Time-dependent density functional theory studies of the photoswitching of the two-photon absorption spectra in stilbene, metacyclophenadiene, and diarylethene chromophores

Photochromophores such as cis-stilbene (1a), metacyclophenadiene (2a), and the diarylethene 3a undergo photoinduced conrotatory opening and closing of a central bond and are currently being sought out as potential candidates for media within 3D optical information storage devices. Strong molecular two-photon absorption (inducing the reversible photoisomerization) is a necessary feature for this application due to the need for high 3D spatial resolution. Here, the one- and two-photon absorption (OPA and TPA) characteristics of the open- and closed-ring isomers of 1-3 have been investigated using time-dependent density functional theory. It was determined that the excited states populated by two-photon absorption were nearly 1 eV higher in energy than the lowest energy excited state populated by one-photon absorption. The electronic structures of the TPA and OPA accessed states were then compared utilizing natural transition orbital analysis. There, it was found that states excited by OPA had pipi* character about the C-C framework associated with the bond formation/scission of the central C-C bond. In contrast, the states populated by TPA have pipi* character along the C-C skeletal periphery, including phenyl excitations. It is postulated that these differences in excited state electronic structure may lead to reaction pathways alternative to photoisomerization about the central C-C bond, impacting the utility of these compounds as 3D information storage media.