Next Article in Journal
Metabolomics Analysis Reveals the Accumulation Patterns of Flavonoids and Volatile Compounds in Camellia oleifera Petals with Different Color
Next Article in Special Issue
Density Functional Theory Investigation of Temperature-Dependent Properties of Cu-Nitrogen-Doped Graphene as a Cathode Material in Fuel Cell Applications
Previous Article in Journal
Arthproliferins A–D, Four New Sesterterpenes from the Mangrove-Sediment-Derived Fungus Arthrinium sp. SCSIO41221
Previous Article in Special Issue
On the Jahn–Teller Effect in Silver Complexes of Dimethyl Amino Phenyl Substituted Phthalocyanine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Time Extension of TAO-DFT

1
Department of Physics, National Taiwan University, Taipei 10617, Taiwan
2
Center for Theoretical Physics and Center for Quantum Science and Engineering, National Taiwan University, Taipei 10617, Taiwan
3
Physics Division, National Center for Theoretical Sciences, Taipei 10617, Taiwan
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(21), 7247; https://doi.org/10.3390/molecules28217247
Submission received: 27 September 2023 / Revised: 19 October 2023 / Accepted: 21 October 2023 / Published: 24 October 2023
(This article belongs to the Special Issue Multiconfigurational and DFT Methods Applied to Chemical Systems)

Abstract

:
Thermally assisted occupation density functional theory (TAO-DFT) has been an efficient electronic structure method for studying the ground-state properties of large electronic systems with multi-reference character over the past few years. To explore the time-dependent (TD) properties of electronic systems (e.g., subject to an intense laser pulse), in this work, we propose a real-time (RT) extension of TAO-DFT, denoted as RT-TAO-DFT. Moreover, we employ RT-TAO-DFT to study the high-order harmonic generation (HHG) spectra and related TD properties of molecular hydrogen H 2 at the equilibrium and stretched geometries, aligned along the polarization of an intense linearly polarized laser pulse. The TD properties obtained with RT-TAO-DFT are compared with those obtained with the widely used time-dependent Kohn–Sham (TDKS) method. In addition, issues related to the possible spin-symmetry breaking effects in the TD properties are discussed.

1. Introduction

Over the last thirty years, Kohn–Sham density functional theory (KS-DFT) [1] has been a popular electronic structure method for the ground-state (GS) properties of physical systems in the presence of static external potentials at zero electronic temperature ( θ e l = 0 ) due to its low computational cost and reasonable accuracy [2,3,4,5]. Conventional time-dependent density functional theory (TD-DFT) [6] (also called the time-dependent Kohn–Sham (TDKS) method, real-time TD-DFT (RT-TD-DFT), or real-time density functional theory (RT-DFT)), which is the time-dependent (TD) extension of KS-DFT, has been recently applied to explore the TD and excited-state properties of electronic systems under the influence of TD external potentials [7,8,9]. Recently, a frequency-domain formulation of linear-response TD-DFT (LR-TD-DFT) [10] has also been adopted to obtain excitation energies (i.e., limited to the weak-field perturbative regime), owing to its computational efficiency and reasonable accuracy [7,8,9]. Nevertheless, for the study of TD phenomena or excitation energies beyond the linear response, conventional TD-DFT [6], which involves propagating the TDKS equation in the time domain without any restriction to the TD external potentials, remains a promising method.
In KS-DFT [1], since the exact exchange-correlation (xc) energy functional E xc [ ρ ] , in terms of the GS density ρ ( r ) , has not been discovered, it remains necessary to adopt density functional approximations (DFAs) for E xc [ ρ ] to perform practical calculations [2,3,4,5]. The xc energy functionals based on the frequently adopted DFAs, such as the LDA (local density approximation) [11,12] and GGAs (generalized gradient approximations) [13], are computationally efficient for the study of large systems. However, the DFA xc energy functionals have a few intrinsic shortcomings [2,3,4,5] and can yield the following qualitative errors: the self-interaction error (SIE), non-covalent interaction error (NCIE), and static correlation error (SCE). Since conventional TD-DFT [6], which usually takes the GS of a physical system as the initial state and often employs the GS xc potential (i.e., the functional derivative of E xc [ ρ ] ) evaluated at the instantaneous density ρ ( r , t ) in the so-called adiabatic approximation [7,8,9], the qualitative errors of E xc [ ρ ] can also degrade the accuracy of conventional TD-DFT results [14,15,16,17].
These qualitative errors can generally be reduced with the modification of the DFA functionals. For example, the SIE can be reduced by mixing the Hartree–Fock (HF) exchange energy into the parent DFA functionals (commonly called hybrid functionals) [18,19,20]. The NCIE can be reduced by combining the parent DFA functionals with the dispersion energy correction (also known as dispersion-corrected functionals) [21,22] or with the second-order Møller–Plesset (MP2) correlation energy (often called double-hybrid functionals) [23]. The SCE can be reduced by incorporating a fully nonlocal correlation energy component, such as the RPA (random phase approximation) correlation energy [24,25], into the parent DFA functionals. Nonetheless, the DFA, dispersion-corrected, hybrid, and double-hybrid functionals fail to resolve the SCE problem, while the RPA and related functionals are very demanding in computational expense and hence are impractical for large systems.
To circumvent the SCE problem at low computational cost, thermally assisted occupation density functional theory (TAO-DFT) [26] (i.e., a density functional theory with fractional orbital occupations) has been recently developed. Note that TAO-DFT is an electronic structure method for the GS properties of physical systems at zero electronic temperature ( θ e l = 0 ), even though it adopts a reference system of noninteracting electrons at some fictitious temperature θ . The xc energy functionals developed in KS-DFT can also be used in TAO-DFT [26,27,28,29]. Nonetheless, in strong contrast to KS-DFT, TAO-DFT, even with the commonly used DFA, dispersion-corrected, and hybrid functionals, can approximately describe strong static correlation effects, especially when an appropriate value of θ is chosen [26,27,28]. Consequently, TAO-DFT is very promising for studying the GS properties of large systems with strong static correlation effects [30,31,32,33,34,35,36,37]. Other TAO-DFT extensions include the schemes that determine the system-independent [38] and system-dependent [39] values of θ , TAO-DFT-based ab initio molecular dynamics (for equilibrium thermodynamic and dynamical properties) [40], and TAO-DFT-based polarizable continuum model (for solvation effects) [41].
Within the framework of TAO-DFT, Yeh et al. have recently proposed a frequency-domain formulation of linear-response time-dependent TAO-DFT [42], denoted as TDTAO-DFT (or, more precisely, LR-TDTAO-DFT by its inherent linear-response (LR) nature), allowing excitation energy calculations in the frequency domain (i.e., using Casida’s formulation [10]). In TDTAO-DFT, the TD effective one-electron potential (see Equations (6) and (B6) of Ref. [42]) is defined with the TD pure state | Ψ TAO ( t ) of a noninteracting reference system (also see Appendix B1 of Ref. [42]). However, the TD density ρ ( r , t ) (see Equation (5) of Ref. [42]) in TDTAO-DFT is generally not associated with a TD noninteracting pure state | Ψ TAO ( t ) but associated with a TD noninteracting ensemble (which should be described by a TD density operator, as will be discussed later). For example, in TDTAO-DFT (with θ 0 ), at the initial time t 0 , the initial density ρ ( r , t 0 ) is simply the TAO-DFT GS density ρ ( r ) (see Equation (1) of Ref. [42]), which should be associated with a thermal ensemble [26] (i.e., not associated with a pure state) of noninteracting electrons at a nonvanishing fictitious temperature ( θ 0 ). Therefore, the underlying assumption of TDTAO-DFT (i.e., that the TD density ρ ( r , t ) is assumed to be associated with the TD pure state | Ψ TAO ( t ) of a noninteracting reference system) is generally incorrect, except only for the θ = 0 case (wherein TDTAO-DFT reduces to conventional TD-DFT [6] or, more precisely, LR-TD-DFT [10] by its inherent LR nature).
To resolve the aforementioned inconsistency of TDTAO-DFT (especially for θ 0 ) [42], in the work, we reformulate the TD extension of TAO-DFT by introducing a new reference system, consisting of an ensemble of noninteracting electrons moving in a TD local potential. This real-time (RT) extension of TAO-DFT is denoted as RT-TAO-DFT. Moreover, since the assumption of a weak perturbation is not required in RT-TAO-DFT, we also employ RT-TAO-DFT to study strong-field electron dynamics in molecules as well as high-order harmonic generation (HHG) [43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60].
The rest of this paper is organized as follows. In Section 2, we review TAO-DFT and discuss closely related electronic structure methods. The formulation of RT-TAO-DFT is presented in Section 3. In Section 4, we describe the details of our RT-TAO-DFT calculations for the HHG spectra and related TD properties of molecular hydrogen H 2 at the equilibrium and stretched geometries, aligned along the polarization of an intense linearly polarized laser pulse. The TD properties computed using RT-TAO-DFT are discussed and compared with the results of conventional TD-DFT [6]. Moreover, issues related to the possible spin-symmetry breaking effects in the TD properties are also discussed. Our conclusions are provided in Section 5.

2. Ground-State Theory: TAO-DFT

2.1. Overview of TAO-DFT

