Search

We model a spin-phase transition in a two-dimensional square array, or a lateral superlattice, of quantum rings in an external perpendicular homogeneous magnetic field. The electron system is placed in a circular cylindrical far-infrared photon cavity with a single circularly symmetric photon mode. Our numerical results reveal that the spin ordering of the two-dimensional electron gas in each quantum ring can be influenced or controlled by the electron-photon coupling strength and the energy of the photons. The Coulomb interaction between the electrons is described by a spin-density functional approach, but the para- and the diamagnetic electron-photon interactions are modeled via a configuration interaction formalism in a truncated many-body Fock-space, which is updated in each iteration step of the density functional approach. In the absence of external electromagnetic pulses this spin-phase transition is replicated in the orbital magnetization of the rings. The spin-phase transition can be suppressed by a strong electron-photon interaction. In addition, fluctuations in the spin configuration are found in dynamical calculations, where the system is excited by a time-dependent scheme specially fit for emphasizing the diamagnetic electron-photon interaction., Comment: RevTeX - pdfLaTeX, 11 pages with 9 included pdf figures

We employ a ``real-time'' excitation scheme to calculate the excitation spectra of a two-dimensional electron system in a square array of quantum dots placed in a circular cylindrical far-infrared photon cavity subjected to a perpendicular homogeneous external magnetic field. The Coulomb interaction of the electrons is handled via spin density functional theory and the para- and the diamagnetic parts of the electron-photon coupling are updated according to a configuration interaction method in each iteration of the density functional calculation. The results show that an excitation scheme built on using the symmetry of the lateral square superlattice of the dots and the cylindrical cavity produces both para- and diamagnetic resonance peaks with oscillator strengths that can be steered by the excitation pulse parameters. The excitation method breaks the conditions for the generalized Kohn theorem and allows for insight into the subband structure of the electron system and can be used both in and outside the linear response regime., Comment: RevTeX - pdfLaTeX, 14 pages with 15 included pdf and png figures

For electronic systems with multi-reference (MR) character, Kohn-Sham density functional theory (KS-DFT) with the conventional exchange-correlation (xc) energy functionals can lead to incorrect spin densities and related properties. For example, for H2 dissociation, the spin-restricted and spin-unrestricted solutions obtained with the same xc energy functional in KS-DFT can be distinctly different, yielding the unphysical spin-symmetry breaking effects in the spin-unrestricted solutions. Recently, thermally-assisted-occupation density functional theory (TAO-DFT) has been shown to resolve the aforementioned spin-symmetry breaking, when the fictitious temperature is properly chosen. In this work, a response theory based on TAO-DFT is developed to demonstrate that TAO-DFT with a sufficiently large fictitious temperature can always resolve the unphysical spin-symmetry breaking in MR systems. To further support this, TAO-DFT calculations with various fictitious temperatures are performed for the dissociation of H2, N2, He2, and Ne2 as well as the twisted ethylene., Comment: accepted for publication in Phys. Rev. A, 23 pages, 5 figures

We model the equilibrium properties of a two-dimensional electron gas in a square lateral superlattice of quantum dots in a GaAs heterostructure subject to an external homogeneous perpendicular magnetic field and a far-infrared circular cylindrical photon cavity with one quantized mode, the TE011 mode. In a truncated linear basis constructed by a tensor product of the single-electron states of the noninteracting system and the eigenstates of the photon number operator, a local spin density approximation of density functional theory is used to compute the electron-photon states of the two-dimensional electron gas in the cavity. The common spatial symmetry of the vector fields for the external magnetic field and the cavity photon field in the long wavelength approximation enhances higher order magnetic single- and multi-photon processes for both the para- and the diamagnetic electron-photon interactions. The electron-photon coupling introduces explicit photon replicas into the bandstructure and all subbands gain a photon content, constant for each subband, that can deviate from an integer value as the coupling is increased or the photon energy is varied. The subbands show a complex Rabi anticrossing behavior when the photon energy and the coupling bring subbands into resonances. The complicated energy subband structure leads to photon density variations in reciprocal space when resonances occur in the spectrum. The electron-photon coupling polarizes the charge density and tends to reduce the Coulomb exchange effects as the coupling strength increases., Comment: RevTeX - pdfLaTeX, 12 pages with 8 included pdf figures

