Source author record

Jeng-Da Chai

Jeng-Da Chai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2020arXiv

Approximate density matrix functionals applied to hetero-atomic bond dissociation

A two-orbital two-electron diatomic model resembling LiH is used to investigate the differences between the exact Löwdin-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.

preprint2020arXiv

Combining density based dynamical correlation with a reduced density matrix strong correlation description

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 Π is employed to correct for double counting a correlation functional of density functional theory (DFT). The proposed ELS-DMΠDFT method with the extended Löwdin-Shull (ELS) 1RDM functional with dispersion and multibond (DM) corrections augmented with the Π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ΠDFT (CASSCF augmented with a Π based scaled DFT correlation correction) method. Furthermore, with the additional M-correction, ELS-DMΠDFT+M reproduces well the benchmark PEC of the N_2 molecule by Lie and Clementi.

preprint2020arXiv

Excitation Energies from Thermally-Assisted-Occupation Density Functional Theory: Theory and Computational Implementation

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Σ_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Σ_g^+$) and the first triplet state ($1^3Σ_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).

preprint2016arXiv

Effect of Li Adsorption on the Electronic and Hydrogen Storage Properties of Acenes: A Dispersion-Corrected TAO-DFT Study

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.

preprint2016arXiv

Electronic and Optical Properties of the Narrowest Armchair Graphene Nanoribbons Studied by Density Functional Methods

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.

preprint2016arXiv

Electronic Properties of Cyclacenes from TAO-DFT

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.

preprint2016arXiv

Role of Kekulé and Non-Kekulé Structures in the Radical Character of Alternant Polycyclic Aromatic Hydrocarbons: A TAO-DFT Study

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.

preprint2016arXiv

SCAN-based hybrid and double-hybrid density functionals from models without fitted parameters

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.

preprint2016arXiv

Short- and long-range corrected hybrid density functionals with the D3 dispersion corrections

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.

preprint2015arXiv

Assessment of Density Functional Methods for Exciton Binding Energies and Related Optoelectronic Properties

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.

preprint2015arXiv

Van der Waals interactions in rare-gas dimers: The role of interparticle interactions

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.

preprint2014arXiv

Assessment of asymptotically corrected model potential scheme for charge-transfer-like excitations in oligoacenes

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.

preprint2014arXiv

Thermally-assisted-occupation density functional theory with generalized-gradient approximations

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.

preprint2013arXiv

Asymptotic Correction Schemes for Semilocal Exchange-Correlation Functionals

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.

preprint2012arXiv

Density functional theory with fractional orbital occupations

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.

preprint2012arXiv

Long-Range Corrected Hybrid Density Functionals with Improved Dispersion Corrections

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.

preprint2012arXiv

Long-range corrected hybrid meta-generalized-gradient approximations with dispersion corrections

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.

preprint2012arXiv

Restoration of the Derivative Discontinuity in Kohn-Sham Density Functional Theory: An Efficient Scheme for Energy Gap Correction

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.

preprint2012arXiv

Seeking for reliable double-hybrid density functionals without fitting parameters: The PBE0-2 functional

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.

preprint2011arXiv

The Quantified NTO Analysis for the Electronic Excitations of Molecular Many-Body Systems

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.

preprint2009arXiv

Orbital-Free Density Functional Theory: Linear Scaling Methods for Kinetic Potentials, and Applications to Solid Al and S

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.