Consider a physical system of N interacting electrons moving in an external potential v ext ( r ) at zero electronic temperature ( θ e l = 0 ). In TAO-DFT [26], the GS density ρ ( r ) of the physical system is represented by the thermal equilibrium density of a reference system (called the thermally assisted occupation (TAO) reference system) of noninteracting electrons in the presence of a local potential v TAO ( r ) (called the TAO potential) at some fictitious temperature θ (i.e., the temperature of the TAO reference system). In other words, ρ ( r ) is represented by the TAO orbitals { ϕ j ( r ) } and TAO orbital occupation numbers (TOONs) { f j } (atomic units (a.u.) are adopted throughout this work):
ρ ( r ) = j f j | ϕ j ( r ) | 2 .
Here, f j is the occupation number of the j-th TAO orbital ϕ j ( r ) , given by the Fermi–Dirac (FD) distribution function
f j = { 1 + exp [ ( ϵ j μ ) / θ ] } 1 ,
where 0 f j 1 , ϵ j is the energy of the j-th TAO orbital ϕ j ( r ) and μ is the chemical potential chosen to conserve N (i.e., the number of electrons):
j { 1 + exp [ ( ϵ j μ ) / θ ] } 1 = N .
On the basis of the Hohenberg–Kohn (HK) theorems [61] for the physical system at θ e l = 0 and the Mermin theorems [62] for the TAO reference system at the fictitious temperature θ , a set of self-consistent equations (i.e., the TAO equations) that determine the TAO orbitals { ϕ j ( r ) } , the TAO orbital energies { ϵ j } , and hence the TOONs { f j } (see Equations (2) and (3)) and the GS density ρ ( r ) (see Equation (1)) are given by [26]
h ^ TAO ( r ) ϕ j ( r ) = ϵ j ϕ j ( r ) .
Here, h ^ TAO ( r ) is the TAO effective one-electron Hamiltonian:
h ^ TAO ( r ) = 1 2 r 2 + v TAO ( r ) ,
with the TAO potential (i.e., the TAO effective one-electron potential)
v TAO ( r ) = v ext ( r ) + δ E H [ ρ ] δ ρ ( r ) + δ E xc θ [ ρ ] δ ρ ( r ) = v ext ( r ) + v H ( r ) + v xc θ ( r ) ,
where v ext ( r ) is the external potential of the physical system, v H ( r ) = δ E H [ ρ ] δ ρ ( r ) = d r ρ ( r ) | r r | is the Hartree potential (i.e., the functional derivative of the Hartree energy functional E H [ ρ ] = 1 2 d r d r ρ ( r ) ρ ( r ) | r r | ), and v xc θ ( r ) = δ E xc θ [ ρ ] δ ρ ( r ) is the xcθ potential, which is the functional derivative of the xcθ energy functional E xc θ [ ρ ] = E xc [ ρ ] + E θ [ ρ ] , with E xc [ ρ ] being the xc energy functional (as defined in KS-DFT [1]) and E θ [ ρ ] being the θ -dependent energy functional (e.g., see Equation (14) of Ref. [26]).
In TAO-DFT [26], to obtain the GS density ρ ( r ) (i.e., represented by Equation (1) with the TAO orbitals { ϕ j ( r ) } and TOONs { f j } ) of the physical system, Equation (1) to (6) should be solved self-consistently. After the self-consistency is achieved, the GS energy E [ ρ ] of the physical system (at θ e l = 0 ) is given by
E [ ρ ] = d r ρ ( r ) v ext ( r ) + A s θ [ { f j , ϕ j } ] + E H [ ρ ] + E xc θ [ ρ ] ,
where the first term is the external potential energy, E H [ ρ ] and E xc θ [ ρ ] are the Hartree and xcθ energy functionals, respectively, and A s θ [ { f j , ϕ j } ] is the noninteracting kinetic free energy at the fictitious temperature θ :
A s θ [ { f j , ϕ j } ] = T s θ [ { f j , ϕ j } ] + E S θ [ { f j } ] ,
i.e., the sum of the kinetic energy
T s θ [ { f j , ϕ j } ] = 1 2 j f j d r ϕ j ( r ) r 2 ϕ j ( r )
and entropy contribution
E S θ [ { f j } ] = θ j [ f j ln ( f j ) + ( 1 f j ) ln ( 1 f j ) ]
of noninteracting electrons at the fictitious temperature θ , which can be exactly computed using the TAO orbitals { ϕ j ( r ) } and TOONs { f j } . Note that, for the special case of θ = 0 , E θ = 0 [ ρ ] = 0 and TAO-DFT (with E xc θ [ ρ ] ) [26] reduce to KS-DFT (with E xc [ ρ ] ) [1].

2.2. Density Representation in TAO-DFT

The GS density ρ ( r ) of a physical system is interacting v-representable (I-VR) as the exact ρ ( r ) belongs to a GS wavefunction of an interacting N-electron Hamiltonian for some external potential v ext ( r ) , which can be exactly computed using the full configuration interaction (FCI) method at the complete basis set limit [63]. Moreover, the exact ρ ( r ) can be represented by the natural orbitals { χ j ( r ) } and natural orbital occupation numbers (NOONs) { n j } [64]:
ρ FCI ( r ) = j n j | χ j ( r ) | 2 ,
where the NOONs { n j } satisfy the following conditions,
0 n j 1 , j n j = N .
Nevertheless, in KS-DFT [1], ρ ( r ) is assumed to be noninteracting pure-state v-representable (NI-PS-VR) as it belongs to a one-determinantal GS wavefunction of a noninteracting N-electron Hamiltonian (i.e., the Kohn–Sham (KS) Hamiltonian) for some local potential (i.e., the KS potential) [65,66,67]. Accordingly, in KS-DFT, ρ ( r ) is represented by the occupied KS orbitals { ϕ i KS ( r ) } :
ρ KS ( r ) = i = 1 N | ϕ i KS ( r ) | 2 .
As has been shown in a number of studies [66,67,68,69,70,71], there are some reasonable GS densities (e.g., the GS densities of some electronic systems with strong static correlation effects) that are not NI-PS-VR. Apparently, these GS densities cannot be obtained with KS-DFT even adopting the exact xc energy functional E xc [ ρ ] .
On the other hand, in TAO-DFT [26], ρ ( r ) (given by Equation (1)) is assumed to be noninteracting thermal ensemble v-representable (NI-TE-VR) as it belongs to a thermal ensemble of a reference system of noninteracting electrons in the presence of a local potential (i.e., the TAO potential) at some fictitious temperature θ . Consequently, in TAO-DFT, ρ ( r ) is represented by the TAO orbitals { ϕ j ( r ) } and TOONs { f j } :
ρ TAO ( r ) = j f j | ϕ j ( r ) | 2 .
where the TOONs { f j } (given by the FD distribution function) satisfy the following conditions,
0 f j 1 , j f j = N .
Owing to the similar expressions of the TAO-DFT GS density ρ TAO ( r ) (see Equation (14)) and the exact GS density ρ FCI ( r ) (see Equation (11)), the fictitious temperature θ in TAO-DFT can be so chosen that the distribution of TOONs is close to the distribution of the exact NOONs, which is closely related to the stability (i.e., the single-reference (SR)/multi-reference (MR) character) of the GS of an electronic system [26]. Accordingly, the exact GS density is more likely to be NI-TE-VR with this choice of θ . In contrast to KS-DFT (i.e., TAO-DFT with θ = 0 ), TAO-DFT has an extra degree of freedom in choosing the θ value to improve the GS density representability.

2.3. Approximate Energy Functionals and Fictitious Temperatures in TAO-DFT

Since the exact xcθ energy functional E xc θ [ ρ ] (i.e., one of the key ingredients in TAO-DFT), in terms of the GS density ρ ( r ) , has not been known, it remains necessary to employ DFAs for E xc θ [ ρ ] to perform practical calculations using TAO-DFT. Conventional DFAs, such as the LDA and GGAs, for E xc θ [ ρ ] (i.e., the DFA xcθ energy functional E xc θ DFA [ ρ ] ) can be adopted [26,27]. In addition to TAO-DFA (i.e., TAO-DFT with the DFA functional E xc θ DFA [ ρ ] ), TAO-DFT with the exact exchange [28] and related hybrid functionals [28,29] may also be employed.
For the GS of an electronic system, the fictitious temperature θ of a given energy functional in TAO-DFT should be so selected that the distribution of TOONs simulates the distribution of the exact NOONs. In this situation, the static correlation associated with the electronic GS can be properly captured by the entropy contribution (see Equation (10)) in TAO-DFT [26]. In other words, the optimal θ should be closely related to the SR/MR character of the electronic GS. For systems with electronic ground states possessing SR character (i.e., SR systems), all the NOONs should be close to either 0 (fully empty) or 1 (fully occupied), and, hence, the optimal θ values in TAO-DFT should be sufficiently small (but nonvanishing for real electronic systems [38]). However, for systems with electronic ground states possessing MR character (i.e., MR systems), the distributions of NOONs (and hence the optimal θ values) can be highly system-dependent. While it remains very challenging to devise a scheme that always yields the best θ of each system for a given energy functional in TAO-DFT, some progress has been achieved in recent years.
For a given energy functional in TAO-DFT, if the optimal θ values of electronic systems can be kept within a narrow range of values, it would be very useful to define an optimal system-independent θ value. Recently, TAO-DFT with the optimal system-independent θ scheme [26,27,28,38], which is as efficient as KS-DFT (i.e., TAO-DFT with θ = 0 ) in computational cost, can be comparable to KS-DFT in performance for various SR systems [26,27,28,38,41], and can outperform KS-DFT for several MR systems [26,27,28,30,31,32,33,34,35,36,37,38,40,41]. To improve the optimal system-independent θ scheme, a self-consistent scheme that determines the optimal θ values of electronic systems has been recently proposed [39].

2.4. Comparison of KS-DFT, TAO-DFT, and FT-DFT

Here, we compare three generally different electronic structure methods (see Table 1): KS-DFT [1], TAO-DFT [26], and FT-DFT (finite-temperature density functional theory, also called the Mermin–Kohn–Sham (MKS) method) [1,62], each of which employs a reference system of noninteracting electrons in the presence of a local potential at some fictitious temperature θ .
Both KS-DFT and TAO-DFT are electronic structure methods for the GS properties of physical systems at zero electronic temperature ( θ e l = 0 ). Note that θ θ e l = 0 is assumed in KS-DFT, while the fictitious temperature ( θ 0 ) can be different from the electronic temperature ( θ e l = 0 ) in TAO-DFT. Accordingly, the GS density ρ ( r ) of a physical system at θ e l = 0 is assumed to be NI-PS-VR in KS-DFT and NI-TE-VR in TAO-DFT. Moreover, in KS-DFT, the HK universal functional (i.e., the sum of the interacting kinetic energy and the electron–electron repulsion energy at θ e l = 0 ) [61], which is a functional of the GS density ρ ( r ) , is given by
F HK [ ρ ] = T s [ { ϕ i KS } ] + E H [ ρ ] + E xc [ ρ ] ,
where T s [ { ϕ i KS } ] (exactly computed using the occupied KS orbitals { ϕ i KS ( r ) } ) is the noninteracting kinetic energy at zero fictitious temperature ( θ = 0 ), and E xc [ ρ ] is the xc energy functional, which needs to be approximated for practical KS-DFT calculations. By contrast, in TAO-DFT, the HK universal functional [61], which is a functional of the GS density ρ ( r ) , is expressed as
F HK [ ρ ] = A s θ [ { f j , ϕ j } ] + E H [ ρ ] + E xc θ [ ρ ] ,
where A s θ [ { f j , ϕ j } ] (exactly computed using the TAO orbitals { ϕ j ( r ) } and TOONs { f j } ) is the noninteracting kinetic free energy at the fictitious temperature θ , and E xc θ [ ρ ] is the xcθ energy functional, which needs to be approximated for practical TAO-DFT calculations. Note that TAO-DFT (with θ = 0 ) reduces to KS-DFT.
On the other hand, FT-DFT (i.e., the MKS method) is an electronic structure method for the thermal equilibrium properties of physical systems at finite electronic temperatures ( θ e l 0 ), and θ θ e l is assumed in FT-DFT. Therefore, the thermal equilibrium density ρ θ e l ( r ) of a physical system at θ e l is assumed to be NI-TE-VR in FT-DFT. Moreover, in FT-DFT, the Mermin (M) universal functional (i.e., the sum of the interacting kinetic free energy and the electron–electron repulsion energy at θ e l ) [62], which is a θ e l -dependent functional of the thermal equilibrium density ρ θ e l ( r ) , is given by
F M θ e l [ ρ θ e l ] = A s θ e l [ { f k MKS , ϕ k MKS } ] + E H [ ρ θ e l ] + F xc θ e l [ ρ θ e l ] ,
where A s θ e l [ { f k MKS , ϕ k MKS } ] (exactly computed using the MKS orbitals { ϕ k MKS ( r ) } and MKS orbital occupation numbers { f k MKS } ) are the noninteracting kinetic free energy at the fictitious temperature θ θ e l , and F xc θ e l [ ρ θ e l ] is the xc free energy functional, which needs to be approximated for practical FT-DFT calculations. Note that FT-DFT (with θ e l = 0 ) reduces to KS-DFT.
Consequently, for the GS properties of physical systems at θ e l = 0 , FT-DFT reduces to KS-DFT, while TAO-DFT (with θ 0 ) can be very different from KS-DFT (especially for MR systems) [26,38,41].