The time-dependent density functional theory (TDDFT) has been broadly used to investigate the excited-state properties of various molecular systems. However, the current TDDFT heavily relies on outcomes from the corresponding ground-state density functional theory (DFT) calculations which may be prone to errors due to the lack of proper treatment in the non-dynamical correlation effects. Recently, thermally-assisted-occupation density functional theory (TAO-DFT) [J.-D. Chai, \textit{J. Chem. Phys.} \textbf{136}, 154104 (2012)], a DFT with fractional orbital occupations, was proposed, explicitly incorporating the non-dynamical correlation effects in the ground-state calculations with low computational complexity. In this work, we develop time-dependent (TD) TAO-DFT, which is a time-dependent, linear-response theory for excited states within the framework of TAO-DFT. With tests on the excited states of H$_{2}$, the first triplet excited state ($1^3\Sigma_u^+$) was described well, with non-imaginary excitation energies. TDTAO-DFT also yields zero singlet-triplet gap in the dissociation limit, for the ground singlet ($1^1\Sigma_g^+$) and the first triplet state ($1^3\Sigma_u^+$). In addition, as compared to traditional TDDFT, the overall excited-state potential energy surfaces obtained from TDTAO-DFT are generally improved and better agree with results from the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD).

A combined density and first-order reduced-density-matrix (1RDM) functional method is proposed for the calculation of potential energy curves (PECs) of molecular multibond dissociation. Its 1RDM functional part, a pair density functional, efficiently approximates the ab initio pair density of the complete active space self-consistent-field (CASSCF) method. The corresponding approximate on top pair density {\Pi} is employed to correct for double counting a correlation functional of density functional theory (DFT). The proposed ELS-DM{\Pi}DFT method with the extended L\"owdin-Shull (ELS) 1RDM functional with dispersion and multibond (DM) corrections augmented with the {\Pi}DFT functional closely reproduces PECs of multibond dissociation in the paradigmatic N_2 , H_2O, and H_2CO molecules calculated with the recently proposed CAS{\Pi}DFT (CASSCF augmented with a {\Pi} based scaled DFT correlation correction) method. Furthermore, with the additional M-correction, ELS-DM{\Pi}DFT+M reproduces well the benchmark PEC of the N_2 molecule by Lie and Clementi.

A two-orbital two-electron diatomic model resembling LiH is used to investigate the differences between the exact L\"owdin-Shull and approximate Hartree-Fock-Bogoliubov and Baerends-Buijse density matrix functionals in the medium- to long-distance dissociation region. In case of homolytic dissociation (one electron on each atom), the approximate functionals fail to generate the correct energy due to a compromise between the Hartree-Fock component (which favors partial charge transfer) and the strong correlation component (which hampers charge transfer). The exact functional is able to generate the physically correct answer by enforcing the equi-charge distribution of the bonding and antibonding orbitals. Besides, the approximate functionals also have issues in correctly describing heterolytic dissociation (two electrons on one atom) due to the strong correlation component hampering charge transfer. In this work, we propose a new scheme in which the homolytic dissociation problem for approximate functionals is avoided by adding a Lagrange multiplier that enforces equi-charge distribution of the bonding and antibonding orbitals. The symmety based nature of the findings implies that they are most likely transferable to other cases in which one uses an approximate one-particle method in conjunction with a symmetrical particle-hole correction factor.

Accurate prediction of the electronic and hydrogen storage properties of linear carbon chains (Cn) and Li-terminated linear carbon chains (Li2Cn), with n carbon atoms (n = 5 - 10), has been very challenging for traditional electronic structure methods, due to the presence of strong static correlation effects. To meet the challenge, we study these properties using our newly developed thermally-assisted-occupation density functional theory (TAO-DFT), a very efficient electronic structure method for the study of large systems with strong static correlation effects. Owing to the alteration of the reactivity of Cn and Li2Cn with n, odd-even oscillations in their electronic properties are found. In contrast to Cn, the binding energies of H2 molecules on Li2Cn are in (or close to) the ideal binding energy range (about 20 to 40 kJ/mol per H2). In addition, the H2 gravimetric storage capacities of Li2Cn are in the range of 10.7 to 17.9 wt%, satisfying the United States Department of Energy (USDOE) ultimate target of 7.5 wt%. On the basis of our results, Li2Cn can be high-capacity hydrogen storage materials for reversible hydrogen uptake and release at near-ambient conditions., Comment: 25 pages, 11 figures. This is an extension of our previous work [e.g., see arXiv:1606.03489]