2.5. TAO-DFT-Related Methods

2.5.1. TAO-DFT with E xc θ [ ρ ] E xc [ ρ ]

Here, we compare two approximate methods (see Table 2) that are closely related to TAO-DFT [26] and FT-DFT (i.e., the MKS method) [1,62].
As mentioned previously, TAO-DFT is an electronic structure method for the GS properties of physical systems at zero electronic temperature ( θ e l = 0 ). In TAO-DFT, the xcθ energy functional E xc θ [ ρ ] = E xc [ ρ ] + E θ [ ρ ] . At zero fictitious temperature ( θ = 0 ), E θ = 0 [ ρ ] = 0 and hence E xc θ [ ρ ] = E xc [ ρ ] . At a sufficiently small fictitious temperature ( θ 0 ), the magnitude of E θ 0 [ ρ ] should remain small compared to that of E xc [ ρ ] , and hence the approximation E xc θ [ ρ ] E xc [ ρ ] can be reasonably justified. Clearly, TAO-DFT with E xc θ [ ρ ] E xc [ ρ ] (also called TAO-DFT without E θ [ ρ ] ) is an approximate TAO-DFT method (good for θ 0 ), which may be adopted to describe the strong static correlation effects of some GS systems (wherever θ 0 can be an appropriate fictitious temperature) at θ e l = 0 .
On the other hand, FT-DFT is an electronic structure method for the thermal equilibrium properties of physical systems at finite electronic temperatures ( θ e l 0 ), wherein θ θ e l is assumed. In FT-DFT, at zero electronic temperature ( θ e l = 0 ), F xc θ e l = 0 [ ρ θ e l = 0 ] = E xc [ ρ ] . At a sufficiently small electronic temperature ( θ e l 0 ), the approximation F xc θ e l [ ρ θ e l ] E xc [ ρ θ e l ] can be reasonably justified. Apparently, FT-DFT with F xc θ e l [ ρ θ e l ] E xc [ ρ θ e l ] is an approximate FT-DFT method (good for θ θ e l 0 ), which may be used to study the temperature effects of thermal equilibrium systems at θ e l 0 [72,73].
According to their mathematical expressions, TAO-DFT with E xc θ [ ρ ] E xc [ ρ ] [26] is strikingly similar to FT-DFT with F xc θ e l [ ρ θ e l ] E xc [ ρ θ e l ] [1,62]. However, owing to their distinctly different physical meanings, one can easily distinguish the two approximate methods simply based on the electronic properties computed. For example, for the GS properties of physical systems at θ e l = 0 , FT-DFT with F xc θ e l [ ρ θ e l ] E xc [ ρ θ e l ] reduces to KS-DFT with E xc [ ρ ] , while TAO-DFT with E xc θ [ ρ ] E xc [ ρ ] at a nonvanishing fictitious temperature ( θ 0 ) can be very different from KS-DFT with E xc [ ρ ] (especially for MR systems) [26,38,41]. Therefore, a number of recent studies on the GS properties of physical systems at θ e l = 0 [74,75,76,77,78,79,80,81] have actually been performed using TAO-DFT with E xc θ [ ρ ] E xc [ ρ ] [26] rather than FT-DFT with F xc θ e l [ ρ θ e l ] E xc [ ρ θ e l ] (which, in fact, should reduce to KS-DFT with E xc [ ρ ] at θ e l = 0 ) [1,62].

2.5.2. KS-DFA with the rTAO Energy Correction

On the basis of the TAO-DFA (with some fictitious temperature θ ) energy expression [26,27], Yeh, Yang, and Hsu have recently proposed a post-KS energy correction, called the rTAO energy correction [82], which is a θ -dependent energy correction evaluated with the KS-DFA (i.e., KS-DFT with the DFA xc energy functional) orbitals. Owing to the post-KS nature, for clarity, we denote this method as the KS-DFA+rTAO method.
For a small fictitious temperature ( θ 40 mhartree), the KS-DFA+rTAO method has been shown to approximately reproduce the TAO-DFA results for some selected electronic properties [82]. This may in part be due to the post-KS nature (i.e., good for small θ ) and a very limited amount of test data. For a considerably large θ (e.g., the optimal θ of TAO-DFT with the exact exchange for the dissociation of molecular hydrogen H 2 [28,38]) or for other electronic properties (e.g., atomization energies), the results obtained with this method can be very different from those obtained with TAO-DFT.
More importantly, the GS density obtained with the KS-DFA+rTAO method (with any value of θ ) is the same as the KS-DFA GS density (i.e., NI-PS-VR). In other words, some reasonable non-NI-PS-VR GS densities (e.g., the GS densities of some MR systems) cannot be obtained with KS-DFA and the KS-DFA+rTAO method (with any value of θ ) [66,67,68,69,70,71].
In particular, whenever the spin-symmetry constraint [3,4,26,83] on the singlet GS density of an electronic system is violated with KS-DFA (which can commonly happen for MR systems), it must also be violated with the KS-DFA+rTAO method (with any value of θ ). In such a situation, the spin-unrestricted KS-DFA/KS-DFA+rTAO results can be very different from the corresponding spin-restricted KS-DFA/KS-DFA+rTAO results, yielding unphysical spin-symmetry breaking effects in the spin-unrestricted KS-DFA/KS-DFA+rTAO calculations. By contrast, the spin-symmetry breaking issues can be greatly resolved by TAO-DFA (with an appropriate θ ) [26,27,28,30,33,38,39,40,41], highlighting the significance of the GS density representation in TAO-DFT.

3. Real-Time Theory: RT-TAO-DFT

3.1. RT-TAO Equation