We propose hybrid schemes incorporating exact exchange into thermally-assisted-occupation density functional theory (TAO-DFT) [J.-D. Chai, J. Chem. Phys. 136, 154104 (2012)] for an improved description of nonlocal exchange effects. With a few simple modifications, global and range-separated hybrid functionals in Kohn-Sham density functional theory (KS-DFT) can be combined seamlessly with TAO-DFT. In comparison with global hybrid functionals in KS-DFT, the resulting global hybrid functionals in TAO-DFT yield promising performance for systems with strong static correlation effects (e.g., the dissociation of H2 and N2, twisted ethylene, and electronic properties of linear acenes), while maintaining similar performance for systems without strong static correlation effects. Besides, a reasonably accurate description of noncovalent interactions can be efficiently achieved through the inclusion of dispersion corrections in hybrid TAO-DFT. Relative to semilocal density functionals in TAO-DFT, global hybrid functionals in TAO-DFT are generally superior in performance for a wide range of applications, such as thermochemistry, kinetics, reaction energies, and optimized geometries., Comment: 49 pages, 17 figures, 4 tables, supplementary material not included. This is an extension of our previous works [e.g., see arXiv:1201.4866, arXiv:1401.0380]

We propose a short- and long-range corrected (SLC) hybrid scheme employing 100% Hartree-Fock (HF) exchange at both zero and infinite interelectronic distances, wherein three SLC hybrid density functionals with the D3 dispersion corrections (SLC-LDA-D3, SLC-PBE-D3, and SLC-B97-D3) are developed. SLC-PBE-D3 and SLC-B97-D3 are shown to be accurate for a very diverse range of applications, such as core ionization and excitation energies, thermochemistry, kinetics, noncovalent interactions, dissociation of symmetric radical cations, vertical ionization potentials, vertical electron affinities, fundamental gaps, and valence, Rydberg, and long-range charge-transfer excitation energies. Relative to omegaB97X-D, SLC-B97-D3 provides significant improvement for core ionization and excitation energies and noticeable improvement for the self-interaction, asymptote, energy-gap, and charge-transfer problems, while performing similarly for thermochemistry, kinetics, and noncovalent interactions., Comment: accepted for publication in J. Chem. Phys., 35 pages, 7 figures, 10 tables, supplementary material not included

Owing to the presence of strong static correlation effects, accurate prediction of the electronic properties (e.g., the singlet-triplet energy gaps, vertical ionization potentials, vertical electron affinities, fundamental gaps, symmetrized von Neumann entropy, active orbital occupation numbers, and real-space representation of active orbitals) of cyclacenes with n fused benzene rings (n = 4-100) has posed a great challenge to traditional electronic structure methods. To meet the challenge, we study these properties using our newly developed thermally-assisted-occupation density functional theory (TAO-DFT), a very efficient method for the study of large systems with strong static correlation effects. Besides, to examine the role of cyclic topology, the electronic properties of cyclacenes are also compared with those of acenes. Similar to acenes, the ground states of cyclacenes are singlets for all the cases studied. In contrast to acenes, the electronic properties of cyclacenes, however, exhibit oscillatory behavior (for n <= 30) in the approach to the corresponding properties of acenes with increasing number of benzene rings. On the basis of the calculated orbitals and their occupation numbers, the larger cyclacenes are shown to exhibit increasing polyradical character in their ground states, with the active orbitals being mainly localized at the peripheral carbon atoms., Comment: accepted for publication in Sci. Rep., 24 pages, 10 figures, supplementary material not included

Due to the presence of strong static correlation effects and noncovalent interactions, accurate prediction of the electronic and hydrogen storage properties of Li-adsorbed acenes with n linearly fused benzene rings (n = 3 - 8) has been very challenging for conventional electronic structure methods. To meet the challenge, we study these properties using our recently developed thermally-assisted-occupation density functional theory (TAO-DFT) with dispersion corrections. In contrast to pure acenes, the binding energies of H2 molecules on Li-adsorbed acenes are in the ideal binding energy range (about 20 to 40 kJ/mol per H2). Besides, the H2 gravimetric storage capacities of Li-adsorbed acenes are in the range of 9.9 to 10.7 wt%, satisfying the United States Department of Energy (USDOE) ultimate target of 7.5 wt%. On the basis of our results, Li-adsorbed acenes can be high-capacity hydrogen storage materials for reversible hydrogen uptake and release at ambient conditions., Comment: 25 pages, 9 figures, 1 table