Consider a physical system of N interacting electrons moving in a TD external potential v ext ( r , t ) . The Hamiltonian operator of the physical system is given by
H ^ ( t ) = T ^ + V ^ ee + v ^ ext ( t ) ,
containing the operators of the kinetic energy
T ^ = 1 2 i = 1 N r i 2 ,
the electron–electron interaction
V ^ ee = i = 1 N j > i N 1 | r i r j | ,
and the TD external potential
v ^ ext ( t ) = i = 1 N v ext ( r i , t ) .
Let | Ψ ( t ) be the TD state of the physical system. For most TD cases of physical interest, in this work, the time propagation is assumed to start from the GS | Ψ GS (i.e., the lowest energy eigenstate of H ^ ( t 0 ) , which is a stationary state) of the unperturbed physical system at time t = t 0 , and | Ψ GS is assumed to be non-degenerate. Accordingly, the TD state | Ψ ( t ) of the physical system is a solution of the TD Schrödinger equation (TDSE):
i t | Ψ ( t ) = H ^ ( t ) | Ψ ( t ) ,
with the initial state
| Ψ ( t 0 ) = | Ψ GS .
From the TD state | Ψ ( t ) , the TD density ρ ( r , t ) of the physical system can be determined by
ρ ( r , t ) = Ψ ( t ) ρ ^ ( r ) Ψ ( t ) ,
where ρ ^ ( r ) = i = 1 N δ ( r r i ) is the number density operator. In particular, the initial density ρ ( r , t 0 ) of the physical system is given by the GS density ρ GS ( r ) of the unperturbed physical system:
ρ ( r , t 0 ) = ρ GS ( r ) = Ψ GS ρ ^ ( r ) Ψ GS .
According to the Runge–Gross (RG) theorems for the physical system (i.e., consisting of a TD pure state) [6], the TD state | Ψ ( t ) is a functional of the TD density ρ ( r , t ) (i.e., formally depending on the density ρ ( r , t ) at all previous times t t ) and the initial state | Ψ ( t 0 ) , i.e., | Ψ ( t ) = | Ψ [ ρ , | Ψ ( t 0 ) ] ( t ) . In this work, the initial state | Ψ ( t 0 ) = | Ψ GS is a functional of the initial density ρ ( r , t 0 ) = ρ GS ( r ) based on the HK theorems [61]. Since the dependence of initial state | Ψ ( t 0 ) is implicitly included in the TD density ρ ( r , t ) , for brevity, | Ψ [ ρ , | Ψ ( t 0 ) ] ( t ) is denoted as | Ψ [ ρ ] ( t ) hereafter. Now, we define the action functional of the physical system:
A [ ρ ] = t 0 t 1 d t Ψ [ ρ ] ( t ) | i t H ^ ( t ) | Ψ [ ρ ] ( t ) = t 0 t 1 d t Ψ [ ρ ] ( t ) | i t T ^ V ^ ee v ^ ext ( t ) | Ψ [ ρ ] ( t ) = B [ ρ ] t 0 t 1 d t d r ρ ( r , t ) v ext ( r , t ) ,
where B [ ρ ] is a universal functional of the TD density ρ ( r , t ) :
B [ ρ ] = t 0 t 1 d t Ψ [ ρ ] ( t ) | i t T ^ V ^ ee | Ψ [ ρ ] ( t ) .
Note that the action functional A [ ρ ] has a stationary point at the exact TD density ρ ( r , t ) of the physical system, given by the Euler equation:
δ A [ ρ ] δ ρ ( r , t ) = δ B [ ρ ] δ ρ ( r , t ) v ext ( r , t ) = 0 .
In order to develop an RT method compatible with TAO-DFT [26], we introduce the RT-TAO reference system, consisting of an ensemble of noninteracting electrons moving in a TD local potential v s ( r , t ) . The RT-TAO reference system can interchange electrons with its environment, and, hence, the electron number N e in the RT-TAO reference system can be varied from 0 to . The Hamiltonian operator of the RT-TAO reference system is given by
H ^ s ( t ) = T ^ s + v ^ s ( t ) ,
containing the operators of the kinetic energy T ^ s and the TD local potential v ^ s ( t ) .
At time t = t 0 , the time propagation starts from the grand canonical ensemble (i.e., a stationary ensemble) of the TAO reference system at some fictitious temperature θ , obtained with TAO-DFT [26]. Therefore, the initial density operator Γ ^ s ( t 0 ) of the RT-TAO reference system is given by the TAO density operator Γ ^ TAO (i.e., the grand canonical density operator [84] of the TAO reference system at the fictitious temperature θ ):
Γ ^ s ( t 0 ) = Γ ^ TAO = N e n w N e , n Φ N e , n 0 Φ N e , n 0
with the equilibrium statistical weights
w N e , n = exp [ ( E N e , n μ N e ) / θ ] N e n exp [ ( E N e , n μ N e ) / θ ] ,
satisfying the following conditions,
0 w N e , n 1 , N e n w N e , n = 1 .
Here, | Φ N e , n 0 denotes the n-th N e -electron eigenstate of H ^ s ( t 0 ) = H ^ TAO (i.e., the Hamiltonian operator of the TAO reference system) and E N e , n the corresponding energy eigenvalue.
For t > t 0 , the noninteracting ensemble of the RT-TAO reference system is generally non-stationary due to the presence of the TD local potential v s ( r , t ) . According to the TDSE
i t | Φ N e , n ( t ) = H ^ s ( t ) | Φ N e , n ( t ) ,
with the initial state | Φ N e , n ( t 0 ) = | Φ N e , n 0 , we know how the TD state | Φ N e , n ( t ) evolves in time. Therefore, the TD noninteracting ensemble of the RT-TAO reference system can be properly defined by the following TD density operator Γ ^ s ( t ) :
Γ ^ s ( t ) = N e n w N e , n Φ N e , n ( t ) Φ N e , n ( t ) ,
where the statistical weights w N e , n , which are assumed to be time-independent, are given by Equation (32), i.e., the initial statistical weights. Note that the TD density operator Γ ^ s ( t ) of the RT-TAO reference system is a solution of the Liouville–von Neumann equation:
i t Γ ^ s ( t ) = [ H ^ s ( t ) , Γ ^ s ( t ) ] ,
with the initial density operator Γ ^ s ( t 0 ) = Γ ^ TAO (given by Equation (31)). From the TD density operator Γ ^ s ( t ) , the TD density ρ s ( r , t ) of the RT-TAO reference system can be determined by
ρ s ( r , t ) = Tr { Γ ^ s ( t ) ρ ^ s ( r ) } = N e n w N e , n Φ N e , n ( t ) | ρ ^ s ( r ) | Φ N e , n ( t ) ,
where ρ ^ s ( r ) is the number density operator [84]. In particular, the initial density ρ s ( r , t 0 ) of the RT-TAO reference system is given by the thermal equilibrium density ρ TAO ( r ) (see Equation (14)) of the TAO reference system [26], which is the same as the GS density ρ GS ( r ) of the unperturbed physical system (note that ρ GS ( r ) is assumed to be NI-TE-VR with this θ ) and hence the initial density ρ ( r , t 0 ) of the physical system (see Equation (26)):
ρ s ( r , t 0 ) = ρ TAO ( r ) = Tr Γ ^ TAO ρ ^ s ( r ) = N e n w N e , n Φ N e , n 0 | ρ ^ s ( r ) | Φ N e , n 0 = ρ GS ( r ) = ρ ( r , t 0 ) .
Here, we seek v s ( r , t ) that yields the same solution ρ s ( r , t ) = ρ ( r , t ) for t t 0 . According to the Li–Tong (LT) theorems for the RT-TAO reference system (i.e., consisting of a TD noninteracting ensemble) [85], the TD density operator Γ ^ s ( t ) is a functional of the TD density ρ s ( r , t ) (i.e., formally depending on the density ρ s ( r , t ) at all previous times t t ) and the initial density operator Γ ^ s ( t 0 ) , i.e., Γ ^ s ( t ) = Γ ^ s [ ρ s , Γ ^ s ( t 0 ) ] ( t ) . In this work, the initial density operator Γ ^ s ( t 0 ) = Γ ^ TAO (see Equation (31)) is a functional of the initial density ρ s ( r , t 0 ) = ρ TAO ( r ) (see Equation (14)) based on the Mermin theorems [62]. Since the dependence of initial density operator Γ ^ s ( t 0 ) is implicitly included in the TD density ρ s ( r , t ) , for brevity, Γ ^ s [ ρ s , Γ ^ s ( t 0 ) ] ( t ) is denoted as Γ ^ s [ ρ s ] ( t ) hereafter. Now, we define the action functional of the RT-TAO reference system:
A s [ ρ s ] = t 0 t 1 d t Tr Γ ^ s [ ρ s ] ( t ) i t H ^ s ( t ) = t 0 t 1 d t Tr Γ ^ s [ ρ s ] ( t ) i t T ^ s v ^ s ( t ) = B s [ ρ s ] t 0 t 1 d t d r ρ s ( r , t ) v s ( r , t ) ,
where B s [ ρ s ] is a universal functional of the TD density ρ s ( r , t ) :
B s [ ρ s ] = t 0 t 1 d t Tr Γ ^ s [ ρ s ] ( t ) i t T ^ s
Note that the action functional A s [ ρ s ] has a stationary point at the exact TD density ρ s ( r , t ) of the RT-TAO reference system, given by the Euler equation:
δ A s [ ρ s ] δ ρ s ( r , t ) = δ B s [ ρ s ] δ ρ s ( r , t ) v s ( r , t ) = 0 .
In RT-TAO-DFT, the universal functional B [ ρ ] (given by Equation (28)) is partitioned as
B [ ρ ] = B s [ ρ ] A H [ ρ ] A xc θ [ ρ ] ,
where the universal functional B s [ ρ ] is given by Equation (40), A H [ ρ ] is the Hartree action functional:
A H [ ρ ] = 1 2 t 0 t 1 d t d r d r ρ ( r , t ) ρ ( r , t ) | r r | ,
and A xc θ [ ρ ] is the xcθ action functional:
A xc θ [ ρ ] B s [ ρ ] B [ ρ ] A H [ ρ ] ,
which is a universal functional of the TD density ρ ( r , t ) . Applying Equation (42) to Equation (29), we obtain
δ B s [ ρ ] δ ρ ( r , t ) δ A H [ ρ ] δ ρ ( r , t ) δ A xc θ [ ρ ] δ ρ ( r , t ) v ext ( r , t ) = 0 .
Comparing Equation (41) with Equation (45) shows that the same solution ρ s ( r , t ) = ρ ( r , t ) can be obtained if we choose the TD effective one-electron potential v s ( r , t ) (up to a purely TD function C ( t ) ) of the RT-TAO reference system as
v s ( r , t ) = v ext ( r , t ) + δ A H [ ρ ] δ ρ ( r , t ) + δ A xc θ [ ρ ] δ ρ ( r , t ) = v ext ( r , t ) + v H ( r , t ) + v xc θ ( r , t ) ,
where v ext ( r , t ) is the TD external potential of the physical system, v H ( r , t ) = δ A H [ ρ ] δ ρ ( r , t ) = d r ρ ( r , t ) | r r | is the TD Hartree potential, and v xc θ ( r , t ) = δ A xc θ [ ρ ] δ ρ ( r , t ) is the TD xcθ potential.
Since the RT-TAO reference system consists of a TD noninteracting ensemble, the Hamiltonian H ^ s ( t ) (see Equation (30)) is separable, with the RT-TAO effective one-electron Hamiltonian
h ^ s ( r , t ) = 1 2 r 2 + v s ( r , t ) ,
where the TD effective one-electron potential v s ( r , t ) (denoted as the RT-TAO potential) is given by Equation (46). The RT-TAO orbitals { ϕ j ( r , t ) } evolve in time according to the effective one-electron TDSE (denoted as the RT-TAO equation):
i t ϕ j ( r , t ) = h ^ s ( r , t ) ϕ j ( r , t ) ,
with
ϕ j ( r , t 0 ) = ϕ j 0 ( r ) ,
i.e., the initial j-th RT-TAO orbital ϕ j ( r , t 0 ) is given by the j-th TAO orbital ϕ j 0 ( r ) (the j-th energy eigenfunction of h ^ s ( r , t 0 ) = h ^ TAO ( r ) (see Equation (5)), associated with the GS of the unperturbed physical system) [26]. The TD density ρ ( r , t ) of the physical system, which is the same as the TD density ρ s ( r , t ) (given by Equation (37)) of the RT-TAO reference system, can be computed using [84]
ρ ( r , t ) = ρ s ( r , t ) = j f j | ϕ j ( r , t ) | 2 ,
where the occupation number f j of the j-th RT-TAO orbital ϕ j ( r , t ) , which is time-independent, is given by Equation (2), i.e., its initial occupation number [26], also satisfying the conditions 0 f j 1 and j f j = N .
For the special case of θ = 0 , RT-TAO-DFT (with the xcθ action functional A xc θ [ ρ ] ) reduces to conventional TD-DFT (with the xc action functional A xc [ ρ ] ) [6], providing that, at time t = t 0 , the initial state of the physical system is the non-degenerate GS of the unperturbed physical system.
Here, we discuss the representation of the TD density ρ ( r , t ) of a physical system in conventional TD-DFT [6] and RT-TAO-DFT. In conventional TD-DFT, ρ ( r , t ) is assumed to be TD noninteracting pure-state v-representable (TD-NI-PS-VR) as it belongs to a TD one-determinantal wavefunction of a noninteracting N-electron Hamiltonian for some TD local potential. By contrast, in RT-TAO-DFT, ρ ( r , t ) (given by Equation (50)) is assumed to be TD noninteracting ensemble v-representable (TD-NI-E-VR) as it belongs to a TD noninteracting ensemble (described by a TD density operator; e.g., see Equation (35)) in the presence of a TD local potential (i.e., the RT-TAO potential).
In RT-TAO-DFT, since we specify the initial state | Ψ ( t 0 ) = | Ψ GS of the physical system and the initial density operator Γ ^ s ( t 0 ) = Γ ^ TAO of the RT-TAO reference system, the two conditions (i) the same initial density ρ ( r , t 0 ) = ρ s ( r , t 0 ) and (ii) the same initial time derivative of the density t ρ ( r , t ) | t = t 0 = t ρ s ( r , t ) | t = t 0 = 0 can be satisfied for the physical and RT-TAO reference systems, providing that the GS density ρ GS ( r ) of the unperturbed physical system is NI-TE-VR with a given value of θ (see Equation (38)). Note that conditions (i) and (ii), which ensure the existence of TD-NI-PS-VR densities [86], may also be the conditions for the existence of TD-NI-E-VR densities [87].
In particular, condition (i) highlights the significance of the initial density representability or the GS density representability (for most TD cases of physical interest, the initial state is chosen as the GS of the unperturbed physical system). For an MR system, condition (i) can be violated with conventional TD-DFT since the corresponding GS theory, i.e., KS-DFT, can suffer from the aforementioned issues related to the GS density representability [66,67,68,69,70,71] and the spin-symmetry constraint [3,4,26,83]. By contrast, these issues can be greatly resolved by TAO-DFT (i.e., the underlying GS theory of RT-TAO-DFT) [26,27,28,30,33,38,39,40,41] when an appropriate θ is chosen.

3.2. Matrix Representation

In RT-TAO-DFT, the j-th RT-TAO orbital ϕ j ( r , t ) can be expanded in the orthonormal one-electron basis { ϕ p 0 ( r ) } , spanned by the GS TAO orbitals (i.e., the TAO orbitals associated with the GS of the unperturbed physical system) [26]:
ϕ j ( r , t ) = p C p j ( t ) ϕ p 0 ( r ) ,
where { C p j ( t ) } are the TD expansion coefficients. Accordingly, the TD density ρ ( r , t ) (see Equation (50)) can be expressed as
ρ ( r , t ) = p q P p q ( t ) ϕ p 0 ( r ) ϕ q 0 ( r ) ,
where P ( t ) is the one-electron density matrix at time t, with matrix elements
P p q ( t ) = j f j C p j ( t ) C q j ( t ) .
Moreover, at time t, F ( t ) is the RT-TAO matrix (commonly known as the Fock matrix), which is the matrix representation of the RT-TAO effective one-electron Hamiltonian h ^ s ( r , t ) (see Equation (47)), with matrix elements
F p q ( t ) = d r ϕ p 0 ( r ) h ^ s ( r , t ) ϕ q 0 ( r ) .
In the orthonormal one-electron basis { ϕ p 0 ( r ) } , the RT-TAO equation (given by Equation (48)) can be reformulated in terms of the TD one-electron density matrix P ( t ) [8]:
i d d t P ( t ) = [ F ( t ) , P ( t ) ] .
As the time propagation is assumed to start from the GS of the unperturbed physical system at time t = 0 (without loss of generality, the initial time t 0 0 is assigned hereafter), the initial one-electron density matrix is given by
P p q ( 0 ) = f p δ p q ,
and the initial RT-TAO matrix is given by
F p q ( 0 ) = ϵ p δ p q ,
where f p and ϵ p are the occupation number and energy, respectively, of the p-th GS TAO orbital ϕ p 0 ( r ) [26].
The formal solution of RT-TAO equation (see Equation (55)) for the TD one-electron density matrix P ( t ) is given by [60,88,89]
P ( t ) = U ( t , 0 ) P ( 0 ) U ( t , 0 ) .
Here, U ( t b , t a ) is a unitary time propagator from t a to t b :
U ( t b , t a ) = T ^ exp i t a t b d t F ( t ) ,
where T ^ denotes time ordering. However, since F ( t a ) does not necessarily commute with F ( t b ) for t a t b (see Equation (59)), it remains challenging to obtain P ( t ) directly from P ( 0 ) (see Equation (58)) for a long time interval [ 0 , t ] in RT-TAO-DFT.
In practical calculations, to reduce the error in the propagation for a long time interval [ 0 , t ] , U ( t , 0 ) is commonly split into a product of multiple time propagators, each corresponding to a small time step Δt:
U ( t , 0 ) = n = 0 m 1 U ( t n + 1 , t n ) ,
where t n = n Δ t denotes the value of t at the n-th time step, noting that t 0 = 0 and t m = m Δ t = t . U ( t n + 1 , t n ) is the time propagator from t n to t n + 1 = t n + Δ t , given by
U ( t n + 1 , t n ) = T ^ exp i t n t n + Δ t d t F ( t ) ,
which takes P ( t n ) to P ( t n + 1 ) :
P ( t n + 1 ) = U ( t n + 1 , t n ) P ( t n ) U ( t n + 1 , t n ) .
We denote U n = U ( t n + 1 , t n ) and P n = P ( t n ) for brevity and apply Equation (60) to Equation (58). Accordingly, the density matrix P m = P ( t m ) = P ( t ) can be obtained from the initial density matrix P 0 = P ( t 0 ) = P ( 0 ) via the following expression:
P m = n = 0 m 1 U n P 0 n = 0 m 1 U n = U m 1 U m 2 U 1 U 0 P 0 U 0 U 1 U m 2 U m 1 .
For a sufficiently small time step Δt, F ( t n ) remains nearly commutative with F ( t n + Δ t ) , and, hence, U n = U ( t n + 1 , t n ) (see Equation (61)) can be computed without considering the time ordering. Note, however, that the exact time-ordered propagator can only be obtained in the limit of an infinitesimal time step (i.e., Δt→0).
Recently, various algorithms [60,88,89] have been developed for the numerical construction of the time propagation of TDKS equation [6], which may also be adopted for the time propagation of RT-TAO equation.
In short, it takes the following key steps to run an RT-TAO-DFT calculation for describing the time evolution of the electron density following a perturbation:
  • Construct the initial one-electron density matrix P ( 0 ) (see Equation (56)) and the initial RT-TAO matrix F ( 0 ) (see Equation (57)) for the GS of the unperturbed physical system at time t = 0 using TAO-DFT (i.e., the respective GS theory).
  • Apply the TD field to the physical system for t > 0 , and propagate the one-electron density matrix P ( t ) and the RT-TAO matrix F ( t ) in the time domain, according to the RT-TAO equation (given by the matrix representation, e.g., see Equation (55)).
  • Post-process the resulting TD observables (electron density, dipole moment, etc.).

4. HHG Spectra from RT-TAO-DFT

HHG from an electronic system (e.g., an atom or molecule) is a nonlinear optical process driven by an intense laser pulse, wherein the laser frequency can be converted into its integer multiples [43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60]. HHG has recently attracted much attention since it can be used to explore the structure and dynamics of electronic systems and chemical reactions on a femtosecond timescale. In addition, HHG can be employed to generate attosecond pulse trains as well as individual attosecond pulses [45,46].
HHG can be qualitatively described by the semiclassical three-step model [43,44]. First, an electron tunnels out from an electronic system in an intense laser field (i.e., tunnel ionization). Second, the electron is driven away from or back to the parent ion by the laser field. Finally, the electron recombines with the parent ion, emitting a high-energy photon.
In this work, we perform RT-TAO-DFT calculations to explicitly obtain the HHG spectra and related TD properties of molecular hydrogen H 2 at the equilibrium and stretched geometries:
  • H 2 with an equilibrium bond length of 1.45 bohr (≈0.767 Å).
  • H 2 with a stretched bond length of 3.78 bohr (≈2.00 Å).
Here, the nuclei of H 2 are positioned along the z-axis (i.e., the laser polarization) with the center of mass being located at the origin.
To obtain the HHG spectrum, the electronic system H 2 , which starts from the GS at time t = 0 , experiences an intense laser pulse for t > 0 . In order to mimic the commonly used Ti:sapphire laser [90], the strong-field interaction is generated by a laser pulse with an oscillating electric field linearly polarized along the z-axis (see Figure 1):
v laser ( r , t ) = z A 0 cos 2 π 2 σ p ( t σ p ) cos [ ω 0 ( t σ p ) ] .
Here, the interaction with the electric field is treated in the dipole approximation and the length gauge [91]. The electric-field amplitude of the laser pulse A 0 = 0.0534 a.u. (corresponding to the peak intensity I 0 1 × 10 14 W/cm2), the laser frequency (also called the fundamental frequency) ω 0 = 1.5498 eV (corresponding to the wavelength λ 0 800 nm), and σ p = 500 a.u. (≈12.1 fs) are adopted.
Figure 1. Electric field of the laser pulse adopted.
Figure 1. Electric field of the laser pulse adopted.
Molecules 28 07247 g001
In the HHG process, the electron released by tunnel ionization can travel far away from the center of the electronic system H 2 . To capture strong-field ionization process and to remove artificial reflections induced by the finite extent of Gaussian basis set (which will be adopted to describe the TAO/RT-TAO orbitals), a complex absorbing potential (CAP) [60], i v CAP ( r ) , is also employed for t > 0 . For an electronic system consisting of N A atoms, the CAP function v CAP ( r ) is defined as the minimum of the values of the atom-centered spherical absorbing potentials:
v CAP ( r ) = min { g 1 ( r ) , . . . , g N A ( r ) } .
Here, g I ( r ) is a spherical absorbing potential around the I-th nucleus:
g I ( r ) = 0 , f o r | r R I | < r 0 η ( | r R I | r 0 ) 2 , f o r r 0 | r R I | < r 0 + V max / η V max , f o r r 0 + V max / η | r R I |
for I = 1, …, N A , where R I is the position of the I-th nucleus. The cutoff radius r 0 should be small enough to interact with the electron density (because the space extended by the Gaussian basis set is finite), while it should be large enough to not overly perturb the original electronic system. Here, we adopt the cutoff radius r 0 = 9.524 bohr (≈5.040 Å), the curvature η = 4.0 hartree/bohr2, and the maximum potential value Vmax = 10 hartree.
To sum up, in the present HHG study, the RT-TAO effective one-electron Hamiltonian h ^ s ( r , t ) (see Equation (47)) is given by
h ^ s ( r , t ) = h ^ s 0 ( r , t ) , f o r t 0 h ^ s 0 ( r , t ) + v laser ( r , t ) i v CAP ( r ) , f o r t > 0
where h ^ s 0 ( r , t ) is the field-free RT-TAO effective one-electron Hamiltonian, v laser ( r , t ) (see Equation (64)) is the strong-field interaction, and i v CAP ( r ) (see Equation (65)) is the CAP.
To propagate the one-electron density matrix P ( t ) in the time domain, we adopt a time step of Δt = 0.02 a.u. (≈ 0.484 as) and a total propagation time of τ = 1000 a.u. (≈24.1 fs), which corresponds to a total of 5 × 10 4 time steps. The modified mid-point unitary transformation (MMUT) algorithm [88,89] is employed for the numerical time propagation of P ( t ) . As the use of the CAP breaks the conservation of the norm of the RT-TAO orbitals, the time propagation is no longer unitary [60].
Following sufficient time propagation, the TD one-electron density matrix P ( t ) is determined, and the TD density ρ ( r , t ) is given by Equation (52), yielding various TD properties. While the RT-TAO orbital occupation numbers, which are the same as the GS TOONs { f j } , are time-independent, the norm of the RT-TAO orbitals can decrease with time t due to the CAP. To describe electron ionization, the number of bound electrons is computed using
N b ( t ) = d r ρ ( r , t ) = d r j f j | ϕ j ( r , t ) | 2 = Tr { P ( t ) } .
Moreover, the induced dipole moment along the laser polarization (i.e., the z-axis) is calculated by
μ ( t ) = d r z ρ ( r , t ) .
Accordingly, the HHG spectrum can be computed using [9,92]
H ( ω ) = 1 2 π 0 τ d t w H ( t ) d 2 μ ( t ) d t 2 e i ω t 2 ,
where the HHG spectrum has been smoothed using the Hamming window function
w H ( t ) = 0.54 0.46 cos ( 2 π t τ )
to reduce the numerical noise. In the HHG spectrum, the harmonic order is defined as ω / ω 0 , with ω 0 being the fundamental frequency (see Equation (64)).
Here, we present the approximations made in the RT-TAO-DFT calculations, the computational details, and the numerical results. As the exact TD xcθ potential v xc θ ( r , t ) = δ A xc θ [ ρ ] δ ρ ( r , t ) (see Equation (46)) remains unknown, approximations to v xc θ ( r , t ) are necessary for practical RT-TAO-DFT calculations. While the exact v xc θ ( r , t ) formally depends on the density ρ ( r , t ) at all previous times t t , in this study, we make the adiabatic approximation:
v xc θ ( r , t ) δ E xc θ [ ρ ] δ ρ ( r ) | ρ ( r ) = ρ ( r , t ) ,
where the TD xcθ potential v xc θ ( r , t ) is approximated by the GS xcθ potential δ E xc θ [ ρ ] δ ρ ( r ) (see Equation (6)) evaluated at the instantaneous density ρ ( r , t ) . In the adiabatic approximation, since the exact xcθ energy functional E xc θ [ ρ ] also remains unknown, a DFA to E xc θ [ ρ ] should be made as well. In this work, we adopt the LDA (i.e., the simplest DFA) xcθ energy functional E xc θ LDA [ ρ ] = E xc LDA [ ρ ] + E θ LDA [ ρ ] , with E xc LDA [ ρ ] being the LDA xc energy functional [11,12] and E θ LDA [ ρ ] being the LDA θ -dependent energy functional [26]. For brevity, RT-TAO-DFT with the adiabatic LDA xcθ potential is denoted as RT-TAO-ALDA. At time t = 0 , the initial state (i.e., the GS of the unperturbed H 2 with a given bond length) is obtained with the underlying GS theory, TAO-LDA (i.e., TAO-DFT with the LDA xcθ energy functional E xc θ LDA [ ρ ] ) [26]. Note that RT-TAO-ALDA (with θ = 0 ) corresponds to TD-ALDA (i.e., conventional TD-DFT with the adiabatic LDA xc potential) [7,8,9], and its underlying GS theory, TAO-LDA (with θ = 0 ), corresponds to KS-LDA (i.e., KS-DFT with the LDA xc energy functional E xc LDA [ ρ ] ) [11,12].
In this study, we also investigate the possible spin-symmetry breaking effects in the TD properties, in analogy to the GS counterparts [3,4,26,83]. At time t = 0 , the initial state (i.e., the GS of the unperturbed H 2 with a given bond length) is a singlet state, and, hence, the initial up-spin and down-spin densities obtained with an exact theory must be the same based on the spin-symmetry constraint [3,4,26,83]. Moreover, for t > 0 , since the TD external potential adopted is spin-independent (see Equation (67)), the up-spin and down-spin densities, which are equally propagated in the time domain, must be the same at any subsequent time t [93]. Therefore, the TD properties (which depend on the TD spin densities) of H 2 obtained with the spin-unrestricted formalism must be identical to those obtained with the spin-restricted formalism.
To examine whether this spin-symmetry constraint can be satisfied by RT-TAO-ALDA, we perform spin-restricted and spin-unrestricted RT-TAO-ALDA (with θ = 0, 7, 20, and 40 mhartree) calculations for the TD properties, such as the number of bound electrons, induced dipole moment, and HHG spectrum, of H 2 at the equilibrium and stretched geometries (aligned along the polarization of an intense linearly polarized laser pulse) using the d-aug-cc-pVTZ basis set and a high-quality grid EML(99,590) containing 99 Euler–Maclaurin radial grid points and 590 Lebedev angular grid points. We note that the choice of basis set can significantly affect the HHG spectrum [94]. For the special case of θ = 0 , RT-TAO-ALDA reduces to TD-ALDA. All numerical results are obtained with a development version of Q-Chem 5.4 [95].
Since the GS of the unperturbed H 2 with an equilibrium bond length of 1.45 bohr exhibits mainly SR character, the spin-symmetry constraint can be satisfied by spin-unrestricted TAO-LDA (with θ = 0, 7, 20, and 40 mhartree) [26], producing the same up-spin and down-spin densities at time t = 0 . In addition, for t > 0 , owing to the use of a TD spin-independent external potential, the up-spin and down-spin densities, which are equally propagated in the time domain, should be the same at any subsequent time t [93]. Therefore, the TD properties, such as the number of bound electrons (see Figure 2), induced dipole moment (see Figure 3), and HHG spectrum (see Figure 4), of H 2 with an equilibrium bond length of 1.45 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with θ = 0, 7, 20, and 40 mhartree), are essentially the same (i.e., within the numerical precision considered).
On the other hand, the GS of the unperturbed H 2 with a stretched bond length of 3.78 bohr exhibits a noticeable MR character [26], and, hence, the spin-symmetry constraint is violated with spin-unrestricted KS-LDA (i.e., TAO-LDA with θ = 0 ), producing symmetry-broken spin densities at time t = 0 . In this situation, even when a TD spin-independent external potential is applied to the initial state (i.e., a spin-symmetry-broken GS) for t > 0 , the TD effective one-electron potentials can be spin-dependent, and, hence, the up-spin and down-spin densities, which are unequally propagated in the time domain, can be very different at any subsequent time t. As shown, the TD properties, such as the number of bound electrons (see Figure 5), induced dipole moment (see Figure 6), and HHG spectrum (see Figure 7), of H2 with a stretched bond length of 3.78 bohr, obtained with spin-restricted and spin-unrestricted TD-ALDA (i.e., RT-TAO-ALDA with θ = 0), are distinctly different, yielding unphysical spin-symmetry breaking effects in all the TD properties examined. Such an unphysical spin-symmetry breaking feature of spin-unrestricted TD-ALDA is apparently undesirable for RT simulations. By contrast, the spin-symmetry breaking effects in the TD properties obtained with RT-TAO-ALDA are shown to be reducible with the increase in θ at essentially no additional computational cost. In particular, the TD properties obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with θ = 40 mhartree) are essentially the same, yielding essentially no unphysical spin-symmetry breaking effects in all the TD properties examined. This desirable feature can be attributed to the satisfaction of spin-symmetry constraint on the singlet GS density of the stretched H2 by spin-unrestricted TAO-LDA (with θ = 40 mhartree) [26].