We investigate the role of Kekule and non-Kekule structures in the radical character of alternant polycyclic aromatic hydrocarbons (PAHs) using thermally-assisted-occupation density functional theory (TAO-DFT), an efficient electronic structure method for the study of large ground-state systems with strong static correlation effects. Our results reveal that the studies of Kekule and non-Kekule structures qualitatively describe the radical character of alternant PAHs, which could be useful when electronic structure calculations are infeasible due to the expensive computational cost. In addition, our results support previous findings on the increase in radical character with increasing system size. For alternant PAHs with the same number of aromatic rings, the geometrical arrangements of aromatic rings are responsible for their radical character., Comment: 37 pages, 12 figures, 5 tables

In the present study, a series of planar poly(p-phenylene) (PPP) oligomers with n phenyl rings (n = 1 - 20), designated as n-PP, are taken as finite-size models of the narrowest armchair graphene nanoribbons with hydrogen passivation. The singlet-triplet energy gap, vertical ionization potential, vertical electron affinity, fundamental gap, optical gap, and exciton binding energy of n-PP are calculated using Kohn-Sham density functional theory and time-dependent density functional theory with various exchange-correlation density functionals. The ground state of n-PP is shown to be singlet for all the chain lengths studied. In contrast to the lowest singlet state (i.e., the ground state), the lowest triplet state and the ground states of the cation and anion of n-PP are found to exhibit some multi-reference character. Overall, the electronic and optical properties of n-PP obtained from the omegaB97 and omegaB97X functionals are in excellent agreement with the available experimental data., Comment: 24 pages, 8 figures, 1 table, supplementary material not included

We examine the performance of our recently developed LFAs-PBE exchange-correlation (XC) potential [C.-R. Pan, P.-T. Fang, and J.-D. Chai, Phys. Rev. A, 2013, 87, 052510] for the high-order harmonic generation (HHG) spectra and related properties of H2+ molecules aligned parallel and perpendicular to the polarization of an intense linearly polarized laser pulse, employing the real-time formulation of time-dependent density functional theory (RT-TDDFT). The results are compared with the exact solutions of the time-dependent Schrodinger equation as well as those obtained with other XC potentials in RT-TDDFT. Owing to its correct (-1/r) asymptote, the LFAs-PBE potential significantly outperforms conventional XC potentials for the HHG spectra and the properties that are sensitive to the XC potential asymptote. Accordingly, the LFAs-PBE potential, which has a computational cost similar to that of the popular Perdew-Burke-Ernzerhof (PBE) potential, can be very promising for the study of the ground-state, excited-state, and time-dependent properties of large electronic systems, extending the applicability of density functional methods for a diverse range of applications., Comment: 23 pages, 9 figures, 1 table

By incorporating the nonempirical SCAN semilocal density functional [Sun, Ruzsinszky, and Perdew, Phys. Rev. Lett. 115, 036402 (2015)] in the underlying expression of four existing hybrid and double-hybrid models, we propose one hybrid (SCAN0) and three double-hybrid (SCAN0-DH, SCAN-QIDH, and SCAN0-2) density functionals, which are free from any fitted parameters. The SCAN-based double-hybrid functionals consistently outperform their parent SCAN semilocal functional for self-interaction problems and noncovalent interactions. In particular, SCAN0-2, which includes about 79% of Hartree-Fock exchange and 50% of second-order Moller-Plesset correlation, is shown to be reliably accurate for a very diverse range of applications, such as thermochemistry, kinetics, noncovalent interactions, and self-interaction problems., Comment: 23 pages, 6 figures, 1 table, supplementary material not included

The exciton binding energy, the energy required to dissociate an excited electron-hole pair into free charge carriers, is one of the key factors to the optoelectronic performance of organic materials. However, it remains unclear whether modern quantum-mechanical calculations, mostly based on Kohn-Sham density functional theory (KS-DFT) and time-dependent density functional theory (TDDFT), are reliably accurate for exciton binding energies. In this study, the exciton binding energies and related optoelectronic properties (e.g., the ionization potentials, electron affinities, fundamental gaps, and optical gaps) of 121 small- to medium-sized molecules are calculated using KS-DFT and TDDFT with various density functionals. Our KS-DFT and TDDFT results are compared with those calculated using highly accurate CCSD and EOM-CCSD methods, respectively. The omegaB97, omegaB97X, and omegaB97X-D functionals are shown to generally outperform (with a mean absolute error of 0.36 eV) other functionals for the properties investigated., Comment: 22 pages, 5 figures, 1 table, supplementary material not included