5. Conclusions

In conclusion, we have developed RT-TAO-DFT (i.e., an RT extension of TAO-DFT [26]), allowing the study of TD properties of both SR and MR systems. By resorting to an ensemble formalism, RT-TAO-DFT has resolved the aforementioned inconsistency of TDTAO-DFT (especially for θ 0 ) [42]. Since the assumption of a weak perturbation is not required in RT-TAO-DFT, spin-restricted and spin-unrestricted RT-TAO-DFT (with various θ ) calculations have been performed to explore the TD properties (e.g., the number of bound electrons, induced dipole moment, and HHG spectrum) of H 2 at the equilibrium and stretched geometries, aligned along the polarization of an intense linearly polarized laser pulse. The TD properties obtained with RT-TAO-DFT (with various θ ) have been compared with those obtained with conventional TD-DFT [6], which corresponds to RT-TAO-DFT (with θ = 0 ). Moreover, issues related to the possible spin-symmetry breaking effects in the TD properties are also discussed.

Author Contributions

Conceptualization, J.-D.C.; Data curation, H.-Y.T.; Formal analysis, H.-Y.T.; Funding acquisition, J.-D.C.; Investigation, H.-Y.T.; Methodology, H.-Y.T. and J.-D.C.; Project administration, J.-D.C.; Resources, J.-D.C.; Software, J.-D.C.; Supervision, J.-D.C.; Validation, H.-Y.T. and J.-D.C.; Visualization, H.-Y.T.; Writing—original draft, H.-Y.T. and J.-D.C.; Writing—review and editing, J.-D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was previously funded by the Ministry of Science and Technology of Taiwan (Grant No. MOST110-2112-M-002-045-MY3).

Data Availability Statement

The numerical data supporting the findings of the present work are available from the authors upon appropriate request.

Acknowledgments