We investigate the potential energy curves of rare-gas dimers with various ranges and strengths of interparticle interactions (nuclear-electron, electron-electron, and nuclear-nuclear interactions). Our investigation is based on the highly accurate coupled-cluster theory associated with those interparticle interactions. For comparison, the performance of the corresponding Hartree-Fock theory, second-order Moller-Plesset perturbation theory, and density functional theory is also investigated. Our results reveal that when the interparticle interactions retain the long-range Coulomb tails, the nature of van der Waals interactions in the rare-gas dimers remains similar. By contrast, when the interparticle interactions are sufficiently short-range, the conventional van der Waals interactions in the rare-gas dimers completely disappear, yielding purely repulsive potential energy curves., Comment: 28 pages, 11 figures

We examine the performance of the asymptotically corrected model potential scheme on the two lowest singlet excitation energies of acenes with different number of linearly fused benzene rings (up to 5), employing both the real-time time-dependent density functional theory and the frequency-domain formulation of linear-response time-dependent density functional theory. The results are compared with the experimental data and those calculated by long-range corrected hybrid functionals and others. The long-range corrected hybrid scheme is shown to outperform the asymptotically corrected model potential scheme for charge-transfer-like excitations., Comment: accepted for publication in Phys. Chem. Chem. Phys., 15 pages, 5 figures

We extend the recently proposed thermally-assisted-occupation density functional theory (TAO-DFT) [J.-D. Chai, J. Chem. Phys. 136, 154104 (2012)] to generalized-gradient approximation (GGA) exchange-correlation density functionals. Relative to our previous TAO-LDA (i.e., the local density approximation to TAO-DFT), the resulting TAO-GGAs are significantly superior for a wide range of applications, such as thermochemistry, kinetics, and reaction energies. For noncovalent interactions, TAO-GGAs with empirical dispersion corrections are shown to yield excellent performance. Due to their computational efficiency for systems with strong static correlation effects, TAO-LDA and TAO-GGAs are applied to study the electronic properties (e.g., the singlet-triplet energy gaps, vertical ionization potentials, vertical electron affinities, fundamental gaps, and symmetrized von Neumann entropy) of acenes with different number of linearly fused benzene rings (up to 100), which is very challenging for conventional electronic structure methods. The ground states of acenes are shown to be singlets for all the chain lengths studied here. With the increase of acene length, the singlet-triplet energy gaps, vertical ionization potentials, and fundamental gaps decrease monotonically, while the vertical electron affinities and symmetrized von Neumann entropy (i.e., a measure of polyradical character) increase monotonically., Comment: 27 pages, 15 figures, 3 tables, supplementary material not included. This is an extension of our previous work [e.g., see arXiv:1201.4866]

Long-range corrected (LC) hybrid functionals and asymptotically corrected (AC) model potentials are two distinct density functional methods with correct asymptotic behavior. They are known to be accurate for properties that are sensitive to the asymptote of the exchange-correlation potential, such as the highest occupied molecular orbital energies and Rydberg excitation energies of molecules. To provide a comprehensive comparison, we investigate the performance of the two schemes and others on a very wide range of applications, including the asymptote problems, self-interaction-error problems, energy-gap problems, charge-transfer problems, and many others. The LC hybrid scheme is shown to consistently outperform the AC model potential scheme. In addition, to be consistent with the molecules collected in the IP131 database [Y.-S. Lin, C.-W. Tsai, G.-D. Li, and J.-D. Chai, J. Chem. Phys., 2012, 136, 154109], we expand the EA115 and FG115 databases to include, respectively, the vertical electron affinities and fundamental gaps of the additional 16 molecules, and develop a new database AE113 (113 atomization energies), consisting of accurate reference values for the atomization energies of the 113 molecules in IP131. These databases will be useful for assessing the accuracy of density functional methods., Comment: accepted for publication in Phys. Chem. Chem. Phys., 46 pages, 4 figures, supplementary material included

From the perspective of perturbation theory, we propose a systematic procedure for the evaluation of the derivative discontinuity (DD) of the exchange-correlation energy functional in Kohn-Sham density functional theory (KS-DFT), wherein the exact DD can in principle be obtained by summing up all the perturbation corrections to infinite order. Truncation of the perturbation series at low order yields an efficient scheme for obtaining the approximate DD. While the zeroth-order theory yields a vanishing DD, the first-order correction to the DD can be expressed as an explicit universal functional of the ground-state density and the KS lowest unoccupied molecular orbital density, allowing the direct evaluation of the DD in the standard KS method without extra computational cost. The fundamental gap can be predicted by adding the estimated DD to the KS gap. This scheme is shown to be accurate in the prediction of the fundamental gaps for a wide variety of atoms and molecules., Comment: accepted for publication in Phys. Rev. Lett., 24 pages, 2 figures, supplemental material included

By incorporating the improved empirical atom-atom dispersion corrections from DFT-D3 [Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. J. Chem. Phys. 2010, 132, 154104], two long-range corrected (LC) hybrid density functionals are proposed. Our resulting LC hybrid functionals, omegaM06-D3 and omegaB97X-D3, are shown to be accurate for a very wide range of applications, such as thermochemistry, kinetics, noncovalent interactions, frontier orbital energies, fundamental gaps, and long-range charge-transfer excitations, when compared with common global and LC hybrid functionals. Relative to omegaB97X-D [Chai, J.-D.; Head-Gordon, M. Phys. Chem. Chem. Phys. 2008, 10, 6615], omegaB97X-D3 (reoptimization of omegaB97X-D with improved dispersion corrections) is shown to be superior for non-bonded interactions, and similar in performance for bonded interactions, while omegaM06-D3 is shown to be superior for general applications., Comment: accepted for publication in J. Chem. Theory Comput., DOI: 10.1021/ct300715s (2012)

Aiming to remedy the incorrect asymptotic behavior of conventional semilocal exchange-correlation (XC) density functionals for finite systems, we propose an asymptotic correction scheme, wherein an exchange density functional whose functional derivative has the correct (-1/r) asymptote can be directly added to any semilocal density functional. In contrast to semilocal approximations, our resulting exchange kernel in reciprocal space exhibits the desirable singularity of the type O(-1/q^2) as q -> 0, which is a necessary feature for describing the excitonic effects in non-metallic solids. By applying this scheme to a popular semilocal density functional, PBE [J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)], the predictions of the properties that are sensitive to the asymptote are significantly improved, while the predictions of the properties that are insensitive to the asymptote remain essentially the same as PBE. Relative to the popular model XC potential scheme, our scheme is significantly superior for ground-state energies and related properties. In addition, without loss of accuracy, two closely related schemes are developed for the efficient treatment of large systems., Comment: 7 pages, 2 figures, 2 tables, supplementary material not included

Due to the severe self-interaction errors associated with some density functional approximations, conventional density functionals often fail to dissociate the hemibonded structure of water dimer radical cation (H2O)2+ into the correct fragments: H2O and H2O+. Consequently, the binding energy of the hemibonded structure (H2O)2+ is not well-defined. For a comprehensive comparison of different functionals for this system, we propose three criteria: (i) The binding energies, (ii) the relative energies between the conformers of the water dimer radical cation, and (iii) the dissociation curves predicted by different functionals. The long-range corrected (LC) double-hybrid functional, omegaB97X-2(LP) [J.-D. Chai and M. Head-Gordon, J. Chem. Phys., 2009, 131, 174105.], is shown to perform reasonably well based on these three criteria. Reasons that LC hybrid functionals generally work better than conventional density functionals for hemibonded systems are also explained in this work., Comment: 10 pages, 5 figures, 4 tables

Without the use of any empirical fitting to experimental or high-level ab initio data, we present a double-hybrid density functional approximation for the exchange-correlation energy, combining the exact Hartree-Fock exchange and second-order Moller-Plesset (MP2) correlation with the Perdew-Burke-Ernzerhof (PBE) functional. This functional, denoted as PBE0-2, is shown to be accurate for a wide range of applications, when compared with other functionals and the ab initio MP2 method. The qualitative failures of conventional density functional approximations, such as self-interaction error and noncovalent interaction error, are significantly reduced by PBE0-2., Comment: accepted for publication in Chem. Phys. Lett., 5 pages, 5 figures, 1 table, supplementary material not included

In contrast to the original Kohn-Sham (KS) formalism, we propose a density functional theory (DFT) with fractional orbital occupations for the study of ground states of many-electron systems, wherein strong static correlation is shown to be described. Even at the simplest level represented by the local density approximation (LDA), our resulting DFT-LDA is shown to improve upon KS-LDA for multi-reference systems, such as dissociation of H2 and N2, and twisted ethylene, while performing similarly to KS-LDA for single-reference systems, such as reaction energies and equilibrium geometries. Because of its computational efficiency (similar to KS-LDA), this DFT-LDA is applied to the study of the singlet-triplet energy gaps (ST gaps) of acenes, which are "challenging problems" for conventional electronic structure methods due to the presence of strong static correlation effects. Our calculated ST gaps are in good agreement with the existing experimental and high-level ab initio data. The ST gaps are shown to decrease monotonically with the increase of chain length, and become vanishingly small (within 0.1 kcal/mol) in the limit of an infinitely large polyacene. In addition, based on our calculated active orbital occupation numbers, the ground states for large acenes are shown to be polyradical singlets., Comment: accepted for publication in J. Chem. Phys., 45 pages, 21 figures, 3 tables, supplementary material not included