This work was supported by the Ministry of Science and Technology of Taiwan (Grant No. MOST110-2112-M-002-045-MY3), National Taiwan University, and the National Center for Theoretical Sciences of Taiwan.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kohn, W.; Sham, L.J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef]
  2. Kümmel, S.; Kronik, L. Orbital-dependent density functionals: Theory and applications. Rev. Mod. Phys. 2008, 80, 3–60. [Google Scholar] [CrossRef]
  3. Cohen, A.J.; Mori-Sánchez, P.; Yang, W. Insights into current limitations of density functional theory. Science 2008, 321, 792–794. [Google Scholar] [CrossRef] [PubMed]
  4. Cohen, A.J.; Mori-Sánchez, P.; Yang, W. Challenges for density functional theory. Chem. Rev. 2012, 112, 289–320. [Google Scholar] [CrossRef] [PubMed]
  5. Engel, E.; Dreizler, R.M. Density Functional Theory: An Advanced Course; Springer: Heidelberg, Germany, 2011. [Google Scholar]
  6. Runge, E.; Gross, E.K.U. Density-functional theory for time-dependent systems. Phys. Rev. Lett. 1984, 52, 997–1000. [Google Scholar] [CrossRef]
  7. Onida, G.; Reining, L.; Rubio, A. Electronic excitations: Density-functional versus many-body Green’s-function approaches. Rev. Mod. Phys. 2002, 74, 601–659. [Google Scholar] [CrossRef]
  8. Dreuw, A.; Head-Gordon, M. Single-reference ab initio methods for the calculation of excited states of large molecules. Chem. Rev. 2005, 105, 4009–4037. [Google Scholar] [CrossRef]
  9. Ullrich, C.A. Time-Dependent Density-Functional Theory: Concepts and Applications; Oxford University: New York, NY, USA, 2012. [Google Scholar]
  10. Casida, M.E. Recent Advances in Density Functional Methods, Part I; World Scientific: Singapore, 1995. [Google Scholar]
  11. Dirac, P.A.M. Note on exchange phenomena in the Thomas atom. Proc. Camb. Philos. Soc. 1930, 26, 376–385. [Google Scholar] [CrossRef]
  12. Perdew, J.P.; Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 1992, 45, 13244–13249. [Google Scholar] [CrossRef]
  13. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef]
  14. Casida, M.E.; Jamorski, C.; Casida, K.C.; Salahub, D.R. Molecular excitation energies to high-lying bound states from time-dependent density-functional response theory: Characterization and correction of the time-dependent local density approximation ionization threshold. J. Chem. Phys. 1998, 108, 4439–4449. [Google Scholar] [CrossRef]
  15. Casida, M.E.; Salahub, D.R. Asymptotic correction approach to improving approximate exchange-correlation potentials: Time-dependent density-functional theory calculations of molecular excitation spectra. J. Chem. Phys. 2000, 113, 8918–8935. [Google Scholar] [CrossRef]
  16. Dreuw, A.; Weisman, J.L.; Head-Gordon, M. Long-range charge-transfer excited states in time-dependent density functional theory require non-local exchange. J. Chem. Phys. 2003, 119, 2943–2946. [Google Scholar] [CrossRef]
  17. Dreuw, A.; Head-Gordon, M. Failure of time-dependent density functional theory for long-range charge-transfer excited states: The zincbacteriochlorin-bacteriochlorin and bacteriochlorophyll-spheroidene complexes. J. Am. Chem. Soc. 2004, 126, 4007–4016. [Google Scholar] [CrossRef] [PubMed]
  18. Becke, A.D. A new mixing of Hartree-Fock and local density-functional theories. J. Chem. Phys. 1993, 98, 1372–1377. [Google Scholar] [CrossRef]
  19. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  20. Stephens, P.J.; Devlin, F.J.; Chabalowski, C.F.; Frisch, M.J. Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem. 1994, 98, 11623–11627. [Google Scholar] [CrossRef]
  21. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef]
  22. Grimme, S.; Hansen, A.; Brandenburg, J.G.; Bannwarth, C. Dispersion-corrected mean-field electronic structure methods. Chem. Rev. 2016, 116, 5105–5154. [Google Scholar] [CrossRef]
  23. Grimme, S. Semiempirical hybrid density functional with perturbative second-order correlation. J. Chem. Phys. 2006, 124, 034108. [Google Scholar] [CrossRef]
  24. Peng, D.; Steinmann, S.N.; van Aggelen, H.; Yang, W. Equivalence of particle-particle random phase approximation correlation energy and ladder-coupled-cluster doubles. J. Chem. Phys. 2013, 139, 104112. [Google Scholar] [CrossRef] [PubMed]
  25. van Aggelen, H.; Yang, Y.; Yang, W. Exchange-correlation energy from pairing matrix fluctuation and the particle-particle random-phase approximation. Phys. Rev. A 2013, 88, 030501(R). [Google Scholar] [CrossRef]
  26. Chai, J.-D. Density functional theory with fractional orbital occupations. J. Chem. Phys. 2012, 136, 154104. [Google Scholar] [CrossRef] [PubMed]
  27. Chai, J.-D. Thermally-assisted-occupation density functional theory with generalized-gradient approximations. J. Chem. Phys. 2014, 140, 18A521. [Google Scholar] [CrossRef] [PubMed]
  28. Chai, J.-D. Role of exact exchange in thermally-assisted-occupation density functional theory: A proposal of new hybrid schemes. J. Chem. Phys. 2017, 146, 044102. [Google Scholar] [CrossRef] [PubMed]
  29. Xuan, F.; Chai, J.-D.; Su, H. Local density approximation for the short-range exchange free energy functional. ACS Omega 2019, 4, 7675–7683. [Google Scholar] [CrossRef] [PubMed]
  30. Wu, C.-S.; Chai, J.-D. Electronic properties of zigzag graphene nanoribbons studied by TAO-DFT. J. Chem. Theory Comput. 2015, 11, 2003–2011. [Google Scholar] [CrossRef] [PubMed]
  31. Yeh, C.-N.; Chai, J.-D. Role of Kekulé and non-Kekulé structures in the radical character of alternant polycyclic aromatic hydrocarbons: A TAO-DFT study. Sci. Rep. 2016, 6, 30562. [Google Scholar] [CrossRef]
  32. Seenithurai, S.; Chai, J.-D. Effect of Li adsorption on the electronic and hydrogen storage properties of acenes: A dispersion-corrected TAO-DFT study. Sci. Rep. 2016, 6, 33081. [Google Scholar] [CrossRef]
  33. Wu, C.-S.; Lee, P.-Y.; Chai, J.-D. Electronic properties of cyclacenes from TAO-DFT. Sci. Rep. 2016, 6, 37249. [Google Scholar] [CrossRef]
  34. Tönshoff, C.; Bettinger, H.F. Pushing the limits of acene chemistry: The recent surge of large acenes. Chem. Eur. J. 2021, 27, 3193–3212. [Google Scholar] [CrossRef] [PubMed]
  35. Gupta, D.; Omont, A.; Bettinger, H.F. Energetics of formation of cyclacenes from 2,3-didehydroacenes and implications for astrochemistry. Chem. Eur. J. 2021, 27, 4605–4616. [Google Scholar] [CrossRef] [PubMed]
  36. Hanson-Heine, M.W.D. Static correlation in vibrational frequencies studied using thermally-assisted-occupation density functional theory. Chem. Phys. Lett. 2020, 739, 137012. [Google Scholar] [CrossRef]
  37. Hanson-Heine, M.W.D. Static electron correlation in anharmonic molecular vibrations: A hybrid TAO-DFT study. J. Phys. Chem. A 2022, 126, 7273–7282. [Google Scholar] [CrossRef] [PubMed]
  38. Chen, B.-J.; Chai, J.-D. TAO-DFT fictitious temperature made simple. RSC Adv. 2022, 12, 12193–12210. [Google Scholar] [CrossRef] [PubMed]
  39. Lin, C.-Y.; Hui, K.; Chung, J.-H.; Chai, J.-D. Self-consistent determination of the fictitious temperature in thermally-assisted-occupation density functional theory. RSC Adv. 2017, 7, 50496–50507. [Google Scholar] [CrossRef]
  40. Li, S.; Chai, J.-D. TAO-DFT-based ab initio molecular dynamics. Front. Chem. 2020, 8, 589432. [Google Scholar] [CrossRef]
  41. Seenithurai, S.; Chai, J.-D. TAO-DFT with the polarizable continuum model. Nanomaterials 2023, 13, 1593. [Google Scholar] [CrossRef]
  42. Yeh, S.-H.; Manjanath, A.; Cheng, Y.-C.; Chai, J.-D.; Hsu, C.-P. Excitation energies from thermally assisted-occupation density functional theory: Theory and computational implementation. J. Chem. Phys. 2020, 153, 084120. [Google Scholar] [CrossRef]
  43. Corkum, P.B. Plasma perspective on strong field multiphoton ionization. Phys. Rev. Lett. 1993, 71, 1994–1997. [Google Scholar] [CrossRef]
  44. Lewenstein, M.; Balcou, P.; Ivanov, M.Y.; L’Huillier, A.; Corkum, P.B. Theory of high-harmonic generation by low-frequency laser fields. Phys. Rev. A 1994, 49, 2117–2132. [Google Scholar] [CrossRef] [PubMed]
  45. Paul, P.M.; Toma, E.S.; Breger, P.; Mullot, G.; Augé, F.; Balcou, P.; Muller, H.G.; Agostini, P. Observation of a train of attosecond pulses from high harmonic generation. Science 2001, 292, 1689–1692. [Google Scholar] [CrossRef] [PubMed]
  46. Kienberger, R.; Goulielmakis, E.; Uiberacker, M.; Baltuska, A.; Yakovlev, V.; Bammer, F.; Scrinzi, A.; Westerwalbesloh, T.; Kleineberg, U.; Heinzmann, U.; et al. Atomic transient recorder. Nature 2004, 427, 817–821. [Google Scholar] [CrossRef] [PubMed]
  47. Chu, X.; Chu, S.-I. Self-interaction-free time-dependent density-functional theory for molecular processes in strong fields: High-order harmonic generation of H2 in intense laser fields. Phys. Rev. A 2001, 63, 023411. [Google Scholar] [CrossRef]
  48. Litvinyuk, I.V.; Lee, K.F.; Dooley, P.W.; Rayner, D.M.; Villeneuve, D.M.; Corkum, P.B. Alignment-dependent strong field ionization of molecules. Phys. Rev. Lett. 2003, 90, 233003. [Google Scholar] [CrossRef]
  49. Kanai, T.; Minemoto, S.; Sakai, H. Quantum interference during high-order harmonic generation from aligned molecules. Nature 2005, 435, 470–474. [Google Scholar] [CrossRef]
  50. Guan, X.X.; Tong, X.M.; Chu, S.-I. Effect of electron correlation on high-order-harmonic generation of helium atoms in intense laser fields: Time-dependent generalized pseudospectral approach in hyperspherical coordinates. Phys. Rev. A 2006, 73, 023403. [Google Scholar] [CrossRef]
  51. Baker, S.; Robinson, J.S.; Haworth, C.A.; Teng, H.; Smith, R.A.; Chirilaa, C.C.; Lein, M.; Tisch, J.W.G.; Marangos, J.P. Probing proton dynamics in molecules on an attosecond time scale. Science 2006, 312, 424–427. [Google Scholar] [CrossRef]
  52. Corkum, P.B.; Krausz, F. Attosecond science. Nat. Phys. 2007, 3, 381–387. [Google Scholar] [CrossRef]
  53. Baker, S.; Robinson, J.S.; Lein, M.; Chirilă, C.C.; Torres, R.; Bandulet, H.C.; Comtois, D.; Kieffer, J.C.; Villeneuve, D.M.; Tisch, J.W.G.; et al. Dynamic two-center interference in high-order harmonic generation from molecules with attosecond nuclear motion. Phys. Rev. Lett. 2008, 101, 053901. [Google Scholar] [CrossRef]
  54. Krausz, F.; Ivanov, M. Attosecond physics. Rev. Mod. Phys. 2009, 81, 163–234. [Google Scholar] [CrossRef]
  55. Mack, M.R.; Whitenack, D.; Wasserman, A. Exchange-correlation asymptotics and high harmonic spectra. Chem. Phys. Lett. 2013, 558, 15–19. [Google Scholar] [CrossRef]
  56. Li, P.-C.; Sheu, Y.-L.; Laughlin, C.; Chu, S.-I. Dynamical origin of near- and below-threshold harmonic generation of Cs in an intense mid-infrared laser field. Nat. Comm. 2015, 6, 7178. [Google Scholar] [CrossRef] [PubMed]
  57. Heslar, J.; Telnov, D.A.; Chu, S.-I. Subcycle dynamics of high-harmonic generation in valence-shell and virtual states of Ar atoms: A self-interaction-free time-dependent density-functional-theory approach. Phys. Rev. A 2015, 91, 023420. [Google Scholar] [CrossRef]
  58. Chou, Y.; Li, P.-C.; Ho, T.-S.; Chu, S.-I. Optimal control of high-order harmonics for the generation of an isolated ultrashort attosecond pulse with two-color midinfrared laser fields. Phys. Rev. A 2015, 91, 063408. [Google Scholar] [CrossRef]
  59. Sun, H.-L.; Peng, W.-T.; Chai, J.-D. Assessment of the LFAs-PBE exchange-correlation potential for high-order harmonic generation of aligned H 2 + molecules. RSC Adv. 2016, 6, 33318–33325. [Google Scholar] [CrossRef]
  60. Zhu, Y.; Herbert, J.M. High harmonic spectra computed using time-dependent Kohn-Sham theory with Gaussian orbitals and a complex absorbing potential. J. Chem. Phys. 2022, 156, 204123. [Google Scholar] [CrossRef]
  61. Hohenberg, P.; Kohn, W. Inhomogeneous electron gas. Phys. Rev. 1964, 136, B864–B871. [Google Scholar] [CrossRef]
  62. Mermin, N.D. Thermal properties of the inhomogeneous electron gas. Phys. Rev. 1965, 137, A1441–A1443. [Google Scholar] [CrossRef]
  63. Helgaker, T.; Jørgensen, P.; Olsen, J. Molecular Electronic-Structure Theory; Wiley: New York, NY, USA, 2000. [Google Scholar]
  64. Löwdin, P.-O.; Shull, H. Natural orbitals in the quantum theory of two-electron systems. Phys. Rev. 1956, 101, 1730–1739. [Google Scholar] [CrossRef]
  65. Levy, M. Universal variational functionals of electron densities, first-order density matrices, and natural spin-orbitals and solution of the v-representability problem. Proc. Natl. Acad. Sci. USA 1979, 76, 6062–6065. [Google Scholar] [CrossRef] [PubMed]
  66. Levy, M. Electron densities in search of Hamiltonians. Phys. Rev. A 1982, 26, 1200–1208. [Google Scholar] [CrossRef]
  67. Lieb, E.H. Density functionals for coulomb systems. Int. J. Quantum Chem. 1983, 24, 243–277. [Google Scholar] [CrossRef]
  68. Schipper, P.R.T.; Gritsenko, O.V.; Baerends, E.J. One-determinantal pure state versus ensemble Kohn-Sham solutions in the case of strong electron correlation: CH2 and C2. Theor. Chem. Acc. 1998, 99, 329–343. [Google Scholar] [CrossRef]
  69. Morrison, R.C. Electron correlation and noninteracting v-representability in density functional theory: The Be isoelectronic series. J. Chem. Phys. 2002, 117, 10506–10511. [Google Scholar] [CrossRef]
  70. Katriel, J.; Roy, S.; Springborg, M. A study of the adiabatic connection for two-electron systems. J. Chem. Phys. 2004, 121, 12179–12190. [Google Scholar] [CrossRef] [PubMed]
  71. Ayers, P.W.; Yang, W. Legendre-transform functionals for spin-density-functional theory. J. Chem. Phys. 2006, 124, 224108. [Google Scholar] [CrossRef] [PubMed]
  72. Kanhere, D.G.; Panat, P.V.; Rajagopal, A.K.; Callaway, J. Exchange-correlation potentials for spin-polarized systems at finite temperatures. Phys. Rev. A 1986, 33, 490–497. [Google Scholar] [CrossRef]
  73. Desjarlais, M.P.; Kress, J.D.; Collins, L.A. Electrical conductivity for warm, dense aluminum plasmas and liquids. Phys. Rev. E 2002, 66, 025401(R). [Google Scholar] [CrossRef]
  74. Grimme, S.; Hansen, A. A practicable real-space measure and visualization of static electron-correlation effects. Angew. Chem. Int. Ed. 2015, 54, 12308–12313. [Google Scholar] [CrossRef]
  75. Bauer, C.A.; Hansen, A.; Grimme, S. The fractional occupation number weighted density as a versatile analysis tool for molecules with a complicated electronic structure. Chem. Eur. J. 2017, 23, 6150–6164. [Google Scholar] [CrossRef] [PubMed]
  76. Pérez-Guardiola, A.; Sandoval-Salinas, M.E.; Casanova, D.; San-Fabián, E.; Pérez-Jiménez, A.J.; Sancho-Garcia, J.C. The role of topology in organic molecules: Origin and comparison of the radical character in linear and cyclic oligoacenes and related oligomers. Phys. Chem. Chem. Phys. 2018, 20, 7112–7124. [Google Scholar] [CrossRef] [PubMed]
  77. Yoshikawa, T.; Doi, T.; Nakai, H. Finite-temperature-based linear-scaling divide-and-conquer self-consistent field method for static electron correlation systems. Chem. Phys. Lett. 2019, 725, 18–23. [Google Scholar] [CrossRef]
  78. Liu, F.; Duan, C.; Kulik, H.J. Rapid detection of strong correlation with machine learning for transition-metal complex high-throughput screening. J. Phys. Chem. Lett. 2020, 11, 8067–8076. [Google Scholar] [CrossRef] [PubMed]
  79. Cho, Y.; Nandy, A.; Duan, C.; Kulik, H.J. DFT-based multireference diagnostics in the solid state: Application to metal-organic frameworks. J. Chem. Theory Comput. 2023, 19, 190–197. [Google Scholar] [CrossRef] [PubMed]
  80. Sandoval-Salinas, M.E.; Bernabeu-Cabanero, R.; Pérez-Jiménez, A.J.; San-Fabián, E.; Sancho-Garcia, J.C. Electronic structure of rhombus-shaped nanographenes: System size evolution from closed- to open-shell ground states. Phys. Chem. Chem. Phys. 2023, 25, 11697–11706. [Google Scholar] [CrossRef] [PubMed]
  81. Nieman, R.; Carvalho, J.R.; Jayee, B.; Hansen, A.; Aquino, A.; Kertesz, M.; Lischka, H. Polyradical character assessment using multireference calculations and comparison with density-functional derived fractional occupation number weighted density analysis. Phys. Chem. Chem. Phys. 2023, 25, 27380–27393. [Google Scholar] [CrossRef]
  82. Yeh, S.-H.; Yang, W.; Hsu, C.-P. Reformulation of thermally assisted-occupation density functional theory in the Kohn-Sham framework. J. Chem. Phys. 2022, 156, 174108. [Google Scholar] [CrossRef]
  83. Teale, A.M.; Helgaker, T.; Savin, A.; Adamo, C.; Aradi, B.; Arbuznikov, A.V.; Ayers, P.W.; Baerends, E.J.; Barone, V.; Calaminici, P.; et al. DFT exchange: Sharing perspectives on the workhorse of quantum chemistry and materials science. Phys. Chem. Chem. Phys. 2022, 24, 28700–28781. [Google Scholar] [CrossRef]
  84. Greiner, M.; Carrier, P.; Görling, A. Extension of exact-exchange density functional theory of solids to finite temperatures. Phys. Rev. B 2010, 81, 155119. [Google Scholar] [CrossRef]
  85. Li, T.-C.; Tong, P.-Q. Hohenberg-Kohn theorem for time-dependent ensembles. Phys. Rev. A 1985, 31, 1950–1951. [Google Scholar] [CrossRef] [PubMed]
  86. van Leeuwen, R. Mapping from densities to potentials in time-dependent density-functional theory. Phys. Rev. Lett. 1999, 82, 3863–3866. [Google Scholar] [CrossRef]
  87. Dufty, J.; Luo, K.; Trickey, S.B. Density response from kinetic theory and time-dependent density-functional theory for matter under extreme conditions. Phys. Rev. E 2018, 98, 033203. [Google Scholar] [CrossRef]
  88. Li, X.; Smith, S.M.; Markevitch, A.N.; Romanov, D.A.; Levis, R.J.; Schlegel, H.B. A time-dependent Hartree-Fock approach for studying the electronic optical response of molecules in intense fields. Phys. Chem. Chem. Phys. 2005, 7, 233–239. [Google Scholar] [CrossRef]
  89. Zhu, Y.; Herbert, J.M. Self-consistent predictor/corrector algorithms for stable and efficient integration of the time-dependent Kohn-Sham equation. J. Chem. Phys. 2018, 148, 044117. [Google Scholar] [CrossRef] [PubMed]
  90. Bedurke, F.; Klamroth, T.; Saalfrank, P. Many-electron dynamics in laser-driven molecules: Wavefunction theory vs. density functional theory. Phys. Chem. Chem. Phys. 2021, 23, 13544–13560. [Google Scholar] [CrossRef] [PubMed]
  91. Cormier, E.; Lambropoulos, P. Optimal gauge and gauge invariance in non-perturbative time-dependent calculation of above-threshold ionization. J. Phys. B At. Mol. Opt. Phys. 1996, 29, 1667–1680. [Google Scholar] [CrossRef]
  92. Bedurke, F.; Klamroth, T.; Krause, P.; Saalfrank, P. Discriminating organic isomers by high harmonic generation: A time-dependent configuration interaction singles study. J. Chem. Phys. 2019, 150, 234114. [Google Scholar] [CrossRef]
  93. Isborn, C.M.; Li, X. Singlet-triplet transitions in real-time time-dependent Hartree-Fock/density functional theory. J. Chem. Theory Comput. 2009, 5, 2415–2419. [Google Scholar] [CrossRef]
  94. White, A.F.; Heide, C.J.; Saalfrank, P.; Head-Gordon, M.; Luppi, E. Computation of high-harmonic generation spectra of the hydrogen molecule using time-dependent configuration-interaction. Mol. Phys. 2016, 114, 947–956. [Google Scholar] [CrossRef]
  95. Epifanovsky, E.; Gilbert, A.T.; Feng, X.; Lee, J.; Mao, Y.; Mardirossian, N.; Pokhilko, P.; White, A.F.; Coons, M.P.; Dempwolff, A.L.; et al. Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package. J. Chem. Phys. 2021, 155, 084801. [Google Scholar] [CrossRef]