We propose a long-range corrected hybrid meta-GGA functional, based on a global hybrid meta-GGA functional, M05 [Y. Zhao, N. E. Schultz, and D. G. Truhlar, J. Chem. Phys. 123, 161103 (2005)], and empirical atom-atom dispersion corrections. Our resulting functional, omegaM05-D, is shown to be accurate for a very wide range of applications, such as thermochemistry, kinetics, noncovalent interactions, equilibrium geometries, frontier orbital energies, fundamental gaps, and excitation energies. In addition, we present three new databases, IP131 (131 ionization potentials), EA115 (115 electron affinities), and FG115 (115 fundamental gaps), consisting of experimental molecular geometries and accurate reference values, which will be useful in the assessment of the accuracy of density functional approximations., Comment: accepted for publication in J. Chem. Phys., 75 pages, 6 figures, supplementary material included

We systematically investigate the possible complex transition origin of electronic excitations of giant molecular systems by using the recently proposed QNTO analysis [J.-H. Li, J.-D. Chai, G. Y. Guo and M. Hayashi, Chem. Phys. Lett., 2011, 514, 362.] combined with long-range corrected TDDFT calculations. Thymine (Thy) related excitations of biomolecule B-DNA are then studied as examples, where the model systems have been constructed extracting from the perfect or a X-ray crystal (PDB code 3BSE) B-DNA structure with at least one Thy included. In the first part, we consider the systems composed of a core molecular segment (e.g. Thy, di-Thy) and a surrounding physical/chemical environment of interest (e.g. backbone, adjacent stacking nucleobases) and examine how the excitation properties of the core vary in response to the environment. We find that the orbitals contributed from DNA backbone and surrounding nucleobases often participate in a transition of Thy-related excitations affecting their composition, absorption energy, and oscillator strength. In the second part, we take into account geometrically induced variation of the excitation properties of various B-DNA segments, e.g. di-Thy, dTpdT etc., obtained from different sources (ideal and 3BSE). It is found that the transition origin of several Thy-related excitations of these segments is sensitive to slight conformational variations, suggesting that DNA with thermal motions in cells may from time to time exhibit very different photo-induced physical and/or chemical processes., Comment: Main Text+Supplemental Material

We show that the origin of electronic transitions of molecular many-body systems can be revealed by a quantified natural transition orbitals (QNTO) analysis and the electronic excitations of the total system can be mapped onto a standard orbitals set of a reference system. We further illustrate QNTO on molecular systems by studying the origin of electronic transitions of DNA moiety, thymine and thymidine. This QNTO analysis also allows us to assess the performance of various functionals used in time-dependent density functional response theory., Comment: Main Text+Supplemental Material; G09 reference corrected

In orbital-free density functional theory the kinetic potential (KP), the functional derivative of the kinetic energy density functional, appears in the Euler equation for the electron density and may be more amenable to simple approximations. We study properties of two solid-state systems, Al and Si, using two nonlocal KPs that gave good results for atoms. Very accurate results are found for Al, but results for Si are much less satisfactory, illustrating the general need for a better treatment of extended covalent systems. A different integration pathway in the KP formalism may prove useful in attacking this fundamental problem., Comment: Submitted to Chemical Physics Letters

In the density functional (DF) theory of Kohn and Sham, the kinetic energy of the ground state of a system of noninteracting electrons in a general external field is calculated using a set of orbitals. Orbital free methods attempt to calculate this directly from the electron density by approximating the universal but unknown kinetic energy density functional. However simple local approximations are inaccurate and it has proved very difficult to devise generally accurate nonlocal approximations. We focus instead on the kinetic potential, the functional derivative of the kinetic energy DF, which appears in the Euler equation for the electron density. We argue that the kinetic potential is more local and more amenable to simple physically motivated approximations in many relevant cases, and describe two pathways by which the value of the kinetic energy can be efficiently calculated. We propose two nonlocal orbital free kinetic potentials that reduce to known exact forms for both slowly varying and rapidly varying perturbations and also reproduce exact results for the linear response of the density of the homogeneous system to small perturbations. A simple and systematic approach for generating accurate and weak ab-initio local pseudopotentials which produce a smooth slowly varying valence component of the electron density is proposed for use in orbital free DF calculations of molecules and solids. The use of these local pseudopotentials further minimizes the possible errors from the kinetic potentials. Our theory yields results for the total energies and ionization energies of atoms, and for the shell structure in the atomic radial density profiles that are in very good agreement with calculations using the full Kohn-Sham theory., Comment: To be published in Phys. Rev. B