Figure 2. Number of bound electrons for H 2 with an equilibrium bond length of 1.45 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Figure 2. Number of bound electrons for H 2 with an equilibrium bond length of 1.45 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Molecules 28 07247 g002
Figure 3. Induced dipole moment for H 2 with an equilibrium bond length of 1.45 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Figure 3. Induced dipole moment for H 2 with an equilibrium bond length of 1.45 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Molecules 28 07247 g003
Figure 4. HHG spectrum for H 2 with an equilibrium bond length of 1.45 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Figure 4. HHG spectrum for H 2 with an equilibrium bond length of 1.45 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Molecules 28 07247 g004
Figure 5. Number of bound electrons for H 2 with a stretched bond length of 3.78 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Figure 5. Number of bound electrons for H 2 with a stretched bond length of 3.78 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Molecules 28 07247 g005
Figure 6. Induced dipole moment for H 2 with a stretched bond length of 3.78 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Figure 6. Induced dipole moment for H 2 with a stretched bond length of 3.78 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Molecules 28 07247 g006aMolecules 28 07247 g006b
Figure 7. HHG spectrum for H 2 with a stretched bond length of 3.78 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Figure 7. HHG spectrum for H 2 with a stretched bond length of 3.78 bohr, obtained with spin-restricted and spin-unrestricted RT-TAO-ALDA (with various θ ). Here, the θ = 0 case corresponds to TD-ALDA.
Molecules 28 07247 g007
Table 1. Comparison of KS-DFT [1], TAO-DFT [26], and FT-DFT [1,62].
Table 1. Comparison of KS-DFT [1], TAO-DFT [26], and FT-DFT [1,62].
KS-DFTTAO-DFTFT-DFT
Electronic Temperature θ e l 00≥0
Fictitious Temperature θ 0≥0≥0
Is θ = θ e l assumed?YesNoYes
Electronic PropertyGSGSThermal Equilibrium
Electron DensityGSGSThermal Equilibrium
Density RepresentationNI-PS-VRNI-TE-VRNI-TE-VR
Universal FunctionalHohenberg–KohnHohenberg–KohnMermin
Approximate Functional E xc [ ρ ] E xc θ [ ρ ] F xc θ e l [ ρ θ e l ]
Table 2. Comparison of TAO-DFT (with E xc θ [ ρ ] E xc [ ρ ] ) [26] and FT-DFT (with F xc θ e l [ ρ θ e l ] E xc [ ρ θ e l ] ) [1,62].
Table 2. Comparison of TAO-DFT (with E xc θ [ ρ ] E xc [ ρ ] ) [26] and FT-DFT (with F xc θ e l [ ρ θ e l ] E xc [ ρ θ e l ] ) [1,62].
TAO-DFT (with E xc θ [ ρ ] E xc [ ρ ] )FT-DFT (with F xc θ el [ ρ θ el ] E xc [ ρ θ el ] )
Electronic Temperature θ e l 0≥0
Fictitious Temperature θ ≥0≥0
Is θ = θ e l assumed?NoYes
Electronic PropertyGSThermal Equilibrium
Electron DensityGSThermal Equilibrium
Density RepresentationNI-TE-VRNI-TE-VR
Approximate Functional E xc [ ρ ] E xc [ ρ θ e l ]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tsai, H.-Y.; Chai, J.-D. Real-Time Extension of TAO-DFT. Molecules 2023, 28, 7247. https://doi.org/10.3390/molecules28217247

AMA Style

Tsai H-Y, Chai J-D. Real-Time Extension of TAO-DFT. Molecules. 2023; 28(21):7247. https://doi.org/10.3390/molecules28217247

Chicago/Turabian Style

Tsai, Hung-Yi, and Jeng-Da Chai. 2023. "Real-Time Extension of TAO-DFT" Molecules 28, no. 21: 7247. https://doi.org/10.3390/molecules28217247

APA Style

Tsai, H. -Y., & Chai, J. -D. (2023). Real-Time Extension of TAO-DFT. Molecules, 28(21), 7247. https://doi.org/10.3390/molecules28217247

Article Metrics

Back to TopTop