We review our previous work on the dynamic structure factor S(k,omega) of liquid Ge (l-Ge) at temperature T = 1250 K, and of amorphous Ge (a-Ge) at T = 300 K, using ab initio molecular dynamics [Phys. Rev. B67, 104205 (2003)]. The electronic energy is computed using density-functional theory, primarily in the generalized gradient approximation, together with a plane wave representation of the wave functions and ultra-soft pseudopotentials. We use a 64-atom cell with periodic boundary conditions, and calculate averages over runs of up to 16 ps. The calculated liquid S(k,omega) agrees qualitatively with that obtained by Hosokawa et al, using inelastic X-ray scattering. In a-Ge, we find that the calculated S(k,omega) is in qualitative agreement with that obtained experimentally by Maley et al. Our results suggest that the ab initio approach is sufficient to allow approximate calculations of S(k,omega) in both liquid and amorphous materials., Comment: Review article, 30 pages Submitted to "Progress in Chemical Physics Research", Nova Science Publishers, New York, 2004

We try to improve the Thomas-Fermi model for the total energy and electron density of atoms and molecules by directly modifying the Euler equation for the electron density, which we argue is less affected by nonlocal corrections. Here we consider the simplest such modification by adding a linear gradient term to the Euler equation. For atoms, the coefficient of the gradient term can be chosen so that the correct exponential decay constant far from the nucleus is obtained. This model then gives a much improved description of the electron density at smaller distances, yielding in particular a finite density at the nucleus that is in good qualitative agreement with exact results. The cusp condition differs from the exact value by a factor of two. Values for the total energy of atomic systems, obtained by coupling parameter integration of the densities given by the Euler equation, are about as accurate as those given by very best Thomas-Fermi-Weizsacker models, and the density is much more accurate. Possible connections to orbital-free methods for the kinetic-energy functional in density functional theory are discussed., Comment: 9 pages, 5 tables, and 6 figures. Accepted for J. Phys. Chem. B

We calculate the dynamic structure factor S(k,omega) of liquid Ge (l-Ge) at temperature T = 1250 K, and of amorphous Ge (a-Ge) at T = 300 K, using ab initio molecular dynamics. The electronic energy is computed using density-functional theory, primarily in the generalized gradient approximation, together with a plane wave representation of the wave functions and ultra-soft pseudopotentials. We use a 64-atom cell with periodic boundary conditions, and calculate averages over runs of up to 16 ps. The calculated liquid S(k,omega) agrees qualitatively with that obtained by Hosokawa et al, using inelastic X-ray scattering. In a-Ge, we find that the calculated S(k,omega) is in qualitative agreement with that obtained experimentally by Maley et al. Our results suggest that the ab initio approach is sufficient to allow approximate calculations of S(k,omega) in both liquid and amorphous materials., Comment: 31 pages and 8 figures. Accepted for Phys. Rev. B

We describe a simple model for calculating the zero-temperature superfluid density of Zn-doped YBa_2Cu_3O_{7-\delta} as a function of the fraction x of in-plane Cu atoms which are replaced by Zn. The basis of the calculation is a ``Swiss cheese'' picture of a single CuO_2 layer, in which a substitutional Zn impurity creates a normal region of area $\pi\xi_{ab}^2$ around it as originally suggested by Nachumi et al. Here $\xi_{ab}$ is the zero-temperature in-plane coherence length at x = 0. We use this picture to calculate the variation of the in-plane superfluid density with x at temperature T = 0, using both a numerical approach and an analytical approximation. For $\delta = 0.37$, if we use the value $\xi_{ab}$ = 18.3 angstrom, we find that the in-plane superfluid decreases with increasing x and vanishes near $x_c = 0.01$ in the analytical approximation, and near $x_c = 0.014$ in the numerical approach. $x_c$ is quite sensitive to $\xi_{ab}$, whose value is not widely agreed upon. The model also predicts a peak in the real part of the conductivity, Re$\sigma_e(\omega, x)$, at concentrations $x \sim x_c$, and low frequencies, and a variation of critical current density with x of the form $J_c(x) \propto n_{S,e}(x)^{7/4}$ near percolation, where $n_{S,e}(x)$ is the in-plane superfluid density., Comment: 19 pages including 6 figures, submitted to Physica C