Hybrid functional pseudopotentials

Science.gov (United States)

Yang, Jing; Tan, Liang Z.; Rappe, Andrew M.

2018-02-01

The consistency between the exchange-correlation functional used in pseudopotential construction and in the actual density functional theory calculation is essential for the accurate prediction of fundamental properties of materials. However, routine hybrid density functional calculations at present still rely on generalized gradient approximation pseudopotentials due to the lack of hybrid functional pseudopotentials. Here, we present a scheme for generating hybrid functional pseudopotentials, and we analyze the importance of pseudopotential density functional consistency for hybrid functionals. For the PBE0 hybrid functional, we benchmark our pseudopotentials for structural parameters and fundamental electronic gaps of the Gaussian-2 (G2) molecular dataset and some simple solids. Our results show that using our PBE0 pseudopotentials in PBE0 calculations improves agreement with respect to all-electron calculations.

An open-source framework for analyzing N-electron dynamics. II. Hybrid density functional theory/configuration interaction methodology.

Science.gov (United States)

Hermann, Gunter; Pohl, Vincent; Tremblay, Jean Christophe

2017-10-30

In this contribution, we extend our framework for analyzing and visualizing correlated many-electron dynamics to non-variational, highly scalable electronic structure method. Specifically, an explicitly time-dependent electronic wave packet is written as a linear combination of N-electron wave functions at the configuration interaction singles (CIS) level, which are obtained from a reference time-dependent density functional theory (TDDFT) calculation. The procedure is implemented in the open-source Python program detCI@ORBKIT, which extends the capabilities of our recently published post-processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). From the output of standard quantum chemistry packages using atom-centered Gaussian-type basis functions, the framework exploits the multideterminental structure of the hybrid TDDFT/CIS wave packet to compute fundamental one-electron quantities such as difference electronic densities, transient electronic flux densities, and transition dipole moments. The hybrid scheme is benchmarked against wave function data for the laser-driven state selective excitation in LiH. It is shown that all features of the electron dynamics are in good quantitative agreement with the higher-level method provided a judicious choice of functional is made. Broadband excitation of a medium-sized organic chromophore further demonstrates the scalability of the method. In addition, the time-dependent flux densities unravel the mechanistic details of the simulated charge migration process at a glance. © 2017 Wiley Periodicals, Inc. © 2017 Wiley Periodicals, Inc.

Recent developments in LIBXC - A comprehensive library of functionals for density functional theory

Science.gov (United States)

Lehtola, Susi; Steigemann, Conrad; Oliveira, Micael J. T.; Marques, Miguel A. L.

2018-01-01

LIBXC is a library of exchange-correlation functionals for density-functional theory. We are concerned with semi-local functionals (or the semi-local part of hybrid functionals), namely local-density approximations, generalized-gradient approximations, and meta-generalized-gradient approximations. Currently we include around 400 functionals for the exchange, correlation, and the kinetic energy, spanning more than 50 years of research. Moreover, LIBXC is by now used by more than 20 codes, not only from the atomic, molecular, and solid-state physics, but also from the quantum chemistry communities.

Correction: General optimization procedure towards the design of a new family of minimal parameter spin-component-scaled double-hybrid density functional theory.

Science.gov (United States)

Roch, Loïc M; Baldridge, Kim K

2018-02-07

Correction for 'General optimization procedure towards the design of a new family of minimal parameter spin-component-scaled double-hybrid density functional theory' by Loïc M. Roch and Kim K. Baldridge, Phys. Chem. Chem. Phys., 2017, 19, 26191-26200.

Explicit polarization (X-Pol) potential using ab initio molecular orbital theory and density functional theory.

Science.gov (United States)

Song, Lingchun; Han, Jaebeom; Lin, Yen-lin; Xie, Wangshen; Gao, Jiali

2009-10-29

The explicit polarization (X-Pol) method has been examined using ab initio molecular orbital theory and density functional theory. The X-Pol potential was designed to provide a novel theoretical framework for developing next-generation force fields for biomolecular simulations. Importantly, the X-Pol potential is a general method, which can be employed with any level of electronic structure theory. The present study illustrates the implementation of the X-Pol method using ab initio Hartree-Fock theory and hybrid density functional theory. The computational results are illustrated by considering a set of bimolecular complexes of small organic molecules and ions with water. The computed interaction energies and hydrogen bond geometries are in good accord with CCSD(T) calculations and B3LYP/aug-cc-pVDZ optimizations.

Density functional theory

International Nuclear Information System (INIS)

Das, M.P.

1984-07-01

The state of the art of the density functional formalism (DFT) is reviewed. The theory is quantum statistical in nature; its simplest version is the well-known Thomas-Fermi theory. The DFT is a powerful formalism in which one can treat the effect of interactions in inhomogeneous systems. After some introductory material, the DFT is outlined from the two basic theorems, and various generalizations of the theorems appropriate to several physical situations are pointed out. Next, various approximations to the density functionals are presented and some practical schemes, discussed; the approximations include an electron gas of almost constant density and an electron gas of slowly varying density. Then applications of DFT in various diverse areas of physics (atomic systems, plasmas, liquids, nuclear matter) are mentioned, and its strengths and weaknesses are pointed out. In conclusion, more recent developments of DFT are indicated

General optimization procedure towards the design of a new family of minimal parameter spin-component-scaled double-hybrid density functional theory.

Science.gov (United States)

Roch, Loïc M; Baldridge, Kim K

2017-10-04

A general optimization procedure towards the development and implementation of a new family of minimal parameter spin-component-scaled double-hybrid (mSD) density functional theory (DFT) is presented. The nature of the proposed exchange-correlation functional establishes a methodology with minimal empiricism. This new family of double-hybrid (DH) density functionals is demonstrated using the PBEPBE functional, illustrating the optimization procedure to the mSD-PBEPBE method, and the performance characteristics shown for a set of non-covalent complexes covering a broad regime of weak interactions. With only two parameters, mSD-PBEPBE and its cost-effective counterpart, RI-mSD-PBEPBE, show a mean absolute error of ca. 0.4 kcal mol -1 averaged over 66 weak interacting systems. Following a successive 2D-grid refinement for a CBS extrapolation of the coefficients, the optimization procedure can be recommended for the design and implementation of a variety of additional DH methods using any of the plethora of currently available functionals.

Density functional theory: Foundations reviewed

Energy Technology Data Exchange (ETDEWEB)

Kryachko, Eugene S., E-mail: eugene.kryachko@ulg.ac.be [Bogolyubov Institute for Theoretical Physics, Kiev, 03680 (Ukraine); Ludeña, Eduardo V., E-mail: popluabe@yahoo.es [Centro de Química, Instituto Venezolano de Investigaciones Científicas, IVIC, Apartado 21827, Caracas 1020-A (Venezuela, Bolivarian Republic of); Prometheus Program, Senescyt (Ecuador); Grupo Ecuatoriano para el Estudio Experimental y Teórico de Nanosistemas, GETNano, USFQ, N104-E, Quito (Ecuador); Escuela Politécnica Superior del Litoral, ESPOL, Guayaquil (Ecuador)

2014-11-10

transformation (local-scaling of coordinates proceeds through density transformation) and so, because these functionals are constructed from prototype N-particle wavefunctions, the ensuing density functionals already have built-in N-representability conditions. This theory is presented in great detail with the purpose of illustrating an alternative way to HKS-DFT which could be used to improve the construction of HKS-DFT functionals. Let us clearly indicate, however, that although appealing from a theoretical point of view, the actual application of LS-DFT to large systems has not taken place mostly because of technical difficulties. Thus, our aim in introducing this theory is to foster a better understanding of its foundations with the hope that it may promote a cross-hybridization with the already existing approaches. Also, to complete our previous discussion on symmetry, in particular, spin-symmetry, we discuss this issue from the perspective of LS-DFT. Finally, in Section  6, we discuss dispersion molecular forces emphasizing their relevance to DFT approaches.

Density-functional theory for internal magnetic fields

Science.gov (United States)

Tellgren, Erik I.

2018-01-01

A density-functional theory is developed based on the Maxwell-Schrödinger equation with an internal magnetic field in addition to the external electromagnetic potentials. The basic variables of this theory are the electron density and the total magnetic field, which can equivalently be represented as a physical current density. Hence, the theory can be regarded as a physical current density-functional theory and an alternative to the paramagnetic current density-functional theory due to Vignale and Rasolt. The energy functional has strong enough convexity properties to allow a formulation that generalizes Lieb's convex analysis formulation of standard density-functional theory. Several variational principles as well as a Hohenberg-Kohn-like mapping between potentials and ground-state densities follow from the underlying convex structure. Moreover, the energy functional can be regarded as the result of a standard approximation technique (Moreau-Yosida regularization) applied to the conventional Schrödinger ground-state energy, which imposes limits on the maximum curvature of the energy (with respect to the magnetic field) and enables construction of a (Fréchet) differentiable universal density functional.

Quantal density functional theory

CERN Document Server

Sahni, Viraht

2016-01-01

This book deals with quantal density functional theory (QDFT) which is a time-dependent local effective potential theory of the electronic structure of matter. The treated time-independent QDFT constitutes a special case. In the 2nd edition, the theory is extended to include the presence of external magnetostatic fields. The theory is a description of matter based on the ‘quantal Newtonian’ first and second laws which is in terms of “classical” fields that pervade all space, and their quantal sources. The fields, which are explicitly defined, are separately representative of electron correlations due to the Pauli exclusion principle, Coulomb repulsion, correlation-kinetic, correlation-current-density, and correlation-magnetic effects. The book further describes Schrödinger theory from the new physical perspective of fields and quantal sources. It also describes traditional Hohenberg-Kohn-Sham DFT, and explains via QDFT the physics underlying the various energy functionals and functional derivatives o...

Multicomponent density functional theory embedding formulation

Energy Technology Data Exchange (ETDEWEB)

Culpitt, Tanner; Brorsen, Kurt R.; Pak, Michael V.; Hammes-Schiffer, Sharon, E-mail: shs3@illinois.edu [Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801 (United States)

2016-07-28

Multicomponent density functional theory (DFT) methods have been developed to treat two types of particles, such as electrons and nuclei, quantum mechanically at the same level. In the nuclear-electronic orbital (NEO) approach, all electrons and select nuclei, typically key protons, are treated quantum mechanically. For multicomponent DFT methods developed within the NEO framework, electron-proton correlation functionals based on explicitly correlated wavefunctions have been designed and used in conjunction with well-established electronic exchange-correlation functionals. Herein a general theory for multicomponent embedded DFT is developed to enable the accurate treatment of larger systems. In the general theory, the total electronic density is separated into two subsystem densities, denoted as regular and special, and different electron-proton correlation functionals are used for these two electronic densities. In the specific implementation, the special electron density is defined in terms of spatially localized Kohn-Sham electronic orbitals, and electron-proton correlation is included only for the special electron density. The electron-proton correlation functional depends on only the special electron density and the proton density, whereas the electronic exchange-correlation functional depends on the total electronic density. This scheme includes the essential electron-proton correlation, which is a relatively local effect, as well as the electronic exchange-correlation for the entire system. This multicomponent DFT-in-DFT embedding theory is applied to the HCN and FHF{sup −} molecules in conjunction with two different electron-proton correlation functionals and three different electronic exchange-correlation functionals. The results illustrate that this approach provides qualitatively accurate nuclear densities in a computationally tractable manner. The general theory is also easily extended to other types of partitioning schemes for multicomponent systems.

Theoretical investigation of pressure-induced structural transitions in americium using GGA+U and hybrid density functional theory methods

DEFF Research Database (Denmark)

Verma, Ashok K.; Modak, P.; Sharma, Surinder M.

2013-01-01

First-principles calculations have been performed for americium (Am) metal using the generalized gradient approximation + orbital-dependent onsite Coulomb repulsion via Hubbard interaction (GGA+U) and hybrid density functional theory (HYB-DFT) methods to investigate various ground state properties......-I to Am-II transition. Good agreement was found between calculated and experimental equations of states for all phases, but the first three phases need larger U (α) parameters (where α represents the fraction of Hartree-Fock exchange energy replacing the DFT exchange energy) than the fourth phase in order...

Multiconfiguration Pair-Density Functional Theory Outperforms Kohn-Sham Density Functional Theory and Multireference Perturbation Theory for Ground-State and Excited-State Charge Transfer.

Science.gov (United States)

Ghosh, Soumen; Sonnenberger, Andrew L; Hoyer, Chad E; Truhlar, Donald G; Gagliardi, Laura

2015-08-11

The correct description of charge transfer in ground and excited states is very important for molecular interactions, photochemistry, electrochemistry, and charge transport, but it is very challenging for Kohn-Sham (KS) density functional theory (DFT). KS-DFT exchange-correlation functionals without nonlocal exchange fail to describe both ground- and excited-state charge transfer properly. We have recently proposed a theory called multiconfiguration pair-density functional theory (MC-PDFT), which is based on a combination of multiconfiguration wave function theory with a new type of density functional called an on-top density functional. Here we have used MC-PDFT to study challenging ground- and excited-state charge-transfer processes by using on-top density functionals obtained by translating KS exchange-correlation functionals. For ground-state charge transfer, MC-PDFT performs better than either the PBE exchange-correlation functional or CASPT2 wave function theory. For excited-state charge transfer, MC-PDFT (unlike KS-DFT) shows qualitatively correct behavior at long-range with great improvement in predicted excitation energies.

Spin theory of the density functional: reduced matrices and density functions

International Nuclear Information System (INIS)

Pavlov, R.; Delchev, Y.; Pavlova, K.; Maruani, J.

1993-01-01

Expressions for the reduced matrices and density functions of N-fermion systems of arbitrary order s (1<=s<=N) are derived within the frame of rigorous spin approach to the density functional theory (DFT). Using the local-scale transformation method and taking into account the particle spin it is shown that the reduced matrices and density functions are functionals of the total one-fermion density. Similar dependence is found for the distribution density of s-particle aggregates. Generalization and applicability of DFT to the case of s-particle ensembles and aggregates is discussed. 14 refs

Magnetic fields and density functional theory

Energy Technology Data Exchange (ETDEWEB)

Salsbury Jr., Freddie [Univ. of California, Berkeley, CA (United States)

1999-02-01

A major focus of this dissertation is the development of functionals for the magnetic susceptibility and the chemical shielding within the context of magnetic field density functional theory (BDFT). These functionals depend on the electron density in the absence of the field, which is unlike any other treatment of these responses. There have been several advances made within this theory. The first of which is the development of local density functionals for chemical shieldings and magnetic susceptibilities. There are the first such functionals ever proposed. These parameters have been studied by constructing functionals for the current density and then using the Biot-Savart equations to obtain the responses. In order to examine the advantages and disadvantages of the local functionals, they were tested numerically on some small molecules.

Magnetic fields and density functional theory

International Nuclear Information System (INIS)

Salsbury, Freddie Jr.

1999-01-01

A major focus of this dissertation is the development of functionals for the magnetic susceptibility and the chemical shielding within the context of magnetic field density functional theory (BDFT). These functionals depend on the electron density in the absence of the field, which is unlike any other treatment of these responses. There have been several advances made within this theory. The first of which is the development of local density functionals for chemical shieldings and magnetic susceptibilities. There are the first such functionals ever proposed. These parameters have been studied by constructing functionals for the current density and then using the Biot-Savart equations to obtain the responses. In order to examine the advantages and disadvantages of the local functionals, they were tested numerically on some small molecules

Orbital functionals in density-matrix- and current-density-functional theory

Energy Technology Data Exchange (ETDEWEB)

Helbig, N

2006-05-15

Density-Functional Theory (DFT), although widely used and very successful in the calculation of several observables, fails to correctly describe strongly correlated materials. In the first part of this work we, therefore, introduce reduced-densitymatrix- functional theory (RDMFT) which is one possible way to treat electron correlation beyond DFT. Within this theory the one-body reduced density matrix (1- RDM) is used as the basic variable. Our main interest is the calculation of the fundamental gap which proves very problematic within DFT. In order to calculate the fundamental gap we generalize RDMFT to fractional particle numbers M by describing the system as an ensemble of an N and an N+1 particle system (with N{<=}M{<=}N+1). For each fixed particle number, M, the total energy is minimized with respect to the natural orbitals and their occupation numbers. This leads to the total energy as a function of M. The derivative of this function with respect to the particle number has a discontinuity at integer particle number which is identical to the gap. In addition, we investigate the necessary and sufficient conditions for the 1- RDM of a system with fractional particle number to be N-representable. Numerical results are presented for alkali atoms, small molecules, and periodic systems. Another problem within DFT is the description of non-relativistic many-electron systems in the presence of magnetic fields. It requires the paramagnetic current density and the spin magnetization to be used as basic variables besides the electron density. However, electron-gas-based functionals of current-spin-density-functional Theory (CSDFT) exhibit derivative discontinuities as a function of the magnetic field whenever a new Landau level is occupied, which makes them difficult to use in practice. Since the appearance of Landau levels is, intrinsically, an orbital effect it is appealing to use orbital-dependent functionals. We have developed a CSDFT version of the optimized

Density functional theory of nuclei

International Nuclear Information System (INIS)

Terasaki, Jun

2008-01-01

The density functional theory of nuclei has come to draw attention of scientists in the field of nuclear structure because the theory is expected to provide reliable numerical data in wide range on the nuclear chart. This article is organized to present an overview of the theory to the people engaged in the theory of other fields as well as those people in the nuclear physics experiments. At first, the outline of the density functional theory widely used in the electronic systems (condensed matter, atoms, and molecules) was described starting from the Kohn-Sham equation derived on the variational principle. Then the theory used in the field of nuclear physics was presented. Hartree-Fock and Hartree-Fock-Bogolyubov approximation by using Skyrme interaction was explained. Comparison of the results of calculations and experiments of binding energies and ground state mean square charge radii of some magic number nuclei were shown. The similarity and dissimilarity between the two streams were summarized. Finally the activities of the international project of Universal Nuclear Energy Density Functional (UNEDF) which was started recently lead by US scientist was reported. This project is programmed for five years. One of the applications of the project is the calculation of the neutron capture cross section of nuclei on the r-process, which is absolutely necessary for the nucleosynthesis research. (S. Funahashi)

Spin-Density Functionals from Current-Density Functional Theory and Vice Versa: A Road towards New Approximations

International Nuclear Information System (INIS)

Capelle, K.; Gross, E.

1997-01-01

It is shown that the exchange-correlation functional of spin-density functional theory is identical, on a certain set of densities, with the exchange-correlation functional of current-density functional theory. This rigorous connection is used to construct new approximations of the exchange-correlation functionals. These include a conceptually new generalized-gradient spin-density functional and a nonlocal current-density functional. copyright 1997 The American Physical Society

Structure factor of blends of solvent-free nanoparticle-organic hybrid materials: density-functional theory and small angle X-ray scattering.

Science.gov (United States)

Yu, Hsiu-Yu; Srivastava, Samanvaya; Archer, Lynden A; Koch, Donald L

2014-12-07

We investigate the static structure factor S(q) of solvent-free nanoparticle-organic hybrid materials consisting of silica nanocores and space-filling polyethylene glycol coronas using a density-functional theory and small angle X-ray scattering measurements. The theory considers a bidisperse suspension of hard spheres with different radii and tethered bead-spring oligomers with different grafting densities to approximate the polydispersity effects in experiments. The experimental systems studied include pure samples with different silica core volume fractions and the associated mean corona grafting densities, and blends with different mixing ratios of the pure samples, in order to introduce varying polydispersity of corona grafting density. Our scattering experiments and theory show that, compared to the hard-sphere suspension with the same core volume fraction, S(q) for pure samples exhibit both substantially smaller values at small q and stronger particle correlations corresponding to a larger effective hard core at large q, indicating that the tethered incompressible oligomers enforce a more uniform particle distribution, and the densely grafted brush gives rise to an additional exclusionary effect between the nanoparticles. According to the theory, polydispersity in the oligomer grafting density controls the deviation of S(q) from the monodisperse system at smaller q, and the interplay of the enhanced effective core size and the entropic attraction among the particles is responsible for complex variations in the particle correlations at larger q. The successful comparison between the predictions and the measurements for the blends further suggests that S(q) can be used to assess the uniformity of grafting density in polymer-grafted nanoparticle materials.

Towards improved local hybrid functionals by calibration of exchange-energy densities

International Nuclear Information System (INIS)

Arbuznikov, Alexei V.; Kaupp, Martin

2014-01-01

A new approach for the calibration of (semi-)local and exact exchange-energy densities in the context of local hybrid functionals is reported. The calibration functions are derived from only the electron density and its spatial derivatives, avoiding spatial derivatives of the exact-exchange energy density or other computationally unfavorable contributions. The calibration functions fulfill the seven more important out of nine known exact constraints. It is shown that calibration improves substantially the definition of a non-dynamical correlation energy term for generalized gradient approximation (GGA)-based local hybrids. Moreover, gauge artifacts in the potential-energy curves of noble-gas dimers may be corrected by calibration. The developed calibration functions are then evaluated for a large range of energy-related properties (atomization energies, reaction barriers, ionization potentials, electron affinities, and total atomic energies) of three sets of local hybrids, using a simple one-parameter local-mixing. The functionals are based on (a) local spin-density approximation (LSDA) or (b) Perdew-Burke-Ernzerhof (PBE) exchange and correlation, and on (c) Becke-88 (B88) exchange and Lee-Yang-Parr (LYP) correlation. While the uncalibrated GGA-based functionals usually provide very poor thermochemical data, calibration allows a dramatic improvement, accompanied by only a small deterioration of reaction barriers. In particular, an optimized BLYP-based local-hybrid functional has been found that is a substantial improvement over the underlying global hybrids, as well as over previously reported LSDA-based local hybrids. It is expected that the present calibration approach will pave the way towards new generations of more accurate hyper-GGA functionals based on a local mixing of exchange-energy densities

Current Density Functional Theory Using Meta-Generalized Gradient Exchange-Correlation Functionals.

Science.gov (United States)

Furness, James W; Verbeke, Joachim; Tellgren, Erik I; Stopkowicz, Stella; Ekström, Ulf; Helgaker, Trygve; Teale, Andrew M

2015-09-08

We present the self-consistent implementation of current-dependent (hybrid) meta-generalized gradient approximation (mGGA) density functionals using London atomic orbitals. A previously proposed generalized kinetic energy density is utilized to implement mGGAs in the framework of Kohn-Sham current density functional theory (KS-CDFT). A unique feature of the nonperturbative implementation of these functionals is the ability to seamlessly explore a wide range of magnetic fields up to 1 au (∼235 kT) in strength. CDFT functionals based on the TPSS and B98 forms are investigated, and their performance is assessed by comparison with accurate coupled-cluster singles, doubles, and perturbative triples (CCSD(T)) data. In the weak field regime, magnetic properties such as magnetizabilities and nuclear magnetic resonance shielding constants show modest but systematic improvements over generalized gradient approximations (GGA). However, in the strong field regime, the mGGA-based forms lead to a significantly improved description of the recently proposed perpendicular paramagnetic bonding mechanism, comparing well with CCSD(T) data. In contrast to functionals based on the vorticity, these forms are found to be numerically stable, and their accuracy at high field suggests that the extension of mGGAs to CDFT via the generalized kinetic energy density should provide a useful starting point for further development of CDFT approximations.

The use of perturbation theory in density-functional theory

International Nuclear Information System (INIS)

Goerling, A.

1996-01-01

Perturbation theory with respect to the electron-electron interaction leads to expressions for the exchange and correlation energies and potentials in terms of Kohn-Sham orbitals and Kohn-Sham eigenvalues. An exact open-quote exchange-only close-quote procedure for solids is introduced. Results for several semiconductors are presented. Perturbation theory expansions for the hardness of molecules and the bad gap of solids are given. Density-functional exchange and correlation energies for excited states are defined and a perturbation theory based Kohn-Sham formalism to treat excited states within density-functional theory is introduced

Dynamical density functional theory for dense atomic liquids

International Nuclear Information System (INIS)

Archer, A J

2006-01-01

Starting from Newton's equations of motion, we derive a dynamical density functional theory (DDFT) applicable to atomic liquids. The theory has the feature that it requires as input the Helmholtz free energy functional from equilibrium density functional theory. This means that, given a reliable equilibrium free energy functional, the correct equilibrium fluid density profile is guaranteed. We show that when the isothermal compressibility is small, the DDFT generates the correct value for the speed of sound in a dense liquid. We also interpret the theory as a dynamical equation for a coarse grained fluid density and show that the theory can be used (making further approximations) to derive the standard mode coupling theory that is used to describe the glass transition. The present theory should provide a useful starting point for describing the dynamics of inhomogeneous atomic fluids

Improving Rydberg Excitations within Time-Dependent Density Functional Theory with Generalized Gradient Approximations: The Exchange-Enhancement-for-Large-Gradient Scheme.

Science.gov (United States)

Li, Shaohong L; Truhlar, Donald G

2015-07-14

Time-dependent density functional theory (TDDFT) with conventional local and hybrid functionals such as the local and hybrid generalized gradient approximations (GGA) seriously underestimates the excitation energies of Rydberg states, which limits its usefulness for applications such as spectroscopy and photochemistry. We present here a scheme that modifies the exchange-enhancement factor to improve GGA functionals for Rydberg excitations within the TDDFT framework while retaining their accuracy for valence excitations and for the thermochemical energetics calculated by ground-state density functional theory. The scheme is applied to a popular hybrid GGA functional and tested on data sets of valence and Rydberg excitations and atomization energies, and the results are encouraging. The scheme is simple and flexible. It can be used to correct existing functionals, and it can also be used as a strategy for the development of new functionals.

Nonlinear electronic excitations in crystalline solids using meta-generalized gradient approximation and hybrid functional in time-dependent density functional theory

Energy Technology Data Exchange (ETDEWEB)

Sato, Shunsuke A. [Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571 (Japan); Taniguchi, Yasutaka [Center for Computational Science, University of Tsukuba, Tsukuba 305-8571 (Japan); Department of Medical and General Sciences, Nihon Institute of Medical Science, 1276 Shimogawara, Moroyama-Machi, Iruma-Gun, Saitama 350-0435 (Japan); Shinohara, Yasushi [Max Planck Institute of Microstructure Physics, 06120 Halle (Germany); Yabana, Kazuhiro [Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571 (Japan); Center for Computational Science, University of Tsukuba, Tsukuba 305-8571 (Japan)

2015-12-14

We develop methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta-generalized gradient approximation was proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional was proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt the predictor-corrector step for a stable time evolution. We have developed a method to evaluate electronic excitation energy without referring to the energy functional which is unknown for the TB-mBJ potential. For the HSE functional, we have developed a method for the operation of the Fock-like term in Fourier space to facilitate efficient use of massive parallel computers equipped with graphic processing units. We compare electronic excitations in silicon and germanium induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: they are close to each other using TB-mBJ and HSE and are much smaller in LDA. At high laser intensities close to the damage threshold, electronic excitation energies do not differ much among the three cases.

The problem of the universal density functional and the density matrix functional theory

International Nuclear Information System (INIS)

Bobrov, V. B.; Trigger, S. A.

2013-01-01

The analysis in this paper shows that the Hohenberg-Kohn theorem is the constellation of two statements: (i) the mathematically rigorous Hohenberg-Kohn lemma, which demonstrates that the same ground-state density cannot correspond to two different potentials of an external field, and (ii) the hypothesis of the existence of the universal density functional. Based on the obtained explicit expression for the nonrel-ativistic particle energy in a local external field, we prove that the energy of the system of more than two non-interacting electrons cannot be a functional of the inhomogeneous density. This result is generalized to the system of interacting electrons. It means that the Hohenberg-Kohn lemma cannot provide justification of the universal density functional for fermions. At the same time, statements of the density functional theory remain valid when considering any number of noninteracting ground-state bosons due to the Bose condensation effect. In the framework of the density matrix functional theory, the hypothesis of the existence of the universal density matrix functional corresponds to the cases of noninteracting particles and to interaction in the Hartree-Fock approximation.

Structure factor of blends of solvent-free nanoparticle–organic hybrid materials: density-functional theory and small angle X-ray scattering

KAUST Repository

Yu, Hsiu-Yu

2014-09-15

© the Partner Organisations 2014. We investigate the static structure factor S(q) of solvent-free nanoparticle-organic hybrid materials consisting of silica nanocores and space-filling polyethylene glycol coronas using a density-functional theory and small angle X-ray scattering measurements. The theory considers a bidisperse suspension of hard spheres with different radii and tethered bead-spring oligomers with different grafting densities to approximate the polydispersity effects in experiments. The experimental systems studied include pure samples with different silica core volume fractions and the associated mean corona grafting densities, and blends with different mixing ratios of the pure samples, in order to introduce varying polydispersity of corona grafting density. Our scattering experiments and theory show that, compared to the hard-sphere suspension with the same core volume fraction, S(q) for pure samples exhibit both substantially smaller values at small q and stronger particle correlations corresponding to a larger effective hard core at large q, indicating that the tethered incompressible oligomers enforce a more uniform particle distribution, and the densely grafted brush gives rise to an additional exclusionary effect between the nanoparticles. According to the theory, polydispersity in the oligomer grafting density controls the deviation of S(q) from the monodisperse system at smaller q, and the interplay of the enhanced effective core size and the entropic attraction among the particles is responsible for complex variations in the particle correlations at larger q. The successful comparison between the predictions and the measurements for the blends further suggests that S(q) can be used to assess the uniformity of grafting density in polymer-grafted nanoparticle materials. This journal is

density functional theory approach

Indian Academy of Sciences (India)

YOGESH ERANDE

2017-07-27

Jul 27, 2017 ... a key role in all optical switching devices, since their optical properties can be .... optimized in the gas phase using Density Functional Theory. (DFT).39 The ...... The Mediation of Electrostatic Effects by Sol- vents J. Am. Chem.

Excited-state density functional theory

International Nuclear Information System (INIS)

Harbola, Manoj K; Hemanadhan, M; Shamim, Md; Samal, P

2012-01-01

Starting with a brief introduction to excited-state density functional theory, we present our method of constructing modified local density approximated (MLDA) energy functionals for the excited states. We show that these functionals give accurate results for kinetic energy and exchange energy compared to the ground state LDA functionals. Further, with the inclusion of GGA correction, highly accurate total energies for excited states are obtained. We conclude with a brief discussion on the further direction of research that include the construction of correlation energy functional and exchange potential for excited states.

Uniform magnetic fields in density-functional theory

Science.gov (United States)

Tellgren, Erik I.; Laestadius, Andre; Helgaker, Trygve; Kvaal, Simen; Teale, Andrew M.

2018-01-01

We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional density functional theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term linear vector potential-DFT (LDFT), the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formulation and a convex formulation in terms of Legendre-Fenchel transformations are constructed. Many theoretical issues in CDFT find simplified analogs in LDFT. We prove results concerning N-representability, Hohenberg-Kohn-like mappings, existence of minimizers in the constrained-search expression, and a restricted analog to gauge invariance. The issue of additivity of the energy over non-interacting subsystems, which is qualitatively different in LDFT and CDFT, is also discussed.

Density functional theory and parallel processing

International Nuclear Information System (INIS)

Ward, R.C.; Geist, G.A.; Butler, W.H.

1987-01-01

The authors demonstrate a method for obtaining the ground state energies and charge densities of a system of atoms described within density functional theory using simulated annealing on a parallel computer

Hole trapping at Al impurities in silica: A challenge for density functional theories

DEFF Research Database (Denmark)

Lægsgaard, Jesper; Stokbro, Kurt

2001-01-01

The atomic geometry and electronic structure around a neutral substitutional Al impurity in silica is investigated using either the unrestricted Hartree-Fock (UHF) approximation, or Beckes three-parameter hybrid functional (B3LYP). It is found that the B3LYP functional fails to describe...... the structural distortions around the Al impurity, while the UHF results are consistent with experimental information. We argue that the failure of the B3LYP functional is caused by the incomplete self-interaction cancellation usually present in density functional theories....

Local electric dipole moments for periodic systems via density functional theory embedding.

Science.gov (United States)

Luber, Sandra

2014-12-21

We describe a novel approach for the calculation of local electric dipole moments for periodic systems. Since the position operator is ill-defined in periodic systems, maximally localized Wannier functions based on the Berry-phase approach are usually employed for the evaluation of local contributions to the total electric dipole moment of the system. We propose an alternative approach: within a subsystem-density functional theory based embedding scheme, subset electric dipole moments are derived without any additional localization procedure, both for hybrid and non-hybrid exchange-correlation functionals. This opens the way to a computationally efficient evaluation of local electric dipole moments in (molecular) periodic systems as well as their rigorous splitting into atomic electric dipole moments. As examples, Infrared spectra of liquid ethylene carbonate and dimethyl carbonate are presented, which are commonly employed as solvents in Lithium ion batteries.

Local electric dipole moments for periodic systems via density functional theory embedding

Energy Technology Data Exchange (ETDEWEB)

Luber, Sandra, E-mail: sandra.luber@chem.uzh.ch [Institut für Chemie, Universität Zürich, Winterthurerstrasse 190, 8057 Zürich (Switzerland)

2014-12-21

We describe a novel approach for the calculation of local electric dipole moments for periodic systems. Since the position operator is ill-defined in periodic systems, maximally localized Wannier functions based on the Berry-phase approach are usually employed for the evaluation of local contributions to the total electric dipole moment of the system. We propose an alternative approach: within a subsystem-density functional theory based embedding scheme, subset electric dipole moments are derived without any additional localization procedure, both for hybrid and non-hybrid exchange–correlation functionals. This opens the way to a computationally efficient evaluation of local electric dipole moments in (molecular) periodic systems as well as their rigorous splitting into atomic electric dipole moments. As examples, Infrared spectra of liquid ethylene carbonate and dimethyl carbonate are presented, which are commonly employed as solvents in Lithium ion batteries.

Characterization of the nitrogen split interstitial defect in wurtzite aluminum nitride using density functional theory

International Nuclear Information System (INIS)

Szállás, A.; Szász, K.; Trinh, X. T.; Son, N. T.; Janzén, E.; Gali, A.

2014-01-01

We carried out Heyd-Scuseria-Ernzerhof hybrid density functional theory plane wave supercell calculations in wurtzite aluminum nitride in order to characterize the geometry, formation energies, transition levels, and hyperfine tensors of the nitrogen split interstitial defect. The calculated hyperfine tensors may provide useful fingerprint of this defect for electron paramagnetic resonance measurement.

Characterization of the nitrogen split interstitial defect in wurtzite aluminum nitride using density functional theory

Energy Technology Data Exchange (ETDEWEB)

Szállás, A., E-mail: szallas.attila@wigner.mta.hu [Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, Hungarian Academy of Sciences, P.O. Box 49, H-1525 Budapest (Hungary); Szász, K. [Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, Hungarian Academy of Sciences, P.O. Box 49, H-1525 Budapest (Hungary); Institute of Physics, Eötvös University, Pázmány Péter sétány 1/A, H-1117 Budapest (Hungary); Trinh, X. T.; Son, N. T.; Janzén, E. [Department of Physics, Chemistry and Biology, Linköping University, SE-581 83 Linköping (Sweden); Gali, A., E-mail: gali.adam@wigner.mta.hu [Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, Hungarian Academy of Sciences, P.O. Box 49, H-1525 Budapest (Hungary); Department of Atomic Physics, Budapest University of Technology and Economics, Budafoki út 8, H-1111 Budapest (Hungary)

2014-09-21

We carried out Heyd-Scuseria-Ernzerhof hybrid density functional theory plane wave supercell calculations in wurtzite aluminum nitride in order to characterize the geometry, formation energies, transition levels, and hyperfine tensors of the nitrogen split interstitial defect. The calculated hyperfine tensors may provide useful fingerprint of this defect for electron paramagnetic resonance measurement.

A Safari Through Density Functional Theory

Science.gov (United States)

Dreizler, Reiner M.; Lüdde, Cora S.

Density functional theory is widely used to treat quantum many body problems in many areas of physics and related fields. A brief survey of this method covering foundations, functionals and applications is presented here.

A density functional theory study of the influence of exchange-correlation functionals on the properties of FeAs.

Science.gov (United States)

Griffin, Sinéad M; Spaldin, Nicola A

2017-06-01

We use density functional theory within the local density approximation (LDA), LDA  +  U, generalised gradient approximation (GGA), GGA  +  U, and hybrid-functional methods to calculate the properties of iron monoarsenide. FeAs, which forms in the MnP structure, is of current interest for potential spintronic applications as well as being the parent compound for the pnictide superconductors. We compare the calculated structural, magnetic and electronic properties obtained using the different functionals to each other and to experiment, and investigate the origin of a recently reported magnetic spiral. Our results indicate the appropriateness or otherwise of the various functionals for describing FeAs and the related Fe-pnictide superconductors.

Long-range corrected density functional theory with accelerated Hartree-Fock exchange integration using a two-Gaussian operator [LC-ωPBE(2Gau)].

Science.gov (United States)

Song, Jong-Won; Hirao, Kimihiko

2015-10-14

Since the advent of hybrid functional in 1993, it has become a main quantum chemical tool for the calculation of energies and properties of molecular systems. Following the introduction of long-range corrected hybrid scheme for density functional theory a decade later, the applicability of the hybrid functional has been further amplified due to the resulting increased performance on orbital energy, excitation energy, non-linear optical property, barrier height, and so on. Nevertheless, the high cost associated with the evaluation of Hartree-Fock (HF) exchange integrals remains a bottleneck for the broader and more active applications of hybrid functionals to large molecular and periodic systems. Here, we propose a very simple yet efficient method for the computation of long-range corrected hybrid scheme. It uses a modified two-Gaussian attenuating operator instead of the error function for the long-range HF exchange integral. As a result, the two-Gaussian HF operator, which mimics the shape of the error function operator, reduces computational time dramatically (e.g., about 14 times acceleration in C diamond calculation using periodic boundary condition) and enables lower scaling with system size, while maintaining the improved features of the long-range corrected density functional theory.

Global hybrids from the semiclassical atom theory satisfying the local density linear response.

Science.gov (United States)

Fabiano, Eduardo; Constantin, Lucian A; Cortona, Pietro; Della Sala, Fabio

2015-01-13

We propose global hybrid approximations of the exchange-correlation (XC) energy functional which reproduce well the modified fourth-order gradient expansion of the exchange energy in the semiclassical limit of many-electron neutral atoms and recover the full local density approximation (LDA) linear response. These XC functionals represent the hybrid versions of the APBE functional [Phys. Rev. Lett. 2011, 106, 186406] yet employing an additional correlation functional which uses the localization concept of the correlation energy density to improve the compatibility with the Hartree-Fock exchange as well as the coupling-constant-resolved XC potential energy. Broad energetic and structural testing, including thermochemistry and geometry, transition metal complexes, noncovalent interactions, gold clusters and small gold-molecule interfaces, as well as an analysis of the hybrid parameters, show that our construction is quite robust. In particular, our testing shows that the resulting hybrid, including 20% of Hartree-Fock exchange and named hAPBE, performs remarkably well for a broad palette of systems and properties, being generally better than popular hybrids (PBE0 and B3LYP). Semiempirical dispersion corrections are also provided.

Multicomponent density-functional theory for time-dependent systems

NARCIS (Netherlands)

Butriy, O.; Ebadi, H.; de Boeij, P. L.; van Leeuwen, R.; Gross, E. K. U.

2007-01-01

We derive the basic formalism of density functional theory for time-dependent electron-nuclear systems. The basic variables of this theory are the electron density in body-fixed frame coordinates and the diagonal of the nuclear N-body density matrix. The body-fixed frame transformation is carried

Density functional theory in quantum chemistry

CERN Document Server

Tsuneda, Takao

2014-01-01

This book examines density functional theory based on the foundation of quantum chemistry. Unconventional in approach, it reviews basic concepts, then describes the physical meanings of state-of-the-art exchange-correlation functionals and their corrections.

Reduced density matrix functional theory via a wave function based approach

Energy Technology Data Exchange (ETDEWEB)

Schade, Robert; Bloechl, Peter [Institute for Theoretical Physics, Clausthal University of Technology, Clausthal (Germany); Pruschke, Thomas [Institute for Theoretical Physics, University of Goettingen, Goettingen (Germany)

2016-07-01

We propose a new method for the calculation of the electronic and atomic structure of correlated electron systems based on reduced density matrix functional theory (rDMFT). The density-matrix functional is evaluated on the fly using Levy's constrained search formalism. The present implementation rests on a local approximation of the interaction reminiscent to that of dynamical mean field theory (DMFT). We focus here on additional approximations to the exact density-matrix functional in the local approximation and evaluate their performance.

Spectral function from Reduced Density Matrix Functional Theory

Science.gov (United States)

Romaniello, Pina; di Sabatino, Stefano; Berger, Jan A.; Reining, Lucia

2015-03-01

In this work we focus on the calculation of the spectral function, which determines, for example, photoemission spectra, from reduced density matrix functional theory. Starting from its definition in terms of the one-body Green's function we derive an expression for the spectral function that depends on the natural occupation numbers and on an effective energy which accounts for all the charged excitations. This effective energy depends on the two-body as well as higher-order density matrices. Various approximations to this expression are explored by using the exactly solvable Hubbard chains.

The force distribution probability function for simple fluids by density functional theory.

Science.gov (United States)

Rickayzen, G; Heyes, D M

2013-02-28

Classical density functional theory (DFT) is used to derive a formula for the probability density distribution function, P(F), and probability distribution function, W(F), for simple fluids, where F is the net force on a particle. The final formula for P(F) ∝ exp(-AF(2)), where A depends on the fluid density, the temperature, and the Fourier transform of the pair potential. The form of the DFT theory used is only applicable to bounded potential fluids. When combined with the hypernetted chain closure of the Ornstein-Zernike equation, the DFT theory for W(F) agrees with molecular dynamics computer simulations for the Gaussian and bounded soft sphere at high density. The Gaussian form for P(F) is still accurate at lower densities (but not too low density) for the two potentials, but with a smaller value for the constant, A, than that predicted by the DFT theory.

Local density approximation for exchange in excited-state density functional theory

OpenAIRE

Harbola, Manoj K.; Samal, Prasanjit

2004-01-01

Local density approximation for the exchange energy is made for treatment of excited-states in density-functional theory. It is shown that taking care of the state-dependence of the LDA exchange energy functional leads to accurate excitation energies.

Hybrid density functional theory study of Cu(In1−xGaxSe2 band structure for solar cell application

Directory of Open Access Journals (Sweden)

Xu-Dong Chen

2014-08-01

Full Text Available Cu(In1−xGaxSe2 (CIGS alloy based thin film photovoltaic solar cells have attracted more and more attention due to its large optical absorption coefficient, long term stability, low cost and high efficiency. However, the previous theoretical investigation of this material with first principle calculation cannot fulfill the requirement of experimental development, especially the accurate description of band structure and density of states. In this work, we use first principle calculation based on hybrid density functional theory to investigate the feature of CIGS, with B3LYP applied in the CuIn1−xGaxSe2 stimulation of the band structure and density of states. We report the simulation of the lattice parameter, band gap and chemical composition. The band gaps of CuGaSe2, CuIn0.25Ga0.75Se2, CuIn0.5Ga0.5Se2, CuIn0.75Ga0.25Se2 and CuInSe2 are obtained as 1.568 eV, 1.445 eV, 1.416 eV, 1.275 eV and 1.205 eV according to our calculation, which agree well with the available experimental values. The band structure of CIGS is also in accordance with the current theory.

Quantal density functional theory. 2. ed.

International Nuclear Information System (INIS)

Sahni, Viraht

2016-01-01

This book is on quantal density functional theory (QDFT) which is a time-dependent local effective potential theory of the electronic structure of matter. The time-independent QDFT constitutes a special case. The 2 nd edition describes the further development of the theory, and extends it to include the presence of an external magnetostatic field. The theory is based on the 'quantal Newtonian' second and first laws for the individual electron. These laws are in terms of 'classical' fields that pervade all space, and their quantal sources. The fields are separately representative of the electron correlations that must be accounted for in local potential theory. Recent developments show that irrespective of the type of external field the electrons are subject to, the only correlations beyond those due to the Pauli exclusion principle and Coulomb repulsion that need be considered are solely of the correlation-kinetic effects. Foundational to QDFT, the book describes Schroedinger theory from the new perspective of the single electron in terms of the 'quantal Newtonian' laws. Hohenberg-Kohn density functional theory (DFT), new understandings of the theory and its extension to the presence of an external uniform magnetostatic field are described. The physical interpretation via QDFT, in terms of electron correlations, of Kohn-Sham DFT, approximations to it and Slater theory are provided.

Quantal density functional theory. 2. ed.

Energy Technology Data Exchange (ETDEWEB)

Sahni, Viraht

2016-07-01

This book is on quantal density functional theory (QDFT) which is a time-dependent local effective potential theory of the electronic structure of matter. The time-independent QDFT constitutes a special case. The 2{sup nd} edition describes the further development of the theory, and extends it to include the presence of an external magnetostatic field. The theory is based on the 'quantal Newtonian' second and first laws for the individual electron. These laws are in terms of 'classical' fields that pervade all space, and their quantal sources. The fields are separately representative of the electron correlations that must be accounted for in local potential theory. Recent developments show that irrespective of the type of external field the electrons are subject to, the only correlations beyond those due to the Pauli exclusion principle and Coulomb repulsion that need be considered are solely of the correlation-kinetic effects. Foundational to QDFT, the book describes Schroedinger theory from the new perspective of the single electron in terms of the 'quantal Newtonian' laws. Hohenberg-Kohn density functional theory (DFT), new understandings of the theory and its extension to the presence of an external uniform magnetostatic field are described. The physical interpretation via QDFT, in terms of electron correlations, of Kohn-Sham DFT, approximations to it and Slater theory are provided.

Four-Component Relativistic Density-Functional Theory Calculations of Nuclear Spin-Rotation Constants: Relativistic Effects in p-Block Hydrides.

Science.gov (United States)

Komorovsky, Stanislav; Repisky, Michal; Malkin, Elena; Demissie, Taye B; Ruud, Kenneth

2015-08-11

We present an implementation of the nuclear spin-rotation (SR) constants based on the relativistic four-component Dirac-Coulomb Hamiltonian. This formalism has been implemented in the framework of the Hartree-Fock and Kohn-Sham theory, allowing assessment of both pure and hybrid exchange-correlation functionals. In the density-functional theory (DFT) implementation of the response equations, a noncollinear generalized gradient approximation (GGA) has been used. The present approach enforces a restricted kinetic balance condition for the small-component basis at the integral level, leading to very efficient calculations of the property. We apply the methodology to study relativistic effects on the spin-rotation constants by performing calculations on XHn (n = 1-4) for all elements X in the p-block of the periodic table and comparing the effects of relativity on the nuclear SR tensors to that observed for the nuclear magnetic shielding tensors. Correlation effects as described by the density-functional theory are shown to be significant for the spin-rotation constants, whereas the differences between the use of GGA and hybrid density functionals are much smaller. Our calculated relativistic spin-rotation constants at the DFT level of theory are only in fair agreement with available experimental data. It is shown that the scaling of the relativistic effects for the spin-rotation constants (varying between Z(3.8) and Z(4.5)) is as strong as for the chemical shieldings but with a much smaller prefactor.

Teaching Density Functional Theory Through Experiential Learning

International Nuclear Information System (INIS)

Narasimhan, Shobhana

2015-01-01

Today, quantum mechanical density functional theory is often the method of choice for performing accurate calculations on atomic, molecular and condensed matter systems. Here, I share some of my experiences in teaching the necessary basics of solid state physics, as well as the theory and practice of density functional theory, in a number of workshops held in developing countries over the past two decades. I discuss the advantages of supplementing the usual mathematically formal teaching methods, characteristic of graduate courses, with the use of visual imagery and analogies. I also describe a successful experiment we carried out, which resulted in a joint publication co-authored by 67 lecturers and students participating in a summer school. (paper)

Density functional theory for polymeric systems in 2D

International Nuclear Information System (INIS)

Słyk, Edyta; Bryk, Paweł; Roth, Roland

2016-01-01

We propose density functional theory for polymeric fluids in two dimensions. The approach is based on Wertheim’s first order thermodynamic perturbation theory (TPT) and closely follows density functional theory for polymers proposed by Yu and Wu (2002 J. Chem. Phys . 117 2368). As a simple application we evaluate the density profiles of tangent hard-disk polymers at hard walls. The theoretical predictions are compared against the results of the Monte Carlo simulations. We find that for short chain lengths the theoretical density profiles are in an excellent agreement with the Monte Carlo data. The agreement is less satisfactory for longer chains. The performance of the theory can be improved by recasting the approach using the self-consistent field theory formalism. When the self-avoiding chain statistics is used, the theory yields a marked improvement in the low density limit. Further improvements for long chains could be reached by going beyond the first order of TPT. (paper)

Density functional theory calculations of the water interactions with ZrO2 nanoparticles Y2O3 doped

Science.gov (United States)

Subhoni, Mekhrdod; Kholmurodov, Kholmirzo; Doroshkevich, Aleksandr; Asgerov, Elmar; Yamamoto, Tomoyuki; Lyubchyk, Andrei; Almasan, Valer; Madadzada, Afag

2018-03-01

Development of a new electricity generation techniques is one of the most relevant tasks, especially nowadays under conditions of extreme growth in energy consumption. The exothermic heterogeneous electrochemical energy conversion to the electric energy through interaction of the ZrO2 based nanopowder system with atmospheric moisture is one of the ways of electric energy obtaining. The questions of conversion into the electric form of the energy of water molecules adsorption in 3 mol% Y2O3 doped ZrO2 nanopowder systems were investigated using the density functional theory calculations. The density functional theory calculations has been realized as in the Kohn-Sham formulation, where the exchange-correlation potential is approximated by a functional of the electronic density. The electronic density, total energy and band structure calculations are carried out using the all-electron, full potential, linear augmented plane wave method of the electronic density and related approximations, i.e. the local density, the generalized gradient and their hybrid approximations.

Open-system Kohn-Sham density functional theory.

Science.gov (United States)

Zhou, Yongxi; Ernzerhof, Matthias

2012-03-07

A simple model for electron transport through molecules is provided by the source-sink potential (SSP) method [F. Goyer, M. Ernzerhof, and M. Zhuang, J. Chem. Phys. 126, 144104 (2007)]. In SSP, the boundary conditions of having an incoming and outgoing electron current are enforced through complex potentials that are added to the Hamiltonian. Depending on the sign of the imaginary part of the potentials, current density is generated or absorbed. In this way, a finite system can be used to model infinite molecular electronic devices. The SSP has originally been developed for the Hückel method and subsequently it has been extended [F. Goyer and M. Ernzerhof, J. Chem. Phys. 134, 174101 (2011)] to the Hubbard model. Here we present a step towards its generalization for first-principles electronic structure theory methods. In particular, drawing on our earlier work, we discuss a new generalized density functional theory for complex non-Hermitian Hamiltonians. This theory enables us to combine SSP and Kohn-Sham theory to obtain a method for the description of open systems that exchange current density with their environment. Similarly, the Hartree-Fock method is extended to the realm of non-Hermitian, SSP containing Hamiltonians. As a proof of principle, we present the first applications of complex-density functional theory (CODFT) as well as non-Hermitian Hartree-Fock theory to electron transport through molecules. © 2012 American Institute of Physics

Communication: Orbital instabilities and triplet states from time-dependent density functional theory and long-range corrected functionals

Science.gov (United States)

Sears, John S.; Koerzdoerfer, Thomas; Zhang, Cai-Rong; Brédas, Jean-Luc

2011-10-01

Long-range corrected hybrids represent an increasingly popular class of functionals for density functional theory (DFT) that have proven to be very successful for a wide range of chemical applications. In this Communication, we examine the performance of these functionals for time-dependent (TD)DFT descriptions of triplet excited states. Our results reveal that the triplet energies are particularly sensitive to the range-separation parameter; this sensitivity can be traced back to triplet instabilities in the ground state coming from the large effective amounts of Hartree-Fock exchange included in these functionals. As such, the use of standard long-range corrected functionals for the description of triplet states at the TDDFT level is not recommended.

Time-dependent quantum fluid density functional theory of hydrogen ...

Indian Academy of Sciences (India)

WINTEC

density functional theory; quantum fluid dynamics. 1. Introduction ... dynamics of strongly non-linear interaction of atoms with intense ... theory and quantum fluid dynamics in real space. .... clear evidence of bond softening since density in the.

Local hybrid functionals with orbital-free mixing functions and balanced elimination of self-interaction error

International Nuclear Information System (INIS)

Silva, Piotr de; Corminboeuf, Clémence

2015-01-01

The recently introduced density overlap regions indicator (DORI) [P. de Silva and C. Corminboeuf, J. Chem. Theory Comput. 10(9), 3745–3756 (2014)] is a density-dependent scalar field revealing regions of high density overlap between shells, atoms, and molecules. In this work, we exploit its properties to construct local hybrid exchange-correlation functionals aiming at balanced reduction of the self-interaction error. We show that DORI can successfully replace the ratio of the von Weizsäcker and exact positive-definite kinetic energy densities, which is commonly used in mixing functions of local hybrids. Additionally, we introduce several semi-empirical parameters to control the local and global admixture of exact exchange. The most promising of our local hybrids clearly outperforms the underlying semi-local functionals as well as their global hybrids

Magnetic Exchange Couplings from Semilocal Functionals Evaluated Nonself-Consistently on Hybrid Densities: Insights on Relative Importance of Exchange, Correlation, and Delocalization.

Science.gov (United States)

Phillips, Jordan J; Peralta, Juan E

2012-09-11

Semilocal functionals generally yield poor magnetic exchange couplings for transition-metal complexes, typically overpredicting in magnitude the experimental values. Here we show that semilocal functionals evaluated nonself-consistently on densities from hybrid functionals can yield magnetic exchange couplings that are greatly improved with respect to their self-consistent semilocal values. Furthermore, when semilocal functionals are evaluated nonself-consistently on densities from a "half-and-half" hybrid, their errors with respect to experimental values can actually be lower than those from self-consistent calculations with standard hybrid functionals such as PBEh or TPSSh. This illustrates that despite their notoriously poor performance for exchange couplings, for many systems semilocal functionals are capable of delivering accurate relative energies for magnetic states provided that their electron delocalization error is corrected. However, while self-consistent calculations with hybrids uniformly improve results for all complexes, evaluating nonself-consistently with semilocal functionals does not give a balanced improvement for both ferro- and antiferromagnetically coupled complexes, indicating that there is more at play with the overestimation problem than simply the delocalization error. Additionally, we show that for some systems the conventional wisdom of choice of exchange functional mattering more than correlation does not hold. This combined with results from the nonself-consistent calculations provide insight on clarifying the relative roles of exchange, correlation, and delocalization in calculating magnetic exchange coupling parameters in Kohn-Sham Density Functional Theory.

Functional renormalization group and Kohn-Sham scheme in density functional theory

Science.gov (United States)

Liang, Haozhao; Niu, Yifei; Hatsuda, Tetsuo

2018-04-01

Deriving accurate energy density functional is one of the central problems in condensed matter physics, nuclear physics, and quantum chemistry. We propose a novel method to deduce the energy density functional by combining the idea of the functional renormalization group and the Kohn-Sham scheme in density functional theory. The key idea is to solve the renormalization group flow for the effective action decomposed into the mean-field part and the correlation part. Also, we propose a simple practical method to quantify the uncertainty associated with the truncation of the correlation part. By taking the φ4 theory in zero dimension as a benchmark, we demonstrate that our method shows extremely fast convergence to the exact result even for the highly strong coupling regime.

Density functional theory, natural bond orbital and quantum theory of ...

Indian Academy of Sciences (India)

Density functional theory, natural bond orbital and quantum theory of atoms in molecule analyses on the hydrogen bonding interactions in tryptophan-water complexes. XIQIAN NIU, ZHENGGUO HUANG. ∗. , LINGLING MA, TINGTING SHEN and LINGFEI GUO. Tianjin Key Laboratory of Structure and Performance for ...

Density functional theory + U modeling of polarons in organohalide lead perovskites

Directory of Open Access Journals (Sweden)

Eric Welch

2016-12-01

Full Text Available We investigate the possible formation of polarons in four organic perovskites (CH3NH3PbI3, CH3NH3PbBr3, CH3NH3PbCl3, and CH3NH3PbI2Cl1 using a density functional theory (DFT calculations with local potentials and hybrid functionals. We show that DFT+U method with U = 8 eV predicts a correct band-gap and matches the forces on ions from hybrid calculations. We then use the DFT + U approach to study the effect of polarons, i.e. to search the configuration space and locate the lowest energy localized band gap state self-trapped hole (STH. STH configurations were found for three pure halides and one mixed halide system. Spin orbit coupling (SOC was also taken into account and the results may be found in the supplementary material. This study focuses on the +U method; however, SOC corrections added to the DFT+U calculations also resulted in STH states in all four systems.

Electronic properties of B and Al doped graphane: A hybrid density functional study

Science.gov (United States)

Mapasha, R. E.; Igumbor, E.; Andriambelaza, N. F.; Chetty, N.

2018-04-01

Using a hybrid density functional theory approach parametrized by Heyd, Scuseria and Ernzerhof (HSE06 hybrid functional), we study the energetics, structural and electronic properties of a graphane monolayer substitutionally doped with the B (BCH) and Al (AlCH) atoms. The BCH defect can be integrated within a graphane monolayer at a relative low formation energy, without major structural distortions and symmetry breaking. The AlCH defect relaxes outward of the monolayer and breaks the symmetry. The density of states plots indicate that BCH doped graphane monolayer is a wide band gap semiconductor, whereas the AlCH defect introduces the spin dependent mid gap states at the vicinity of the Fermi level, revealing a metallic character with the pronounced magnetic features. We further examine the response of the Al dependent spin states on the multiple charge states doping. We find that the defect formation energy, structural and electronic properties can be altered via charge state modulation. The +1 charge doping opens an energy band gap of 1.75 eV. This value corresponds to the wavelength in the visible spectrum, suggesting an ideal material for solar cell absorbers. Our study fine tunes the graphane band gap through the foreign atom doping as well as via defect charge state modulation.

Bayesian error estimation in density-functional theory

DEFF Research Database (Denmark)

Mortensen, Jens Jørgen; Kaasbjerg, Kristen; Frederiksen, Søren Lund

2005-01-01

We present a practical scheme for performing error estimates for density-functional theory calculations. The approach, which is based on ideas from Bayesian statistics, involves creating an ensemble of exchange-correlation functionals by comparing with an experimental database of binding energies...

Rational Density Functional Selection Using Game Theory.

Science.gov (United States)

McAnanama-Brereton, Suzanne; Waller, Mark P

2018-01-22

Theoretical chemistry has a paradox of choice due to the availability of a myriad of density functionals and basis sets. Traditionally, a particular density functional is chosen on the basis of the level of user expertise (i.e., subjective experiences). Herein we circumvent the user-centric selection procedure by describing a novel approach for objectively selecting a particular functional for a given application. We achieve this by employing game theory to identify optimal functional/basis set combinations. A three-player (accuracy, complexity, and similarity) game is devised, through which Nash equilibrium solutions can be obtained. This approach has the advantage that results can be systematically improved by enlarging the underlying knowledge base, and the deterministic selection procedure mathematically justifies the density functional and basis set selections.

Quantal density functional theory II. Approximation methods and applications

International Nuclear Information System (INIS)

Sahni, Viraht

2010-01-01

This book is on approximation methods and applications of Quantal Density Functional Theory (QDFT), a new local effective-potential-energy theory of electronic structure. What distinguishes the theory from traditional density functional theory is that the electron correlations due to the Pauli exclusion principle, Coulomb repulsion, and the correlation contribution to the kinetic energy -- the Correlation-Kinetic effects -- are separately and explicitly defined. As such it is possible to study each property of interest as a function of the different electron correlations. Approximations methods based on the incorporation of different electron correlations, as well as a many-body perturbation theory within the context of QDFT, are developed. The applications are to the few-electron inhomogeneous electron gas systems in atoms and molecules, as well as to the many-electron inhomogeneity at metallic surfaces. (orig.)

Density functional theory a practical introduction

CERN Document Server

Sholl, David

2009-01-01

Demonstrates how anyone in math, science, and engineering can master DFT calculations Density functional theory (DFT) is one of the most frequently used computational tools for studying and predicting the properties of isolated molecules, bulk solids, and material interfaces, including surfaces. Although the theoretical underpinnings of DFT are quite complicated, this book demonstrates that the basic concepts underlying the calculations are simple enough to be understood by anyone with a background in chemistry, physics, engineering, or mathematics. The authors show how the widespread availability of powerful DFT codes makes it possible for students and researchers to apply this important computational technique to a broad range of fundamental and applied problems. Density Functional Theory: A Practical Introduction offers a concise, easy-to-follow introduction to the key concepts and practical applications of DFT, focusing on plane-wave DFT. The authors have many years of experience introducing DFT to studen...

Chemical hardness and density functional theory

Indian Academy of Sciences (India)

Unknown

RALPH G PEARSON. Chemistry Department, University of California, Santa Barbara, CA 93106, USA. Abstract. The concept of chemical hardness is reviewed from a personal point of view. Keywords. Hardness; softness; hard & soft acids bases (HSAB); principle of maximum hardness. (PMH) density functional theory (DFT) ...

Periodic subsystem density-functional theory

International Nuclear Information System (INIS)

Genova, Alessandro; Pavanello, Michele; Ceresoli, Davide

2014-01-01

By partitioning the electron density into subsystem contributions, the Frozen Density Embedding (FDE) formulation of subsystem Density Functional Theory (DFT) has recently emerged as a powerful tool for reducing the computational scaling of Kohn–Sham DFT. To date, however, FDE has been employed to molecular systems only. Periodic systems, such as metals, semiconductors, and other crystalline solids have been outside the applicability of FDE, mostly because of the lack of a periodic FDE implementation. To fill this gap, in this work we aim at extending FDE to treat subsystems of molecular and periodic character. This goal is achieved by a dual approach. On one side, the development of a theoretical framework for periodic subsystem DFT. On the other, the realization of the method into a parallel computer code. We find that periodic FDE is capable of reproducing total electron densities and (to a lesser extent) also interaction energies of molecular systems weakly interacting with metallic surfaces. In the pilot calculations considered, we find that FDE fails in those cases where there is appreciable density overlap between the subsystems. Conversely, we find FDE to be in semiquantitative agreement with Kohn–Sham DFT when the inter-subsystem density overlap is low. We also conclude that to make FDE a suitable method for describing molecular adsorption at surfaces, kinetic energy density functionals that go beyond the GGA level must be employed

Periodic subsystem density-functional theory

Science.gov (United States)

Genova, Alessandro; Ceresoli, Davide; Pavanello, Michele

2014-11-01

By partitioning the electron density into subsystem contributions, the Frozen Density Embedding (FDE) formulation of subsystem Density Functional Theory (DFT) has recently emerged as a powerful tool for reducing the computational scaling of Kohn-Sham DFT. To date, however, FDE has been employed to molecular systems only. Periodic systems, such as metals, semiconductors, and other crystalline solids have been outside the applicability of FDE, mostly because of the lack of a periodic FDE implementation. To fill this gap, in this work we aim at extending FDE to treat subsystems of molecular and periodic character. This goal is achieved by a dual approach. On one side, the development of a theoretical framework for periodic subsystem DFT. On the other, the realization of the method into a parallel computer code. We find that periodic FDE is capable of reproducing total electron densities and (to a lesser extent) also interaction energies of molecular systems weakly interacting with metallic surfaces. In the pilot calculations considered, we find that FDE fails in those cases where there is appreciable density overlap between the subsystems. Conversely, we find FDE to be in semiquantitative agreement with Kohn-Sham DFT when the inter-subsystem density overlap is low. We also conclude that to make FDE a suitable method for describing molecular adsorption at surfaces, kinetic energy density functionals that go beyond the GGA level must be employed.

Extending density functional embedding theory for covalently bonded systems.

Science.gov (United States)

Yu, Kuang; Carter, Emily A

2017-12-19

Quantum embedding theory aims to provide an efficient solution to obtain accurate electronic energies for systems too large for full-scale, high-level quantum calculations. It adopts a hierarchical approach that divides the total system into a small embedded region and a larger environment, using different levels of theory to describe each part. Previously, we developed a density-based quantum embedding theory called density functional embedding theory (DFET), which achieved considerable success in metals and semiconductors. In this work, we extend DFET into a density-matrix-based nonlocal form, enabling DFET to study the stronger quantum couplings between covalently bonded subsystems. We name this theory density-matrix functional embedding theory (DMFET), and we demonstrate its performance in several test examples that resemble various real applications in both chemistry and biochemistry. DMFET gives excellent results in all cases tested thus far, including predicting isomerization energies, proton transfer energies, and highest occupied molecular orbital-lowest unoccupied molecular orbital gaps for local chromophores. Here, we show that DMFET systematically improves the quality of the results compared with the widely used state-of-the-art methods, such as the simple capped cluster model or the widely used ONIOM method.

JDFTx: Software for joint density-functional theory

Directory of Open Access Journals (Sweden)

Ravishankar Sundararaman

2017-01-01

Full Text Available Density-functional theory (DFT has revolutionized computational prediction of atomic-scale properties from first principles in physics, chemistry and materials science. Continuing development of new methods is necessary for accurate predictions of new classes of materials and properties, and for connecting to nano- and mesoscale properties using coarse-grained theories. JDFTx is a fully-featured open-source electronic DFT software designed specifically to facilitate rapid development of new theories, models and algorithms. Using an algebraic formulation as an abstraction layer, compact C++11 code automatically performs well on diverse hardware including GPUs (Graphics Processing Units. This code hosts the development of joint density-functional theory (JDFT that combines electronic DFT with classical DFT and continuum models of liquids for first-principles calculations of solvated and electrochemical systems. In addition, the modular nature of the code makes it easy to extend and interface with, facilitating the development of multi-scale toolkits that connect to ab initio calculations, e.g. photo-excited carrier dynamics combining electron and phonon calculations with electromagnetic simulations.

Angle-dependent strong-field molecular ionization rates with tuned range-separated time-dependent density functional theory

Energy Technology Data Exchange (ETDEWEB)

Sissay, Adonay [Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803 (United States); Abanador, Paul; Mauger, François; Gaarde, Mette; Schafer, Kenneth J. [Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana 70803 (United States); Lopata, Kenneth, E-mail: klopata@lsu.edu [Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803 (United States); Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803 (United States)

2016-09-07

Strong-field ionization and the resulting electronic dynamics are important for a range of processes such as high harmonic generation, photodamage, charge resonance enhanced ionization, and ionization-triggered charge migration. Modeling ionization dynamics in molecular systems from first-principles can be challenging due to the large spatial extent of the wavefunction which stresses the accuracy of basis sets, and the intense fields which require non-perturbative time-dependent electronic structure methods. In this paper, we develop a time-dependent density functional theory approach which uses a Gaussian-type orbital (GTO) basis set to capture strong-field ionization rates and dynamics in atoms and small molecules. This involves propagating the electronic density matrix in time with a time-dependent laser potential and a spatial non-Hermitian complex absorbing potential which is projected onto an atom-centered basis set to remove ionized charge from the simulation. For the density functional theory (DFT) functional we use a tuned range-separated functional LC-PBE*, which has the correct asymptotic 1/r form of the potential and a reduced delocalization error compared to traditional DFT functionals. Ionization rates are computed for hydrogen, molecular nitrogen, and iodoacetylene under various field frequencies, intensities, and polarizations (angle-dependent ionization), and the results are shown to quantitatively agree with time-dependent Schrödinger equation and strong-field approximation calculations. This tuned DFT with GTO method opens the door to predictive all-electron time-dependent density functional theory simulations of ionization and ionization-triggered dynamics in molecular systems using tuned range-separated hybrid functionals.

Angle-dependent strong-field molecular ionization rates with tuned range-separated time-dependent density functional theory

International Nuclear Information System (INIS)

Sissay, Adonay; Abanador, Paul; Mauger, François; Gaarde, Mette; Schafer, Kenneth J.; Lopata, Kenneth

2016-01-01

Strong-field ionization and the resulting electronic dynamics are important for a range of processes such as high harmonic generation, photodamage, charge resonance enhanced ionization, and ionization-triggered charge migration. Modeling ionization dynamics in molecular systems from first-principles can be challenging due to the large spatial extent of the wavefunction which stresses the accuracy of basis sets, and the intense fields which require non-perturbative time-dependent electronic structure methods. In this paper, we develop a time-dependent density functional theory approach which uses a Gaussian-type orbital (GTO) basis set to capture strong-field ionization rates and dynamics in atoms and small molecules. This involves propagating the electronic density matrix in time with a time-dependent laser potential and a spatial non-Hermitian complex absorbing potential which is projected onto an atom-centered basis set to remove ionized charge from the simulation. For the density functional theory (DFT) functional we use a tuned range-separated functional LC-PBE*, which has the correct asymptotic 1/r form of the potential and a reduced delocalization error compared to traditional DFT functionals. Ionization rates are computed for hydrogen, molecular nitrogen, and iodoacetylene under various field frequencies, intensities, and polarizations (angle-dependent ionization), and the results are shown to quantitatively agree with time-dependent Schrödinger equation and strong-field approximation calculations. This tuned DFT with GTO method opens the door to predictive all-electron time-dependent density functional theory simulations of ionization and ionization-triggered dynamics in molecular systems using tuned range-separated hybrid functionals.

On the universality of the long-/short-range separation in multiconfigurational density-functional theory

Science.gov (United States)

Fromager, Emmanuel; Toulouse, Julien; Jensen, Hans Jørgen Aa.

2007-02-01

In many cases, the dynamic correlation can be calculated quite accurately and at a fairly low computational cost in Kohn-Sham density-functional theory (KS-DFT), using current standard approximate functionals. However, in general, KS-DFT does not treat static correlation effects (near degeneracy) adequately which, on the other hand, can be described in wave-function theory (WFT), for example, with a multiconfigurational self-consistent field (MCSCF) model. It is therefore of high interest to develop a hybrid model which combines the best of both WFT and DFT approaches. The merge of WFT and DFT can be achieved by splitting the two-electron interaction into long-range and short-range parts. The long-range part is then treated by WFT and the short-range part by DFT. In this work the authors consider the so-called "erf" long-range interaction erf(μr12)/r12, which is based on the standard error function, and where μ is a free parameter which controls the range of the long-/short-range decomposition. In order to formulate a general method, they propose a recipe for the definition of an optimal μopt parameter, which is independent of the approximate short-range functional and the approximate wave function, and they discuss its universality. Calculations on a test set consisting of He, Be, Ne, Mg, H2, N2, and H2O yield μopt≈0.4a.u.. A similar analysis on other types of test systems such as actinide compounds is currently in progress. Using the value of 0.4a.u. for μ, encouraging results are obtained with the hybrid MCSCF-DFT method for the dissociation energies of H2, N2, and H2O, with both short-range local-density approximation and PBE-type functionals.

Density-functional theory in one dimension for contact-interacting fermions

International Nuclear Information System (INIS)

Magyar, R.J.; Burke, K.

2004-01-01

A density-functional theory is developed for fermions in one dimension, interacting via a δ function. Such systems provide a natural testing ground for questions of principle, as the local-density approximation should be highly accurate since for this interaction type the exchange contribution to the local-density approximation is intrinsically self-interaction-free. The exact-exchange contribution to the total energy is a local functional of the density. A local-density approximation for correlation is obtained using perturbation theory and Bethe ansatz results for the one-dimensional contact-interacting uniform Fermi gas. The ground-state energies are calculated for two finite systems, the analogs of helium and of Hooke's atom. The local-density approximation is shown to be excellent as expected

Benchmark density functional theory calculations for nanoscale conductance

DEFF Research Database (Denmark)

Strange, Mikkel; Bækgaard, Iben Sig Buur; Thygesen, Kristian Sommer

2008-01-01

We present a set of benchmark calculations for the Kohn-Sham elastic transmission function of five representative single-molecule junctions. The transmission functions are calculated using two different density functional theory methods, namely an ultrasoft pseudopotential plane-wave code...

Benchmarks for electronically excited states: Time-dependent density functional theory and density functional theory based multireference configuration interaction

DEFF Research Database (Denmark)

Silva-Junior, Mario R.; Schreiber, Marko; Sauer, Stephan P. A.

2008-01-01

Time-dependent density functional theory (TD-DFT) and DFT-based multireference configuration interaction (DFT/MRCI) calculations are reported for a recently proposed benchmark set of 28 medium-sized organic molecules. Vertical excitation energies, oscillator strengths, and excited-state dipole...

Applications of Density Functional Theory in Soft Condensed Matter

Science.gov (United States)

Löwen, Hartmut

Applications of classical density functional theory (DFT) to soft matter systems like colloids, liquid crystals and polymer solutions are discussed with a focus on the freezing transition and on nonequilibrium Brownian dynamics. First, after a brief reminder of equilibrium density functional theory, DFT is applied to the freezing transition of liquids into crystalline lattices. In particular, spherical particles with radially symmetric pair potentials will be treated (like hard spheres, the classical one-component plasma or Gaussian-core particles). Second, the DFT will be generalized towards Brownian dynamics in order to tackle nonequilibrium problems. After a general introduction to Brownian dynamics using the complementary Smoluchowski and Langevin pictures appropriate for the dynamics of colloidal suspensions, the dynamical density functional theory (DDFT) will be derived from the Smoluchowski equation. This will be done first for spherical particles (e.g. hard spheres or Gaussian-cores) without hydrodynamic interactions. Then we show how to incorporate hydrodynamic interactions between the colloidal particles into the DDFT framework and compare to Brownian dynamics computer simulations. Third orientational degrees of freedom (rod-like particles) will be considered as well. In the latter case, the stability of intermediate liquid crystalline phases (isotropic, nematic, smectic-A, plastic crystals etc) can be predicted. Finally, the corresponding dynamical extension of density functional theory towards orientational degrees of freedom is proposed and the collective behaviour of "active" (self-propelled) Brownian particles is briefly discussed.

Joint density-functional theory and its application to systems in solution

Science.gov (United States)

Petrosyan, Sahak A.

The physics of solvation, the interaction of water with solutes, plays a central role in chemistry and biochemistry, and it is essential for the very existence of life. Despite the central importance of water and the advent of the quantum theory early in the twentieth century, the link between the fundamental laws of physics and the observable properties of water remain poorly understood to this day. The central goal of this thesis is to develop a new formalism and framework to make the study of systems (solutes or surfaces) in contact with liquid water as practical and accurate as standard electronic structure calculations without the need for explicit averaging over large ensembles of configurations of water molecules. The thesis introduces a new form of density functional theory for the ab initio description of electronic systems in contact with a molecular liquid environment. This theory rigorously joins an electron density-functional for the electrons of a solute with a classical density-functional theory for the liquid into a single variational principle for the free energy of the combined system. Using the new form of density-functional theory for the ab initio description of electronic systems in contact with a molecular liquid environment, the thesis then presents the first detailed study of the impact of a solvent on the surface chemistry of Cr2O3, the passivating layer of stainless steel alloys. In comparison to a vacuum, we predict that the presence of water has little impact on the adsorption of chloride ions to the oxygen-terminated surface but has a dramatic effect on the binding of hydrogen to that surface. A key ingredient of a successful joint density functional theory is a good approximate functional for describing the solvent. We explore how the simplest examples of the best known class of approximate forms for the classical density functional fail when applied directly to water. The thesis then presents a computationally efficient density-functional

Ground-state densities from the Rayleigh-Ritz variation principle and from density-functional theory.

Science.gov (United States)

Kvaal, Simen; Helgaker, Trygve

2015-11-14

The relationship between the densities of ground-state wave functions (i.e., the minimizers of the Rayleigh-Ritz variation principle) and the ground-state densities in density-functional theory (i.e., the minimizers of the Hohenberg-Kohn variation principle) is studied within the framework of convex conjugation, in a generic setting covering molecular systems, solid-state systems, and more. Having introduced admissible density functionals as functionals that produce the exact ground-state energy for a given external potential by minimizing over densities in the Hohenberg-Kohn variation principle, necessary and sufficient conditions on such functionals are established to ensure that the Rayleigh-Ritz ground-state densities and the Hohenberg-Kohn ground-state densities are identical. We apply the results to molecular systems in the Born-Oppenheimer approximation. For any given potential v ∈ L(3/2)(ℝ(3)) + L(∞)(ℝ(3)), we establish a one-to-one correspondence between the mixed ground-state densities of the Rayleigh-Ritz variation principle and the mixed ground-state densities of the Hohenberg-Kohn variation principle when the Lieb density-matrix constrained-search universal density functional is taken as the admissible functional. A similar one-to-one correspondence is established between the pure ground-state densities of the Rayleigh-Ritz variation principle and the pure ground-state densities obtained using the Hohenberg-Kohn variation principle with the Levy-Lieb pure-state constrained-search functional. In other words, all physical ground-state densities (pure or mixed) are recovered with these functionals and no false densities (i.e., minimizing densities that are not physical) exist. The importance of topology (i.e., choice of Banach space of densities and potentials) is emphasized and illustrated. The relevance of these results for current-density-functional theory is examined.

Exchange-correlation energies of atoms from efficient density functionals: influence of the electron density

Science.gov (United States)

Tao, Jianmin; Ye, Lin-Hui; Duan, Yuhua

2017-12-01

The primary goal of Kohn-Sham density functional theory is to evaluate the exchange-correlation contribution to electronic properties. However, the accuracy of a density functional can be affected by the electron density. Here we apply the nonempirical Tao-Mo (TM) semilocal functional to study the influence of the electron density on the exchange and correlation energies of atoms and ions, and compare the results with the commonly used nonempirical semilocal functionals local spin-density approximation (LSDA), Perdew-Burke-Ernzerhof (PBE), Tao-Perdew-Staroverov-Scuseria (TPSS), and hybrid functional PBE0. We find that the spin-restricted Hartree-Fock density yields the exchange and correlation energies in good agreement with the Optimized Effective Potential method, particularly for spherical atoms and ions. However, the errors of these semilocal and hybrid functionals become larger for self-consistent densities. We further find that the quality of the electron density have greater effect on the exchange-correlation energies of kinetic energy density-dependent meta-GGA functionals TPSS and TM than on those of the LSDA and GGA, and therefore, should have greater influence on the performance of meta-GGA functionals. Finally, we show that the influence of the density quality on PBE0 is slightly reduced, compared to that of PBE, due to the exact mixing.

Rydberg energies using excited state density functional theory

International Nuclear Information System (INIS)

Cheng, C.-L.; Wu Qin; Van Voorhis, Troy

2008-01-01

We utilize excited state density functional theory (eDFT) to study Rydberg states in atoms. We show both analytically and numerically that semilocal functionals can give quite reasonable Rydberg energies from eDFT, even in cases where time dependent density functional theory (TDDFT) fails catastrophically. We trace these findings to the fact that in eDFT the Kohn-Sham potential for each state is computed using the appropriate excited state density. Unlike the ground state potential, which typically falls off exponentially, the sequence of excited state potentials has a component that falls off polynomially with distance, leading to a Rydberg-type series. We also address the rigorous basis of eDFT for these systems. Perdew and Levy have shown using the constrained search formalism that every stationary density corresponds, in principle, to an exact stationary state of the full many-body Hamiltonian. In the present context, this means that the excited state DFT solutions are rigorous as long as they deliver the minimum noninteracting kinetic energy for the given density. We use optimized effective potential techniques to show that, in some cases, the eDFT Rydberg solutions appear to deliver the minimum kinetic energy because the associated density is not pure state v-representable. We thus find that eDFT plays a complementary role to constrained DFT: The former works only if the excited state density is not the ground state of some potential while the latter applies only when the density is a ground state density.

Single-particle energies and density of states in density functional theory

Science.gov (United States)

van Aggelen, H.; Chan, G. K.-L.

2015-07-01

Time-dependent density functional theory (TD-DFT) is commonly used as the foundation to obtain neutral excited states and transition weights in DFT, but does not allow direct access to density of states and single-particle energies, i.e. ionisation energies and electron affinities. Here we show that by extending TD-DFT to a superfluid formulation, which involves operators that break particle-number symmetry, we can obtain the density of states and single-particle energies from the poles of an appropriate superfluid response function. The standard Kohn- Sham eigenvalues emerge as the adiabatic limit of the superfluid response under the assumption that the exchange- correlation functional has no dependence on the superfluid density. The Kohn- Sham eigenvalues can thus be interpreted as approximations to the ionisation energies and electron affinities. Beyond this approximation, the formalism provides an incentive for creating a new class of density functionals specifically targeted at accurate single-particle eigenvalues and bandgaps.

Multiconfiguration Pair-Density Functional Theory Is Free From Delocalization Error.

Science.gov (United States)

Bao, Junwei Lucas; Wang, Ying; He, Xiao; Gagliardi, Laura; Truhlar, Donald G

2017-11-16

Delocalization error has been singled out by Yang and co-workers as the dominant error in Kohn-Sham density functional theory (KS-DFT) with conventional approximate functionals. In this Letter, by computing the vertical first ionization energy for well separated He clusters, we show that multiconfiguration pair-density functional theory (MC-PDFT) is free from delocalization error. To put MC-PDFT in perspective, we also compare it with some Kohn-Sham density functionals, including both traditional and modern functionals. Whereas large delocalization errors are almost universal in KS-DFT (the only exception being the very recent corrected functionals of Yang and co-workers), delocalization error is removed by MC-PDFT, which bodes well for its future as a step forward from KS-DFT.

A density functional theory-based chemical potential equalisation

Indian Academy of Sciences (India)

A chemical potential equalisation scheme is proposed for the calculation of these quantities and hence the dipole polarizability within the framework of density functional theory based linear response theory. The resulting polarizability is expressed in terms of the contributions from individual atoms in the molecule. A few ...

Analysis of self-consistency effects in range-separated density-functional theory with Møller-Plesset perturbation theory

DEFF Research Database (Denmark)

Fromager, Emmanuel; Jensen, Hans Jørgen Aagaard

2011-01-01

Range-separated density-functional theory combines wave function theory for the long-range part of the two-electron interaction with density-functional theory for the short-range part. When describing the long-range interaction with non-variational methods, such as perturbation or coupled......-cluster theories, self-consistency effects are introduced in the density functional part, which for an exact solution requires iterations. They are generally assumed to be small but no detailed study has been performed so far. Here, the authors analyze self-consistency when using Møller-Plesset-type (MP......) perturbation theory for the long range interaction. The lowest-order self-consistency corrections to the wave function and the energy, that enter the perturbation expansions at the second and fourth order, respectively, are both expressed in terms of the one-electron reduced density matrix. The computational...

Size-dependent error of the density functional theory ionization potential in vacuum and solution.

Science.gov (United States)

Sosa Vazquez, Xochitl A; Isborn, Christine M

2015-12-28

Density functional theory is often the method of choice for modeling the energetics of large molecules and including explicit solvation effects. It is preferable to use a method that treats systems of different sizes and with different amounts of explicit solvent on equal footing. However, recent work suggests that approximate density functional theory has a size-dependent error in the computation of the ionization potential. We here investigate the lack of size-intensivity of the ionization potential computed with approximate density functionals in vacuum and solution. We show that local and semi-local approximations to exchange do not yield a constant ionization potential for an increasing number of identical isolated molecules in vacuum. Instead, as the number of molecules increases, the total energy required to ionize the system decreases. Rather surprisingly, we find that this is still the case in solution, whether using a polarizable continuum model or with explicit solvent that breaks the degeneracy of each solute, and we find that explicit solvent in the calculation can exacerbate the size-dependent delocalization error. We demonstrate that increasing the amount of exact exchange changes the character of the polarization of the solvent molecules; for small amounts of exact exchange the solvent molecules contribute a fraction of their electron density to the ionized electron, but for larger amounts of exact exchange they properly polarize in response to the cationic solute. In vacuum and explicit solvent, the ionization potential can be made size-intensive by optimally tuning a long-range corrected hybrid functional.

Functional development in density functional theory for superconductors

Energy Technology Data Exchange (ETDEWEB)

Sanna, Antonio; Gross, E.K.U.; Essenberger, Frank [Max Planck Institute of Microstructure Physics, Halle (Saale) (Germany)

2015-07-01

Density functional theory for superconductors (SCDFT) is a fully parameter-free approach to superconductivity that allows for accurate predictions of critical temperature and properties of superconductors. We report on the most recent extensions of the method, in particular the development of new functionals to: (1) incorporate in a correct fashion Migdal's theorem; (2) compute the excitation spectrum; (3) include spin-fluctuation mediated pairing Applications and predictions are shown for a set of materials, including conventional and unconventional superconductors.

A classical density functional theory of ionic liquids.

Science.gov (United States)

Forsman, Jan; Woodward, Clifford E; Trulsson, Martin

2011-04-28

We present a simple, classical density functional approach to the study of simple models of room temperature ionic liquids. Dispersion attractions as well as ion correlation effects and excluded volume packing are taken into account. The oligomeric structure, common to many ionic liquid molecules, is handled by a polymer density functional treatment. The theory is evaluated by comparisons with simulations, with an emphasis on the differential capacitance, an experimentally measurable quantity of significant practical interest.

Density functional theory calculations of charge transport properties ...

Indian Academy of Sciences (India)

ZIRAN CHEN

2017-08-04

Aug 4, 2017 ... properties of 'plate-like' coronene topological structures ... Keywords. Organic semiconductors; density functional theory; charge carrier mobility; ambipolar transport; ..... nology Department of Sichuan Province (Grant Number.

Decay of autoionizing states in time-dependent density functional and reduced density matrix functional theory

Energy Technology Data Exchange (ETDEWEB)

Kapoor, Varun; Brics, Martins; Bauer, Dieter [Institut fuer Physik, Universitaet Rostock, 18051 Rostock (Germany)

2013-07-01

Autoionizing states are inaccessible to time-dependent density functional theory (TDDFT) using known, adiabatic Kohn-Sham (KS) potentials. We determine the exact KS potential for a numerically exactly solvable model Helium atom interacting with a laser field that is populating an autoionizing state. The exact single-particle density of the population in the autoionizing state corresponds to that of the energetically lowest quasi-stationary state in the exact KS potential. We describe how this exact potential controls the decay by a barrier whose height and width allows for the density to tunnel out and decay with the same rate as in the ab initio time-dependent Schroedinger calculation. However, devising a useful exchange-correlation potential that is capable of governing such a scenario in general and in more complex systems is hopeless. As an improvement over TDDFT, time-dependent reduced density matrix functional theory has been proposed. We are able to obtain for the above described autoionization process the exact time-dependent natural orbitals (i.e., the eigenfunctions of the exact, time-dependent one-body reduced density matrix) and study the potentials that appear in the equations of motion for the natural orbitals and the structure of the two-body density matrix expanded in them.

Approximate spin projected spin-unrestricted density functional theory method: Application to diradical character dependences of second hyperpolarizabilities

Energy Technology Data Exchange (ETDEWEB)

Nakano, Masayoshi, E-mail: mnaka@cheng.es.osaka-u.ac.jp; Minami, Takuya, E-mail: mnaka@cheng.es.osaka-u.ac.jp; Fukui, Hitoshi, E-mail: mnaka@cheng.es.osaka-u.ac.jp; Yoneda, Kyohei, E-mail: mnaka@cheng.es.osaka-u.ac.jp; Shigeta, Yasuteru, E-mail: mnaka@cheng.es.osaka-u.ac.jp; Kishi, Ryohei, E-mail: mnaka@cheng.es.osaka-u.ac.jp [Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531 (Japan); Champagne, Benoît; Botek, Edith [Laboratoire de Chimie Théorique, Facultés Universitaires Notre-Dame de la Paix (FUNDP), rue de Bruxelles, 61, 5000 Namur (Belgium)

2015-01-22

We develop a novel method for the calculation and the analysis of the one-electron reduced densities in open-shell molecular systems using the natural orbitals and approximate spin projected occupation numbers obtained from broken symmetry (BS), i.e., spin-unrestricted (U), density functional theory (DFT) calculations. The performance of this approximate spin projection (ASP) scheme is examined for the diradical character dependence of the second hyperpolarizability (γ) using several exchange-correlation functionals, i.e., hybrid and long-range corrected UDFT schemes. It is found that the ASP-LC-UBLYP method with a range separating parameter μ = 0.47 reproduces semi-quantitatively the strongly-correlated [UCCSD(T)] result for p-quinodimethane, i.e., the γ variation as a function of the diradical character.

Source-Free Exchange-Correlation Magnetic Fields in Density Functional Theory.

Science.gov (United States)

Sharma, S; Gross, E K U; Sanna, A; Dewhurst, J K

2018-03-13

Spin-dependent exchange-correlation energy functionals in use today depend on the charge density and the magnetization density: E xc [ρ, m]. However, it is also correct to define the functional in terms of the curl of m for physical external fields: E xc [ρ,∇ × m]. The exchange-correlation magnetic field, B xc , then becomes source-free. We study this variation of the theory by uniquely removing the source term from local and generalized gradient approximations to the functional. By doing so, the total Kohn-Sham moments are improved for a wide range of materials for both functionals. Significantly, the moments for the pnictides are now in good agreement with experiment. This source-free method is simple to implement in all existing density functional theory codes.

Density functional theory

International Nuclear Information System (INIS)

Freyss, M.

2015-01-01

This chapter gives an introduction to first-principles electronic structure calculations based on the density functional theory (DFT). Electronic structure calculations have a crucial importance in the multi-scale modelling scheme of materials: not only do they enable one to accurately determine physical and chemical properties of materials, they also provide data for the adjustment of parameters (or potentials) in higher-scale methods such as classical molecular dynamics, kinetic Monte Carlo, cluster dynamics, etc. Most of the properties of a solid depend on the behaviour of its electrons, and in order to model or predict them it is necessary to have an accurate method to compute the electronic structure. DFT is based on quantum theory and does not make use of any adjustable or empirical parameter: the only input data are the atomic number of the constituent atoms and some initial structural information. The complicated many-body problem of interacting electrons is replaced by an equivalent single electron problem, in which each electron is moving in an effective potential. DFT has been successfully applied to the determination of structural or dynamical properties (lattice structure, charge density, magnetisation, phonon spectra, etc.) of a wide variety of solids. Its efficiency was acknowledged by the attribution of the Nobel Prize in Chemistry in 1998 to one of its authors, Walter Kohn. A particular attention is given in this chapter to the ability of DFT to model the physical properties of nuclear materials such as actinide compounds. The specificities of the 5f electrons of actinides will be presented, i.e., their more or less high degree of localisation around the nuclei and correlations. The limitations of the DFT to treat the strong 5f correlations are one of the main issues for the DFT modelling of nuclear fuels. Various methods that exist to better treat strongly correlated materials will finally be presented. (author)

Plato: A localised orbital based density functional theory code

Science.gov (United States)

Kenny, S. D.; Horsfield, A. P.

2009-12-01

The Plato package allows both orthogonal and non-orthogonal tight-binding as well as density functional theory (DFT) calculations to be performed within a single framework. The package also provides extensive tools for analysing the results of simulations as well as a number of tools for creating input files. The code is based upon the ideas first discussed in Sankey and Niklewski (1989) [1] with extensions to allow high-quality DFT calculations to be performed. DFT calculations can utilise either the local density approximation or the generalised gradient approximation. Basis sets from minimal basis through to ones containing multiple radial functions per angular momenta and polarisation functions can be used. Illustrations of how the package has been employed are given along with instructions for its utilisation. Program summaryProgram title: Plato Catalogue identifier: AEFC_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFC_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 219 974 No. of bytes in distributed program, including test data, etc.: 1 821 493 Distribution format: tar.gz Programming language: C/MPI and PERL Computer: Apple Macintosh, PC, Unix machines Operating system: Unix, Linux and Mac OS X Has the code been vectorised or parallelised?: Yes, up to 256 processors tested RAM: Up to 2 Gbytes per processor Classification: 7.3 External routines: LAPACK, BLAS and optionally ScaLAPACK, BLACS, PBLAS, FFTW Nature of problem: Density functional theory study of electronic structure and total energies of molecules, crystals and surfaces. Solution method: Localised orbital based density functional theory. Restrictions: Tight-binding and density functional theory only, no exact exchange. Unusual features: Both atom centred and uniform meshes available

What Density Functional Theory could do for Quantum Information

Science.gov (United States)

Mattsson, Ann

2015-03-01

The Hohenberg-Kohn theorem of Density Functional Theory (DFT), and extensions thereof, tells us that all properties of a system of electrons can be determined through their density, which uniquely determines the many-body wave-function. Given access to the appropriate, universal, functionals of the density we would, in theory, be able to determine all observables of any electronic system, without explicit reference to the wave-function. On the other hand, the wave-function is at the core of Quantum Information (QI), with the wave-function of a set of qubits being the central computational resource in a quantum computer. While there is seemingly little overlap between DFT and QI, reliance upon observables form a key connection. Though the time-evolution of the wave-function and associated phase information is fundamental to quantum computation, the initial and final states of a quantum computer are characterized by observables of the system. While observables can be extracted directly from a system's wave-function, DFT tells us that we may be able to intuit a method for extracting them from its density. In this talk, I will review the fundamentals of DFT and how these principles connect to the world of QI. This will range from DFT's utility in the engineering of physical qubits, to the possibility of using it to efficiently (but approximately) simulate Hamiltonians at the logical level. The apparent paradox of describing algorithms based on the quantum mechanical many-body wave-function with a DFT-like theory based on observables will remain a focus throughout. The ultimate goal of this talk is to initiate a dialog about what DFT could do for QI, in theory and in practice. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Tailoring the electronic structure of β-Ga2O3 by non-metal doping from hybrid density functional theory calculations.

Science.gov (United States)

Guo, Weiyan; Guo, Yating; Dong, Hao; Zhou, Xin

2015-02-28

A systematic study using density functional theory has been performed for β-Ga2O3 doped with non-metal elements X (X = C, N, F, Si, P, S, Cl, Se, Br, and I) to evaluate the effect of doping on the band edges and photocatalytic activity of β-Ga2O3. The utilization of a more reliable hybrid density functional, as prescribed by Heyd, Scuseria and Ernzerhof, is found to be effective in predicting the band gap of β-Ga2O3 (4.5 eV), in agreement with the experimental result (4.59 eV). Based on the relaxed structures of X-doped systems, the defect formation energies and the plots of density of states have been calculated to analyze the band edges, the band gap states and the preferred doping sites. Our results show that the doping is energetically favored under Ga-rich growth conditions with respect to O-rich growth conditions. It is easier to replace the threefold coordinated O atom with non-metal elements compared to the fourfold coordinated O atom. X-doped systems (X = C, Si, P) show no change in the band gap, with the presence of discrete midgap states, which have adverse effect on the photocatalytic properties. The photocatalytic redox ability can be improved to a certain extent by doping with N, S, Cl, Se, Br, and I. The band alignments for Se-doped and I-doped β-Ga2O3 are well positioned for the feasibility of both photo-oxidation and photo-reduction of water, which are promising photocatalysts for water splitting in the visible region.

Perspective: Fundamental aspects of time-dependent density functional theory

Energy Technology Data Exchange (ETDEWEB)

Maitra, Neepa T. [Department of Physics and Astronomy, Hunter College and the Physics Program at the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065 (United States)

2016-06-14

In the thirty-two years since the birth of the foundational theorems, time-dependent density functional theory has had a tremendous impact on calculations of electronic spectra and dynamics in chemistry, biology, solid-state physics, and materials science. Alongside the wide-ranging applications, there has been much progress in understanding fundamental aspects of the functionals and the theory itself. This Perspective looks back to some of these developments, reports on some recent progress and current challenges for functionals, and speculates on future directions to improve the accuracy of approximations used in this relatively young theory.

Difficulties in applying pure Kohn-Sham density functional theory electronic structure methods to protein molecules

Science.gov (United States)

Rudberg, Elias

2012-02-01

Self-consistency-based Kohn-Sham density functional theory (KS-DFT) electronic structure calculations with Gaussian basis sets are reported for a set of 17 protein-like molecules with geometries obtained from the Protein Data Bank. It is found that in many cases such calculations do not converge due to vanishing HOMO-LUMO gaps. A sequence of polyproline I helix molecules is also studied and it is found that self-consistency calculations using pure functionals fail to converge for helices longer than six proline units. Since the computed gap is strongly correlated to the fraction of Hartree-Fock exchange, test calculations using both pure and hybrid density functionals are reported. The tested methods include the pure functionals BLYP, PBE and LDA, as well as Hartree-Fock and the hybrid functionals BHandHLYP, B3LYP and PBE0. The effect of including solvent molecules in the calculations is studied, and it is found that the inclusion of explicit solvent molecules around the protein fragment in many cases gives a larger gap, but that convergence problems due to vanishing gaps still occur in calculations with pure functionals. In order to achieve converged results, some modeling of the charge distribution of solvent water molecules outside the electronic structure calculation is needed. Representing solvent water molecules by a simple point charge distribution is found to give non-vanishing HOMO-LUMO gaps for the tested protein-like systems also for pure functionals.

Difficulties in applying pure Kohn-Sham density functional theory electronic structure methods to protein molecules

International Nuclear Information System (INIS)

Rudberg, Elias

2012-01-01

Self-consistency-based Kohn-Sham density functional theory (KS-DFT) electronic structure calculations with Gaussian basis sets are reported for a set of 17 protein-like molecules with geometries obtained from the Protein Data Bank. It is found that in many cases such calculations do not converge due to vanishing HOMO-LUMO gaps. A sequence of polyproline I helix molecules is also studied and it is found that self-consistency calculations using pure functionals fail to converge for helices longer than six proline units. Since the computed gap is strongly correlated to the fraction of Hartree-Fock exchange, test calculations using both pure and hybrid density functionals are reported. The tested methods include the pure functionals BLYP, PBE and LDA, as well as Hartree-Fock and the hybrid functionals BHandHLYP, B3LYP and PBE0. The effect of including solvent molecules in the calculations is studied, and it is found that the inclusion of explicit solvent molecules around the protein fragment in many cases gives a larger gap, but that convergence problems due to vanishing gaps still occur in calculations with pure functionals. In order to achieve converged results, some modeling of the charge distribution of solvent water molecules outside the electronic structure calculation is needed. Representing solvent water molecules by a simple point charge distribution is found to give non-vanishing HOMO-LUMO gaps for the tested protein-like systems also for pure functionals. (fast track communication)

Fundamentals of time-dependent density functional theory

International Nuclear Information System (INIS)

Marques, Miguel A.L.; Rubio, Angel

2012-01-01

There have been many significant advances in time-dependent density functional theory over recent years, both in enlightening the fundamental theoretical basis of the theory, as well as in computational algorithms and applications. This book, as successor to the highly successful volume Time-Dependent Density Functional Theory (Lect. Notes Phys. 706, 2006) brings together for the first time all recent developments in a systematic and coherent way. First, a thorough pedagogical presentation of the fundamental theory is given, clarifying aspects of the original proofs and theorems, as well as presenting fresh developments that extend the theory into new realms such as alternative proofs of the original Runge-Gross theorem, open quantum systems, and dispersion forces to name but a few. Next, all of the basic concepts are introduced sequentially and building in complexity, eventually reaching the level of open problems of interest. Contemporary applications of the theory are discussed, from real-time coupled-electron-ion dynamics, to excited-state dynamics and molecular transport. Last but not least, the authors introduce and review recent advances in computational implementation, including massively parallel architectures and graphical processing units. Special care has been taken in editing this volume as a multi-author textbook, following a coherent line of thought, and making all the relevant connections between chapters and concepts consistent throughout. As such it will prove to be the text of reference in this field, both for beginners as well as expert researchers and lecturers teaching advanced quantum mechanical methods to model complex physical systems, from molecules to nanostructures, from biocomplexes to surfaces, solids and liquids. (orig.)

Screened exchange hybrid density functional for accurate and efficient structures and interaction energies.

Science.gov (United States)

Brandenburg, Jan Gerit; Caldeweyher, Eike; Grimme, Stefan

2016-06-21

We extend the recently introduced PBEh-3c global hybrid density functional [S. Grimme et al., J. Chem. Phys., 2015, 143, 054107] by a screened Fock exchange variant based on the Henderson-Janesko-Scuseria exchange hole model. While the excellent performance of the global hybrid is maintained for small covalently bound molecules, its performance for computed condensed phase mass densities is further improved. Most importantly, a speed up of 30 to 50% can be achieved and especially for small orbital energy gap cases, the method is numerically much more robust. The latter point is important for many applications, e.g., for metal-organic frameworks, organic semiconductors, or protein structures. This enables an accurate density functional based electronic structure calculation of a full DNA helix structure on a single core desktop computer which is presented as an example in addition to comprehensive benchmark results.

Density-functional theory of atoms and molecules

CERN Document Server

Parr, Robert G

1995-01-01

Provides an account of the fundamental principles of the density-functional theory of the electronic structure of matter and its applications to atoms and molecules. This book contains a discussion of the chemical potential and its derivatives. It is intended for physicists, chemists, and advanced students in chemistry.

Assessment of oscillator strengths with multiconfigurational short-range density functional theory for electronic excitations in organic molecules

DEFF Research Database (Denmark)

Hedegård, Erik Donovan

2017-01-01

considered the large collection of organic molecules whose excited states were investigated with a range of electronic structure methods by Thiel et al. As a by-product of our calculations of oscillator strengths, we also obtain electronic excitation energies, which enable us to compare the performance......We have in a series of recent papers investigated electronic excited states with a hybrid between a complete active space self-consistent field (CASSCF) wave function and density functional theory (DFT). This method has been dubbed the CAS short-range DFT method (CAS–srDFT). The previous papers...

Density Functional Theory Study of Leaching Performance of Different Acids on Pyrochlore (100) Surface

Science.gov (United States)

Yang, Xiuli; Fang, Qing; Ouyang, Hui

2018-06-01

Pyrochlore leaching using hydrofluoric, sulfuric, and hydrochloric acids has been studied via experimental methods for years, but the interactions between niobium atoms on the pyrochlore surface and different acids have not been investigated. In this work, first-principles calculations based on density functional theory were used to elucidate the leaching performance of these three acids from the viewpoint of geometrical and electronic structures. The calculation results indicate that sulfate, chloride, and fluoride anions influence the geometric structure of pyrochlore (100) to different extents, decreasing in the order: sulfate, fluoride, chloride. Orbitals of O1 and O2 atoms of sulfate hybridized with those of surface niobium atom. Fluorine orbitals hybridized with those of surface niobium atoms. However, no obvious overlap exists between any orbitals of chlorine and surface niobium, revealing that chlorine does not interact chemically with surface niobium atoms.

Prediction of the equilibrium structures and photomagnetic properties of the Prussian blue analogue RbMn[Fe(CN)(6)] by density functional theory

NARCIS (Netherlands)

Luzon, Javier; Castro, Miguel; Vertelman, Esther J.M.; Gengler, Régis Y.N.; van Koningsbruggen, Petra J.; Molodtsova, Olga; Knupfer, Martin; Rudolf, Petra; Loosdrecht, Paul H.M. van; Broer, Ria

2008-01-01

A periodic density functional theory method using the B3LYP hybrid exchange-correlation potential is applied to the Prussian blue analogue RbMn[Fe(CN)(6)] to evaluate the suitability of the method for studying, and predicting, the photomagnetic behavior of Prussian blue analogues and related

Multiconfiguration pair-density functional theory: barrier heights and main group and transition metal energetics.

Science.gov (United States)

Carlson, Rebecca K; Li Manni, Giovanni; Sonnenberger, Andrew L; Truhlar, Donald G; Gagliardi, Laura

2015-01-13

Kohn-Sham density functional theory, resting on the representation of the electronic density and kinetic energy by a single Slater determinant, has revolutionized chemistry, but for open-shell systems, the Kohn-Sham Slater determinant has the wrong symmetry properties as compared to an accurate wave function. We have recently proposed a theory, called multiconfiguration pair-density functional theory (MC-PDFT), in which the electronic kinetic energy and classical Coulomb energy are calculated from a multiconfiguration wave function with the correct symmetry properties, and the rest of the energy is calculated from a density functional, called the on-top density functional, that depends on the density and the on-top pair density calculated from this wave function. We also proposed a simple way to approximate the on-top density functional by translation of Kohn-Sham exchange-correlation functionals. The method is much less expensive than other post-SCF methods for calculating the dynamical correlation energy starting with a multiconfiguration self-consistent-field wave function as the reference wave function, and initial tests of the theory were quite encouraging. Here, we provide a broader test of the theory by applying it to bond energies of main-group molecules and transition metal complexes, barrier heights and reaction energies for diverse chemical reactions, proton affinities, and the water dimerization energy. Averaged over 56 data points, the mean unsigned error is 3.2 kcal/mol for MC-PDFT, as compared to 6.9 kcal/mol for Kohn-Sham theory with a comparable density functional. MC-PDFT is more accurate on average than complete active space second-order perturbation theory (CASPT2) for main-group small-molecule bond energies, alkyl bond dissociation energies, transition-metal-ligand bond energies, proton affinities, and the water dimerization energy.

On the accuracy of density-functional theory exchange-correlation functionals for H bonds in small water clusters: Benchmarks approaching the complete basis set limit

Science.gov (United States)

Santra, Biswajit; Michaelides, Angelos; Scheffler, Matthias

2007-11-01

The ability of several density-functional theory (DFT) exchange-correlation functionals to describe hydrogen bonds in small water clusters (dimer to pentamer) in their global minimum energy structures is evaluated with reference to second order Møller-Plesset perturbation theory (MP2). Errors from basis set incompleteness have been minimized in both the MP2 reference data and the DFT calculations, thus enabling a consistent systematic evaluation of the true performance of the tested functionals. Among all the functionals considered, the hybrid X3LYP and PBE0 functionals offer the best performance and among the nonhybrid generalized gradient approximation functionals, mPWLYP and PBE1W perform best. The popular BLYP and B3LYP functionals consistently underbind and PBE and PW91 display rather variable performance with cluster size.

Subsystem density functional theory with meta-generalized gradient approximation exchange-correlation functionals.

Science.gov (United States)

Śmiga, Szymon; Fabiano, Eduardo; Laricchia, Savio; Constantin, Lucian A; Della Sala, Fabio

2015-04-21

We analyze the methodology and the performance of subsystem density functional theory (DFT) with meta-generalized gradient approximation (meta-GGA) exchange-correlation functionals for non-bonded molecular systems. Meta-GGA functionals depend on the Kohn-Sham kinetic energy density (KED), which is not known as an explicit functional of the density. Therefore, they cannot be directly applied in subsystem DFT calculations. We propose a Laplacian-level approximation to the KED which overcomes this limitation and provides a simple and accurate way to apply meta-GGA exchange-correlation functionals in subsystem DFT calculations. The so obtained density and energy errors, with respect to the corresponding supermolecular calculations, are comparable with conventional approaches, depending almost exclusively on the approximations in the non-additive kinetic embedding term. An embedding energy error decomposition explains the accuracy of our method.

Density functional theory study of bulk and single-layer magnetic semiconductor CrPS4

Science.gov (United States)

Zhuang, Houlong L.; Zhou, Jia

2016-11-01

Searching for two-dimensional (2D) materials with multifunctionality is one of the main goals of current research in 2D materials. Magnetism and semiconducting are certainly two desirable functional properties for a single 2D material. In line with this goal, here we report a density functional theory (DFT) study of bulk and single-layer magnetic semiconductor CrPS4. We find that the ground-state magnetic structure of bulk CrPS4 exhibits the A-type antiferromagnetic ordering, which transforms to ferromagnetic (FM) ordering in single-layer CrPS4. The calculated formation energy and phonon spectrum confirm the stability of single-layer CrPS4. The band gaps of FM single-layer CrPS4 calculated with a hybrid density functional are within the visible-light range. We also study the effects of FM ordering on the optical absorption spectra and band alignments for water splitting, indicating that single-layer CrPS4 could be a potential photocatalyst. Our work opens up ample opportunities of energy-related applications of single-layer CrPS4.

Relativistic Adiabatic Time-Dependent Density Functional Theory Using Hybrid Functionals and Noncollinear Spin Magnetization

DEFF Research Database (Denmark)

Bast, Radovan; Jensen, Hans Jørgen Aagaard; Saue, Trond

2009-01-01

into reduction of algebra from quaternion to complex or real. For hybrid GGAs with noncollinear spin magnetization we derive a new computationally advantageous equation for the full second variational derivatives of such exchange-correlation functionals. We apply our implementation to calculations on the ns2...... → ns1np1 excitation energies in the Zn, Cd, and Hg atoms (n = 4-6) and (vertical) excitation energies of UO2+ 2 ; and we test the performance of various functionals by comparison with experimental data (group 12 atoms) or higher-level computational results (UO2+2 ). The results indicate...

Long-Range Corrected Hybrid Density Functionals with Damped Atom-Atom Dispersion Corrections

Energy Technology Data Exchange (ETDEWEB)

Chai, Jeng-Da; Head-Gordon, Martin

2008-06-14

We report re-optimization of a recently proposed long-range corrected (LC) hybrid density functionals [J.-D. Chai and M. Head-Gordon, J. Chem. Phys. 128, 084106 (2008)] to include empirical atom-atom dispersion corrections. The resulting functional, {omega}B97X-D yields satisfactory accuracy for thermochemistry, kinetics, and non-covalent interactions. Tests show that for non-covalent systems, {omega}B97X-D shows slight improvement over other empirical dispersion-corrected density functionals, while for covalent systems and kinetics, it performs noticeably better. Relative to our previous functionals, such as {omega}B97X, the new functional is significantly superior for non-bonded interactions, and very similar in performance for bonded interactions.

Density functional theory of polydisperse fluid interfaces

International Nuclear Information System (INIS)

Baus, M.; Bellier-Castella, L.; Xu, H.

2002-01-01

Most colloids usually exhibit one or several polydispersities. A natural framework for the theoretical description of polydisperse systems is provided by the extension of density functional theory to 'continuous' mixtures. This will be illustrated here by the study of both the bulk and interfacial properties of a simple van der Waals model for a polydisperse colloidal fluid. (author)

Optical Absorption in Molecular Crystals from Time-Dependent Density Functional Theory

Science.gov (United States)

2017-04-23

Our approach represents a full solid-state calculation, allowing for polarization ef- fects while still capable of capturing inter-molecular dis...AFRL-AFOSR-UK-TR-2017-0030 Optical absorption in molecular crystals from time-dependent density functional theory Leeor Kronik WEIZMANN INSTITUTE OF...from time-dependent density functional theory 5a.  CONTRACT NUMBER 5b.  GRANT NUMBER FA9550-15-1-0290 5c.  PROGRAM ELEMENT NUMBER 61102F 6. AUTHOR(S

Orbital-dependent exchange-correlation functionals in density-functional theory realized by the FLAPW method

Energy Technology Data Exchange (ETDEWEB)

Betzinger, Markus

2011-12-14

In this thesis, we extended the applicability of the full-potential linearized augmented-plane-wave (FLAPW) method, one of the most precise, versatile and generally applicable electronic structure methods for solids working within the framework of density-functional theory (DFT), to orbital-dependent functionals for the exchange-correlation (xc) energy. Two different schemes that deal with orbital-dependent functionals, the Kohn-Sham (KS) and the generalized Kohn-Sham (gKS) formalism, have been realized. Hybrid functionals, combining some amount of the orbital-dependent exact exchange energy with local or semi-local functionals of the density, are implemented within the gKS scheme. We work in particular with the PBE0 hybrid of Perdew, Burke, and Ernzerhof. Our implementation relies on a representation of the non-local exact exchange potential - its calculation constitutes the most time consuming step in a practical calculation - by an auxiliary mixed product basis (MPB). In this way, the matrix elements of the Hamiltonian corresponding to the non-local potential become a Brillouin-zone (BZ) sum over vector-matrix-vector products. Several techniques are developed and explored to further accelerate our numerical scheme. We show PBE0 results for a variety of semiconductors and insulators. In comparison with experiment, the PBE0 functional leads to improved band gaps and an improved description of localized states. Even for the ferromagnetic semiconductor EuO with localized 4f electrons, the electronic and magnetic properties are correctly described by the PBE0 functional. Subsequently, we discuss the construction of the local, multiplicative exact exchange (EXX) potential from the non-local, orbital-dependent exact exchange energy. For this purpose we employ the optimized effective potential (OEP) method. Central ingredients of the OEP equation are the KS wave-function response and the single-particle density response function. We show that a balance between the LAPW

Orbital-dependent exchange-correlation functionals in density-functional theory realized by the FLAPW method

International Nuclear Information System (INIS)

Betzinger, Markus

2011-01-01

In this thesis, we extended the applicability of the full-potential linearized augmented-plane-wave (FLAPW) method, one of the most precise, versatile and generally applicable electronic structure methods for solids working within the framework of density-functional theory (DFT), to orbital-dependent functionals for the exchange-correlation (xc) energy. Two different schemes that deal with orbital-dependent functionals, the Kohn-Sham (KS) and the generalized Kohn-Sham (gKS) formalism, have been realized. Hybrid functionals, combining some amount of the orbital-dependent exact exchange energy with local or semi-local functionals of the density, are implemented within the gKS scheme. We work in particular with the PBE0 hybrid of Perdew, Burke, and Ernzerhof. Our implementation relies on a representation of the non-local exact exchange potential - its calculation constitutes the most time consuming step in a practical calculation - by an auxiliary mixed product basis (MPB). In this way, the matrix elements of the Hamiltonian corresponding to the non-local potential become a Brillouin-zone (BZ) sum over vector-matrix-vector products. Several techniques are developed and explored to further accelerate our numerical scheme. We show PBE0 results for a variety of semiconductors and insulators. In comparison with experiment, the PBE0 functional leads to improved band gaps and an improved description of localized states. Even for the ferromagnetic semiconductor EuO with localized 4f electrons, the electronic and magnetic properties are correctly described by the PBE0 functional. Subsequently, we discuss the construction of the local, multiplicative exact exchange (EXX) potential from the non-local, orbital-dependent exact exchange energy. For this purpose we employ the optimized effective potential (OEP) method. Central ingredients of the OEP equation are the KS wave-function response and the single-particle density response function. We show that a balance between the LAPW

Reduced density matrix functional theory at finite temperature

Energy Technology Data Exchange (ETDEWEB)

Baldsiefen, Tim

2012-10-15

Density functional theory (DFT) is highly successful in many fields of research. There are, however, areas in which its performance is rather limited. An important example is the description of thermodynamical variables of a quantum system in thermodynamical equilibrium. Although the finite-temperature version of DFT (FT-DFT) rests on a firm theoretical basis and is only one year younger than its brother, groundstate DFT, it has been successfully applied to only a few problems. Because FT-DFT, like DFT, is in principle exact, these shortcomings can be attributed to the difficulties of deriving valuable functionals for FT-DFT. In this thesis, we are going to present an alternative theoretical description of quantum systems in thermal equilibrium. It is based on the 1-reduced density matrix (1RDM) of the system, rather than on its density and will rather cumbersomly be called finite-temperature reduced density matrix functional theory (FT-RDMFT). Its zero-temperature counterpart (RDMFT) proved to be successful in several fields, formerly difficult to address via DFT. These fields include, for example, the calculation of dissociation energies or the calculation of the fundamental gap, also for Mott insulators. This success is mainly due to the fact that the 1RDM carries more directly accessible ''manybody'' information than the density alone, leading for example to an exact description of the kinetic energy functional. This sparks the hope that a description of thermodynamical systems employing the 1RDM via FT-RDMFT can yield an improvement over FT-DFT. Giving a short review of RDMFT and pointing out difficulties when describing spin-polarized systems initiates our work. We then lay the theoretical framework for FT-RDMFT by proving the required Hohenberg-Kohn-like theorems, investigating and determining the domain of FT-RDMFT functionals and by deriving several properties of the exact functional. Subsequently, we present a perturbative method to

Reduced density matrix functional theory at finite temperature

International Nuclear Information System (INIS)

Baldsiefen, Tim

2012-10-01

Density functional theory (DFT) is highly successful in many fields of research. There are, however, areas in which its performance is rather limited. An important example is the description of thermodynamical variables of a quantum system in thermodynamical equilibrium. Although the finite-temperature version of DFT (FT-DFT) rests on a firm theoretical basis and is only one year younger than its brother, groundstate DFT, it has been successfully applied to only a few problems. Because FT-DFT, like DFT, is in principle exact, these shortcomings can be attributed to the difficulties of deriving valuable functionals for FT-DFT. In this thesis, we are going to present an alternative theoretical description of quantum systems in thermal equilibrium. It is based on the 1-reduced density matrix (1RDM) of the system, rather than on its density and will rather cumbersomly be called finite-temperature reduced density matrix functional theory (FT-RDMFT). Its zero-temperature counterpart (RDMFT) proved to be successful in several fields, formerly difficult to address via DFT. These fields include, for example, the calculation of dissociation energies or the calculation of the fundamental gap, also for Mott insulators. This success is mainly due to the fact that the 1RDM carries more directly accessible ''manybody'' information than the density alone, leading for example to an exact description of the kinetic energy functional. This sparks the hope that a description of thermodynamical systems employing the 1RDM via FT-RDMFT can yield an improvement over FT-DFT. Giving a short review of RDMFT and pointing out difficulties when describing spin-polarized systems initiates our work. We then lay the theoretical framework for FT-RDMFT by proving the required Hohenberg-Kohn-like theorems, investigating and determining the domain of FT-RDMFT functionals and by deriving several properties of the exact functional. Subsequently, we present a perturbative method to iteratively construct

Comparison of exact-exchange calculations for solids in current-spin-density- and spin-density-functional theory

DEFF Research Database (Denmark)

Sharma, S.; Pittalis, S.; Kurth, S.

2007-01-01

The relative merits of current-spin-density- and spin-density-functional theory are investigated for solids treated within the exact-exchange-only approximation. Spin-orbit splittings and orbital magnetic moments are determined at zero external magnetic field. We find that for magnetic (Fe, Co......, and Ni) and nonmagnetic (Si and Ge) solids, the exact-exchange current-spin-density functional approach does not significantly improve the accuracy of the corresponding spin-density functional results....

Can Unrestricted Density-Functional Theory Describe Open Shell Singlet Biradicals?

Directory of Open Access Journals (Sweden)

Dieter Cremer

2002-04-01

Full Text Available Abstract: Unrestricted density functional theory (UDFT can be used for the description of open-shell singlet (OSS biradicals provided a number of precautions are considered. Biradicals that require a two-determinantal wave function (e.g. OSS state of carbenes cannot be described by UDFT for principal reasons. However, if the overlap between the open-shell orbitals is small (the single electrons are located at different atomic centers errors become small and, then, the principal failure of UDFT in these cases is not apparent and may even be disguised by the fact that UDFT has the advantage of describing spin polarization better than any restricted open shell DFT method. In the case of OSS biradicals with two- or multiconfigurational character (but a onedeterminantal form of the leading configuration, reasonable results can be obtained by broken-symmetry (BS-UDFT, however in each case this has to be checked. In no case is it reasonable to lower the symmetry of a molecule to get a suitable UDFT description. Hybrid functionals such as B3LYP perform better than pure DFT functionals in BS-UDFT calculations because the former reduce the self-interaction error of DFT exchange functionals, which mimics unspecified static electron correlation effects, so that the inclusion of specific static electron correlation effects via the form of the wavefunction becomes more effective.

Application of gradient-corrected density functional theory to the structures and thermochemistries of ScF3, TiF4, VF5, and CrF6

International Nuclear Information System (INIS)

Russo, T.V.; Martin, R.L.; Jeffrey Hay, P.

1995-01-01

Density functional theory (DFT) and Hartree--Fock (HF) calculations are reported for the family of transition metal fluorides ScF 3 , TiF 4 , VF 5 , and CrF 6 . Both HF and the local-density approximation (LDA) yield excellent agreement with experimental bond lengths, while the B-LYP gradient-corrected density functional gives bond lengths 0.04-0.05 A too long. An investigation of various combinations of exchange and correlation functionals shows that, for this series, the origin of this behavior lies in the Becke exchange functional. Much improved bond distances are found using the hybrid HF/DFT functional advocated by Becke. This approximation also leads to much improved thermochemistries. The LDA overestimates average bond energies in this series by 30-40 kcal/mol, whereas the B-LYP functional overbinds by only ∼8-12 kcal/mol, and the hybrid HF/DFT method overbinds by only ∼2 kcal/mol. The hybrid method predicts the octahedral isomer of CrF 6 to be more stable than the trigonal prismatic form by 14 kcal/mol. Comparison of theoretical vibrational frequencies with experiment supports the assignment of an octahedral geometry

Local hybrid functionals: An assessment for thermochemical kinetics

International Nuclear Information System (INIS)

Kaupp, Martin; Bahmann, Hilke; Arbuznikov, Alexei V.

2007-01-01

Local hybrid functionals with position-dependent exact-exchange admixture are a new class of exchange-correlation functionals in density functional theory that promise to advance the available accuracy in many areas of application. Local hybrids with different local mixing functions (LMFs) governing the position dependence are validated for the heats of formation of the extended G3/99 set, and for two sets of barriers of hydrogen-transfer and heavy-atom transfer reactions (HTBH38 and NHTBH38 databases). A simple local hybrid Lh-SVWN with only Slater and exact exchange plus local correlation and a one-parameter LMF, g(r)=b(τ W (r)/τ(r)), performs best and provides overall mean absolute errors for thermochemistry and kinetics that are a significant improvement over standard state-of-the-art global hybrid functionals. In particular, this local hybrid functional does not suffer from the systematic deterioration that standard functionals exhibit for larger molecules. In contrast, local hybrids based on generalized gradient approximation exchange tend to give rise to nonintuitive LMFs, and no improved functionals have been obtained along this route. The LMF is a real-space function and thus can be analyzed in detail. We use, in particular, graphical analyses to rationalize the performance of different local hybrids for thermochemistry and reaction barriers

Time-dependent density functional theory for multi-component systems

International Nuclear Information System (INIS)

Tiecheng Li; Peiqing Tong

1985-10-01

The Runge-Gross version of Hohenberg-Kohn-Sham's density functional theory is generalized to multi-component systems, both for arbitrary time-dependent pure states and for arbitrary time-dependent ensembles. (author)

Consistent Atomic Geometries and Electronic Structure of Five Phases of Potassium Niobate from Density-Functional Theory

Directory of Open Access Journals (Sweden)

Falko Schmidt

2017-01-01

Full Text Available We perform a comprehensive theoretical study of the structural and electronic properties of potassium niobate (KNbO3 in the cubic, tetragonal, orthorhombic, monoclinic, and rhombohedral phase, based on density-functional theory. The influence of different parametrizations of the exchange-correlation functional on the investigated properties is analyzed in detail, and the results are compared to available experimental data. We argue that the PBEsol and AM05 generalized gradient approximations as well as the RTPSS meta-generalized gradient approximation yield consistently accurate structural data for both the external and internal degrees of freedom and are overall superior to the local-density approximation or other conventional generalized gradient approximations for the structural characterization of KNbO3. Band-structure calculations using a HSE-type hybrid functional further indicate significant near degeneracies of band-edge states in all phases which are expected to be relevant for the optical response of the material.

Time-dependent density functional theory for many-electron systems interacting with cavity photons.

Science.gov (United States)

Tokatly, I V

2013-06-07

Time-dependent (current) density functional theory for many-electron systems strongly coupled to quantized electromagnetic modes of a microcavity is proposed. It is shown that the electron-photon wave function is a unique functional of the electronic (current) density and the expectation values of photonic coordinates. The Kohn-Sham system is constructed, which allows us to calculate the above basic variables by solving self-consistent equations for noninteracting particles. We suggest possible approximations for the exchange-correlation potentials and discuss implications of this approach for the theory of open quantum systems. In particular we show that it naturally leads to time-dependent density functional theory for systems coupled to the Caldeira-Leggett bath.

Approximate self-consistent potentials for density-functional-theory exchange-correlation functionals

International Nuclear Information System (INIS)

Cafiero, Mauricio; Gonzalez, Carlos

2005-01-01

We show that potentials for exchange-correlation functionals within the Kohn-Sham density-functional-theory framework may be written as potentials for simpler functionals multiplied by a factor close to unity, and in a self-consistent field calculation, these effective potentials find the correct self-consistent solutions. This simple theory is demonstrated with self-consistent exchange-only calculations of the atomization energies of some small molecules using the Perdew-Kurth-Zupan-Blaha (PKZB) meta-generalized-gradient-approximation (meta-GGA) exchange functional. The atomization energies obtained with our method agree with or surpass previous meta-GGA calculations performed in a non-self-consistent manner. The results of this work suggest the utility of this simple theory to approximate exchange-correlation potentials corresponding to energy functionals too complicated to generate closed forms for their potentials. We hope that this method will encourage the development of complex functionals which have correct boundary conditions and are free of self-interaction errors without the worry that the functionals are too complex to differentiate to obtain potentials

Solvation in atomic liquids: connection between Gaussian field theory and density functional theory

Directory of Open Access Journals (Sweden)

V. Sergiievskyi

2017-12-01

Full Text Available For the problem of molecular solvation, formulated as a liquid submitted to the external potential field created by a molecular solute of arbitrary shape dissolved in that solvent, we draw a connection between the Gaussian field theory derived by David Chandler [Phys. Rev. E, 1993, 48, 2898] and classical density functional theory. We show that Chandler's results concerning the solvation of a hard core of arbitrary shape can be recovered by either minimising a linearised HNC functional using an auxiliary Lagrange multiplier field to impose a vanishing density inside the core, or by minimising this functional directly outside the core — indeed a simpler procedure. Those equivalent approaches are compared to two other variants of DFT, either in the HNC, or partially linearised HNC approximation, for the solvation of a Lennard-Jones solute of increasing size in a Lennard-Jones solvent. Compared to Monte-Carlo simulations, all those theories give acceptable results for the inhomogeneous solvent structure, but are completely out-of-range for the solvation free-energies. This can be fixed in DFT by adding a hard-sphere bridge correction to the HNC functional.

Covariant density functional theory: The role of the pion

International Nuclear Information System (INIS)

Lalazissis, G. A.; Karatzikos, S.; Serra, M.; Otsuka, T.; Ring, P.

2009-01-01

We investigate the role of the pion in covariant density functional theory. Starting from conventional relativistic mean field (RMF) theory with a nonlinear coupling of the σ meson and without exchange terms we add pions with a pseudovector coupling to the nucleons in relativistic Hartree-Fock approximation. In order to take into account the change of the pion field in the nuclear medium the effective coupling constant of the pion is treated as a free parameter. It is found that the inclusion of the pion to this sort of density functionals does not destroy the overall description of the bulk properties by RMF. On the other hand, the noncentral contribution of the pion (tensor coupling) does have effects on single particle energies and on binding energies of certain nuclei.

Ground-State Gas-Phase Structures of Inorganic Molecules Predicted by Density Functional Theory Methods

KAUST Repository

Minenkov, Yury

2017-11-29

We tested a battery of density functional theory (DFT) methods ranging from generalized gradient approximation (GGA) via meta-GGA to hybrid meta-GGA schemes as well as Møller–Plesset perturbation theory of the second order and a single and double excitation coupled-cluster (CCSD) theory for their ability to reproduce accurate gas-phase structures of di- and triatomic molecules derived from microwave spectroscopy. We obtained the most accurate molecular structures using the hybrid and hybrid meta-GGA approximations with B3PW91, APF, TPSSh, mPW1PW91, PBE0, mPW1PBE, B972, and B98 functionals, resulting in lowest errors. We recommend using these methods to predict accurate three-dimensional structures of inorganic molecules when intramolecular dispersion interactions play an insignificant role. The structures that the CCSD method predicts are of similar quality although at considerably larger computational cost. The structures that GGA and meta-GGA schemes predict are less accurate with the largest absolute errors detected with BLYP and M11-L, suggesting that these methods should not be used if accurate three-dimensional molecular structures are required. Because of numerical problems related to the integration of the exchange–correlation part of the functional and large scattering of errors, most of the Minnesota models tested, particularly MN12-L, M11, M06-L, SOGGA11, and VSXC, are also not recommended for geometry optimization. When maintaining a low computational budget is essential, the nonseparable gradient functional N12 might work within an acceptable range of error. As expected, the DFT-D3 dispersion correction had a negligible effect on the internuclear distances when combined with the functionals tested on nonweakly bonded di- and triatomic inorganic molecules. By contrast, the dispersion correction for the APF-D functional has been found to shorten the bonds significantly, up to 0.064 Å (AgI), in Ag halides, BaO, BaS, BaF, BaCl, Cu halides, and Li and

Evaluation of exchange-correlation functionals for time-dependent density functional theory calculations on metal complexes.

Science.gov (United States)

Holland, Jason P; Green, Jennifer C

2010-04-15

The electronic absorption spectra of a range of copper and zinc complexes have been simulated by using time-dependent density functional theory (TD-DFT) calculations implemented in Gaussian03. In total, 41 exchange-correlation (XC) functionals including first-, second-, and third-generation (meta-generalized gradient approximation) DFT methods were compared in their ability to predict the experimental electronic absorption spectra. Both pure and hybrid DFT methods were tested and differences between restricted and unrestricted calculations were also investigated by comparison of analogous neutral zinc(II) and copper(II) complexes. TD-DFT calculated spectra were optimized with respect to the experimental electronic absorption spectra by use of a Matlab script. Direct comparison of the performance of each XC functional was achieved both qualitatively and quantitatively by comparison of optimized half-band widths, root-mean-squared errors (RMSE), energy scaling factors (epsilon(SF)), and overall quality-of-fit (Q(F)) parameters. Hybrid DFT methods were found to outperform all pure DFT functionals with B1LYP, B97-2, B97-1, X3LYP, and B98 functionals providing the highest quantitative and qualitative accuracy in both restricted and unrestricted systems. Of the functionals tested, B1LYP gave the most accurate results with both average RMSE and overall Q(F) 0.990) for the copper complexes. The XC functional performance in spin-restricted TD-DFT calculations on the zinc complexes was found to be slightly worse. PBE1PBE, mPW1PW91 and B1LYP gave the most accurate results with typical RMSE and Q(F) values between 5.3 and 7.3%, and epsilon(SF) around 0.930. These studies illustrate the power of modern TD-DFT calculations for exploring excited state transitions of metal complexes. 2009 Wiley Periodicals, Inc.

Density functional theory of the electrical double layer: the RFD functional

International Nuclear Information System (INIS)

Gillespie, Dirk; Valisko, Monika; Boda, Dezso

2005-01-01

Density functional theory (DFT) of electrolytes is applied to the electrical double layer under a wide range of conditions. The ions are charged, hard spheres of different size and valence, and the wall creating the double layer is uncharged, weakly charged, and strongly charged. Under all conditions, the density and electrostatic potential profiles calculated using the recently proposed RFD electrostatic functional (Gillespie et al 2002 J. Phys.: Condens. Matter 14 12129; 2003 Phys. Rev. E 68 031503) compare well to Monte Carlo simulations. When the wall is strongly charged, the RFD functional results agree with the results of a simpler perturbative electrostatic DFT, but the two functionals' results qualitatively disagree when the wall is uncharged or weakly charged. The RFD functional reproduces these phenomena of weakly charged double layers. It also reproduces bulk thermodynamic quantities calculated from pair correlation functions

Dispersion correction derived from first principles for density functional theory and Hartree-Fock theory.

Science.gov (United States)

Guidez, Emilie B; Gordon, Mark S

2015-03-12

The modeling of dispersion interactions in density functional theory (DFT) is commonly performed using an energy correction that involves empirically fitted parameters for all atom pairs of the system investigated. In this study, the first-principles-derived dispersion energy from the effective fragment potential (EFP) method is implemented for the density functional theory (DFT-D(EFP)) and Hartree-Fock (HF-D(EFP)) energies. Overall, DFT-D(EFP) performs similarly to the semiempirical DFT-D corrections for the test cases investigated in this work. HF-D(EFP) tends to underestimate binding energies and overestimate intermolecular equilibrium distances, relative to coupled cluster theory, most likely due to incomplete accounting for electron correlation. Overall, this first-principles dispersion correction yields results that are in good agreement with coupled-cluster calculations at a low computational cost.

Comments on the locality in density-functional theory

International Nuclear Information System (INIS)

Lindgren, Ingvar; Salomonson, Sten

2003-01-01

The 'locality hypothesis' in density-functional theory (DFT), implying that the functional derivative is equivalent to a multiplicative local function, forms the basis of models of Kohn-Sham type. This has been generally accepted by the community since the advent of the model, and has later been formally proved for a large class of functionals. The hypothesis has recently been questioned by Nesbet [Phys. Rev. A 58, R12 (1998) and Phys. Rev. A 65, 010502 (2001)], who claims that it fails for the kinetic-energy functional for a system with more than two noninteracting electrons with a nondegenerate ground state. This conclusion has been questioned by Gal [Phys. Rev. A 62, 044501 (2000)] and by Holas and March [Phys. Rev. A 64, 016501 (2001)]. We claim that the arguments of Nesbet are incorrect, since the orbital functional used for the kinetic energy is not a unique functional of the total density in the domain of unnormalized orbitals. We have demonstrated that with a proper definition of the kinetic energy, which is a unique density functional also in the unnormalized region, the derivative can be represented by a single local multiplicative function for all v-representable densities. Therefore, we consider the controversy connected with the issue raised by Nesbet as resolved. We believe that the proof of the differentiability given here can be extended to larger groups of DFT functionals, and works along these lines are in progress

Conversion of localized lower hybrid oscillations and fast magnetosonic waves at a plasma density cavity

International Nuclear Information System (INIS)

Hall, J.O.

2004-01-01

Analytic expressions are presented for conversion of localized lower hybrid oscillations and magnetosonic waves by scattering off a small scale density cavity. The governing equations are solved in slab geometry with wave vectors perpendicular to both the ambient magnetic field and the density gradient associated with density cavity using a scale length separation method. The theory predicts strong excitation of localized lower hybrid oscillations for a set of frequencies between the lower hybrid frequency of the ambient plasma and the minimum lower hybrid frequency inside the cavity. The theory is relevant for the lower hybrid solitary structures observed in space plasmas

Statistical theory of electron densities

International Nuclear Information System (INIS)

Pratt, L.R.; Hoffman, G.G.; Harris, R.A.

1988-01-01

An optimized Thomas--Fermi theory is proposed which retains the simplicity of the original theory and is a suitable reference theory for Monte Carlo density functional treatments of condensed materials. The key ingredient of the optimized theory is a neighborhood sampled potential which contains effects of the inhomogeneities in the one-electron potential. In contrast to the traditional Thomas--Fermi approach, the optimized theory predicts a finite electron density in the vicinity of a nucleus. Consideration of the example of an ideal electron gas subject to a central Coulomb field indicates that implementation of the approach is straightforward. The optimized theory is found to fail completely when a classically forbidden region is approached. However, these circumstances are not of primary interest for calculations of interatomic forces. It is shown how the energy functional of the density may be constructed by integration of a generalized Hellmann--Feynman relation. This generalized Hellmann--Feynman relation proves to be equivalent to the variational principle of density functional quantum mechanics, and, therefore, the present density theory can be viewed as a variational consequence of the constructed energy functional

Triplet state photochemistry and the three-state crossing of acetophenone within time-dependent density-functional theory

Energy Technology Data Exchange (ETDEWEB)

Huix-Rotllant, Miquel, E-mail: miquel.huix@gmail.com; Ferré, Nicolas, E-mail: nicolas.ferre@univ-amu.fr [Institut de Chimie Radicalaire (UMR-7273), Aix-Marseille Université, CNRS, 13397 Marseille Cedex 20 (France)

2014-04-07

Even though time-dependent density-functional theory (TDDFT) works generally well for describing excited states energies and properties in the Franck-Condon region, it can dramatically fail in predicting photochemistry, notably when electronic state crossings occur. Here, we assess the ability of TDDFT to describe the photochemistry of an important class of triplet sensitizers, namely, aromatic ketones. We take acetophenone as a test molecule, for which accurate ab initio results exist in the literature. Triplet acetophenone is generated thanks to an exotic three-state crossing involving one singlet and two triplets states (i.e., a simultaneous intersystem crossing and triplet conical intersection), thus being a stringent test for approximate TDDFT. We show that most exchange-correlation functionals can only give a semi-qualitative picture of the overall photochemistry, in which the three-state crossing is rather represented as a triplet conical intersection separated from the intersystem crossing. The best result overall is given by the double hybrid functional mPW2PLYP, which is even able to reproduce quantitatively the three-state crossing region. We rationalize this results by noting that double hybrid functionals include a larger portion of double excitation character to the excited states.

Quantum Crystallography: Density Matrix-Density Functional Theory and the X-Ray Diffraction Experiment

Science.gov (United States)

Soirat, Arnaud J. A.

Density Matrix Theory is a Quantum Mechanical formalism in which the wavefunction is eliminated and its role taken over by reduced density matrices. The interest of this is that, it allows one, in principle, to calculate any electronic property of a physical system, without having to solve the Schrodinger equation, using only two entities much simpler than an N-body wavefunction: first and second -order reduced density matrices. In practice, though, this very promising possibility faces the tremendous theoretical problem of N-representability, which has been solved for the former, but, until now, voids any hope of theoretically determining the latter. However, it has been shown that single determinant reduced density matrices of any order may be recovered from coherent X-ray diffraction data, if one provides a proper Quantum Mechanical description of the Crystallography experiment. A deeper investigation of this method is the purpose of this work, where we, first, further study the calculation of X-ray reduced density matrices N-representable by a single Slater determinant. In this context, we independently derive necessary and sufficient conditions for the uniqueness of the method. We then show how to account for electron correlation in this model. For the first time, indeed, we derive highly accurate, yet practical, density matrices approximately N-representable by correlated-determinant wavefunctions. The interest of such a result lies in the Quantum Mechanical validity of these density matrices, their property of being entirely obtainable from X-ray coherent diffraction data, their very high accuracy conferred by this known property of the N-representing wavefunction, as well as their definition as explicit functionals of the density. All of these properties are finally used in both a theoretical and a numerical application: in the former, we show that these density matrices may be used in the context of Density Functional Theory to highly accurately determine

Multireference Density Functional Theory with Generalized Auxiliary Systems for Ground and Excited States.

Science.gov (United States)

Chen, Zehua; Zhang, Du; Jin, Ye; Yang, Yang; Su, Neil Qiang; Yang, Weitao

2017-09-21

To describe static correlation, we develop a new approach to density functional theory (DFT), which uses a generalized auxiliary system that is of a different symmetry, such as particle number or spin, from that of the physical system. The total energy of the physical system consists of two parts: the energy of the auxiliary system, which is determined with a chosen density functional approximation (DFA), and the excitation energy from an approximate linear response theory that restores the symmetry to that of the physical system, thus rigorously leading to a multideterminant description of the physical system. The electron density of the physical system is different from that of the auxiliary system and is uniquely determined from the functional derivative of the total energy with respect to the external potential. Our energy functional is thus an implicit functional of the physical system density, but an explicit functional of the auxiliary system density. We show that the total energy minimum and stationary states, describing the ground and excited states of the physical system, can be obtained by a self-consistent optimization with respect to the explicit variable, the generalized Kohn-Sham noninteracting density matrix. We have developed the generalized optimized effective potential method for the self-consistent optimization. Among options of the auxiliary system and the associated linear response theory, reformulated versions of the particle-particle random phase approximation (pp-RPA) and the spin-flip time-dependent density functional theory (SF-TDDFT) are selected for illustration of principle. Numerical results show that our multireference DFT successfully describes static correlation in bond dissociation and double bond rotation.

Efficient molecular density functional theory using generalized spherical harmonics expansions.

Science.gov (United States)

Ding, Lu; Levesque, Maximilien; Borgis, Daniel; Belloni, Luc

2017-09-07

We show that generalized spherical harmonics are well suited for representing the space and orientation molecular density in the resolution of the molecular density functional theory. We consider the common system made of a rigid solute of arbitrary complexity immersed in a molecular solvent, both represented by molecules with interacting atomic sites and classical force fields. The molecular solvent density ρ(r,Ω) around the solute is a function of the position r≡(x,y,z) and of the three Euler angles Ω≡(θ,ϕ,ψ) describing the solvent orientation. The standard density functional, equivalent to the hypernetted-chain closure for the solute-solvent correlations in the liquid theory, is minimized with respect to ρ(r,Ω). The up-to-now very expensive angular convolution products are advantageously replaced by simple products between projections onto generalized spherical harmonics. The dramatic gain in speed of resolution enables to explore in a systematic way molecular solutes of up to nanometric sizes in arbitrary solvents and to calculate their solvation free energy and associated microscopic solvent structure in at most a few minutes. We finally illustrate the formalism by tackling the solvation of molecules of various complexities in water.

Study of homogeneous bubble nucleation in liquid carbon dioxide by a hybrid approach combining molecular dynamics simulation and density gradient theory

Science.gov (United States)

Langenbach, K.; Heilig, M.; Horsch, M.; Hasse, H.

2018-03-01

A new method for predicting homogeneous bubble nucleation rates of pure compounds from vapor-liquid equilibrium (VLE) data is presented. It combines molecular dynamics simulation on the one side with density gradient theory using an equation of state (EOS) on the other. The new method is applied here to predict bubble nucleation rates in metastable liquid carbon dioxide (CO2). The molecular model of CO2 is taken from previous work of our group. PC-SAFT is used as an EOS. The consistency between the molecular model and the EOS is achieved by adjusting the PC-SAFT parameters to VLE data obtained from the molecular model. The influence parameter of density gradient theory is fitted to the surface tension of the molecular model. Massively parallel molecular dynamics simulations are performed close to the spinodal to compute bubble nucleation rates. From these simulations, the kinetic prefactor of the hybrid nucleation theory is estimated, whereas the nucleation barrier is calculated from density gradient theory. This enables the extrapolation of molecular simulation data to the whole metastable range including technically relevant densities. The results are tested against available experimental data and found to be in good agreement. The new method does not suffer from typical deficiencies of classical nucleation theory concerning the thermodynamic barrier at the spinodal and the bubble size dependence of surface tension, which is typically neglected in classical nucleation theory. In addition, the density in the center of critical bubbles and their surface tension is determined as a function of their radius. The usual linear Tolman correction to the capillarity approximation is found to be invalid.

Study of homogeneous bubble nucleation in liquid carbon dioxide by a hybrid approach combining molecular dynamics simulation and density gradient theory.

Science.gov (United States)

Langenbach, K; Heilig, M; Horsch, M; Hasse, H

2018-03-28

A new method for predicting homogeneous bubble nucleation rates of pure compounds from vapor-liquid equilibrium (VLE) data is presented. It combines molecular dynamics simulation on the one side with density gradient theory using an equation of state (EOS) on the other. The new method is applied here to predict bubble nucleation rates in metastable liquid carbon dioxide (CO 2 ). The molecular model of CO 2 is taken from previous work of our group. PC-SAFT is used as an EOS. The consistency between the molecular model and the EOS is achieved by adjusting the PC-SAFT parameters to VLE data obtained from the molecular model. The influence parameter of density gradient theory is fitted to the surface tension of the molecular model. Massively parallel molecular dynamics simulations are performed close to the spinodal to compute bubble nucleation rates. From these simulations, the kinetic prefactor of the hybrid nucleation theory is estimated, whereas the nucleation barrier is calculated from density gradient theory. This enables the extrapolation of molecular simulation data to the whole metastable range including technically relevant densities. The results are tested against available experimental data and found to be in good agreement. The new method does not suffer from typical deficiencies of classical nucleation theory concerning the thermodynamic barrier at the spinodal and the bubble size dependence of surface tension, which is typically neglected in classical nucleation theory. In addition, the density in the center of critical bubbles and their surface tension is determined as a function of their radius. The usual linear Tolman correction to the capillarity approximation is found to be invalid.

General framework for fluctuating dynamic density functional theory

Science.gov (United States)

Durán-Olivencia, Miguel A.; Yatsyshin, Peter; Goddard, Benjamin D.; Kalliadasis, Serafim

2017-12-01

We introduce a versatile bottom-up derivation of a formal theoretical framework to describe (passive) soft-matter systems out of equilibrium subject to fluctuations. We provide a unique connection between the constituent-particle dynamics of real systems and the time evolution equation of their measurable (coarse-grained) quantities, such as local density and velocity. The starting point is the full Hamiltonian description of a system of colloidal particles immersed in a fluid of identical bath particles. Then, we average out the bath via Zwanzig’s projection-operator techniques and obtain the stochastic Langevin equations governing the colloidal-particle dynamics. Introducing the appropriate definition of the local number and momentum density fields yields a generalisation of the Dean-Kawasaki (DK) model, which resembles the stochastic Navier-Stokes description of a fluid. Nevertheless, the DK equation still contains all the microscopic information and, for that reason, does not represent the dynamical law of observable quantities. We address this controversial feature of the DK description by carrying out a nonequilibrium ensemble average. Adopting a natural decomposition into local-equilibrium and nonequilibrium contribution, where the former is related to a generalised version of the canonical distribution, we finally obtain the fluctuating-hydrodynamic equation governing the time-evolution of the mesoscopic density and momentum fields. Along the way, we outline the connection between the ad hoc energy functional introduced in previous DK derivations and the free-energy functional from classical density-functional theory. The resultant equation has the structure of a dynamical density-functional theory (DDFT) with an additional fluctuating force coming from the random interactions with the bath. We show that our fluctuating DDFT formalism corresponds to a particular version of the fluctuating Navier-Stokes equations, originally derived by Landau and Lifshitz

Derivation of the density functional theory from the cluster expansion.

Science.gov (United States)

Hsu, J Y

2003-09-26

The density functional theory is derived from a cluster expansion by truncating the higher-order correlations in one and only one term in the kinetic energy. The formulation allows self-consistent calculation of the exchange correlation effect without imposing additional assumptions to generalize the local density approximation. The pair correlation is described as a two-body collision of bound-state electrons, and modifies the electron- electron interaction energy as well as the kinetic energy. The theory admits excited states, and has no self-interaction energy.

Density functional theory in surface science and heterogeneous catalysis

DEFF Research Database (Denmark)

Nørskov, Jens Kehlet; Scheffler, M.; Toulhoat, H.

2006-01-01

Solid surfaces are used extensively as catalysts throughout the chemical industry, in the energy sector, and in environmental protection. Recently, density functional theory has started providing new insight into the atomic-scale mechanisms of heterogeneous catalysis, helping to interpret the large...

Density-functional theory simulation of large quantum dots

Science.gov (United States)

Jiang, Hong; Baranger, Harold U.; Yang, Weitao

2003-10-01

Kohn-Sham spin-density functional theory provides an efficient and accurate model to study electron-electron interaction effects in quantum dots, but its application to large systems is a challenge. Here an efficient method for the simulation of quantum dots using density-function theory is developed; it includes the particle-in-the-box representation of the Kohn-Sham orbitals, an efficient conjugate-gradient method to directly minimize the total energy, a Fourier convolution approach for the calculation of the Hartree potential, and a simplified multigrid technique to accelerate the convergence. We test the methodology in a two-dimensional model system and show that numerical studies of large quantum dots with several hundred electrons become computationally affordable. In the noninteracting limit, the classical dynamics of the system we study can be continuously varied from integrable to fully chaotic. The qualitative difference in the noninteracting classical dynamics has an effect on the quantum properties of the interacting system: integrable classical dynamics leads to higher-spin states and a broader distribution of spacing between Coulomb blockade peaks.

Constructive definition of functional derivatives in density-functional theory

International Nuclear Information System (INIS)

Luo Ji

2006-01-01

It is shown that the functional derivatives in density-functional theory (DFT) can be explicitly defined within the domain of electron densities restricted by the electron number, and a constructive definition of such restricted derivatives is suggested. With this definition, Kohn-Sham (KS) equations can be established for an N-electron system without extending the functional domain and introducing a Lagrange multiplier. This may clarify some of the fundamental questions raised by Nesbet (1998 Phys. Rev. A 58 R12). The definition naturally leads to the fact that the KS effective potential is determined only to within an additive constant, thus the KS levels can shift freely and the relation between the highest occupied molecular orbital (HOMO) energy and the ionization potential of the system depends on the choice of the constant. On the other hand, if the domain of functionals is indeed extended beyond the electron number restriction, conclusions depend on whether the extended functionals have unrestricted derivatives or not. It is shown that the ensemble extension of DFT to open systems of mixed states (Perdew et al 1982 Phys. Rev. Lett. 49 1691) leads to an energy functional which has no unrestricted derivative at integer electron numbers. Hence after this extension, the relation between the HOMO energy and the ionization potential for an N-electron system is still uncertain. Besides, there are different extensions of the energy functional to a domain of densities unrestricted by the integer electron number, resulting in different unrestricted derivatives and electron systems with different chemical potentials. Even for the exact exchange-correlation potential, there is still an undetermined constant, whether it is a restricted or unrestricted derivative

Covariant density functional theory for nuclear matter

Energy Technology Data Exchange (ETDEWEB)

Badarch, U.

2007-07-01

The present thesis is organized as follows. In Chapter 2 we study the Nucleon-Nucleon (NN) interaction in Dirac-Brueckner (DB) approach. We start by considering the NN interaction in free-space in terms of the Bethe-Salpeter (BS) equation to the meson exchange potential model. Then we present the DB approach for nuclear matter by extending the BS equation for the in-medium NN interaction. From the solution of the three-dimensional in-medium BS equation, we derive the DB self-energies and total binding energy which are the main results of the DB approach, which we later incorporate in the field theoretical calculation of the nuclear equation of state. In Chapter 3, we introduce the basic concepts of density functional theory in the context of Quantum Hadrodynamics (QHD-I). We reach the main point of this work in Chapter 4 where we introduce the DDRH approach. In the DDRH theory, the medium dependence of the meson-nucleon vertices is expressed as functionals of the baryon field operators. Because of the complexities of the operator-valued functionals we decide to use the mean-field approximation. In Chapter 5, we contrast microscopic and phenomenological approaches to extracting density dependent meson-baryon vertices. Chapter 6 gives the results of our studies of the EOS of infinite nuclear matter in detail. Using formulas derived in Chapters 4 and 5 we calculate the properties of symmetric and asymmetric nuclear matter and pure neutron matter. (orig.)

Covariant density functional theory for nuclear matter

International Nuclear Information System (INIS)

Badarch, U.

2007-01-01

The present thesis is organized as follows. In Chapter 2 we study the Nucleon-Nucleon (NN) interaction in Dirac-Brueckner (DB) approach. We start by considering the NN interaction in free-space in terms of the Bethe-Salpeter (BS) equation to the meson exchange potential model. Then we present the DB approach for nuclear matter by extending the BS equation for the in-medium NN interaction. From the solution of the three-dimensional in-medium BS equation, we derive the DB self-energies and total binding energy which are the main results of the DB approach, which we later incorporate in the field theoretical calculation of the nuclear equation of state. In Chapter 3, we introduce the basic concepts of density functional theory in the context of Quantum Hadrodynamics (QHD-I). We reach the main point of this work in Chapter 4 where we introduce the DDRH approach. In the DDRH theory, the medium dependence of the meson-nucleon vertices is expressed as functionals of the baryon field operators. Because of the complexities of the operator-valued functionals we decide to use the mean-field approximation. In Chapter 5, we contrast microscopic and phenomenological approaches to extracting density dependent meson-baryon vertices. Chapter 6 gives the results of our studies of the EOS of infinite nuclear matter in detail. Using formulas derived in Chapters 4 and 5 we calculate the properties of symmetric and asymmetric nuclear matter and pure neutron matter. (orig.)

Description of plasmon-like band in silver clusters: the importance of the long-range Hartree-Fock exchange in time-dependent density-functional theory simulations.

Science.gov (United States)

Rabilloud, Franck

2014-10-14

Absorption spectra of Ag20 and Ag55(q) (q = +1, -3) nanoclusters are investigated in the framework of the time-dependent density functional theory in order to analyse the role of the d electrons in plasmon-like band of silver clusters. The description of the plasmon-like band from calculations using density functionals containing an amount of Hartree-Fock exchange at long range, namely, hybrid and range-separated hybrid (RSH) density functionals, is in good agreement with the classical interpretation of the plasmon-like structure as a collective excitation of valence s-electrons. In contrast, using local or semi-local exchange functionals (generalized gradient approximations (GGAs) or meta-GGAs) leads to a strong overestimation of the role of d electrons in the plasmon-like band. The semi-local asymptotically corrected model potentials also describe the plasmon as mainly associated to d electrons, though calculated spectra are in fairly good agreement with those calculated using the RSH scheme. Our analysis shows that a portion of non-local exchange modifies the description of the plasmon-like band.

Local and linear chemical reactivity response functions at finite temperature in density functional theory

International Nuclear Information System (INIS)

Franco-Pérez, Marco; Ayers, Paul W.; Gázquez, José L.; Vela, Alberto

2015-01-01

We explore the local and nonlocal response functions of the grand canonical potential density functional at nonzero temperature. In analogy to the zero-temperature treatment, local (e.g., the average electron density and the local softness) and nonlocal (e.g., the softness kernel) intrinsic response functions are defined as partial derivatives of the grand canonical potential with respect to its thermodynamic variables (i.e., the chemical potential of the electron reservoir and the external potential generated by the atomic nuclei). To define the local and nonlocal response functions of the electron density (e.g., the Fukui function, the linear density response function, and the dual descriptor), we differentiate with respect to the average electron number and the external potential. The well-known mathematical relationships between the intrinsic response functions and the electron-density responses are generalized to nonzero temperature, and we prove that in the zero-temperature limit, our results recover well-known identities from the density functional theory of chemical reactivity. Specific working equations and numerical results are provided for the 3-state ensemble model

Inclusion of Dispersion Effects in Density Functional Theory

DEFF Research Database (Denmark)

Møgelhøj, Andreas

on fitting to high-level ab initio and experimental results. The fitting scheme, based on Baysian theory, focuses on the three aspects: a) model space, b) datasets, and c) model selection. The model space consists of a flexible expansion of the exchange enhancement factor in the generalized gradient......In this thesis, applications and development will be presented within the field of van der Waals interactions in density functional theory. The thesis is based on the three projects: i) van der Waals interactions effect on the structure of liquid water at ambient conditions, ii) development......-range van der Waals interactions is essential to describe the adsorption/desorption process and commonly used generalized gradient approximation functionals are seen to be incapable of this....

Increasing the applicability of density functional theory. IV. Consequences of ionization-potential improved exchange-correlation potentials.

Science.gov (United States)

Verma, Prakash; Bartlett, Rodney J

2014-05-14

This paper's objective is to create a "consistent" mean-field based Kohn-Sham (KS) density functional theory (DFT) meaning the functional should not only provide good total energy properties, but also the corresponding KS eigenvalues should be accurate approximations to the vertical ionization potentials (VIPs) of the molecule, as the latter condition attests to the viability of the exchange-correlation potential (VXC). None of the prominently used DFT approaches show these properties: the optimized effective potential VXC based ab initio dft does. A local, range-separated hybrid potential cam-QTP-00 is introduced as the basis for a "consistent" KS DFT approach. The computed VIPs as the negative of KS eigenvalue have a mean absolute error of 0.8 eV for an extensive set of molecule's electron ionizations, including the core. Barrier heights, equilibrium geometries, and magnetic properties obtained from the potential are in good agreement with experiment. A similar accuracy with less computational efforts can be achieved by using a non-variational global hybrid variant of the QTP-00 approach.

Describing long-range charge-separation processes with subsystem density-functional theory

Energy Technology Data Exchange (ETDEWEB)

Solovyeva, Alisa; Neugebauer, Johannes, E-mail: j.neugebauer@uni-muenster.de [Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster (Germany); Pavanello, Michele, E-mail: m.pavanello@rutgers.edu [Department of Chemistry, Rutgers University, 73 Warren St., Newark, New Jersey 07102 (United States)

2014-04-28

Long-range charge-transfer processes in extended systems are difficult to describe with quantum chemical methods. In particular, cost-effective (non-hybrid) approximations within time-dependent density functional theory (DFT) are not applicable unless special precautions are taken. Here, we show that the efficient subsystem DFT can be employed as a constrained DFT variant to describe the energetics of long-range charge-separation processes. A formal analysis of the energy components in subsystem DFT for such excitation energies is presented, which demonstrates that both the distance dependence and the long-range limit are correctly described. In addition, electronic couplings for these processes as needed for rate constants in Marcus theory can be obtained from this method. It is shown that the electronic structure of charge-separated states constructed by a positively charged subsystem interacting with a negatively charged one is difficult to converge — charge leaking from the negative subsystem to the positive one can occur. This problem is related to the delocalization error in DFT and can be overcome with asymptotically correct exchange–correlation (XC) potentials or XC potentials including a sufficiently large amount of exact exchange. We also outline an approximate way to obtain charge-transfer couplings between locally excited and charge-separated states.

Describing long-range charge-separation processes with subsystem density-functional theory

International Nuclear Information System (INIS)

Solovyeva, Alisa; Neugebauer, Johannes; Pavanello, Michele

2014-01-01

Long-range charge-transfer processes in extended systems are difficult to describe with quantum chemical methods. In particular, cost-effective (non-hybrid) approximations within time-dependent density functional theory (DFT) are not applicable unless special precautions are taken. Here, we show that the efficient subsystem DFT can be employed as a constrained DFT variant to describe the energetics of long-range charge-separation processes. A formal analysis of the energy components in subsystem DFT for such excitation energies is presented, which demonstrates that both the distance dependence and the long-range limit are correctly described. In addition, electronic couplings for these processes as needed for rate constants in Marcus theory can be obtained from this method. It is shown that the electronic structure of charge-separated states constructed by a positively charged subsystem interacting with a negatively charged one is difficult to converge — charge leaking from the negative subsystem to the positive one can occur. This problem is related to the delocalization error in DFT and can be overcome with asymptotically correct exchange–correlation (XC) potentials or XC potentials including a sufficiently large amount of exact exchange. We also outline an approximate way to obtain charge-transfer couplings between locally excited and charge-separated states

International Workshop on Electronic Density Functional Theory : Recent Progress and New Directions

CERN Document Server

Vignale, Giovanni; Das, Mukunda

1998-01-01

This book is an outcome of the International Workshop on Electronic Density Functional Theory, held at Griffith University in Brisbane, Australia, in July 1996. Density functional theory, standing as it does at the boundary between the disciplines of physics, chemistry, and materials science, is a great mixer. Invited experts from North America, Europe, and Australia mingled with students from several disciplines, rapidly taking up the informal style for which Australia is famous. A list of participants is given at the end of the book. Density functional theory (DFT) is a subtle approach to the very difficult problem of predicting the behavior of many interacting particles. A major application is the study of many-electron systems. This was the workshop theme, embracing inter alia computational chemistry and condensed matter physics. DFT circumvents the more conceptually straightforward (but more computationally intensive) approach in which one solves the many-body Schrodinger equation. It relies instead on r...

Density Functional Theory An Advanced Course

CERN Document Server

Dreizler, Reiner M

2011-01-01

Density Functional Theory (DFT) has firmly established itself as the workhorse for the atomic-level simulation of condensed matter phases, pure or composite materials and quantum chemical systems. The present book is a rigorous and detailed introduction to the foundations up to and including such advanced topics as orbital-dependent functionals and both time-dependent and relativistic DFT. Given the many ramifications of contemporary DFT, this text concentrates on the self-contained presentation of the basics of the most widely used DFT variants. This implies a thorough discussion of the corresponding existence theorems and effective single particle equations, as well as of key approximations utilized in implementations. The formal results are complemented by selected quantitative results, which primarily aim at illustrating strengths and weaknesses of a particular approach or functional. DFT for superconducting or nuclear and hadronic systems are not addressed in this work. The structure and material contain...

Grand canonical electronic density-functional theory: Algorithms and applications to electrochemistry

International Nuclear Information System (INIS)

Sundararaman, Ravishankar; Goddard, William A. III; Arias, Tomas A.

2017-01-01

First-principles calculations combining density-functional theory and continuum solvation models enable realistic theoretical modeling and design of electrochemical systems. When a reaction proceeds in such systems, the number of electrons in the portion of the system treated quantum mechanically changes continuously, with a balancing charge appearing in the continuum electrolyte. A grand-canonical ensemble of electrons at a chemical potential set by the electrode potential is therefore the ideal description of such systems that directly mimics the experimental condition. We present two distinct algorithms: a self-consistent field method and a direct variational free energy minimization method using auxiliary Hamiltonians (GC-AuxH), to solve the Kohn-Sham equations of electronic density-functional theory directly in the grand canonical ensemble at fixed potential. Both methods substantially improve performance compared to a sequence of conventional fixed-number calculations targeting the desired potential, with the GC-AuxH method additionally exhibiting reliable and smooth exponential convergence of the grand free energy. Lastly, we apply grand-canonical density-functional theory to the under-potential deposition of copper on platinum from chloride-containing electrolytes and show that chloride desorption, not partial copper monolayer formation, is responsible for the second voltammetric peak.

Grand canonical electronic density-functional theory: Algorithms and applications to electrochemistry

Science.gov (United States)

Sundararaman, Ravishankar; Goddard, William A.; Arias, Tomas A.

2017-03-01

First-principles calculations combining density-functional theory and continuum solvation models enable realistic theoretical modeling and design of electrochemical systems. When a reaction proceeds in such systems, the number of electrons in the portion of the system treated quantum mechanically changes continuously, with a balancing charge appearing in the continuum electrolyte. A grand-canonical ensemble of electrons at a chemical potential set by the electrode potential is therefore the ideal description of such systems that directly mimics the experimental condition. We present two distinct algorithms: a self-consistent field method and a direct variational free energy minimization method using auxiliary Hamiltonians (GC-AuxH), to solve the Kohn-Sham equations of electronic density-functional theory directly in the grand canonical ensemble at fixed potential. Both methods substantially improve performance compared to a sequence of conventional fixed-number calculations targeting the desired potential, with the GC-AuxH method additionally exhibiting reliable and smooth exponential convergence of the grand free energy. Finally, we apply grand-canonical density-functional theory to the under-potential deposition of copper on platinum from chloride-containing electrolytes and show that chloride desorption, not partial copper monolayer formation, is responsible for the second voltammetric peak.

Time-dependent quantum fluid density functional theory of hydrogen ...

Indian Academy of Sciences (India)

A time-dependent generalized non-linear Schrödinger equation (GNLSE) of motion was earlier derived in our laboratory by combining density functional theory and quantum fluid dynamics in threedimensional space. In continuation of the work reported previously, the GNLSE is applied to provide additional knowledge on ...

Extended screened exchange functional derived from transcorrelated density functional theory.

Science.gov (United States)

Umezawa, Naoto

2017-09-14

We propose a new formulation of the correlation energy functional derived from the transcorrelated method in use in density functional theory (TC-DFT). An effective Hamiltonian, H TC , is introduced by a similarity transformation of a many-body Hamiltonian, H, with respect to a complex function F: H TC =1FHF. It is proved that an expectation value of H TC for a normalized single Slater determinant, D n , corresponds to the total energy: E[n] = ⟨Ψ n |H|Ψ n ⟩/⟨Ψ n |Ψ n ⟩ = ⟨D n |H TC |D n ⟩ under the two assumptions: (1) The electron density nr associated with a trial wave function Ψ n = D n F is v-representable and (2) Ψ n and D n give rise to the same electron density nr. This formulation, therefore, provides an alternative expression of the total energy that is useful for the development of novel correlation energy functionals. By substituting a specific function for F, we successfully derived a model correlation energy functional, which resembles the functional form of the screened exchange method. The proposed functional, named the extended screened exchange (ESX) functional, is described within two-body integrals and is parametrized for a numerically exact correlation energy of the homogeneous electron gas. The ESX functional does not contain any ingredients of (semi-)local functionals and thus is totally free from self-interactions. The computational cost for solving the self-consistent-field equation is comparable to that of the Hartree-Fock method. We apply the ESX functional to electronic structure calculations for a solid silicon, H - ion, and small atoms. The results demonstrate that the TC-DFT formulation is promising for the systematic improvement of the correlation energy functional.

Density functional and many-body theories of Hydrogen plasmas

International Nuclear Information System (INIS)

Perrot, F.; Dharma-Wardana, M.W.C.

1983-11-01

This work is an attempt to go beyond the standard description of hot condensed matter using the well-known ''average atom model''. The first part describes a static model using ''Density functional theory'' to calculate self-consistent coupled electron and ion density profiles of the plasma not restricted to a single average atomic sphere. In a second part, the results are used as ingredients for a many-body approach to electronic properties: the one-particle Green-function self-energy is calculated, from which shifted levels, populations and level-widths are deduced. Results for the Hydrogen plasma are reported, with emphasis on the 1s bound state

Time-dependent density-functional theory concepts and applications

CERN Document Server

Ullrich, Carsten A

2011-01-01

Time-dependent density-functional theory (TDDFT) describes the quantum dynamics of interacting electronic many-body systems formally exactly and in a practical and efficient manner. TDDFT has become the leading method for calculating excitation energies and optical properties of large molecules, with accuracies that rival traditional wave-function based methods, but at a fraction of the computational cost.This book is the first graduate-level text on the concepts and applications of TDDFT, including many examples and exercises, and extensive coverage of the literature. The book begins with a s

Optical rotation calculated with time-dependent density functional theory: the OR45 benchmark.

Science.gov (United States)

Srebro, Monika; Govind, Niranjan; de Jong, Wibe A; Autschbach, Jochen

2011-10-13

Time-dependent density functional theory (TDDFT) computations are performed for 42 organic molecules and three transition metal complexes, with experimental molar optical rotations ranging from 2 to 2 × 10(4) deg cm(2) dmol(-1). The performances of the global hybrid functionals B3LYP, PBE0, and BHLYP, and of the range-separated functionals CAM-B3LYP and LC-PBE0 (the latter being fully long-range corrected), are investigated. The performance of different basis sets is studied. When compared to liquid-phase experimental data, the range-separated functionals do, on average, not perform better than B3LYP and PBE0. Median relative deviations between calculations and experiment range from 25 to 29%. A basis set recently proposed for optical rotation calculations (LPol-ds) on average does not give improved results compared to aug-cc-pVDZ in TDDFT calculations with B3LYP. Individual cases are discussed in some detail, among them norbornenone for which the LC-PBE0 functional produced an optical rotation that is close to available data from coupled-cluster calculations, but significantly smaller in magnitude than the liquid-phase experimental value. Range-separated functionals and BHLYP perform well for helicenes and helicene derivatives. Metal complexes pose a challenge to first-principles calculations of optical rotation.

Geometry-based density functional theory an overview

CERN Document Server

Schmidt, M

2003-01-01

An overview of recent developments and applications of a specific density functional approach that originates from Rosenfeld's fundamental measure theory for hard spheres is given. Model systems that were treated include penetrable spheres that interact with a step function pair potential, the Widom-Rowlinson model, the Asakura-Oosawa colloid-polymer mixture, ternary mixtures of spheres, needles, and globular polymers, hard-body amphiphilic mixtures, fluids in porous media, and random sequential adsorption that describes non-equilibrium processes such as colloidal deposition and random car parking. In these systems various physical phenomena were studied, such as correlations in liquids, freezing and demixing phase behaviour, the properties of fluid interfaces with and without orientational order, and wetting and layering phenomena at walls.

Geometry-based density functional theory: an overview

Science.gov (United States)

Schmidt, Matthias

2003-01-01

An overview of recent developments and applications of a specific density functional approach that originates from Rosenfeld's fundamental measure theory for hard spheres is given. Model systems that were treated include penetrable spheres that interact with a step function pair potential, the Widom-Rowlinson model, the Asakura-Oosawa colloid-polymer mixture, ternary mixtures of spheres, needles, and globular polymers, hard-body amphiphilic mixtures, fluids in porous media, and random sequential adsorption that describes non-equilibrium processes such as colloidal deposition and random car parking. In these systems various physical phenomena were studied, such as correlations in liquids, freezing and demixing phase behaviour, the properties of fluid interfaces with and without orientational order, and wetting and layering phenomena at walls.

Geometry-based density functional theory: an overview

International Nuclear Information System (INIS)

Schmidt, Matthias

2003-01-01

An overview of recent developments and applications of a specific density functional approach that originates from Rosenfeld's fundamental measure theory for hard spheres is given. Model systems that were treated include penetrable spheres that interact with a step function pair potential, the Widom-Rowlinson model, the Asakura-Oosawa colloid-polymer mixture, ternary mixtures of spheres, needles, and globular polymers, hard-body amphiphilic mixtures, fluids in porous media, and random sequential adsorption that describes non-equilibrium processes such as colloidal deposition and random car parking. In these systems various physical phenomena were studied, such as correlations in liquids, freezing and demixing phase behaviour, the properties of fluid interfaces with and without orientational order, and wetting and layering phenomena at walls

Reducing Systematic Errors in Oxide Species with Density Functional Theory Calculations

DEFF Research Database (Denmark)

Christensen, Rune; Hummelshøj, Jens S.; Hansen, Heine Anton

2015-01-01

Density functional theory calculations can be used to gain valuable insight into the fundamental reaction processes in metal−oxygen systems, e.g., metal−oxygen batteries. Here, the ability of a range of different exchange-correlation functionals to reproduce experimental enthalpies of formation...

2nd derivatives of the electronic energy in density functional theory

Energy Technology Data Exchange (ETDEWEB)

Dam, H. van

2001-08-01

This document details the equations needed to implement the calculation of vibrational frequencies within the density functional formalism of electronic structure theory. This functionality has been incorporated into the CCP1 DFT module and the required changes to the application programmers interface are outlined. Throughout it is assumed that an implementation of Hartree-Fock vibrational frequencies is available that can be modified to incorporate the density functional formalism. Employing GAMESS-UK as an example the required changes to the Hartree-Fock code are outlined. (author)

Exact-exchange time-dependent density-functional theory for static and dynamic polarizabilities

International Nuclear Information System (INIS)

Hirata, So; Ivanov, Stanislav; Bartlett, Rodney J.; Grabowski, Ireneusz

2005-01-01

Time-dependent density-functional theory (TDDFT) employing the exact-exchange functional has been formulated on the basis of the optimized-effective-potential (OEP) method of Talman and Shadwick for second-order molecular properties and implemented into a Gaussian-basis-set, trial-vector algorithm. The only approximation involved, apart from the lack of correlation effects and the use of Gaussian-type basis functions, was the consistent use of the adiabatic approximation in the exchange kernel and in the linear response function. The static and dynamic polarizabilities and their anisotropy predicted by the TDDFT with exact exchange (TDOEP) agree accurately with the corresponding values from time-dependent Hartree-Fock theory, the exact-exchange counterpart in the wave function theory. The TDOEP is free from the nonphysical asymptotic decay of the exchange potential of most conventional density functionals or from any other manifestations of the incomplete cancellation of the self-interaction energy. The systematic overestimation of the absolute values and dispersion of polarizabilities that plagues most conventional TDDFT cannot be seen in the TDOEP

Temperature-dependent study of isotropic-nematic transition for a Gay-Berne fluid using density-functional theory

International Nuclear Information System (INIS)

Singh, Ram Chandra

2007-01-01

We have used the density-functional theory to study the effect of varying temperature on the isotropic-nematic transition of a fluid of molecules interacting via the Gay-Berne intermolecular potential. The nematic phase is found to be stable with respect to isotropic phase in the temperature range 0.80≤T*≤1.25. Pair correlation functions needed as input information in density-functional theory is calculated using the Percus-Yevick integral equation theory. We find that the density-functional theory is good for studying the isotropic-nematic transition in molecular fluids if the values of the pair-correlation functions in the isotropic phase are known accurately. We have also compared our results with computer simulation results wherever they are available

Survival of the most transferable at the top of Jacob's ladder: Defining and testing the ωB97M(2) double hybrid density functional

Science.gov (United States)

Mardirossian, Narbe; Head-Gordon, Martin

2018-06-01

A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resulting functional, ωB97M(2), is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 data points clearly surpasses that of all of the tested density functionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.

Pressure and surface tension of soild-liquid interface using Tarazona density functional theory

Directory of Open Access Journals (Sweden)

M. M.

2000-12-01

Full Text Available   The weighted density functional theory proposed by Tarazona is applied to study the solid-liquid interface. In the last two decades the weighted density functional became a useful tool to consider the properties of inhomogeneous liquids. In this theory, the role of the size of molecules or the particles of which the matter is composed, was found to be important. In this resarch we study a hard sphere fluid beside a hard wall. For this study the liquid is an inhomogeneous system. We use the definition of the direct correlation function as a second derivative of free energy with respect to the density. We use this definition and the definition of the weighting function, then we minimize the grand potential with respect to the density to get the Euler Lagrange equation and we obtain an integral equation to find the inhomogeneous density profile. The obtained density profile as a function of the distance from the wall, for different bulk density is plotted in three dimensions. We also calculate the pressure and compare it with the Carnahan-starling results, and finally we obtained the surface tension at liquid-solid interface and compared it with the results of Monte Carlo simulation.

Density Functional Theory and Materials Modeling at Atomistic Length Scales

Directory of Open Access Journals (Sweden)

Swapan K. Ghosh

2002-04-01

Full Text Available Abstract: We discuss the basic concepts of density functional theory (DFT as applied to materials modeling in the microscopic, mesoscopic and macroscopic length scales. The picture that emerges is that of a single unified framework for the study of both quantum and classical systems. While for quantum DFT, the central equation is a one-particle Schrodinger-like Kohn-Sham equation, the classical DFT consists of Boltzmann type distributions, both corresponding to a system of noninteracting particles in the field of a density-dependent effective potential, the exact functional form of which is unknown. One therefore approximates the exchange-correlation potential for quantum systems and the excess free energy density functional or the direct correlation functions for classical systems. Illustrative applications of quantum DFT to microscopic modeling of molecular interaction and that of classical DFT to a mesoscopic modeling of soft condensed matter systems are highlighted.

Theoretical investigation of cyromazine tautomerism using density functional theory and Møller–Plesset perturbation theory methods

Science.gov (United States)

A computational chemistry analysis of six unique tautomers of cyromazine, a pesticide used for fly control, was performed with density functional theory (DFT) and canonical second order Møller–Plesset perturbation theory (MP2) methods to gain insight into the contributions of molecular structure to ...

Hybrid density functional study on the mechanism for the enhanced photocatalytic properties of the ultrathin hybrid layered nanocomposite g-C3N4/BiOCl

Science.gov (United States)

Yao, Wenzhi; Zhang, Jihua; Wang, Yuanxu; Ren, Fengzhu

2018-03-01

To investigate the origin of the high photocatalytic performance of experimentally synthesized g-C3N4/ BiOCl, we studied its geometry structure, electronic structure, and photocatalytic properties by means of hybrid density-functional theory (DFT). The calculated band alignment of g-C3N4 and few-layer BiOCl sheets clearly shows that g-C3N4/ BiOCl is a standard type-II nanocomposite. The density of states, Bader charge, partial charge density, charge density difference, and the effective masses show that electron-hole pair can be effectively separated in the g-C3N4/BiOCl interface. The calculated absorption coefficients indicate an obvious redshift of the absorption edge. The band gap of g-C3N4/BiOCl can be modulated by external electric field, and a semiconductor-semimetal transition is observed. The type-II vdW heterostructure is still maintained during the changes of external electric field. Especially, when the electric field reaches to +0.7 V/Å, the impurity states have been eliminated with the band gap of 2.3 eV. An analysis of optical properties shows that the absorption coefficient in the visible-light region is enhanced considerably as the electric-field strength increases. Our calculation results suggest that the ultrathin hybrid layered g-C3N4/BiOCl nanocomposite may have significant advantages for visible-light photocatalysis.

Existence of time-dependent density-functional theory for open electronic systems: time-dependent holographic electron density theorem.

Science.gov (United States)

Zheng, Xiao; Yam, ChiYung; Wang, Fan; Chen, GuanHua

2011-08-28

We present the time-dependent holographic electron density theorem (TD-HEDT), which lays the foundation of time-dependent density-functional theory (TDDFT) for open electronic systems. For any finite electronic system, the TD-HEDT formally establishes a one-to-one correspondence between the electron density inside any finite subsystem and the time-dependent external potential. As a result, any electronic property of an open system in principle can be determined uniquely by the electron density function inside the open region. Implications of the TD-HEDT on the practicality of TDDFT are also discussed.

Self-consistent hybrid functionals for solids: a fully-automated implementation

Science.gov (United States)

Erba, A.

2017-08-01

A fully-automated algorithm for the determination of the system-specific optimal fraction of exact exchange in self-consistent hybrid functionals of the density-functional-theory is illustrated, as implemented into the public Crystal program. The exchange fraction of this new class of functionals is self-consistently updated proportionally to the inverse of the dielectric response of the system within an iterative procedure (Skone et al 2014 Phys. Rev. B 89, 195112). Each iteration of the present scheme, in turn, implies convergence of a self-consistent-field (SCF) and a coupled-perturbed-Hartree-Fock/Kohn-Sham (CPHF/KS) procedure. The present implementation, beside improving the user-friendliness of self-consistent hybrids, exploits the unperturbed and electric-field perturbed density matrices from previous iterations as guesses for subsequent SCF and CPHF/KS iterations, which is documented to reduce the overall computational cost of the whole process by a factor of 2.

Structure of solvent-free grafted nanoparticles: Molecular dynamics and density-functional theory

KAUST Repository

Chremos, Alexandros

2011-01-01

The structure of solvent-free oligomer-grafted nanoparticles has been investigated using molecular dynamics simulations and density-functional theory. At low temperatures and moderate to high oligomer lengths, the qualitative features of the core particle pair probability, structure factor, and the oligomer brush configuration obtained from the simulations can be explained by a density-functional theory that incorporates the configurational entropy of the space-filling oligomers. In particular, the structure factor at small wave numbers attains a value much smaller than the corresponding hard-sphere suspension, the first peak of the pair distribution function is enhanced due to entropic attractions among the particles, and the oligomer brush expands with decreasing particle volume fraction to fill the interstitial space. At higher temperatures, the simulations reveal effects that differ from the theory and are likely caused by steric repulsions of the expanded corona chains. © 2011 American Institute of Physics.

Many-body perturbation theory using the density-functional concept: beyond the GW approximation.

Science.gov (United States)

Bruneval, Fabien; Sottile, Francesco; Olevano, Valerio; Del Sole, Rodolfo; Reining, Lucia

2005-05-13

We propose an alternative formulation of many-body perturbation theory that uses the density-functional concept. Instead of the usual four-point integral equation for the polarizability, we obtain a two-point one, which leads to excellent optical absorption and energy-loss spectra. The corresponding three-point vertex function and self-energy are then simply calculated via an integration, for any level of approximation. Moreover, we show the direct impact of this formulation on the time-dependent density-functional theory. Numerical results for the band gap of bulk silicon and solid argon illustrate corrections beyond the GW approximation for the self-energy.

Affordable and accurate large-scale hybrid-functional calculations on GPU-accelerated supercomputers

Science.gov (United States)

Ratcliff, Laura E.; Degomme, A.; Flores-Livas, José A.; Goedecker, Stefan; Genovese, Luigi

2018-03-01

Performing high accuracy hybrid functional calculations for condensed matter systems containing a large number of atoms is at present computationally very demanding or even out of reach if high quality basis sets are used. We present a highly optimized multiple graphics processing unit implementation of the exact exchange operator which allows one to perform fast hybrid functional density-functional theory (DFT) calculations with systematic basis sets without additional approximations for up to a thousand atoms. With this method hybrid DFT calculations of high quality become accessible on state-of-the-art supercomputers within a time-to-solution that is of the same order of magnitude as traditional semilocal-GGA functionals. The method is implemented in a portable open-source library.

FDE-vdW: A van der Waals inclusive subsystem density-functional theory

Energy Technology Data Exchange (ETDEWEB)

Kevorkyants, Ruslan; Pavanello, Michele, E-mail: m.pavanello@rutgers.edu [Department of Chemistry, Rutgers University, Newark, New Jersey 07102 (United States); Eshuis, Henk [Department of Chemistry and Biochemistry, Montclair State University, Montclair, New Jersey 07043 (United States)

2014-07-28

We present a formally exact van der Waals inclusive electronic structure theory, called FDE-vdW, based on the Frozen Density Embedding formulation of subsystem Density-Functional Theory. In subsystem DFT, the energy functional is composed of subsystem additive and non-additive terms. We show that an appropriate definition of the long-range correlation energy is given by the value of the non-additive correlation functional. This functional is evaluated using the fluctuation–dissipation theorem aided by a formally exact decomposition of the response functions into subsystem contributions. FDE-vdW is derived in detail and several approximate schemes are proposed, which lead to practical implementations of the method. We show that FDE-vdW is Casimir-Polder consistent, i.e., it reduces to the generalized Casimir-Polder formula for asymptotic inter-subsystems separations. Pilot calculations of binding energies of 13 weakly bound complexes singled out from the S22 set show a dramatic improvement upon semilocal subsystem DFT, provided that an appropriate exchange functional is employed. The convergence of FDE-vdW with basis set size is discussed, as well as its dependence on the choice of associated density functional approximant.

Ground-state properties of rare-earth metals: an evaluation of density-functional theory

International Nuclear Information System (INIS)

Söderlind, Per; Turchi, P E A; Landa, A; Lordi, V

2014-01-01

The rare-earth metals have important technological applications due to their magnetic properties, but are scarce and expensive. Development of high-performance magnetic materials with less rare-earth content is desired, but theoretical modeling is hampered by complexities of the rare earths electronic structure. The existence of correlated (atomic-like) 4f electrons in the vicinity of the valence band makes any first-principles theory challenging. Here, we apply and evaluate the efficacy of density-functional theory for the series of lanthanides (rare earths), investigating the influence of the electron exchange and correlation functional, spin-orbit interaction, and orbital polarization. As a reference, the results are compared with those of the so-called ‘standard model’ of the lanthanides in which electrons are constrained to occupy 4f core states with no hybridization with the valence electrons. Some comparisons are also made with models designed for strong electron correlations. Our results suggest that spin–orbit coupling and orbital polarization are important, particularly for the magnitude of the magnetic moments, and that calculated equilibrium volumes, bulk moduli, and magnetic moments show correct trends overall. However, the precision of the calculated properties is not at the level of that found for simpler metals in the Periodic Table of Elements, and the electronic structures do not accurately reproduce x-ray photoemission spectra. (paper)

Molecular Excitation Energies from Time-Dependent Density Functional Theory Employing Random-Phase Approximation Hessians with Exact Exchange.

Science.gov (United States)

Heßelmann, Andreas

2015-04-14

Molecular excitation energies have been calculated with time-dependent density-functional theory (TDDFT) using random-phase approximation Hessians augmented with exact exchange contributions in various orders. It has been observed that this approach yields fairly accurate local valence excitations if combined with accurate asymptotically corrected exchange-correlation potentials used in the ground-state Kohn-Sham calculations. The inclusion of long-range particle-particle with hole-hole interactions in the kernel leads to errors of 0.14 eV only for the lowest excitations of a selection of three alkene, three carbonyl, and five azabenzene molecules, thus surpassing the accuracy of a number of common TDDFT and even some wave function correlation methods. In the case of long-range charge-transfer excitations, the method typically underestimates accurate reference excitation energies by 8% on average, which is better than with standard hybrid-GGA functionals but worse compared to range-separated functional approximations.

Computational complexity of time-dependent density functional theory

International Nuclear Information System (INIS)

Whitfield, J D; Yung, M-H; Tempel, D G; Aspuru-Guzik, A; Boixo, S

2014-01-01

Time-dependent density functional theory (TDDFT) is rapidly emerging as a premier method for solving dynamical many-body problems in physics and chemistry. The mathematical foundations of TDDFT are established through the formal existence of a fictitious non-interacting system (known as the Kohn–Sham system), which can reproduce the one-electron reduced probability density of the actual system. We build upon these works and show that on the interior of the domain of existence, the Kohn–Sham system can be efficiently obtained given the time-dependent density. We introduce a V-representability parameter which diverges at the boundary of the existence domain and serves to quantify the numerical difficulty of constructing the Kohn-Sham potential. For bounded values of V-representability, we present a polynomial time quantum algorithm to generate the time-dependent Kohn–Sham potential with controllable error bounds. (paper)

Melting slope of MgO from molecular dynamics and density functional theory

Science.gov (United States)

Tangney, Paul; Scandolo, Sandro

2009-09-01

We combine density functional theory (DFT) with molecular dynamics simulations based on an accurate atomistic force field to calculate the pressure derivative of the melting temperature of magnesium oxide at ambient pressure—a quantity for which a serious disagreement between theory and experiment has existed for almost 15 years. We find reasonable agreement with previous DFT results and with a very recent experimental determination of the slope. We pay particular attention to areas of possible weakness in theoretical calculations and conclude that the long-standing discrepancy with experiment could only be explained by a dramatic failure of existing density functionals or by flaws in the original experiment.

Semiclassical neutral atom as a reference system in density functional theory.

Science.gov (United States)

Constantin, Lucian A; Fabiano, E; Laricchia, S; Della Sala, F

2011-05-06

We use the asymptotic expansions of the semiclassical neutral atom as a reference system in density functional theory to construct accurate generalized gradient approximations (GGAs) for the exchange-correlation and kinetic energies without any empiricism. These asymptotic functionals are among the most accurate GGAs for molecular systems, perform well for solid state, and overcome current GGA state of the art in frozen density embedding calculations. Our results also provide evidence for the conjointness conjecture between exchange and kinetic energies of atomic systems.

Optical properties of Al nanostructures from time dependent density functional theory

KAUST Repository

Mokkath, Junais Habeeb; Schwingenschlö gl, Udo

2016-01-01

The optical properties of Al nanostructures are investigated by means of time dependent density functional theory, considering chains of varying length and ladders/stripes of varying aspect ratio. The absorption spectra show redshifting

Thermodynamics as a Foundation for Density Functional Theory

International Nuclear Information System (INIS)

Argaman, Nathan

2014-01-01

Density Functional Theory (DFT) is the method of choice for an ever increasing number of electronic structure computations (recently reaching 30,000 publications per year). It was founded in the sixties on the basis of the Hohenberg-Kohn theorem and the Kohn-Sham equations, which were originally proved and derived for electronic ground states. Alternatively, one may use thermodynamics to derive DFT for finite-temperature ensembles, with the ground-state theory recovered in the zero temperature limit. Specifically, the transformation from chemical potential µ to electron number N as a free variable may be directly generalized to clarify how DFT uses the density distribution n(r), rather than the external potential v(r), to specify a particular inhomogeneous electronic system. Relating interacting and non-interacting systems with the same n(r) distribution, one recovers not only the Kohn-Sham formulation, but also the so-called adiabatic connection theorem, which gives an explicit expression for the exchange-correlation energy in terms of the 'exchangecorrelation hole.' This derivation has the advantage of being constructive, rather than being based on a reductio ad absurdum argument. It thus serves as an excellent basis for a discussion of the approximations which are inevitably introduced, including the Local Density Approximation (LDA) and the Generalized Gradient Approximation (GGA)

A density functional theory investigation of the electronic structure and spin moments of magnetite

KAUST Repository

Noh, Junghyun

2014-08-01

We present the results of density functional theory (DFT) calculations on magnetite, Fe3O4, which has been recently considered as electrode in the emerging field of organic spintronics. Given the nature of the potential applications, we evaluated the magnetite room-temperature cubic phase in terms of structural, electronic, and magnetic properties. We considered GGA (PBE), GGA + U (PBE + U), and range-separated hybrid (HSE06 and HSE(15%)) functionals. Calculations using HSE06 and HSE(15%) functionals underline the impact that inclusion of exact exchange has on the electronic structure. While the modulation of the band gap with exact exchange has been seen in numerous situations, the dramatic change in the valence band nature and states near the Fermi level has major implications for even a qualitative interpretation of the DFT results. We find that HSE06 leads to highly localized states below the Fermi level while HSE(15%) and PBE + U result in delocalized states around the Fermi level. The significant differences in local magnetic moments and atomic charges indicate that describing room-temperature bulk materials, surfaces and interfaces may require different functionals than their low-temperature counterparts.

A density functional theory investigation of the electronic structure and spin moments of magnetite

KAUST Repository

Noh, Junghyun; Osman, Osman I; Aziz, Saadullah G; Winget, Paul; Bredas, Jean-Luc

2014-01-01

We present the results of density functional theory (DFT) calculations on magnetite, Fe3O4, which has been recently considered as electrode in the emerging field of organic spintronics. Given the nature of the potential applications, we evaluated the magnetite room-temperature cubic phase in terms of structural, electronic, and magnetic properties. We considered GGA (PBE), GGA + U (PBE + U), and range-separated hybrid (HSE06 and HSE(15%)) functionals. Calculations using HSE06 and HSE(15%) functionals underline the impact that inclusion of exact exchange has on the electronic structure. While the modulation of the band gap with exact exchange has been seen in numerous situations, the dramatic change in the valence band nature and states near the Fermi level has major implications for even a qualitative interpretation of the DFT results. We find that HSE06 leads to highly localized states below the Fermi level while HSE(15%) and PBE + U result in delocalized states around the Fermi level. The significant differences in local magnetic moments and atomic charges indicate that describing room-temperature bulk materials, surfaces and interfaces may require different functionals than their low-temperature counterparts.

Basic concepts of Density Functional Theory: Electronic structure calculation

International Nuclear Information System (INIS)

Sharma, B. Indrajit

2016-01-01

We are looking for a material which possesses the required properties as demanded for technological applications. For this we have to repeat the preparation of the appropriate materials and its characterizations. So, before proceeding to experiments, one can study on computer generated structure and predict the properties of the desired material. To do this, a concept of Density Functional Theory comes out. (paper)

A Hybrid Density Functional Theory/Molecular Mechanics Approach for Linear Response Properties in Heterogeneous Environments.

Science.gov (United States)

Rinkevicius, Zilvinas; Li, Xin; Sandberg, Jaime A R; Mikkelsen, Kurt V; Ågren, Hans

2014-03-11

We introduce a density functional theory/molecular mechanical approach for computation of linear response properties of molecules in heterogeneous environments, such as metal surfaces or nanoparticles embedded in solvents. The heterogeneous embedding environment, consisting from metallic and nonmetallic parts, is described by combined force fields, where conventional force fields are used for the nonmetallic part and capacitance-polarization-based force fields are used for the metallic part. The presented approach enables studies of properties and spectra of systems embedded in or placed at arbitrary shaped metallic surfaces, clusters, or nanoparticles. The capability and performance of the proposed approach is illustrated by sample calculations of optical absorption spectra of thymidine absorbed on gold surfaces in an aqueous environment, where we study how different organizations of the gold surface and how the combined, nonadditive effect of the two environments is reflected in the optical absorption spectrum.

The QTP family of consistent functionals and potentials in Kohn-Sham density functional theory

Energy Technology Data Exchange (ETDEWEB)

Jin, Yifan; Bartlett, Rodney J., E-mail: bartlett@qtp.ufl.edu [Quantum Theory Project and Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611 (United States)

2016-07-21

This manuscript presents the second, consistent density functional in the QTP (Quantum Theory Project) family, that is, the CAM-QTP(01). It is a new range-separated exchange-correlation functional in which the non-local exchange contribution is 100% at large separation. It follows the same basic principles of this family that the Kohn-Sham eigenvalues of the occupied orbitals approximately equal the vertical ionization energies, which is not fulfilled by most of the traditional density functional methods. This new CAM-QTP(01) functional significantly improves the accuracy of the vertical excitation energies especially for the Rydberg states in the test set. It also reproduces many other properties such as geometries, reaction barrier heights, and atomization energies.

Differentiability in density-functional theory: Further study of the locality theorem

International Nuclear Information System (INIS)

Lindgren, Ingvar; Salomonson, Sten

2004-01-01

The locality theorem in density-functional theory (DFT) states that the functional derivative of the Hohenberg-Kohn universal functional can be expressed as a local multiplicative potential function, and this is the basis of DFT and of the successful Kohn-Sham model. Nesbet has in several papers [Phys. Rev. A 58, R12 (1998); ibid.65, 010502 (2001); Adv. Quant. Chem, 43, 1 (2003)] claimed that this theorem is in conflict with fundamental quantum physics, and as a consequence that the Hohenberg-Kohn theory cannot be generally valid. We have commented upon these works [Comment, Phys. Rev. A 67, 056501 (2003)] and recently extended the arguments [Adv. Quantum Chem. 43, 95 (2003)]. We have shown that there is no such conflict and that the locality theorem is inherently exact. In the present work we have furthermore verified this numerically by constructing a local Kohn-Sham potential for the 1s2s 3 S state of helium that generates the many-body electron density and shown that the corresponding 2s Kohn-Sham orbital eigenvalue agrees with the ionization energy to nine digits. Similar result is obtained with the Hartree-Fock density. Therefore, in addition to verifying the locality theorem, this result also confirms the so-called ionization-potential theorem

Correlation functional in screened-exchange density functional theory procedures.

Science.gov (United States)

Chan, Bun; Kawashima, Yukio; Hirao, Kimihiko

2017-10-15

In the present study, we have explored several prospects for the further development of screened-exchange density functional theory (SX-DFT) procedures. Using the performance of HSE06 as our measure, we find that the use of alternative correlation functionals (as oppose to PBEc in HSE06) also yields adequate results for a diverse set of thermochemical properties. We have further examined the performance of new SX-DFT procedures (termed HSEB-type methods) that comprise the HSEx exchange and a (near-optimal) reparametrized B97c (c OS,0  = c SS,0  = 1, c OS,1  = -1.5, c OS,2  = -0.644, c SS,1  = -0.5, and c SS,2  = 1.10) correlation functionals. The different variants of HSEB all perform comparably to or slightly better than the original HSE-type procedures. These results, together with our fundamental analysis of correlation functionals, point toward various directions for advancing SX-DFT methods. © 2017 Wiley Periodicals, Inc. © 2017 Wiley Periodicals, Inc.

Microhartree precision in density functional theory calculations

Science.gov (United States)

Gulans, Andris; Kozhevnikov, Anton; Draxl, Claudia

2018-04-01

To address ultimate precision in density functional theory calculations we employ the full-potential linearized augmented plane-wave + local-orbital (LAPW + lo) method and justify its usage as a benchmark method. LAPW + lo and two completely unrelated numerical approaches, the multiresolution analysis (MRA) and the linear combination of atomic orbitals, yield total energies of atoms with mean deviations of 0.9 and 0.2 μ Ha , respectively. Spectacular agreement with the MRA is reached also for total and atomization energies of the G2-1 set consisting of 55 molecules. With the example of α iron we demonstrate the capability of LAPW + lo to reach μ Ha /atom precision also for periodic systems, which allows also for the distinction between the numerical precision and the accuracy of a given functional.

Spin-adapted open-shell time-dependent density functional theory. II. Theory and pilot application.

Science.gov (United States)

Li, Zhendong; Liu, Wenjian; Zhang, Yong; Suo, Bingbing

2011-04-07

The excited states of open-shell systems calculated by unrestricted Kohn-Sham-based time-dependent density functional theory (U-TD-DFT) are often heavily spin-contaminated and hence meaningless. This is solved ultimately by the recently proposed spin-adapted time-dependent density functional theory (TD-DFT) (S-TD-DFT) [J. Chem. Phys. 133, 064106 (2010)]. Unlike the standard restricted open-shell Kohn-Sham-based TD-DFT (R-TD-DFT) which can only access the singlet-coupled single excitations, the S-TD-DFT can capture both the singlet- and triplet-coupled single excitations with the same computational effort as the U-TD-DFT. The performances of the three approaches (U-TD-DFT, R-TD-DFT, and S-TD-DFT) are compared for both the spin-conserving and spin-flip excitations of prototypical open-shell systems, the nitrogen (N(2)(+)) and naphthalene (C(10)H(8)(+)) cations. The results show that the S-TD-DFT gives rise to balanced descriptions of excited states of open-shell systems.

Radical scavenging activity of some natural tropolones by density functional theory

Directory of Open Access Journals (Sweden)

A. G. Al-Sehemi

2017-07-01

Full Text Available The ground state neutral geometries of some natural tropolones, i.e. stipitatonic acid (AF1, stipitalide (AF2, stipitaldehydic acid (AF3 and methyl stipitate (AF4 have been optimized by using Density Functional Theory (DFT at B3LYP/6-31G*, B3LYP/6-31G**, B3LYP/6-31+G* and B3LYP/6-31+G** levels of theory. The excited state geometries of AF1-AF4 were optimized by adopting the Time Dependent Density Functional Theory (TDDFT at the same levels of theory. The frequencies and cation species of AF1-AF4 were also computed at all the above mentioned levels of theory. We shed light on the electro-optical and molecular properties, e.g. energy gaps, highest occupied molecular orbitals, lowest unoccupied molecular orbitals, absorption wavelengths, electronegativity (χ, hardness (η, electrophilicity (ω, softness (S, electrophilicity index (ωi and the radical scavenging activity (RSA. Hydrogen atom transfer (HAT and one-electron transfer mechanisms have been discussed to shed light on the RSA. The smallest ionization potential and bond dissociation energy of AF4 are revealing that this compound would have more RSA than those of other counterparts.

Structures, electronic properties and magnetisms of FeBN (N ≤ 15) clusters: density functional theory investigations

International Nuclear Information System (INIS)

Liu Huoyan; Lel Xueling; Chen Hang; Liu Zhifeng; Liu Liren; Zhu Hengjiang

2011-01-01

The equilibrium structures, electronic properties and magnetisms of FeB N (N ≤ 15) clusters have been investigated by generalized gradient approximation (GGA) of density functional theory at different spin multiplicities. The average atomic binding energies, second-order energy differences and gaps of ground-state structures are calculated and discussed. The results show that FeB 3 , FeB 8 , FeB 12 and FeB 14 possess relatively higher stabilities. Moreover, there is a distinct hybridization between the d orbital of Fe and the p orbital of B for the ground-state cluster. The total magnetic moment for groundstate cluster is mainly provided by 3 d orbital of Fe atom, and exhibits the odd-even oscillation tendency with the increasing of cluster size. (authors)

A J matrix engine for density functional theory calculations

International Nuclear Information System (INIS)

White, C.A.; Head-Gordon, M.

1996-01-01

We introduce a new method for the formation of the J matrix (Coulomb interaction matrix) within a basis of Cartesian Gaussian functions, as needed in density functional theory and Hartree endash Fock calculations. By summing the density matrix into the underlying Gaussian integral formulas, we have developed a J matrix open-quote open-quote engine close-quote close-quote which forms the exact J matrix without explicitly forming the full set of two electron integral intermediates. Several precomputable quantities have been identified, substantially reducing the number of floating point operations and memory accesses needed in a J matrix calculation. Initial timings indicate a speedup of greater than four times for the (pp parallel pp) class of integrals with speedups increasing to over ten times for (ff parallel ff) integrals. copyright 1996 American Institute of Physics

Multi-configuration time-dependent density-functional theory based on range separation

DEFF Research Database (Denmark)

Fromager, E.; Knecht, S.; Jensen, Hans Jørgen Aagaard

2013-01-01

Multi-configuration range-separated density-functional theory is extended to the time-dependent regime. An exact variational formulation is derived. The approximation, which consists in combining a long-range Multi-Configuration- Self-Consistent Field (MCSCF) treatment with an adiabatic short...... (srGGA) approximations. As expected, when modeling long-range interactions with the MCSCF model instead of the adiabatic Buijse-Baerends density-matrix functional as recently proposed by Pernal [J. Chem. Phys. 136, 184105 (2012)10.1063/1.4712019], the description of both the 1D doubly-excited state...

Perspective: Fifty years of density-functional theory in chemical physics

International Nuclear Information System (INIS)

Becke, Axel D.

2014-01-01

Since its formal inception in 1964–1965, Kohn-Sham density-functional theory (KS-DFT) has become the most popular electronic structure method in computational physics and chemistry. Its popularity stems from its beautifully simple conceptual framework and computational elegance. The rise of KS-DFT in chemical physics began in earnest in the mid 1980s, when crucial developments in its exchange-correlation term gave the theory predictive power competitive with well-developed wave-function methods. Today KS-DFT finds itself under increasing pressure to deliver higher and higher accuracy and to adapt to ever more challenging problems. If we are not mindful, however, these pressures may submerge the theory in the wave-function sea. KS-DFT might be lost. I am hopeful the Kohn-Sham philosophical, theoretical, and computational framework can be preserved. This Perspective outlines the history, basic concepts, and present status of KS-DFT in chemical physics, and offers suggestions for its future development

Perspective: Fifty years of density-functional theory in chemical physics

Energy Technology Data Exchange (ETDEWEB)

Becke, Axel D., E-mail: axel.becke@dal.ca [Department of Chemistry, Dalhousie University, 6274 Coburg Rd., P.O. Box 15000, Halifax, Nova Scotia B3H 4R2 (Canada)

2014-05-14

Since its formal inception in 1964–1965, Kohn-Sham density-functional theory (KS-DFT) has become the most popular electronic structure method in computational physics and chemistry. Its popularity stems from its beautifully simple conceptual framework and computational elegance. The rise of KS-DFT in chemical physics began in earnest in the mid 1980s, when crucial developments in its exchange-correlation term gave the theory predictive power competitive with well-developed wave-function methods. Today KS-DFT finds itself under increasing pressure to deliver higher and higher accuracy and to adapt to ever more challenging problems. If we are not mindful, however, these pressures may submerge the theory in the wave-function sea. KS-DFT might be lost. I am hopeful the Kohn-Sham philosophical, theoretical, and computational framework can be preserved. This Perspective outlines the history, basic concepts, and present status of KS-DFT in chemical physics, and offers suggestions for its future development.

Density-functional theory for fluid-solid and solid-solid phase transitions.

Science.gov (United States)

Bharadwaj, Atul S; Singh, Yashwant

2017-03-01

We develop a theory to describe solid-solid phase transitions. The density functional formalism of classical statistical mechanics is used to find an exact expression for the difference in the grand thermodynamic potentials of the two coexisting phases. The expression involves both the symmetry conserving and the symmetry broken parts of the direct pair correlation function. The theory is used to calculate phase diagram of systems of soft spheres interacting via inverse power potentials u(r)=ε(σ/r)^{n}, where parameter n measures softness of the potential. We find that for 1/nfcc) structure while for 1/n≥0.154 the body-centred-cubic (bcc) structure is preferred. The bcc structure transforms into the fcc structure upon increasing the density. The calculated phase diagram is in good agreement with the one found from molecular simulations.

Critical Assessment of Time-Dependent Density Functional Theory for Excited States of Open-Shell Systems: II. Doublet-Quartet Transitions.

Science.gov (United States)

Li, Zhendong; Liu, Wenjian

2016-06-14

Compared with closed-shell systems, open-shell systems place three additional challenges to time-dependent density functional theory (TD-DFT) for electronically excited states: (a) the spin-contamination problem is a serious issue; (b) the exchange-correlation (XC) kernel may be numerically instable; and (c) the single-determinant description of open-shell ground states readily becomes energetically instable. Confined to flip-up single excitations, the spin-contamination problem can largely be avoided by using the spin-flip TD-DFT (SF-TD-DFT) formalism, provided that a noncollinear XC kernel is employed. As for the numerical instabilities associated with such a kernel, only an ad hoc scheme has been proposed so far, viz., the ALDA0 kernel, which amounts to setting the divergent components (arising from density gradients and kinetic energy density) simply to zero. The ground-state instability problem can effectively be avoided by introducing the Tamm-Dancoff approximation (TDA) to TD-DFT. Therefore, on a general basis, the SF-TDA/ALDA0 Ansatz is so far the only promising means within the TD-DFT framework for flip-up single excitations of open-shell systems. To assess systematically the performance of SF-TDA/ALDA0, in total 61 low-lying quartet excited states of the benchmark set of 11 small radicals [J. Chem. Theory Comput. 2016, 12, 238] are investigated with various XC functionals. Taking the MRCISD+Q (multireference configuration interaction with singles and doubles plus the Davidson correction) results as benchmark, it is found that the mean absolute errors of SF-TDA/ALDA0 with the SAOP (statistical averaging of model orbital potentials), global hybrid, and range-separated hybrid functionals are in the range of 0.2-0.4 eV. This is in line not only with the typical accuracy of TD-DFT for singlet and triplet excited states of closed-shell systems but also with the gross accuracy of spin-adapted TD-DFT for spin-conserving excited states of open-shell systems.

Pressure and surface tension of solid-liquid interface using Tara zona density functional theory

International Nuclear Information System (INIS)

Moradi, M.; Kavosh Tehrani, M.

2001-01-01

The weighted density functional theory proposed by Tara zona is applied to study the solid-liquid interface. In the last two decades the weighted density functional became a useful tool to consider the properties of inhomogeneous liquids. In this theory, the role of the size of molecules or the particles of which the matter is composed, was found to be important. In this research we study a hard sphere fluid beside a hard wall. For this study the liquid is an inhomogeneous system. We use the definition of the direct correlation function as a second derivative of free energy with respect to the density. We use this definition and the definition of the weighting function, then we minimize the grand potential with respect to the density to get the Euler Lagrange equation and we obtain an integral equation to find the inhomogeneous density profile. The obtained density profile as a function of the distance from the wall, for different bulk density is pitted in three dimensions. We also calculate the pressure and compare it with the Carnahan-Starling results, and finally we obtained the surface tension at liquid-solid interface and compared it with the results of Monte Carlo simulation

Strontium, nickel, cadmium, and lead substitution into calcite, studied by density functional theory

DEFF Research Database (Denmark)

Andersson, Martin Peter; Sakuma, Hiroshi; Stipp, Susan Louise Svane

2014-01-01

We have used density functional theory to predict the ion exchange energies for divalent cations Ni(2+), Sr(2+), Cd(2+), and Pb(2+) into a calcite {10.4} surface in equilibrium with water. Exchange energies were calculated for substitution into the topmost surface layer, at the mineral-fluid inte......We have used density functional theory to predict the ion exchange energies for divalent cations Ni(2+), Sr(2+), Cd(2+), and Pb(2+) into a calcite {10.4} surface in equilibrium with water. Exchange energies were calculated for substitution into the topmost surface layer, at the mineral...

Calculation of exchange coupling constants in triply-bridged dinuclear Cu(II) compounds based on spin-flip constricted variational density functional theory.

Science.gov (United States)

Seidu, Issaka; Zhekova, Hristina R; Seth, Michael; Ziegler, Tom

2012-03-08

The performance of the second-order spin-flip constricted variational density functional theory (SF-CV(2)-DFT) for the calculation of the exchange coupling constant (J) is assessed by application to a series of triply bridged Cu(II) dinuclear complexes. A comparison of the J values based on SF-CV(2)-DFT with those obtained by the broken symmetry (BS) DFT method and experiment is provided. It is demonstrated that our methodology constitutes a viable alternative to the BS-DFT method. The strong dependence of the calculated exchange coupling constants on the applied functionals is demonstrated. Both SF-CV(2)-DFT and BS-DFT affords the best agreement with experiment for hybrid functionals.

Analysis of the photophysical properties of zearalenone using density functional theory

Science.gov (United States)

The intrinsic photophysical properties of the resorcylic acid moiety of zearalenone offer a convenient label free method to determine zearalenone levels in contaminated agricultural products. Density functional theory and steady-state fluorescence methods were applied to investigate the role of stru...

Stochastic density functional theory at finite temperatures

Science.gov (United States)

Cytter, Yael; Rabani, Eran; Neuhauser, Daniel; Baer, Roi

2018-03-01

Simulations in the warm dense matter regime using finite temperature Kohn-Sham density functional theory (FT-KS-DFT), while frequently used, are computationally expensive due to the partial occupation of a very large number of high-energy KS eigenstates which are obtained from subspace diagonalization. We have developed a stochastic method for applying FT-KS-DFT, that overcomes the bottleneck of calculating the occupied KS orbitals by directly obtaining the density from the KS Hamiltonian. The proposed algorithm scales as O (" close=")N3T3)">N T-1 and is compared with the high-temperature limit scaling O density approximation (LDA); we demonstrate its efficiency, statistical errors, and bias in the estimation of the free energy per electron for a diamond structure silicon. The bias is small compared to the fluctuations and is independent of system size. In addition to calculating the free energy itself, one can also use the method to calculate its derivatives and obtain the equations of state.

Magnetic behavior study of samarium nitride using density functional theory

Science.gov (United States)

Som, Narayan N.; Mankad, Venu H.; Dabhi, Shweta D.; Patel, Anjali; Jha, Prafulla K.

2018-02-01

In this work, the state-of-art density functional theory is employed to study the structural, electronic and magnetic properties of samarium nitride (SmN). We have performed calculation for both ferromagnetic and antiferromagnetic states in rock-salt phase. The calculated results of optimized lattice parameter and magnetic moment agree well with the available experimental and theoretical values. From energy band diagram and electronic density of states, we observe a half-metallic behaviour in FM phase of rock salt SmN in while metallicity in AFM I and AFM III phases. We present and discuss our current understanding of the possible half-metallicity together with the magnetic ordering in SmN. The calculated phonon dispersion curves shows dynamical stability of the considered structures. The phonon density of states and Eliashberg functional have also been analysed to understand the superconductivity in SmN.

Regarding the use and misuse of retinal protonated Schiff base photochemistry as a test case for time-dependent density-functional theory

Energy Technology Data Exchange (ETDEWEB)

Valsson, Omar [Department of Chemistry and Applied Biosciences, ETH Zurich and Facoltà di Informatica, Instituto di Scienze Computationali, Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6900 Lugano (Switzerland); Filippi, Claudia, E-mail: c.filippi@utwente.nl [MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede (Netherlands); Casida, Mark E., E-mail: mark.casida@ujf-grenoble.fr [Laboratoire de Chimie Théorique, Département de Chimie Moléculaire (DCM), Institut de Chimie Moléculaire de Grenoble (ICMG), Université Joseph Fourier, Grenoble I, F-3801 Grenoble (France)

2015-04-14

The excited-state relaxation of retinal protonated Schiff bases (PSBs) is an important test case for biological applications of time-dependent (TD) density-functional theory (DFT). While well-known shortcomings of approximate TD-DFT might seem discouraging for application to PSB relaxation, progress continues to be made in the development of new functionals and of criteria allowing problematic excitations to be identified within the framework of TD-DFT itself. Furthermore, experimental and theoretical ab initio advances have recently lead to a revised understanding of retinal PSB photochemistry, calling for a reappraisal of the performance of TD-DFT in describing this prototypical photoactive system. Here, we re-investigate the performance of functionals in (TD-)DFT calculations in light of these new benchmark results, which we extend to larger PSB models. We focus on the ability of the functionals to describe primarily the early skeletal relaxation of the chromophore and investigate how far along the out-of-plane pathways these functionals are able to describe the subsequent rotation around formal single and double bonds. Conventional global hybrid and range-separated hybrid functionals are investigated as the presence of Hartree-Fock exchange reduces problems with charge-transfer excitations as determined by the Peach-Benfield-Helgaker-Tozer Λ criterion and by comparison with multi-reference perturbation theory results. While we confirm that most functionals cannot render the complex photobehavior of the retinal PSB, do we also observe that LC-BLYP gives the best description of the initial part of the photoreaction.

Pressure-Dependent Electronic and Transport Properties of Bulk Platinum Oxide by Density Functional Theory

Science.gov (United States)

Kansara, Shivam; Gupta, Sanjeev K.; Sonvane, Yogesh; Nekrasov, Kirill A.; Kichigina, Natalia V.

2018-02-01

The structural, electronic, and vibrational properties of bulk platinum oxide (PtO) at compressive pressures in the interval from 0 GPa to 35 GPa are investigated using the density functional theory. The calculated electronic band structure of PtO shows poor metallicity at very low density of states on the Fermi level. However, the hybrid pseudopotential calculation yielded 0.78 eV and 1.30 eV direct band and indirect gap, respectively. Importantly, our results predict that PtO has a direct band gap within the framework of HSE06, and it prefers equally zero magnetic order at different pressures. In the Raman spectra, peaks are slightly shifted towards higher frequency with the decrease in pressure. We have also calculated the thermoelectric properties, namely the electronic thermal conductivity and electrical conductivity, with respect to temperature and thermodynamic properties such as entropy, specific heat at constant volume, enthalpy and Gibbs free energy with respect to pressure. The result shows that PtO is a promising candidate for use as a catalyst, in sensors, as a photo-cathode in water electrolysis, for thermal decomposition of inorganic salt and fuel cells.

Density functional theory for field emission from carbon nano-structures.

Science.gov (United States)

Li, Zhibing

2015-12-01

Electron field emission is understood as a quantum mechanical many-body problem in which an electronic quasi-particle of the emitter is converted into an electron in vacuum. Fundamental concepts of field emission, such as the field enhancement factor, work-function, edge barrier and emission current density, will be investigated, using carbon nanotubes and graphene as examples. A multi-scale algorithm basing on density functional theory is introduced. We will argue that such a first principle approach is necessary and appropriate for field emission of nano-structures, not only for a more accurate quantitative description, but, more importantly, for deeper insight into field emission. Copyright © 2015 The Author. Published by Elsevier B.V. All rights reserved.

The X3LYP extended density functional for accurate descriptions of nonbond interactions, spin states, and thermochemical properties

OpenAIRE

Xu, Xin; Goddard, William A., III

2004-01-01

We derive the form for an exact exchange energy density for a density decaying with Gaussian-like behavior at long range. Based on this, we develop the X3LYP (extended hybrid functional combined with Lee-Yang-Parr correlation functional) extended functional for density functional theory to significantly improve the accuracy for hydrogen-bonded and van der Waals complexes while also improving the accuracy in heats of formation, ionization potentials, electron affinities, and total atomic energ...

Time-dependent current-density functional theory for generalized open quantum systems.

Science.gov (United States)

Yuen-Zhou, Joel; Rodríguez-Rosario, César; Aspuru-Guzik, Alán

2009-06-14

In this article, we prove the one-to-one correspondence between vector potentials and particle and current densities in the context of master equations with arbitrary memory kernels, therefore extending time-dependent current-density functional theory (TD-CDFT) to the domain of generalized many-body open quantum systems (OQS). We also analyse the issue of A-representability for the Kohn-Sham (KS) scheme proposed by D'Agosta and Di Ventra for Markovian OQS [Phys. Rev. Lett. 2007, 98, 226403] and discuss its domain of validity. We suggest ways to expand their scheme, but also propose a novel KS scheme where the auxiliary system is both closed and non-interacting. This scheme is tested numerically with a model system, and several considerations for the future development of functionals are indicated. Our results formalize the possibility of practising TD-CDFT in OQS, hence expanding the applicability of the theory to non-Hamiltonian evolutions.

Density functional theory studies of screw dislocation core structures in bcc metals

DEFF Research Database (Denmark)

Frederiksen, Søren Lund; Jacobsen, Karsten Wedel

2003-01-01

The core structures of (I 11) screw dislocations in bee metals are studied using density functional theory in the local-density approximation. For Mo and Fe, direct calculations of the core structures show the cores to be symmetric with respect to 180degrees rotations around an axis perpendicular...... to symmetric core structures for all the studied metals....

Alternative separation of exchange and correlation energies in range-separated density-functional perturbation theory

DEFF Research Database (Denmark)

Cornaton, Y.; Stoyanova, A.; Jensen, Hans Jørgen Aagaard

2013-01-01

of the noninteracting Kohn-Sham one. When second-order corrections to the density are neglected, the energy expression reduces to a range-separated double-hybrid (RSDH) type of functional, RSDHf, where "f" stands for "full-range integrals" as the regular full-range interaction appears explicitly in the energy...

Introduction to Density Functional Theory: Calculations by Hand on the Helium Atom

Science.gov (United States)

Baseden, Kyle A.; Tye, Jesse W.

2014-01-01

Density functional theory (DFT) is a type of electronic structure calculation that has rapidly gained popularity. In this article, we provide a step-by-step demonstration of a DFT calculation by hand on the helium atom using Slater's X-Alpha exchange functional on a single Gaussian-type orbital to represent the atomic wave function. This DFT…

Natural excitation orbitals from linear response theories : Time-dependent density functional theory, time-dependent Hartree-Fock, and time-dependent natural orbital functional theory

NARCIS (Netherlands)

Van Meer, R.; Gritsenko, O. V.; Baerends, E. J.

2017-01-01

Straightforward interpretation of excitations is possible if they can be described as simple single orbital-to-orbital (or double, etc.) transitions. In linear response time-dependent density functional theory (LR-TDDFT), the (ground state) Kohn-Sham orbitals prove to be such an orbital basis. In

Time-dependent density-functional theory in the projector augmented-wave method

DEFF Research Database (Denmark)

Walter, Michael; Häkkinen, Hannu; Lehtovaara, Lauri

2008-01-01

We present the implementation of the time-dependent density-functional theory both in linear-response and in time-propagation formalisms using the projector augmented-wave method in real-space grids. The two technically very different methods are compared in the linear-response regime where we...

On the accuracy of van der Waals inclusive density-functional theory exchange-correlation functionals for ice at ambient and high pressures.

Science.gov (United States)

Santra, Biswajit; Klimes, Jirí; Tkatchenko, Alexandre; Alfè, Dario; Slater, Ben; Michaelides, Angelos; Car, Roberto; Scheffler, Matthias

2013-10-21

Density-functional theory (DFT) has been widely used to study water and ice for at least 20 years. However, the reliability of different DFT exchange-correlation (xc) functionals for water remains a matter of considerable debate. This is particularly true in light of the recent development of DFT based methods that account for van der Waals (vdW) dispersion forces. Here, we report a detailed study with several xc functionals (semi-local, hybrid, and vdW inclusive approaches) on ice Ih and six proton ordered phases of ice. Consistent with our previous study [B. Santra, J. Klimeš, D. Alfè, A. Tkatchenko, B. Slater, A. Michaelides, R. Car, and M. Scheffler, Phys. Rev. Lett. 107, 185701 (2011)] which showed that vdW forces become increasingly important at high pressures, we find here that all vdW inclusive methods considered improve the relative energies and transition pressures of the high-pressure ice phases compared to those obtained with semi-local or hybrid xc functionals. However, we also find that significant discrepancies between experiment and the vdW inclusive approaches remain in the cohesive properties of the various phases, causing certain phases to be absent from the phase diagram. Therefore, room for improvement in the description of water at ambient and high pressures remains and we suggest that because of the stern test the high pressure ice phases pose they should be used in future benchmark studies of simulation methods for water.

Exact density functional and wave function embedding schemes based on orbital localization

International Nuclear Information System (INIS)

Hégely, Bence; Nagy, Péter R.; Kállay, Mihály; Ferenczy, György G.

2016-01-01

Exact schemes for the embedding of density functional theory (DFT) and wave function theory (WFT) methods into lower-level DFT or WFT approaches are introduced utilizing orbital localization. First, a simple modification of the projector-based embedding scheme of Manby and co-workers [J. Chem. Phys. 140, 18A507 (2014)] is proposed. We also use localized orbitals to partition the system, but instead of augmenting the Fock operator with a somewhat arbitrary level-shift projector we solve the Huzinaga-equation, which strictly enforces the Pauli exclusion principle. Second, the embedding of WFT methods in local correlation approaches is studied. Since the latter methods split up the system into local domains, very simple embedding theories can be defined if the domains of the active subsystem and the environment are treated at a different level. The considered embedding schemes are benchmarked for reaction energies and compared to quantum mechanics (QM)/molecular mechanics (MM) and vacuum embedding. We conclude that for DFT-in-DFT embedding, the Huzinaga-equation-based scheme is more efficient than the other approaches, but QM/MM or even simple vacuum embedding is still competitive in particular cases. Concerning the embedding of wave function methods, the clear winner is the embedding of WFT into low-level local correlation approaches, and WFT-in-DFT embedding can only be more advantageous if a non-hybrid density functional is employed.

Exact density functional and wave function embedding schemes based on orbital localization

Science.gov (United States)

Hégely, Bence; Nagy, Péter R.; Ferenczy, György G.; Kállay, Mihály

2016-08-01

Exact schemes for the embedding of density functional theory (DFT) and wave function theory (WFT) methods into lower-level DFT or WFT approaches are introduced utilizing orbital localization. First, a simple modification of the projector-based embedding scheme of Manby and co-workers [J. Chem. Phys. 140, 18A507 (2014)] is proposed. We also use localized orbitals to partition the system, but instead of augmenting the Fock operator with a somewhat arbitrary level-shift projector we solve the Huzinaga-equation, which strictly enforces the Pauli exclusion principle. Second, the embedding of WFT methods in local correlation approaches is studied. Since the latter methods split up the system into local domains, very simple embedding theories can be defined if the domains of the active subsystem and the environment are treated at a different level. The considered embedding schemes are benchmarked for reaction energies and compared to quantum mechanics (QM)/molecular mechanics (MM) and vacuum embedding. We conclude that for DFT-in-DFT embedding, the Huzinaga-equation-based scheme is more efficient than the other approaches, but QM/MM or even simple vacuum embedding is still competitive in particular cases. Concerning the embedding of wave function methods, the clear winner is the embedding of WFT into low-level local correlation approaches, and WFT-in-DFT embedding can only be more advantageous if a non-hybrid density functional is employed.

Exact density functional and wave function embedding schemes based on orbital localization

Energy Technology Data Exchange (ETDEWEB)

Hégely, Bence; Nagy, Péter R.; Kállay, Mihály, E-mail: kallay@mail.bme.hu [MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest (Hungary); Ferenczy, György G. [Medicinal Chemistry Research Group, Research Centre for Natural Sciences, Hungarian Academy of Sciences, Magyar tudósok körútja 2, H-1117 Budapest (Hungary); Department of Biophysics and Radiation Biology, Semmelweis University, Tűzoltó u. 37-47, H-1094 Budapest (Hungary)

2016-08-14

Exact schemes for the embedding of density functional theory (DFT) and wave function theory (WFT) methods into lower-level DFT or WFT approaches are introduced utilizing orbital localization. First, a simple modification of the projector-based embedding scheme of Manby and co-workers [J. Chem. Phys. 140, 18A507 (2014)] is proposed. We also use localized orbitals to partition the system, but instead of augmenting the Fock operator with a somewhat arbitrary level-shift projector we solve the Huzinaga-equation, which strictly enforces the Pauli exclusion principle. Second, the embedding of WFT methods in local correlation approaches is studied. Since the latter methods split up the system into local domains, very simple embedding theories can be defined if the domains of the active subsystem and the environment are treated at a different level. The considered embedding schemes are benchmarked for reaction energies and compared to quantum mechanics (QM)/molecular mechanics (MM) and vacuum embedding. We conclude that for DFT-in-DFT embedding, the Huzinaga-equation-based scheme is more efficient than the other approaches, but QM/MM or even simple vacuum embedding is still competitive in particular cases. Concerning the embedding of wave function methods, the clear winner is the embedding of WFT into low-level local correlation approaches, and WFT-in-DFT embedding can only be more advantageous if a non-hybrid density functional is employed.

Classical density functional theory & simulations on a coarse-grained model of aromatic ionic liquids.

Science.gov (United States)

Turesson, Martin; Szparaga, Ryan; Ma, Ke; Woodward, Clifford E; Forsman, Jan

2014-05-14

A new classical density functional approach is developed to accurately treat a coarse-grained model of room temperature aromatic ionic liquids. Our major innovation is the introduction of charge-charge correlations, which are treated in a simple phenomenological way. We test this theory on a generic coarse-grained model for aromatic RTILs with oligomeric forms for both cations and anions, approximating 1-alkyl-3-methyl imidazoliums and BF₄⁻, respectively. We find that predictions by the new density functional theory for fluid structures at charged surfaces are very accurate, as compared with molecular dynamics simulations, across a range of surface charge densities and lengths of the alkyl chain. Predictions of interactions between charged surfaces are also presented.

A comprehensive study of g-factors, elastic, structural and electronic properties of III-V semiconductors using hybrid-density functional theory

Science.gov (United States)

Bastos, Carlos M. O.; Sabino, Fernando P.; Sipahi, Guilherme M.; Da Silva, Juarez L. F.

2018-02-01

Despite the large number of theoretical III-V semiconductor studies reported every year, our atomistic understanding is still limited. The limitations of the theoretical approaches to yield accurate structural and electronic properties on an equal footing, is due to the unphysical self-interaction problem that mainly affects the band gap and spin-orbit splitting (SOC) in semiconductors and, in particular, III-V systems with similar magnitude of the band gap and SOC. In this work, we report a consistent study of the structural and electronic properties of the III-V semiconductors by using the screening hybrid-density functional theory framework, by fitting the α parameters for 12 different III-V compounds, namely, AlN, AlP, AlAs, AlSb, GaN, GaP, GaAs, GaSb, InN, InP, InAs, and InSb, to minimize the deviation between the theoretical and experimental values of the band gap and SOC. Structural relaxation effects were also included. Except for AlP, whose α = 0.127, we obtained α values that ranged from 0.209 to 0.343, which deviate by less than 0.1 from the universal value of 0.25. Our results for the lattice parameter and elastic constants indicate that the fitting of α does not affect those structural parameters when compared with the HSE06 functional, where α = 0.25. Our analysis of the band structure based on the k ṡ p method shows that the effective masses are in agreement with the experimental values, which can be attributed to the simultaneous fitting of the band gap and SOC. Also, we estimate the values of g-factors, extracted directly from the band structure, which are close to experimental results, which indicate that the obtained band structure produced a realistic set of k ṡ p parameters.

Hartree and Exchange in Ensemble Density Functional Theory: Avoiding the Nonuniqueness Disaster.

Science.gov (United States)

Gould, Tim; Pittalis, Stefano

2017-12-15

Ensemble density functional theory is a promising method for the efficient and accurate calculation of excitations of quantum systems, at least if useful functionals can be developed to broaden its domain of practical applicability. Here, we introduce a guaranteed single-valued "Hartree-exchange" ensemble density functional, E_{Hx}[n], in terms of the right derivative of the universal ensemble density functional with respect to the coupling constant at vanishing interaction. We show that E_{Hx}[n] is straightforwardly expressible using block eigenvalues of a simple matrix [Eq. (14)]. Specialized expressions for E_{Hx}[n] from the literature, including those involving superpositions of Slater determinants, can now be regarded as originating from the unifying picture presented here. We thus establish a clear and practical description for Hartree and exchange in ensemble systems.

Analytic calculations of hyper-Raman spectra from density functional theory hyperpolarizability gradients

Energy Technology Data Exchange (ETDEWEB)

Ringholm, Magnus; Ruud, Kenneth [Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø – The Arctic University of Norway, 9037 Tromsø (Norway); Bast, Radovan [Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, AlbaNova University Center, S-10691 Stockholm (Sweden); PDC Center for High Performance Computing, Royal Institute of Technology, S-10044 Stockholm (Sweden); Oggioni, Luca [Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø – The Arctic University of Norway, 9037 Tromsø (Norway); Department of Physics G. Occhialini, University of Milano Bicocca, Piazza della scienza 3, 20126 Milan (Italy); Ekström, Ulf [Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, 0315 Oslo (Norway)

2014-10-07

We present the first analytic calculations of the geometrical gradients of the first hyperpolarizability tensors at the density-functional theory (DFT) level. We use the analytically calculated hyperpolarizability gradients to explore the importance of electron correlation effects, as described by DFT, on hyper-Raman spectra. In particular, we calculate the hyper-Raman spectra of the all-trans and 11-cis isomers of retinal at the Hartree-Fock (HF) and density-functional levels of theory, also allowing us to explore the sensitivity of the hyper-Raman spectra on the geometrical characteristics of these structurally related molecules. We show that the HF results, using B3LYP-calculated vibrational frequencies and force fields, reproduce the experimental data for all-trans-retinal well, and that electron correlation effects are of minor importance for the hyper-Raman intensities.

A density functional study of backbone structures of polydiacetylene: destabilization of butatriene structure

International Nuclear Information System (INIS)

Katagiri, Hideki; Shimoi, Yukihiro; Abe, Shuji

2004-01-01

Backbone structures of polydiacetylene are studied with first-principles electronic structure method using plane-waves within generalized gradient approximation (GGA) of density functional theory. In spin-restricted calculations a coarse k-point sampling gives a potential energy curve with two local minima corresponding to acetylene and butatriene structures. However, the potential barrier between the two structures rapidly decreases with increasing number of k-points, which results in destabilization of the butatriene structure. Spin polarization effects also destabilize the butatriene structure, inducing atom-centered spin-density-wave state. These potential energies were compared with those obtained by Hartree-Fock, density functional within local density approximation (LDA) and GGA, and hybrid density functional methods using a gaussian basis set. The comparison shows that the density functional methods within LDA and GGA favor the destabilization of the butatriene structure in contrast to the Hartree-Fock method

Time-dependent density functional theory for open quantum systems with unitary propagation.

Science.gov (United States)

Yuen-Zhou, Joel; Tempel, David G; Rodríguez-Rosario, César A; Aspuru-Guzik, Alán

2010-01-29

We extend the Runge-Gross theorem for a very general class of open quantum systems under weak assumptions about the nature of the bath and its coupling to the system. We show that for Kohn-Sham (KS) time-dependent density functional theory, it is possible to rigorously include the effects of the environment within a bath functional in the KS potential. A Markovian bath functional inspired by the theory of nonlinear Schrödinger equations is suggested, which can be readily implemented in currently existing real-time codes. Finally, calculations on a helium model system are presented.

Ionization Potentials of Chemical Warfare Agents and Related Compounds Determined with Density Functional Theory

National Research Council Canada - National Science Library

Wright, J

2000-01-01

...) agents at contaminated sites. Reported herein are theoretical ionization potentials for CW agents and their related compounds calculated using density functional theory at the B3LYP/6-311+G(2d,p) level of theory...

Multicomponent Time-Dependent Density Functional Theory: Proton and Electron Excitation Energies.

Science.gov (United States)

Yang, Yang; Culpitt, Tanner; Hammes-Schiffer, Sharon

2018-04-05

The quantum mechanical treatment of both electrons and protons in the calculation of excited state properties is critical for describing nonadiabatic processes such as photoinduced proton-coupled electron transfer. Multicomponent density functional theory enables the consistent quantum mechanical treatment of more than one type of particle and has been implemented previously for studying ground state molecular properties within the nuclear-electronic orbital (NEO) framework, where all electrons and specified protons are treated quantum mechanically. To enable the study of excited state molecular properties, herein the linear response multicomponent time-dependent density functional theory (TDDFT) is derived and implemented within the NEO framework. Initial applications to FHF - and HCN illustrate that NEO-TDDFT provides accurate proton and electron excitation energies within a single calculation. As its computational cost is similar to that of conventional electronic TDDFT, the NEO-TDDFT approach is promising for diverse applications, particularly nonadiabatic proton transfer reactions, which may exhibit mixed electron-proton vibronic excitations.

Magnetic exchange couplings from constrained density functional theory: an efficient approach utilizing analytic derivatives.

Science.gov (United States)

Phillips, Jordan J; Peralta, Juan E

2011-11-14

We introduce a method for evaluating magnetic exchange couplings based on the constrained density functional theory (C-DFT) approach of Rudra, Wu, and Van Voorhis [J. Chem. Phys. 124, 024103 (2006)]. Our method shares the same physical principles as C-DFT but makes use of the fact that the electronic energy changes quadratically and bilinearly with respect to the constraints in the range of interest. This allows us to use coupled perturbed Kohn-Sham spin density functional theory to determine approximately the corrections to the energy of the different spin configurations and construct a priori the relevant energy-landscapes obtained by constrained spin density functional theory. We assess this methodology in a set of binuclear transition-metal complexes and show that it reproduces very closely the results of C-DFT. This demonstrates a proof-of-concept for this method as a potential tool for studying a number of other molecular phenomena. Additionally, routes to improving upon the limitations of this method are discussed. © 2011 American Institute of Physics

Electronic structure of bulk AnO{sub 2} (An = U, Np, Pu) and water adsorption on the (111) and (110) surfaces of UO{sub 2} and PuO{sub 2} from hybrid density functional theory within the periodic electrostatic embedded cluster method

Energy Technology Data Exchange (ETDEWEB)

Wellington, Joseph P.W. [Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ (United Kingdom); Kerridge, Andrew [Department of Chemistry, Lancaster University, Bailrigg, Lancaster LA1 4YP (United Kingdom); Austin, Jonathan [National Nuclear Laboratory, Chadwick House, Warrington Road, Birchwood Park, Warrington WA3 6AE (United Kingdom); Kaltsoyannis, Nikolas, E-mail: nikolas.kaltsoyannis@manchester.ac.uk [School of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL (United Kingdom)

2016-12-15

Generalised gradient approximation (PBE) and hybrid (PBE0) density functional theory (DFT) within the periodic electrostatic embedded cluster method have been used to study AnO{sub 2} bulk and surfaces (An = U, Np, Pu). The electronic structure has been investigated by examining the projected density of states (PDOS). While PBE incorrectly predicts these systems to be metallic, PBE0 finds them to be insulators, with the composition of the valence and conduction levels agreeing well with experiment. Molecular and dissociative water adsorption on the (111) and (110) surfaces of UO{sub 2} and PuO{sub 2} has been investigated, with that on the (110) surface being stronger than on the (111). Similar energies are found for molecular and dissociative adsorption on the (111) surfaces, while on the (110) there is a clear preference for dissociative adsorption. Adsorption energies and geometries on the (111) surface of UO{sub 2} are in good agreement with recent periodic DFT studies using the GGA+U approach, and our data for dissociative adsorption on the (110) surface of PuO{sub 2} match experiment rather well, especially when dispersion corrections are included. - Graphical abstract: The electronic structures of AnO{sub 2} (An = U, Np, Pu) are studied computationally with hybrid density functional theory, and the geometries and energetics of water adsorption on the low index surfaces are presented.

Gaussian-3 theory using density functional geometries and zero-point energies

International Nuclear Information System (INIS)

Baboul, A.G.; Curtiss, L.A.; Redfern, P.C.; Raghavachari, K.

1999-01-01

A variation of Gaussian-3 (G3) theory is presented in which the geometries and zero-point energies are obtained from B3LYP density functional theory [B3LYP/6-31G(d)] instead of geometries from second-order perturbation theory [MP2(FU)/6-31G(d)] and zero-point energies from Hartree - Fock theory [HF/6-31G(d)]. This variation, referred to as G3//B3LYP, is assessed on 299 energies (enthalpies of formation, ionization potentials, electron affinities, proton affinities) from the G2/97 test set [J. Chem. Phys. 109, 42 (1998)]. The G3//B3LYP average absolute deviation from experiment for the 299 energies is 0.99 kcal/mol compared to 1.01 kcal/mol for G3 theory. Generally, the results from the two methods are similar, with some exceptions. G3//B3LYP theory gives significantly improved results for several cases for which MP2 theory is deficient for optimized geometries, such as CN and O 2 + . However, G3//B3LYP does poorly for ionization potentials that involve a Jahn - Teller distortion in the cation (CH 4 + , BF 3 + , BCl 3 + ) because of the B3LYP/6-31G(d) geometries. The G3(MP2) method is also modified to use B3LYP/6-31G(d) geometries and zero-point energies. This variation, referred to as G3(MP2)//B3LYP, has an average absolute deviation of 1.25 kcal/mol compared to 1.30 kcal/mol for G3(MP2) theory. Thus, use of density functional geometries and zero-point energies in G3 and G3(MP2) theories is a useful alternative to MP2 geometries and HF zero-point energies. copyright 1999 American Institute of Physics

Recent progress in orbital-free density functional theory (recent advances in computational chemistry)

CERN Document Server

Wesolowski, Tomasz A

2013-01-01

This is a comprehensive overview of state-of-the-art computational methods based on orbital-free formulation of density functional theory completed by the most recent developments concerning the exact properties, approximations, and interpretations of the relevant quantities in density functional theory. The book is a compilation of contributions stemming from a series of workshops which had been taking place since 2002. It not only chronicles many of the latest developments but also summarises some of the more significant ones. The chapters are mainly reviews of sub-domains but also include original research. Readership: Graduate students, academics and researchers in computational chemistry. Atomic & molecular physicists, theoretical physicists, theoretical chemists, physical chemists and chemical physicists.

Charge Transfer Enhancement in the D-π-A Type Porphyrin Dyes: A Density Functional Theory (DFT and Time-Dependent Density Functional Theory (TD-DFT Study

Directory of Open Access Journals (Sweden)

Guo-Jun Kang

2016-11-01

Full Text Available The electronic geometries and optical properties of two D-π-A type zinc porphyrin dyes (NCH3-YD2 and TPhe-YD were systematically investigated by density functional theory (DFT and time-dependent density functional theory (TD-DFT to reveal the origin of significantly altered charge transfer enhancement by changing the electron donor of the famous porphyrin-based sensitizer YD2-o-C8. The molecular geometries and photophysical properties of dyes before and after binding to the TiO2 cluster were fully investigated. From the analyses of natural bond orbital (NBO, extended charge decomposition analysis (ECDA, and electron density variations (Δρ between the excited state and ground state, it was found that the introduction of N(CH32 and 1,1,2-triphenylethene groups enhanced the intramolecular charge-transfer (ICT character compared to YD2-o-C8. The absorption wavelength and transition possess character were significantly influenced by N(CH32 and 1,1,2-triphenylethene groups. NCH3-YD2 with N(CH32 groups in the donor part is an effective way to improve the interactions between the dyes and TiO2 surface, light having efficiency (LHE, and free energy change (ΔGinject, which is expected to be an efficient dye for use in dye-sensitized solar cells (DSSCs.

Conjugation-promoted reaction of open-cage fullerene: A density functional theory study

KAUST Repository

Guo, Yong; Yan, Jingjing; Khashab, Niveen M.

2012-01-01

Density functional theory calculations are performed to study the addition mechanism of e-rich moieties such as triethyl phosphite to a carbonyl group on the rim of a fullerene orifice. Three possible reaction channels have been investigated

van der Waals forces in density functional theory: Perturbational long-range electron-interaction corrections

International Nuclear Information System (INIS)

Angyan, Janos G.; Gerber, Iann C.; Savin, Andreas; Toulouse, Julien

2005-01-01

Long-range exchange and correlation effects, responsible for the failure of currently used approximate density functionals in describing van der Waals forces, are taken into account explicitly after a separation of the electron-electron interaction in the Hamiltonian into short- and long-range components. We propose a 'range-separated hybrid' functional based on a local density approximation for the short-range exchange-correlation energy, combined with a long-range exact exchange energy. Long-range correlation effects are added by a second-order perturbational treatment. The resulting scheme is general and is particularly well adapted to describe van der Waals complexes, such as rare gas dimers

Discontinuities of Green functions in field theory at finite temperature and density

International Nuclear Information System (INIS)

Kobes, R.L.; Semenoff, G.W.

1985-01-01

We derive systematic rules for calculating the imaginary parts of Minkowski space Green functions in quantum field theory at finite temperature and density. Self-energy corrections are used as an example of the application of these rules. (orig.)

Range-separated density-functional theory for molecular excitation energies

International Nuclear Information System (INIS)

Rebolini, E.

2014-01-01

Linear-response time-dependent density-functional theory (TDDFT) is nowadays a method of choice to compute molecular excitation energies. However, within the usual adiabatic semi-local approximations, it is not able to describe properly Rydberg, charge-transfer or multiple excitations. Range separation of the electronic interaction allows one to mix rigorously density-functional methods at short range and wave function or Green's function methods at long range. When applied to the exchange functional, it already corrects most of these deficiencies but multiple excitations remain absent as they need a frequency-dependent kernel. In this thesis, the effects of range separation are first assessed on the excitation energies of a partially-interacting system in an analytic and numerical study in order to provide guidelines for future developments of range-separated methods for excitation energy calculations. It is then applied on the exchange and correlation TDDFT kernels in a single-determinant approximation in which the long-range part of the correlation kernel vanishes. A long-range frequency-dependent second-order correlation kernel is then derived from the Bethe-Salpeter equation and added perturbatively to the range-separated TDDFT kernel in order to take into account the effects of double excitations. (author)

Buckled graphene: A model study based on density functional theory

KAUST Repository

Khan, Yasser

2010-09-01

We make use of ab initio calculations within density functional theory to investigate the influence of buckling on the electronic structure of single layer graphene. Our systematic study addresses a wide range of bond length and bond angle variations in order to obtain insights into the energy scale associated with the formation of ripples in a graphene sheet. © 2010 Elsevier B.V. All rights reserved.

Buckled graphene: A model study based on density functional theory

KAUST Repository

Khan, Yasser; Mukaddam, Mohsin Ahmed; Schwingenschlö gl, Udo

2010-01-01

We make use of ab initio calculations within density functional theory to investigate the influence of buckling on the electronic structure of single layer graphene. Our systematic study addresses a wide range of bond length and bond angle variations in order to obtain insights into the energy scale associated with the formation of ripples in a graphene sheet. © 2010 Elsevier B.V. All rights reserved.

ONETEP: linear-scaling density-functional theory with plane-waves

International Nuclear Information System (INIS)

Haynes, P D; Mostof, A A; Skylaris, C-K; Payne, M C

2006-01-01

This paper provides a general overview of the methodology implemented in onetep (Order-N Electronic Total Energy Package), a parallel density-functional theory code for largescale first-principles quantum-mechanical calculations. The distinctive features of onetep are linear-scaling in both computational effort and resources, obtained by making well-controlled approximations which enable simulations to be performed with plane-wave accuracy. Titanium dioxide clusters of increasing size designed to mimic surfaces are studied to demonstrate the accuracy and scaling of onetep

Hybrid functional calculation of electronic and phonon structure of BaSnO3

International Nuclear Information System (INIS)

Kim, Bog G.; Jo, J.Y.; Cheong, S.W.

2013-01-01

Barium stannate, BaSnO 3 (BSO), with a cubic perovskite structure, has been highlighted as a promising host material for the next generation transparent oxide electrodes. This study examined theoretically the electronic structure and phonon structure of BSO using hybrid density functional theory based on the HSE06 functional. The electronic structure results of BSO were corrected by extending the phonon calculations based on the hybrid density functional. The fundamental thermal properties were also predicted based on a hybrid functional calculation. Overall, a detailed understanding of the electronic structure, phonon modes and phonon dispersion of BSO will provide a theoretical starting-point for engineering applications of this material. - Graphical Abstract: (a) Crystal structure of BaSnO 3 . The center ball is Ba and small (red) ball on edge is oxygen and SnO 6 octahedrons are plotted as polyhedron. (b) Electronic band structure along the high symmetry point in the Brillouin zone using the HSE06 hybrid functional. (c) The phonon dispersion curve calculated using the HSE06 hybrid functional (d) Zone center lowest energy F 1u phonon mode. Highlights: ► We report the full hybrid functional calculation of not only the electronic structure but also the phonon structure for BaSnO 3 . ► The band gap calculation of HSE06 revealed an indirect gap with 2.48 eV. ► The effective mass at the conduction band minimum and valence band maximum was calculated. ► In addition, the phonon structure of BSO was calculated using the HSE06 functional. ► Finally, the heat capacity was calculated and compared with the recent experimental result.

Anatase TiO{sub 2} nanowires functionalized by organic sensitizers for solar cells: A screened Coulomb hybrid density functional study

Energy Technology Data Exchange (ETDEWEB)

Ünal, Hatice; Mete, Ersen, E-mail: emete@balikesir.edu.tr [Deparment of Physics, Balikesir University, Balikesir 10145 (Turkey); Gunceler, Deniz [Deparment of Physics, Cornell University, Ithaca, New York 14853 (United States); Gülseren, Oğuz [Department of Physics, Bilkent University, Ankara 06800 (Turkey); Ellialtioğlu, Şinasi [Basic Sciences, TED University, Ankara 06420 (Turkey)

2015-11-21

The adsorption of two different organic molecules cyanidin glucoside (C{sub 21}O{sub 11}H{sub 20}) and TA-St-CA on anatase (101) and (001) nanowires has been investigated using the standard and the range separated hybrid density functional theory calculations. The electronic structures and optical spectra of resulting dye–nanowire combined systems show distinct features for these types of photochromophores. The lowest unoccupied molecular orbital of the natural dye cyanidin glucoside is located below the conduction band of the semiconductor while, in the case of TA-St-CA, it resonates with the states inside the conduction band. The wide-bandgap anatase nanowires can be functionalized for solar cells through electron-hole generation and subsequent charge injection by these dye sensitizers. The intermolecular charge transfer character of Donor-π-Acceptor type dye TA-St-CA is substantially modified by its adsorption on TiO{sub 2} surfaces. Cyanidin glucoside exhibits relatively stronger anchoring on the nanowires through its hydroxyl groups. The atomic structures of dye–nanowire systems re-optimized with the inclusion of nonlinear solvation effects showed that the binding strengths of both dyes remain moderate even in ionic solutions.

A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions.

Science.gov (United States)

Goerigk, Lars; Hansen, Andreas; Bauer, Christoph; Ehrlich, Stephan; Najibi, Asim; Grimme, Stefan

2017-12-13

We present the GMTKN55 benchmark database for general main group thermochemistry, kinetics and noncovalent interactions. Compared to its popular predecessor GMTKN30 [Goerigk and Grimme J. Chem. Theory Comput., 2011, 7, 291], it allows assessment across a larger variety of chemical problems-with 13 new benchmark sets being presented for the first time-and it also provides reference values of significantly higher quality for most sets. GMTKN55 comprises 1505 relative energies based on 2462 single-point calculations and it is accessible to the user community via a dedicated website. Herein, we demonstrate the importance of better reference values, and we re-emphasise the need for London-dispersion corrections in density functional theory (DFT) treatments of thermochemical problems, including Minnesota methods. We assessed 217 variations of dispersion-corrected and -uncorrected density functional approximations, and carried out a detailed analysis of 83 of them to identify robust and reliable approaches. Double-hybrid functionals are the most reliable approaches for thermochemistry and noncovalent interactions, and they should be used whenever technically feasible. These are, in particular, DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), and B2GPPLYP-D3(BJ). The best hybrids are ωB97X-V, M052X-D3(0), and ωB97X-D3, but we also recommend PW6B95-D3(BJ) as the best conventional global hybrid. At the meta-generalised-gradient (meta-GGA) level, the SCAN-D3(BJ) method can be recommended. Other meta-GGAs are outperformed by the GGA functionals revPBE-D3(BJ), B97-D3(BJ), and OLYP-D3(BJ). We note that many popular methods, such as B3LYP, are not part of our recommendations. In fact, with our results we hope to inspire a change in the user community's perception of common DFT methods. We also encourage method developers to use GMTKN55 for cross-validation studies of new methodologies.

Critique of the foundations of time-dependent density-functional theory

International Nuclear Information System (INIS)

Schirmer, J.; Dreuw, A.

2007-01-01

The general expectation that, in principle, the time-dependent density-functional theory (TDDFT) is an exact formulation of the time evolution of an interacting N-electron system is critically reexamined. It is demonstrated that the previous TDDFT foundation, resting on four theorems by Runge and Gross (RG) [Phys. Rev. Lett. 52, 997 (1984)], is invalid because undefined phase factors corrupt the RG action integral functionals. Our finding confirms much of a previous analysis by van Leeuwen [Int. J. Mod. Phys. B 15, 1969 (2001)]. To analyze the RG theorems and other aspects of TDDFT, an utmost simplification of the Kohn-Sham (KS) concept has been introduced, in which the ground-state density is obtained from a single KS equation for one spatial (spinless) orbital. The time-dependent (TD) form of this radical Kohn-Sham (rKS) scheme, which has the same validity status as the ordinary KS version, has proved to be a valuable tool for analysis. The rKS concept is used to clarify also the alternative nonvariational formulation of TD KS theory. We argue that it is just a formal theory, allowing one to reproduce but not predict the time development of the exact density of the interacting N-electron system. Besides the issue of the formal exactness of TDDFT, it is shown that both the static and time-dependent KS linear response equations neglect the particle-particle (p-p) and hole-hole (h-h) matrix elements of the perturbing operator. For a local (multiplicative) operator this does not lead to a loss of information due to a remarkable general property of local operators. Accordingly, no logical inconsistency arises with respect to DFT, because DFT requires any external potential to be local. For a general nonlocal operator the error resulting from the neglected matrix elements is of second order in the electronic repulsion

PREFACE: Classical density functional theory methods in soft and hard matter Classical density functional theory methods in soft and hard matter

Science.gov (United States)

Haataja, Mikko; Gránásy, László; Löwen, Hartmut

2010-08-01

Herein we provide a brief summary of the background, events and results/outcome of the CECAM workshop 'Classical density functional theory methods in soft and hard matter held in Lausanne between October 21 and October 23 2009, which brought together two largely separately working communities, both of whom employ classical density functional techniques: the soft-matter community and the theoretical materials science community with interests in phase transformations and evolving microstructures in engineering materials. After outlining the motivation for the workshop, we first provide a brief overview of the articles submitted by the invited speakers for this special issue of Journal of Physics: Condensed Matter, followed by a collection of outstanding problems identified and discussed during the workshop. 1. Introduction Classical density functional theory (DFT) is a theoretical framework, which has been extensively employed in the past to study inhomogeneous complex fluids (CF) [1-4] and freezing transitions for simple fluids, amongst other things. Furthermore, classical DFT has been extended to include dynamics of the density field, thereby opening a new avenue to study phase transformation kinetics in colloidal systems via dynamical DFT (DDFT) [5]. While DDFT is highly accurate, the computations are numerically rather demanding, and cannot easily access the mesoscopic temporal and spatial scales where diffusional instabilities lead to complex solidification morphologies. Adaptation of more efficient numerical methods would extend the domain of DDFT towards this regime of particular interest to materials scientists. In recent years, DFT has re-emerged in the form of the so-called 'phase-field crystal' (PFC) method for solid-state systems [6, 7], and it has been successfully employed to study a broad variety of interesting materials phenomena in both atomic and colloidal systems, including elastic and plastic deformations, grain growth, thin film growth, solid

Optical excitation and electron relaxation dynamics at semiconductor surfaces: a combined approach of density functional and density matrix theory applied to the silicon (001) surface

Energy Technology Data Exchange (ETDEWEB)

Buecking, N

2007-11-05

In this work a new theoretical formalism is introduced in order to simulate the phononinduced relaxation of a non-equilibrium distribution to equilibrium at a semiconductor surface numerically. The non-equilibrium distribution is effected by an optical excitation. The approach in this thesis is to link two conventional, but approved methods to a new, more global description: while semiconductor surfaces can be investigated accurately by density-functional theory, the dynamical processes in semiconductor heterostructures are successfully described by density matrix theory. In this work, the parameters for density-matrix theory are determined from the results of density-functional calculations. This work is organized in two parts. In Part I, the general fundamentals of the theory are elaborated, covering the fundamentals of canonical quantizations as well as the theory of density-functional and density-matrix theory in 2{sup nd} order Born approximation. While the formalism of density functional theory for structure investigation has been established for a long time and many different codes exist, the requirements for density matrix formalism concerning the geometry and the number of implemented bands exceed the usual possibilities of the existing code in this field. A special attention is therefore attributed to the development of extensions to existing formulations of this theory, where geometrical and fundamental symmetries of the structure and the equations are used. In Part II, the newly developed formalism is applied to a silicon (001)surface in a 2 x 1 reconstruction. As first step, density-functional calculations using the LDA functional are completed, from which the Kohn-Sham-wave functions and eigenvalues are used to calculate interaction matrix elements for the electron-phonon-coupling an the optical excitation. These matrix elements are determined for the optical transitions from valence to conduction bands and for electron-phonon processes inside the

Metal-insulator transition at the LaAlO3/SrTiO3 interface revisited: A hybrid functional study

KAUST Repository

Cossu, Fabrizio; Eyert, V.; Schwingenschlö gl, Udo

2013-01-01

We investigate the electronic properties of the LaAlO3/SrTiO3 interface using density functional theory. In contrast to previous studies, which relied on (semi-)local functionals and the GGA+U method, we here use a recently developed hybrid

Describing a Strongly Correlated Model System with Density Functional Theory.

Science.gov (United States)

Kong, Jing; Proynov, Emil; Yu, Jianguo; Pachter, Ruth

2017-07-06

The linear chain of hydrogen atoms, a basic prototype for the transition from a metal to Mott insulator, is studied with a recent density functional theory model functional for nondynamic and strong correlation. The computed cohesive energy curve for the transition agrees well with accurate literature results. The variation of the electronic structure in this transition is characterized with a density functional descriptor that yields the atomic population of effectively localized electrons. These new methods are also applied to the study of the Peierls dimerization of the stretched even-spaced Mott insulator to a chain of H 2 molecules, a different insulator. The transitions among the two insulating states and the metallic state of the hydrogen chain system are depicted in a semiquantitative phase diagram. Overall, we demonstrate the capability of studying strongly correlated materials with a mean-field model at the fundamental level, in contrast to the general pessimistic view on such a feasibility.

Density functional theory and beyond-opportunities for quantum methods in materials modeling semiconductor technology

International Nuclear Information System (INIS)

Shankar, Sadasivan; Simka, Harsono; Haverty, Michael

2008-01-01

In the semiconductor industry, the use of new materials has been increasing with the advent of nanotechnology. As critical dimensions decrease, and the number of materials increases, the interactions between heterogeneous materials themselves and processing increase in complexity. Traditionally, applications of ab initio techniques are confined to electronic structure and band gap calculations of bulk materials, which are then used in coarse-grained models such as mesoscopic and continuum models. Density functional theory is the most widely used ab initio technique that was successfully extended to several applications. This paper illustrates applications of density functional theory to semiconductor processes and proposes further opportunities for use of such techniques in process development

Locality of correlation in density functional theory

Energy Technology Data Exchange (ETDEWEB)

Burke, Kieron [Department of Chemistry, University of California, Irvine, California 92697 (United States); Cancio, Antonio [Department of Physics and Astronomy, Ball State University, Muncie, Indiana 47306 (United States); Gould, Tim [Qld Micro- and Nanotechnology Centre, Griffith University, Nathan, Qld 4111 (Australia); Pittalis, Stefano [CNR-Istituto di Nanoscienze, Via Campi 213A, I-41125 Modena (Italy)

2016-08-07

The Hohenberg-Kohn density functional was long ago shown to reduce to the Thomas-Fermi (TF) approximation in the non-relativistic semiclassical (or large-Z) limit for all matter, i.e., the kinetic energy becomes local. Exchange also becomes local in this limit. Numerical data on the correlation energy of atoms support the conjecture that this is also true for correlation, but much less relevant to atoms. We illustrate how expansions around a large particle number are equivalent to local density approximations and their strong relevance to density functional approximations. Analyzing highly accurate atomic correlation energies, we show that E{sub C} → −A{sub C} ZlnZ + B{sub C}Z as Z → ∞, where Z is the atomic number, A{sub C} is known, and we estimate B{sub C} to be about 37 mhartree. The local density approximation yields A{sub C} exactly, but a very incorrect value for B{sub C}, showing that the local approximation is less relevant for the correlation alone. This limit is a benchmark for the non-empirical construction of density functional approximations. We conjecture that, beyond atoms, the leading correction to the local density approximation in the large-Z limit generally takes this form, but with B{sub C} a functional of the TF density for the system. The implications for the construction of approximate density functionals are discussed.

Performance of Density Functional Theory Procedures for the Calculation of Proton-Exchange Barriers: Unusual Behavior of M06-Type Functionals.

Science.gov (United States)

Chan, Bun; Gilbert, Andrew T B; Gill, Peter M W; Radom, Leo

2014-09-09

We have examined the performance of a variety of density functional theory procedures for the calculation of complexation energies and proton-exchange barriers, with a focus on the Minnesota-class of functionals that are generally highly robust and generally show good accuracy. A curious observation is that M05-type and M06-type methods show an atypical decrease in calculated barriers with increasing proportion of Hartree-Fock exchange. To obtain a clearer picture of the performance of the underlying components of M05-type and M06-type functionals, we have investigated the combination of MPW-type and PBE-type exchange and B95-type and PBE-type correlation procedures. We find that, for the extensive E3 test set, the general performance of the various hybrid-DFT procedures improves in the following order: PBE1-B95 → PBE1-PBE → MPW1-PBE → PW6-B95. As M05-type and M06-type procedures are related to PBE1-B95, it would be of interest to formulate and examine the general performance of an alternative Minnesota DFT method related to PW6-B95.

Surface effects on mean inner potentials studied using density functional theory

Energy Technology Data Exchange (ETDEWEB)

Pennington, Robert S., E-mail: robert.pennington@uni-ulm.de [Institute for Experimental Physics, Ulm University, Albert-Einstein-Allee 11, 89081 Ulm (Germany); Boothroyd, Chris B.; Dunin-Borkowski, Rafal E. [Ernst Ruska-Centre and Peter Grüneberg Institute, Forschungzentrum Jülich, 52425 Jülich (Germany)

2015-12-15

Quantitative materials characterization using electron holography frequently requires knowledge of the mean inner potential, but reported experimental mean inner potential measurements can vary widely. Using density functional theory, we have simulated the mean inner potential for materials with a range of different surface conditions and geometries. We use both “thin-film” and “nanowire” specimen geometries. We consider clean bulk-terminated surfaces with different facets and surface reconstructions using atom positions from both structural optimization and experimental data and we also consider surfaces both with and without adsorbates. We find that the mean inner potential is surface-dependent, with the strongest dependency on surface adsorbates. We discuss the outlook and perspective for future mean inner potential measurements. - Highlights: • Density functional theory (DFT) is used to simulate mean inner potentials (MIP). • Applications for MIP electron holography measurements are considered. • MIPs are found to be surface-dependent, for thin-film and nanowire geometries. • The DFT simulation precision is extensively tested for multiple materials. • Surface adsorbates can create a strong positive or negative effect.

Density-density functionals and effective potentials in many-body electronic structure calculations

International Nuclear Information System (INIS)

Reboredo, Fernando A.; Kent, Paul R.

2008-01-01

We demonstrate the existence of different density-density functionals designed to retain selected properties of the many-body ground state in a non-interacting solution starting from the standard density functional theory ground state. We focus on diffusion quantum Monte Carlo applications that require trial wave functions with optimal Fermion nodes. The theory is extensible and can be used to understand current practices in several electronic structure methods within a generalized density functional framework. The theory justifies and stimulates the search of optimal empirical density functionals and effective potentials for accurate calculations of the properties of real materials, but also cautions on the limits of their applicability. The concepts are tested and validated with a near-analytic model.

Wavelet-based linear-response time-dependent density-functional theory

Science.gov (United States)

Natarajan, Bhaarathi; Genovese, Luigi; Casida, Mark E.; Deutsch, Thierry; Burchak, Olga N.; Philouze, Christian; Balakirev, Maxim Y.

2012-06-01

Linear-response time-dependent (TD) density-functional theory (DFT) has been implemented in the pseudopotential wavelet-based electronic structure program BIGDFT and results are compared against those obtained with the all-electron Gaussian-type orbital program DEMON2K for the calculation of electronic absorption spectra of N2 using the TD local density approximation (LDA). The two programs give comparable excitation energies and absorption spectra once suitably extensive basis sets are used. Convergence of LDA density orbitals and orbital energies to the basis-set limit is significantly faster for BIGDFT than for DEMON2K. However the number of virtual orbitals used in TD-DFT calculations is a parameter in BIGDFT, while all virtual orbitals are included in TD-DFT calculations in DEMON2K. As a reality check, we report the X-ray crystal structure and the measured and calculated absorption spectrum (excitation energies and oscillator strengths) of the small organic molecule N-cyclohexyl-2-(4-methoxyphenyl)imidazo[1, 2-a]pyridin-3-amine.

External field as the functional of inhomogeneous density and the density matrix functional approach

NARCIS (Netherlands)

Bobrov, V.B.; Trigger, S.A.; Vlasov, Y.P.

2012-01-01

Based on the Hohenberg-Kohn lemma and the hypotheses of the density functional existence for the external-field potential, it is shown that the strict result of the density functional theory is the equation of the external-field potential as the density functional. This result leads to the

The X3LYP extended density functional for accurate descriptions of nonbond interactions, spin states, and thermochemical properties

Science.gov (United States)

Xu, Xin; Goddard, William A.

2004-01-01

We derive the form for an exact exchange energy density for a density decaying with Gaussian-like behavior at long range. Based on this, we develop the X3LYP (extended hybrid functional combined with Lee–Yang–Parr correlation functional) extended functional for density functional theory to significantly improve the accuracy for hydrogen-bonded and van der Waals complexes while also improving the accuracy in heats of formation, ionization potentials, electron affinities, and total atomic energies [over the most popular and accurate method, B3LYP (Becke three-parameter hybrid functional combined with Lee–Yang–Parr correlation functional)]. X3LYP also leads to a good description of dipole moments, polarizabilities, and accurate excitation energies from s to d orbitals for transition metal atoms and ions. We suggest that X3LYP will be useful for predicting ligand binding in proteins and DNA. PMID:14981235

Restoring the consistency with the contact density theorem of a classical density functional theory of ions at a planar electrical double layer.

Science.gov (United States)

Gillespie, Dirk

2014-11-01

Classical density functional theory (DFT) of fluids is a fast and efficient theory to compute the structure of the electrical double layer in the primitive model of ions where ions are modeled as charged, hard spheres in a background dielectric. While the hard-core repulsive component of this ion-ion interaction can be accurately computed using well-established DFTs, the electrostatic component is less accurate. Moreover, many electrostatic functionals fail to satisfy a basic theorem, the contact density theorem, that relates the bulk pressure, surface charge, and ion densities at their distances of closest approach for ions in equilibrium at a smooth, hard, planar wall. One popular electrostatic functional that fails to satisfy the contact density theorem is a perturbation approach developed by Kierlik and Rosinberg [Phys. Rev. A 44, 5025 (1991)PLRAAN1050-294710.1103/PhysRevA.44.5025] and Rosenfeld [J. Chem. Phys. 98, 8126 (1993)JCPSA60021-960610.1063/1.464569], where the full free-energy functional is Taylor-expanded around a bulk (homogeneous) reference fluid. Here, it is shown that this functional fails to satisfy the contact density theorem because it also fails to satisfy the known low-density limit. When the functional is corrected to satisfy this limit, a corrected bulk pressure is derived and it is shown that with this pressure both the contact density theorem and the Gibbs adsorption theorem are satisfied.

Hydrogen plasmas beyond density-functional theory: dynamic correlations and the onset of localization

International Nuclear Information System (INIS)

Perrot, F.; Dharma-Wardana, M.W.C.

1984-01-01

The density-functional theory (DFT) equations - previously considered in their application to the study of a system of ions and electrons in thermodynamic equilibrium at arbitrary temperatures and pressure - are reviewed with attention given to extending their validity in obtaining the one-electron excitation spectrum. The DFT model developed here provides structure factors and Kohn-Sham eigenstates which are then used to calculate the self-energy of the one-electron Green function, thus transcending the local-density approximations and the well-known limitations of DFT, especially with regard to the excitation spectrum. The one-particle formalism used makes contact with the multiple-scattering theories of disordered materials, liquid metals, etc., and is a necessary first step to a future calculation of two-particle propagators and related properties. 28 references

Use of density functional theory in drug metabolism studies

DEFF Research Database (Denmark)

Rydberg, Patrik; Jørgensen, Flemming Steen; Olsen, Lars

2014-01-01

INTRODUCTION: The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool...... isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account...

Self-interaction corrections in density functional theory

International Nuclear Information System (INIS)

Tsuneda, Takao; Hirao, Kimihiko

2014-01-01

Self-interaction corrections for Kohn-Sham density functional theory are reviewed for their physical meanings, formulations, and applications. The self-interaction corrections get rid of the self-interaction error, which is the sum of the Coulomb and exchange self-interactions that remains because of the use of an approximate exchange functional. The most frequently used self-interaction correction is the Perdew-Zunger correction. However, this correction leads to instabilities in the electronic state calculations of molecules. To avoid these instabilities, several self-interaction corrections have been developed on the basis of the characteristic behaviors of self-interacting electrons, which have no two-electron interactions. These include the von Weizsäcker kinetic energy and long-range (far-from-nucleus) asymptotic correction. Applications of self-interaction corrections have shown that the self-interaction error has a serious effect on the states of core electrons, but it has a smaller than expected effect on valence electrons. This finding is supported by the fact that the distribution of self-interacting electrons indicates that they are near atomic nuclei rather than in chemical bonds

Current Issues in Finite-T Density-Functional Theory and Warm-Correlated Matter †

Directory of Open Access Journals (Sweden)

M. W. C. Dharma-wardana

2016-03-01

Full Text Available Finite-temperature density functional theory (DFT has become of topical interest, partly due to the increasing ability to create novel states of warm-correlated matter (WCM.Warm-dense matter (WDM, ultra-fast matter (UFM, and high-energy density matter (HEDM may all be regarded as subclasses of WCM. Strong electron-electron, ion-ion and electron-ion correlation effects and partial degeneracies are found in these systems where the electron temperature Te is comparable to the electron Fermi energy EF. Thus, many electrons are in continuum states which are partially occupied. The ion subsystem may be solid, liquid or plasma, with many states of ionization with ionic charge Zj. Quasi-equilibria with the ion temperature Ti ≠ Te are common. The ion subsystem in WCM can no longer be treated as a passive “external potential”, as is customary in T = 0 DFT dominated by solid-state theory or quantum chemistry. Many basic questions arise in trying to implement DFT for WCM. Hohenberg-Kohn-Mermin theory can be adapted for treating these systems if suitable finite-T exchange-correlation (XC functionals can be constructed. They are functionals of both the one-body electron density ne and the one-body ion densities ρj. Here, j counts many species of nuclei or charge states. A method of approximately but accurately mapping the quantum electrons to a classical Coulomb gas enables one to treat electron-ion systems entirely classically at any temperature and arbitrary spin polarization, using exchange-correlation effects calculated in situ, directly from the pair-distribution functions. This eliminates the need for any XC-functionals. This classical map has been used to calculate the equation of state of WDM systems, and construct a finite-T XC functional that is found to be in close agreement with recent quantum path-integral simulation data. In this review, current developments and concerns in finite-T DFT, especially in the context of non-relativistic warm

Formation energies of rutile metal dioxides using density functional theory

DEFF Research Database (Denmark)

Martinez, Jose Ignacio; Hansen, Heine Anton; Rossmeisl, Jan

2009-01-01

We apply standard density functional theory at the generalized gradient approximation (GGA) level to study the stability of rutile metal oxides. It is well known that standard GGA exchange and correlation in some cases is not sufficient to address reduction and oxidation reactions. Especially...... and due to a more accurate description of exchange for this particular GGA functional compared to PBE. Furthermore, we would expect the self-interaction problem to be largest for the most localized d orbitals; that means the late 3d metals and since Co, Fe, Ni, and Cu do not form rutile oxides...

Density-functional theory based on the electron distribution on the energy coordinate

Science.gov (United States)

Takahashi, Hideaki

2018-03-01

We developed an electronic density functional theory utilizing a novel electron distribution n(ɛ) as a basic variable to compute ground state energy of a system. n(ɛ) is obtained by projecting the electron density n({\\boldsymbol{r}}) defined on the space coordinate {\\boldsymbol{r}} onto the energy coordinate ɛ specified with the external potential {\\upsilon }ext}({\\boldsymbol{r}}) of interest. It was demonstrated that the Kohn-Sham equation can also be formulated with the exchange-correlation functional E xc[n(ɛ)] that employs the density n(ɛ) as an argument. It turned out an exchange functional proposed in our preliminary development suffices to describe properly the potential energies of several types of chemical bonds with comparable accuracies to the corresponding functional based on local density approximation. As a remarkable feature of the distribution n(ɛ) it inherently involves the spatially non-local information of the exchange hole at the bond dissociation limit in contrast to conventional approximate functionals. By taking advantage of this property we also developed a prototype of the static correlation functional E sc including no empirical parameters, which showed marked improvements in describing the dissociations of covalent bonds in {{{H}}}2,{{{C}}}2{{{H}}}4 and {CH}}4 molecules.

Reduced density-matrix functional theory: Correlation and spectroscopy.

Science.gov (United States)

Di Sabatino, S; Berger, J A; Reining, L; Romaniello, P

2015-07-14

In this work, we explore the performance of approximations to electron correlation in reduced density-matrix functional theory (RDMFT) and of approximations to the observables calculated within this theory. Our analysis focuses on the calculation of total energies, occupation numbers, removal/addition energies, and spectral functions. We use the exactly solvable Hubbard dimer at 1/4 and 1/2 fillings as test systems. This allows us to analyze the underlying physics and to elucidate the origin of the observed trends. For comparison, we also report the results of the GW approximation, where the self-energy functional is approximated, but no further hypothesis is made concerning the approximations of the observables. In particular, we focus on the atomic limit, where the two sites of the dimer are pulled apart and electrons localize on either site with equal probability, unless a small perturbation is present: this is the regime of strong electron correlation. In this limit, using the Hubbard dimer at 1/2 filling with or without a spin-symmetry-broken ground state allows us to explore how degeneracies and spin-symmetry breaking are treated in RDMFT. We find that, within the used approximations, neither in RDMFT nor in GW, the signature of strong correlation is present, when looking at the removal/addition energies and spectral function from the spin-singlet ground state, whereas both give the exact result for the spin-symmetry broken case. Moreover, we show how the spectroscopic properties change from one spin structure to the other.

Hybrid Fluid/Kinetic Modeling Of Magnetized High Energy Density Plasmas

Science.gov (United States)

Hansen, David; Held, Eric; King, Jacob; Stoltz, Peter; Masti, Robert; Srinivasan, Bhuvana

2017-10-01

MHD modeling with an equation of state (EOS) of the Rayleigh-Taylor (RT) instabily in Z indicates that it is seeded by the electro-thermal instability. Large thermodynamic drives associated with gradients at the interface between the liner and the coronal regions distort distribution functions and likely lead to non-local transport effects in a plasma which varies from weakly to strongly coupled. In this work, we discuss using effective potential theory along with a Chapman-Ensksog-like (CEL) formalism to develop hybrid fluid/kinetic modeling capabilities for these plasmas. Effective potential theory addresses the role of Coulomb collisions on transport across coupling regimes and the CEL approach bridges the gap between full-blow kinetic simulations and the EOS tables, which only depend locally on density and temperature. Quantitative results on the Spitzer problem across coupling coupling regimes will be presented as a first step. DOE Grant No. DE-SC0016525.

No need for external orthogonality in subsystem density-functional theory.

Science.gov (United States)

Unsleber, Jan P; Neugebauer, Johannes; Jacob, Christoph R

2016-08-03

Recent reports on the necessity of using externally orthogonal orbitals in subsystem density-functional theory (SDFT) [Annu. Rep. Comput. Chem., 8, 2012, 53; J. Phys. Chem. A, 118, 2014, 9182] are re-investigated. We show that in the basis-set limit, supermolecular Kohn-Sham-DFT (KS-DFT) densities can exactly be represented as a sum of subsystem densities, even if the subsystem orbitals are not externally orthogonal. This is illustrated using both an analytical example and in basis-set free numerical calculations for an atomic test case. We further show that even with finite basis sets, SDFT calculations using accurate reconstructed potentials can closely approach the supermolecular KS-DFT density, and that the deviations between SDFT and KS-DFT decrease as the basis-set limit is approached. Our results demonstrate that formally, there is no need to enforce external orthogonality in SDFT, even though this might be a useful strategy when developing projection-based DFT embedding schemes.

Isobaric-Isothermal Molecular Dynamics Utilizing Density Functional Theory: An Assessment of the Structure and Density of Water at Near-Ambient Conditions

International Nuclear Information System (INIS)

Schmidt, J.; VandeVondele, J.; Kuo, I.W.; Sebastiani, D.; Siepmann, J.I.; Hutter, J.; Mundy, C.J.

2009-01-01

We present herein a comprehensive density functional theory study toward assessing the accuracy of two popular gradient-corrected exchange correlation functionals on the structure and density of liquid water at near ambient conditions in the isobaric-isothermal ensemble. Our results indicate that both PBE and BLYP functionals under predict the density and over structure the liquid. Adding the dispersion correction due to Grimme(1, 2) improves the predicted densities for both BLYP and PBE in a significant manner. Moreover, the addition of the dispersion correction for BLYP yields an oxygen-oxygen radial distribution function in excellent agreement with experiment. Thus, we conclude that one can obtain a very satisfactory model for water using BLYP and a correction for dispersion.

Locating the rate-limiting step for the interaction of hydrogen with Mg(0001) using density-functional theory calculations and rate theory

DEFF Research Database (Denmark)

Vegge, Tejs

2004-01-01

The dissociation of molecular hydrogen on a Mgs0001d surface and the subsequent diffusion of atomic hydrogen into the magnesium substrate is investigated using Density Functional Theory (DFT) calculations and rate theory. The minimum energy path and corresponding transition states are located usi...... to be rate-limiting for the ab- and desorption of hydrogen, respectively. Zero-point energy contributions are found to be substantial for the diffusion of atomic hydrogen, but classical rates are still found to be within an order of magnitude at room temperature.......The dissociation of molecular hydrogen on a Mgs0001d surface and the subsequent diffusion of atomic hydrogen into the magnesium substrate is investigated using Density Functional Theory (DFT) calculations and rate theory. The minimum energy path and corresponding transition states are located using...

Dynamic density functional theory of solid tumor growth: Preliminary models

Directory of Open Access Journals (Sweden)

Arnaud Chauviere

2012-03-01

Full Text Available Cancer is a disease that can be seen as a complex system whose dynamics and growth result from nonlinear processes coupled across wide ranges of spatio-temporal scales. The current mathematical modeling literature addresses issues at various scales but the development of theoretical methodologies capable of bridging gaps across scales needs further study. We present a new theoretical framework based on Dynamic Density Functional Theory (DDFT extended, for the first time, to the dynamics of living tissues by accounting for cell density correlations, different cell types, phenotypes and cell birth/death processes, in order to provide a biophysically consistent description of processes across the scales. We present an application of this approach to tumor growth.

Descriptions of carbon isotopes within the energy density functional theory

International Nuclear Information System (INIS)

Ismail, Atef; Cheong, Lee Yen; Yahya, Noorhana; Tammam, M.

2014-01-01

Within the energy density functional (EDF) theory, the structure properties of Carbon isotopes are systematically studied. The shell model calculations are done for both even-A and odd-A nuclei, to study the structure of rich-neutron Carbon isotopes. The EDF theory indicates the single-neutron halo structures in 15 C, 17 C and 19 C, and the two-neutron halo structures in 16 C and 22 C nuclei. It is also found that close to the neutron drip-line, there exist amazing increase in the neutron radii and decrease on the binding energies BE, which are tightly related with the blocking effect and correspondingly the blocking effect plays a significant role in the shell model configurations

Magnetism in Sc-doped ZnO with zinc vacancies: A hybrid density functional and GGA + U approaches

KAUST Repository

Kanoun, Mohammed; Goumri-Said, Souraya; Schwingenschlö gl, Udo; Manchon, Aurelien

2012-01-01

We investigate the zinc vacancy effects on the electronic structures and magnetic properties of Sc-doped ZnO, by performing first-principles calculations within both GGA + U and Heyd-Scuseria-Ernzerhof hybrid functional methods. We find that Sc impurities stabilize considerably Zn vacancies. The electronic and magnetic analysis shows a half metallic ferromagnetic character with a total magnetic moment of 2.01 μ B. The magnetism mainly stems from the O 2p states around the Zn vacancies. Calculations with the hybrid density functional agree with the GGA + U results but give an accurate description of the electronic structure for pure ZnO and Sc-doped ZnO with Zn vacancies. © 2012 Elsevier B.V. All rights reserved.

Magnetism in Sc-doped ZnO with zinc vacancies: A hybrid density functional and GGA + U approaches

KAUST Repository

Kanoun, Mohammed

2012-04-01

We investigate the zinc vacancy effects on the electronic structures and magnetic properties of Sc-doped ZnO, by performing first-principles calculations within both GGA + U and Heyd-Scuseria-Ernzerhof hybrid functional methods. We find that Sc impurities stabilize considerably Zn vacancies. The electronic and magnetic analysis shows a half metallic ferromagnetic character with a total magnetic moment of 2.01 μ B. The magnetism mainly stems from the O 2p states around the Zn vacancies. Calculations with the hybrid density functional agree with the GGA + U results but give an accurate description of the electronic structure for pure ZnO and Sc-doped ZnO with Zn vacancies. © 2012 Elsevier B.V. All rights reserved.

Central depression in nucleonic densities: Trend analysis in the nuclear density functional theory approach

Science.gov (United States)

Schuetrumpf, B.; Nazarewicz, W.; Reinhard, P.-G.

2017-08-01

Background: The central depression of nucleonic density, i.e., a reduction of density in the nuclear interior, has been attributed to many factors. For instance, bubble structures in superheavy nuclei are believed to be due to the electrostatic repulsion. In light nuclei, the mechanism behind the density reduction in the interior has been discussed in terms of shell effects associated with occupations of s orbits. Purpose: The main objective of this work is to reveal mechanisms behind the formation of central depression in nucleonic densities in light and heavy nuclei. To this end, we introduce several measures of the internal nucleonic density. Through the statistical analysis, we study the information content of these measures with respect to nuclear matter properties. Method: We apply nuclear density functional theory with Skyrme functionals. Using the statistical tools of linear least square regression, we inspect correlations between various measures of central depression and model parameters, including nuclear matter properties. We study bivariate correlations with selected quantities as well as multiple correlations with groups of parameters. Detailed correlation analysis is carried out for 34Si for which a bubble structure has been reported recently, 48Ca, and N =82 , 126, and 184 isotonic chains. Results: We show that the central depression in medium-mass nuclei is very sensitive to shell effects, whereas for superheavy systems it is firmly driven by the electrostatic repulsion. An appreciable semibubble structure in proton density is predicted for 294Og, which is currently the heaviest nucleus known experimentally. Conclusion: Our correlation analysis reveals that the central density indicators in nuclei below 208Pb carry little information on parameters of nuclear matter; they are predominantly driven by shell structure. On the other hand, in the superheavy nuclei there exists a clear relationship between the central nucleonic density and symmetry energy.

Relativistic density functional for nuclear structure

CERN Document Server

2016-01-01

This book aims to provide a detailed introduction to the state-of-the-art covariant density functional theory, which follows the Lorentz invariance from the very beginning and is able to describe nuclear many-body quantum systems microscopically and self-consistently. Covariant density functional theory was introduced in nuclear physics in the 1970s and has since been developed and used to describe the diversity of nuclear properties and phenomena with great success. In order to provide an advanced and updated textbook of covariant density functional theory for graduate students and nuclear physics researchers, this book summarizes the enormous amount of material that has accumulated in the field of covariant density functional theory over the last few decades as well as the latest developments in this area. Moreover, the book contains enough details for readers to follow the formalism and theoretical results, and provides exhaustive references to explore the research literature.

On the universality of the long-/short-range separation in multiconfigurational density-functional theory

DEFF Research Database (Denmark)

Fromager, Emmanuel; Toulouse, Julien; Jensen, Hans Jørgen Aagaard

2007-01-01

In many cases, the dynamic correlation can be calculated quite accurately and at a fairly low computational cost in Kohn-Sham density-functional theory (KS-DFT), using current standard approximate functionals. However, in general, KS-DFT does not treat static correlation effects (near degeneracy...

Molecular dynamics simulations of the amino acid-ZnO (10-10) interface: A comparison between density functional theory and density functional tight binding results

Energy Technology Data Exchange (ETDEWEB)

Holthaus, Svea große; Köppen, Susan, E-mail: koeppen@hmi.uni-bremen.de; Frauenheim, Thomas; Ciacchi, Lucio Colombi [Bremen Centre for Computational Materials Science, University of Bremen, 28359 Bremen (Germany)

2014-06-21

We investigate the adsorption behavior of four different amino acids (glutamine, glutamate, serine, cysteine) on the zinc oxide (101{sup ¯}0) surface, comparing the geometry and energy associated with a number of different adsorption configurations. In doing this, we highlight the benefits and limits of using density-functional tight-binding (DFTB) with respect to standard density functional theory (DFT). The DFTB method is found to reliably reproduce the DFT adsorption geometries. Analysis of the adsorption configurations emphasizes the fundamental role of the first hydration layer in mediating the interactions between the amino acids and the surface. Direct surface-molecule bonds are found to form predominantly via the carboxylate groups of the studied amino acids. No surface-mediated chemical reactions are observed, with the notable exception of a proton transfer from the thiol group of cysteine to a hydroxyl group of the surface hydration layer. The adsorption energies are found to be dominated both by the formation of direct or indirect surface-molecule hydrogen bonds, but also by the rearrangement of the hydrogen-bond network in surface proximity in a non-intuitive way. Energetic comparisons between DFTB and DFT are made difficult on one side by the long time necessary to achieve convergence of potential energy values in MD simulations and on the other side by the necessity of including higher-order corrections to DFTB to obtain a good description of the hydrogen bond energetics. Overall, our results suggest that DFTB is a good reference method to set the correct chemical states and the initial geometries of hybrid biomolecule/ZnO systems to be simulated with non-reactive force fields.

Molecular dynamics simulations of the amino acid-ZnO (10-10) interface: a comparison between density functional theory and density functional tight binding results.

Science.gov (United States)

grosse Holthaus, Svea; Köppen, Susan; Frauenheim, Thomas; Ciacchi, Lucio Colombi

2014-06-21

We investigate the adsorption behavior of four different amino acids (glutamine, glutamate, serine, cysteine) on the zinc oxide (101̄0) surface, comparing the geometry and energy associated with a number of different adsorption configurations. In doing this, we highlight the benefits and limits of using density-functional tight-binding (DFTB) with respect to standard density functional theory (DFT). The DFTB method is found to reliably reproduce the DFT adsorption geometries. Analysis of the adsorption configurations emphasizes the fundamental role of the first hydration layer in mediating the interactions between the amino acids and the surface. Direct surface-molecule bonds are found to form predominantly via the carboxylate groups of the studied amino acids. No surface-mediated chemical reactions are observed, with the notable exception of a proton transfer from the thiol group of cysteine to a hydroxyl group of the surface hydration layer. The adsorption energies are found to be dominated both by the formation of direct or indirect surface-molecule hydrogen bonds, but also by the rearrangement of the hydrogen-bond network in surface proximity in a non-intuitive way. Energetic comparisons between DFTB and DFT are made difficult on one side by the long time necessary to achieve convergence of potential energy values in MD simulations and on the other side by the necessity of including higher-order corrections to DFTB to obtain a good description of the hydrogen bond energetics. Overall, our results suggest that DFTB is a good reference method to set the correct chemical states and the initial geometries of hybrid biomolecule/ZnO systems to be simulated with non-reactive force fields.

An assessment of pure, hybrid, meta, and hybrid-meta GGA density functional theory methods for open-shell systems: the case of the nonheme iron enzyme 8R-LOX.

Science.gov (United States)

Bushnell, Eric A C; Gauld, James W

2013-01-15

The performance of a range density functional theory functionals combined in a quantum mechanical (QM)/molecular mechanical (MM) approach was investigated in their ability to reliably provide geometries, electronic distributions, and relative energies of a multicentered open-shell mechanistic intermediate in the mechanism 8R-Lipoxygenase. With the use of large QM/MM active site chemical models, the smallest average differences in geometries between the catalytically relevant quartet and sextet complexes were obtained with the B3LYP(*) functional. Moreover, in the case of the relative energies between (4) II and (6) II, the use of the B3LYP(*) functional provided a difference of 0.0 kcal mol(-1). However, B3LYP(±) and B3LYP also predicted differences in energies of less than 1 kcal mol(-1). In the case of describing the electronic distribution (i.e., spin density), the B3LYP(*), B3LYP, or M06-L functionals appeared to be the most suitable. Overall, the results obtained suggest that for systems with multiple centers having unpaired electrons, the B3LYP(*) appears most well rounded to provide reliable geometries, electronic structures, and relative energies. Copyright © 2012 Wiley Periodicals, Inc.

Application of time-dependent current-density-functional theory to nonlocal exchange-correlation effects in polymers

NARCIS (Netherlands)

van Faassen, M; de Boeij, PL; van Leeuwen, R; Berger, JA; Snijders, JG

2003-01-01

We provide a successful approach towards the solution of the longstanding problem of the large overestimation of the static polarizability of conjugated oligomers obtained using the local density approximation within density-functional theory. The local approximation is unable to describe the highly

Excitonic effects in solids : time-dependent density functional theory versus the Bethe-Salpeter equation

International Nuclear Information System (INIS)

Sagmeister, S.

2009-01-01

The aim of this work is to compare two state-of-the-art methods for the investigation of excitonic effects in solids, namely Time-Dependent Density Functional Theory (TDDFT) and Many-Body Perturbation Theory (MBPT), for selected simple gap systems as well as semiconducting polymers. Within TDDFT, the linear response framework is used and the Dyson equation for the density-density response function is solved, whereas within MBPT, the Bethe-Salpeter equation (BSE) for the electron-hole correlation function is solved. The dielectric function is obtained as a last step. Both techniques take into account the excitonic effects caused by the interaction of electron-hole pairs. In the former these effects are included in the exchange-correlation (xc) kernel, whereas in the latter they are located in the interaction kernel of the BSE. Kohn-Sham single-particle wave functions obtained from Density Functional Theory within the linearized augmented planewave (LAPW) method are used to calculate all relevant quantities of the formalism. For the simple systems GaAs, Si and LiF are chosen. The role of several approximations to the xc kernel is studied and it is found that for GaAs and Si simple semi-empirical models provide a dielectric function in accordance with the BSE. For the case of LiF, being a system with a weak screening and a strongly bound exciton, only an xc kernel derived from MBPT yields reasonable results but still a slight discrepancy to the BSE is observed. Finally, the semiconducting polymers poly-acetylene and poly(phenylene-vinylene) (PPV) are studied. For both materials the concept of semi-empirical approximations to the xc kernel turns out to be ambiguous due to their low-dimensional character. In the case of poly-acetylene, the xc kernel derived from MBPT yields a dielectric function which is in close but not exact agreement with the one obtained from the BSE. (author) [de

Hybrid Density Functional Study of the Local Structures and Energy Levels of CaAl2O4:Ce3.

Science.gov (United States)

Lou, Bibo; Jing, Weiguo; Lou, Liren; Zhang, Yongfan; Yin, Min; Duan, Chang-Kui

2018-05-03

First-principles calculations were carried out for the electronic structures of Ce 3+ in calcium aluminate phosphors, CaAl 2 O 4 , and their effects on luminescence properties. Hybrid density functional approaches were used to overcome the well-known underestimation of band gaps of conventional density functional approaches and to calculate the energy levels of Ce 3+ ions more accurately. The obtained 4f-5d excitation and emission energies show good consistency with measured values. A detailed energy diagram of all three sites is obtained, which explains qualitatively all of the luminescent phenomena. With the results of energy levels calculated by combining the hybrid functional of Heyd, Scuseria, and Ernzerhof (HSE06) and the constraint occupancy approach, we are able to construct a configurational coordinate diagram to analyze the processes of capture of a hole or an electron and luminescence. This approach can be applied for systematic high-throughput calculations in predicting Ce 3+ activated luminescent materials with a moderate computing requirement.

Optical properties of Al nanostructures from time dependent density functional theory

KAUST Repository

Mokkath, Junais Habeeb

2016-04-05

The optical properties of Al nanostructures are investigated by means of time dependent density functional theory, considering chains of varying length and ladders/stripes of varying aspect ratio. The absorption spectra show redshifting for increasing length and aspect ratio. For the chains the absorption is dominated by HOMO → LUMO transitions, whereas ladders and stripes reveal more complex spectra of plasmonic nature above a specific aspect ratio.

Decay of hollow states in time-dependent density functional theory

Energy Technology Data Exchange (ETDEWEB)

Kapoor, Varun; Bauer, Dieter [Institut fuer Physik, Wismarsche Str. 43-45, Universitaet Rostock, Rostock-18051 (Germany)

2012-07-01

Hollow or multiply excited states are inaccessible in time dependent density functional theory (TDDFT) using adiabatic Kohn-Sham potentials. We determine the exact Kohn Sham (KS) potential for doubly excited states in an exactly solvable model Helium atom. The exact single-particle density corresponds to the energetically lowest quasi-stationary state in the exact KS potential. We describe how this exact potential controls the decay by a barrier whose origin is traced back to phase of the exact KS orbital. The potential controls the barrier height and width in order for the density to tunnel out and decay with the same rate as the doubly excited state in the ab initio time-dependent Schroedinger calculation. Instead, adiabatic KS potentials only show direct photoionization but no autoionization. A frequency-dependent linear response kernel would be necessary in order to capture the decay of autoionizing states.

Time-Dependent Density Functional Theory for Open Systems and Its Applications.

Science.gov (United States)

Chen, Shuguang; Kwok, YanHo; Chen, GuanHua

2018-02-20

Photovoltaic devices, electrochemical cells, catalysis processes, light emitting diodes, scanning tunneling microscopes, molecular electronics, and related devices have one thing in common: open quantum systems where energy and matter are not conserved. Traditionally quantum chemistry is confined to isolated and closed systems, while quantum dissipation theory studies open quantum systems. The key quantity in quantum dissipation theory is the reduced system density matrix. As the reduced system density matrix is an O(M! × M!) matrix, where M is the number of the particles of the system of interest, quantum dissipation theory can only be employed to simulate systems of a few particles or degrees of freedom. It is thus important to combine quantum chemistry and quantum dissipation theory so that realistic open quantum systems can be simulated from first-principles. We have developed a first-principles method to simulate the dynamics of open electronic systems, the time-dependent density functional theory for open systems (TDDFT-OS). Instead of the reduced system density matrix, the key quantity is the reduced single-electron density matrix, which is an N × N matrix where N is the number of the atomic bases of the system of interest. As the dimension of the key quantity is drastically reduced, the TDDFT-OS can thus be used to simulate the dynamics of realistic open electronic systems and efficient numerical algorithms have been developed. As an application, we apply the method to study how quantum interference develops in a molecular transistor in time domain. We include electron-phonon interaction in our simulation and show that quantum interference in the given system is robust against nuclear vibration not only in the steady state but also in the transient dynamics. As another application, by combining TDDFT-OS with Ehrenfest dynamics, we study current-induced dissociation of water molecules under scanning tunneling microscopy and follow its time dependent

The density functional theory and the charged fluid molecular dynamics

International Nuclear Information System (INIS)

Hansen, J.P.; Zerah, G.

1993-01-01

Car and Parrinello had the idea of combining the density functional theory (Hohenberg, Kohn and Sham) to the 'molecular dynamics' numerical modelling method, in order to simulate metallic or co-valent solids and liquids from the first principles. The objective of this paper is to present a simplified version of this method ab initio, applicable to classical and quantal charged systems. The method is illustrated with recent results on charged colloidal suspensions and highly correlated electron-proton plasmas. 1 fig., 21 refs

Validation of missed space-group symmetry in X-ray powder diffraction structures with dispersion-corrected density functional theory

DEFF Research Database (Denmark)

Hempler, Daniela; Schmidt, Martin U.; Van De Streek, Jacco

2017-01-01

More than 600 molecular crystal structures with correct, incorrect and uncertain space-group symmetry were energy-minimized with dispersion-corrected density functional theory (DFT-D, PBE-D3). For the purpose of determining the correct space-group symmetry the required tolerance on the atomic...... with missed symmetry were investigated by dispersion-corrected density functional theory. In 98.5% of the cases the correct space group is found....

Descriptions of carbon isotopes within the energy density functional theory

Energy Technology Data Exchange (ETDEWEB)

Ismail, Atef [Fundamental and Applied Sciences Department, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 31750 Tronoh, Perak, Malaysia and Department of Physics, Al-Azhar University, 71524 Assiut (Egypt); Cheong, Lee Yen; Yahya, Noorhana [Fundamental and Applied Sciences Department, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 31750 Tronoh, Perak (Malaysia); Tammam, M. [Department of Physics, Al-Azhar University, 71524 Assiut (Egypt)

2014-10-24

Within the energy density functional (EDF) theory, the structure properties of Carbon isotopes are systematically studied. The shell model calculations are done for both even-A and odd-A nuclei, to study the structure of rich-neutron Carbon isotopes. The EDF theory indicates the single-neutron halo structures in {sup 15}C, {sup 17}C and {sup 19}C, and the two-neutron halo structures in {sup 16}C and {sup 22}C nuclei. It is also found that close to the neutron drip-line, there exist amazing increase in the neutron radii and decrease on the binding energies BE, which are tightly related with the blocking effect and correspondingly the blocking effect plays a significant role in the shell model configurations.

From The Cover: The X3LYP extended density functional for accurate descriptions of nonbond interactions, spin states, and thermochemical properties

Science.gov (United States)

Xu, Xin; Goddard, William A., III

2004-03-01

We derive the form for an exact exchange energy density for a density decaying with Gaussian-like behavior at long range. Based on this, we develop the X3LYP (extended hybrid functional combined with Lee-Yang-Parr correlation functional) extended functional for density functional theory to significantly improve the accuracy for hydrogen-bonded and van der Waals complexes while also improving the accuracy in heats of formation, ionization potentials, electron affinities, and total atomic energies [over the most popular and accurate method, B3LYP (Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional)]. X3LYP also leads to a good description of dipole moments, polarizabilities, and accurate excitation energies from s to d orbitals for transition metal atoms and ions. We suggest that X3LYP will be useful for predicting ligand binding in proteins and DNA.

Solving large nonlinear generalized eigenvalue problems from Density Functional Theory calculations in parallel

DEFF Research Database (Denmark)

Bendtsen, Claus; Nielsen, Ole Holm; Hansen, Lars Bruno

2001-01-01

The quantum mechanical ground state of electrons is described by Density Functional Theory, which leads to large minimization problems. An efficient minimization method uses a self-consistent field (SCF) solution of large eigenvalue problems. The iterative Davidson algorithm is often used, and we...

An efficient method for hybrid density functional calculation with spin-orbit coupling

Science.gov (United States)

Wang, Maoyuan; Liu, Gui-Bin; Guo, Hong; Yao, Yugui

2018-03-01

In first-principles calculations, hybrid functional is often used to improve accuracy from local exchange correlation functionals. A drawback is that evaluating the hybrid functional needs significantly more computing effort. When spin-orbit coupling (SOC) is taken into account, the non-collinear spin structure increases computing effort by at least eight times. As a result, hybrid functional calculations with SOC are intractable in most cases. In this paper, we present an approximate solution to this problem by developing an efficient method based on a mixed linear combination of atomic orbital (LCAO) scheme. We demonstrate the power of this method using several examples and we show that the results compare very well with those of direct hybrid functional calculations with SOC, yet the method only requires a computing effort similar to that without SOC. The presented technique provides a good balance between computing efficiency and accuracy, and it can be extended to magnetic materials.

Benchmarking density-functional-theory calculations of rotational g tensors and magnetizabilities using accurate coupled-cluster calculations.

Science.gov (United States)

Lutnaes, Ola B; Teale, Andrew M; Helgaker, Trygve; Tozer, David J; Ruud, Kenneth; Gauss, Jürgen

2009-10-14

An accurate set of benchmark rotational g tensors and magnetizabilities are calculated using coupled-cluster singles-doubles (CCSD) theory and coupled-cluster single-doubles-perturbative-triples [CCSD(T)] theory, in a variety of basis sets consisting of (rotational) London atomic orbitals. The accuracy of the results obtained is established for the rotational g tensors by careful comparison with experimental data, taking into account zero-point vibrational corrections. After an analysis of the basis sets employed, extrapolation techniques are used to provide estimates of the basis-set-limit quantities, thereby establishing an accurate benchmark data set. The utility of the data set is demonstrated by examining a wide variety of density functionals for the calculation of these properties. None of the density-functional methods are competitive with the CCSD or CCSD(T) methods. The need for a careful consideration of vibrational effects is clearly illustrated. Finally, the pure coupled-cluster results are compared with the results of density-functional calculations constrained to give the same electronic density. The importance of current dependence in exchange-correlation functionals is discussed in light of this comparison.

Unsupervised Learning and Pattern Recognition of Biological Data Structures with Density Functional Theory and Machine Learning.

Science.gov (United States)

Chen, Chien-Chang; Juan, Hung-Hui; Tsai, Meng-Yuan; Lu, Henry Horng-Shing

2018-01-11

By introducing the methods of machine learning into the density functional theory, we made a detour for the construction of the most probable density function, which can be estimated by learning relevant features from the system of interest. Using the properties of universal functional, the vital core of density functional theory, the most probable cluster numbers and the corresponding cluster boundaries in a studying system can be simultaneously and automatically determined and the plausibility is erected on the Hohenberg-Kohn theorems. For the method validation and pragmatic applications, interdisciplinary problems from physical to biological systems were enumerated. The amalgamation of uncharged atomic clusters validated the unsupervised searching process of the cluster numbers and the corresponding cluster boundaries were exhibited likewise. High accurate clustering results of the Fisher's iris dataset showed the feasibility and the flexibility of the proposed scheme. Brain tumor detections from low-dimensional magnetic resonance imaging datasets and segmentations of high-dimensional neural network imageries in the Brainbow system were also used to inspect the method practicality. The experimental results exhibit the successful connection between the physical theory and the machine learning methods and will benefit the clinical diagnoses.

Characterization of Phase Transition in Heisenberg Fluids from Density Functional Theory

International Nuclear Information System (INIS)

Li Liangsheng; Li Li; Chen Xiaosong

2009-01-01

The phase transition of Heisenberg fluid has been investigated with the density functional theory in mean-field approximation (MF). The matrix of the second derivatives of the grand canonical potential Ω with respect to the particle density fluctuations and the magnetization fluctuations has been investigated and diagonalized. The smallest eigenvalue being 0 signalizes the phase instability and the related eigenvector characterizes this phase transition. We find a Curie line where the order parameter is pure magnetization and a spinodal where the order parameter is a mixture of particle density and magnetization. Along the spinodal, the character of phase instability changes continuously from predominant condensation to predominant ferromagnetic phase transition with the decrease of total density. The spinodal meets the Curie line at the critical endpoint with the reduced density ρ* = ρσ 3 = 0.224 and the reduced temperature T* = kT/ element of = 1.87 (σ is the diameter of Heisenberg hard sphere and element of is the coupling constant).

Prediction of half-metallic properties in TlCrS{sub 2} and TlCrSe{sub 2} based on density functional theory

Energy Technology Data Exchange (ETDEWEB)

Hashimzade, F.M.; Huseinova, D.A. [Institute of Physics, National Academy of Sciences of Azerbaijan, AZ 1143 Baku (Azerbaijan); Jahangirli, Z.A. [Institute of Physics, National Academy of Sciences of Azerbaijan, AZ 1143 Baku (Azerbaijan); Institute of Radiation Problems, National Academy of Sciences of Azerbaijan, AZ 1143 Baku (Azerbaijan); Mehdiyev, B.H., E-mail: bachschi@yahoo.de [Institute of Physics, National Academy of Sciences of Azerbaijan, AZ 1143 Baku (Azerbaijan)

2017-08-01

Highlights: • Half-metallic properties of TlCrS2, TlCrSe2 and hypothetical TlCrSSe have been investigated by first-principles all-electron full-potential linearized augmented plane wave plus local orbital (FP-LAPW+lo) method based on density functional theory (DFT). • Total magnetic moment keeps its integer value on a relatively wide range of changes in volume (−10% ÷ 10%) for TlCrS2 and TlCrSSe, while total magnetic moment TlCrSe2 decreases with increasing volume, approaching to integer value 3 μB. • The states at the Fermi level in the case of spin-up channel consist of a hybridization of p-states of the atom S(Se) with d-states of Cr. - Abstract: Half-metallic properties of TlCrS{sub 2}, TlCrSe{sub 2} and hypothetical TlCrSSe have been investigated by first-principles all-electron full-potential linearized augmented plane wave plus local orbital (FP-LAPW+lo) method based on density functional theory (DFT). The results of calculations show that TlCrS{sub 2} and TlCrSSe are half-metals with energy gap (E{sub g}) ∼0.12 eV for spin-down channel. Strong hybridization of p-state of chalchogen and d-state of Cr leads to bonding and antibonding states and subsequently to the appearance of a gap in spin-down channel of TlCrS{sub 2} and TlCrSSe. In the case of TlCrSe{sub 2}, there is a partial hybridization and p-state is partially present in the DOS at Fermi level making this compound nearly half-metallic. The present calculations revealed that total magnetic moment keeps its integer value on a relatively wide range of changes in volume (−10% ÷ 10%) for TlCrS{sub 2} and TlCrSSe, while total magnetic moment of TlCrSe{sub 2} decreases with increasing volume approaching to integer value 3 μB.

Density functional theory study of the concerted pyrolysis mechanism for lignin models

Science.gov (United States)

Thomas Elder; Ariana Beste

2014-01-01

ABSTRACT: Studies on the pyrolysis mechanisms of lignin model compounds have largely focused on initial homolytic cleavage reactions. It has been noted, however, that concerted mechanisms may also account for observed product formation. In the current work, the latter processes are examined and compared to the former, by the application of density functional theory...

Density functional theory of simple polymers in a slit pore. III. Surface tension

International Nuclear Information System (INIS)

Hooper, Justin B.; McCoy, John D.; Curro, John G.; Swol, Frank van

2000-01-01

In a previous study of tangent hard-site chains near a surface, the inhomogeneous density profiles were found through density functional theory. In the current study, the surface tensions of these systems are found from the results of the previous study through a thermodynamic integration. The calculated surface tensions are then compared to those found directly through computer simulation. Both the surface tension and surface excess for polymeric systems are shown to differ qualitatively from those of atomic systems, although certain similarities are seen at high densities. (c) 2000 American Institute of Physics

Hot-electron-assisted femtochemistry at surfaces: A time-dependent density functional theory approach

DEFF Research Database (Denmark)

Gavnholt, Jeppe; Rubio, Angel; Olsen, Thomas

2009-01-01

Using time-evolution time-dependent density functional theory (TDDFT) within the adiabatic local-density approximation, we study the interactions between single electrons and molecular resonances at surfaces. Our system is a nitrogen molecule adsorbed on a ruthenium surface. The surface is modele...... resonance and the lowering of the resonance energy due to an image charge effect. Finally we apply the TDDFT procedure to only consider the decay of molecular excitations and find that it agrees quite well with the width of the projected density of Kohn-Sham states....

Nematic-smectic A and nematic-solid transitions of parallel hard spherocylinders from density functional theory

NARCIS (Netherlands)

University Utrecht

1992-01-01

A simple density functional theory for the various liquid-crystalline phases of parallel hard spherocylinders is formulated on the basis of Pynn's ansatz for the direct correlation function of the spherocylinders. Fair agreement with the computer simulations is found.

Improved modification for the density-functional theory calculation of thermodynamic properties for C-H-O composite compounds.

Science.gov (United States)

Liu, Min Hsien; Chen, Cheng; Hong, Yaw Shun

2005-02-08

A three-parametric modification equation and the least-squares approach are adopted to calibrating hybrid density-functional theory energies of C(1)-C(10) straight-chain aldehydes, alcohols, and alkoxides to accurate enthalpies of formation DeltaH(f) and Gibbs free energies of formation DeltaG(f), respectively. All calculated energies of the C-H-O composite compounds were obtained based on B3LYP6-311++G(3df,2pd) single-point energies and the related thermal corrections of B3LYP6-31G(d,p) optimized geometries. This investigation revealed that all compounds had 0.05% average absolute relative error (ARE) for the atomization energies, with mean value of absolute error (MAE) of just 2.1 kJ/mol (0.5 kcal/mol) for the DeltaH(f) and 2.4 kJ/mol (0.6 kcal/mol) for the DeltaG(f) of formation.

Tuning electronic properties in graphene quantum dots by chemical functionalization: Density functional theory calculations

Science.gov (United States)

Abdelsalam, Hazem; Elhaes, Hanan; Ibrahim, Medhat A.

2018-03-01

The energy gap and dipole moment of chemically functionalized graphene quantum dots are investigated by density functional theory. The energy gap can be tuned through edge passivation by different elements or groups. Edge passivation by oxygen considerably decreases the energy gap in hexagonal nanodots. Edge states in triangular quantum dots can also be manipulated by passivation with fluorine. The dipole moment depends on: (a) shape and edge termination of the quantum dot, (b) attached group, and (c) position to which the groups are attached. Depending on the position of attached groups, the total dipole can be increased, decreased, or eliminated.

Electronic states of aryl radical functionalized graphenes: Density functional theory study

Science.gov (United States)

Tachikawa, Hiroto; Kawabata, Hiroshi

2016-06-01

Functionalized graphenes are known as a high-performance molecular device. In the present study, the structures and electronic states of the aryl radical functionalized graphene have been investigated by the density functional theory (DFT) method to elucidate the effects of functionalization on the electronic states of graphene (GR). Also, the mechanism of aryl radical reaction with GR was investigated. The benzene, biphenyl, p-terphenyl, and p-quaterphenyl radicals [denoted by (Bz) n (n = 1-4), where n means numbers of benzene rings in aryl radical] were examined as aryl radicals. The DFT calculation of GR-(Bz) n (n = 1-4) showed that the aryl radical binds to the carbon atom of GR, and a C-C single bond was formed. The binding energies of aryl radicals to GR were calculated to be ca. 6.0 kcal mol-1 at the CAM-B3LYP/6-311G(d,p) level. It was found that the activation barrier exists in the aryl radical addition: the barrier heights were calculated to be 10.0 kcal mol-1. The electronic states of GR-(Bz) n were examined on the basis of theoretical results.

Nuclear reactivity indices in the context of spin polarized density functional theory

International Nuclear Information System (INIS)

Cardenas, Carlos; Lamsabhi, Al Mokhtar; Fuentealba, Patricio

2006-01-01

In this work, the nuclear reactivity indices of density functional theory have been generalized to the spin polarized case and their relationship to electron spin polarized indices has been established. In particular, the spin polarized version of the nuclear Fukui function has been proposed and a finite difference approximation has been used to evaluate it. Applications to a series of triatomic molecules demonstrate the ability of the new functions to predict the geometrical changes due to a change in the spin multiplicity. The main equations in the different ensembles have also been presented

The study of structural characteristic of yttrium(3+) aquaion by density functional methods

International Nuclear Information System (INIS)

Buz'ko, V.Yu.; Sukhno, I.V.; Polushin, A.A.; Panyushkin, V.T.

2006-01-01

Structural characteristics of the yttrium aquaion Y(H 2 O) 8 3+ are calculated by DFT technique using density functionals SVWN5, B3LYP, B3P86, O3LYP, B3PW91, B1LYP, B971, MPW1PW9, PBE1PBE, BHandH and BNandHLYP. All calculations are carried out by means of GAUSSIAN-03 code. Structural characteristics of the Y(H 2 O) 8 3+ aquaion obtained by the use of the density functional method agree satisfactorily with the experiment. Hybrid nonlocal GGA-functionals foretell worse the structural characteristics of the Y(H 2 O) 8 3+ aquaion as compared with the reasonable simple combined functional BHandH and the simplest functional SVWN5 of the LSDA theory [ru

Electronic zero-point oscillations in the strong-interaction limit of density functional theory

NARCIS (Netherlands)

Gori Giorgi, P.; Vignale, G.; Seidl, M.

2009-01-01

The exchange-correlation energy in Kohn-Sham density functional theory can be expressed exactly in terms of the change in the expectation of the electron-electron repulsion operator when, in the many-electron Hamiltonian, this same operator is multiplied by a real parameter λ varying between 0

Basic Density-Functional Theory—an Overview

Science.gov (United States)

von Barth, U.

In these notes I have given a personally flavored exposA~© of static density-functional theory (DFT). I have started from standard many-body physics at a very elementary level and then gradually introduced the basic concepts of DFT. Successively more advanced topics are added and at the end I even discuss a few not yet published theories. The discussion represents many of the personal views of the author and there is no attempt at being comprehensive. I fully realize that I am often ‘unfair’ in treating the achievements of other researchers. Many topics of standard DFT are deliberately left out like, e.g., time-dependence, excitations, and magnetic or relativistic effects. These notes represent a compilation of a series of lectures given at at the EXC!TING Summer School DFT beyond the ground state at RiksgrA~¤nsen, Sweden in June of 2003.

A physically motivated sparse cubature scheme with applications to molecular density-functional theory

International Nuclear Information System (INIS)

Rodriguez, Juan I; Thompson, David C; Anderson, James S M; Thomson, Jordan W; Ayers, Paul W

2008-01-01

We present a novel approach for performing multi-dimensional integration of arbitrary functions. The method starts with Smolyak-type sparse grids as cubature formulae on the unit cube and uses a transformation of coordinates based on the conditional distribution method to adapt those formulae to real space. Our method is tested on integrals in one, two, three and six dimensions. The three dimensional integration formulae are used to evaluate atomic interaction energies via the Gordon-Kim model. The six dimensional integration formulae are tested in conjunction with the nonlocal exchange-correlation energy functional proposed by Lee and Parr. This methodology is versatile and powerful; we contemplate application to frozen-density embedding, next-generation molecular-mechanics force fields, 'kernel-type' exchange-correlation energy functionals and pair-density functional theory

Strong Correlation in Kohn-Sham Density Functional Theory

NARCIS (Netherlands)

Malet, F.; Gori Giorgi, P.

2012-01-01

We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approximation corresponds to a highly nonlocal density functional whose functional derivative can be easily constructed,

Applications of the Conceptual Density Functional Theory Indices to Organic Chemistry Reactivity.

Science.gov (United States)

Domingo, Luis R; Ríos-Gutiérrez, Mar; Pérez, Patricia

2016-06-09

Theoretical reactivity indices based on the conceptual Density Functional Theory (DFT) have become a powerful tool for the semiquantitative study of organic reactivity. A large number of reactivity indices have been proposed in the literature. Herein, global quantities like the electronic chemical potential μ, the electrophilicity ω and the nucleophilicity N indices, and local condensed indices like the electrophilic P k + and nucleophilic P k - Parr functions, as the most relevant indices for the study of organic reactivity, are discussed.

Simple descriptors for proton-conducting perovskites from density functional theory

DEFF Research Database (Denmark)

Bork, Nicolai Christian; Bonanos, Nikolaos; Rossmeisl, Jan

2010-01-01

series of (pseudo)cubic perovskites, ABO3, have been investigated using density functional theory calculations. The structures have been optimized and thermodynamic properties and activation energies for the relevant steps of the hydrogen/proton diffusion mechanism have been calculated using...... the nudged elastic band path technique. We find a strong correlation between the O-H binding energy for hydrogen/proton uptake in perovskites and the energy barriers involved in the observed Grotthuss-type diffusion process. We demonstrate the possibility of estimating diffusion rates based on O-H binding...

Functionalized carbon nanotube based hybrid electrochemical capacitors using neutral bromide redox-active electrolyte for enhancing energy density

Science.gov (United States)

Tang, Xiaohui; Lui, Yu Hui; Chen, Bolin; Hu, Shan

2017-06-01

A hybrid electrochemical capacitor (EC) with enhanced energy density is realized by integrating functionalized carbon nanotube (FCNT) electrodes with redox-active electrolyte that has a neutral pH value (1 M Na2SO4 and 0.5 M KBr mixed aqueous solution). The negative electrode shows an electric double layer capacitor-type behavior. On the positive electrode, highly reversible Br-/Br3- redox reactions take place, presenting a battery-type behavior, which contributes to increase the capacitance of the hybrid cell. The voltage window of the whole cell is extended up to 1.5 V because of the high over-potentials of oxygen and hydrogen evolution reactions in the neutral electrolyte. Compared with raw CNT, the FCNT has better wettability in the aqueous electrolyte and contributes to increase the electric double layer capacitance of the cell. As a result, the maximum energy density of 28.3 Wh kg-1 is obtained from the hybrid EC at 0.5 A g-1 without sacrificing its power density, which is around 4 times larger than that of the electrical double layer capacitor constructed by FCNT electrodes and 1 M Na2SO4 electrolyte. Moreover, the discharge capacity retained 86.3% of its initial performance after 10000 cycles of galvanostatic charge and discharge test (10 A/g), suggesting its long life cycle even at high current loading.

Simulation of electron energy loss spectra of nanomaterials with linear-scaling density functional theory

International Nuclear Information System (INIS)

Tait, E W; Payne, M C; Ratcliff, L E; Haynes, P D; Hine, N D M

2016-01-01

Experimental techniques for electron energy loss spectroscopy (EELS) combine high energy resolution with high spatial resolution. They are therefore powerful tools for investigating the local electronic structure of complex systems such as nanostructures, interfaces and even individual defects. Interpretation of experimental electron energy loss spectra is often challenging and can require theoretical modelling of candidate structures, which themselves may be large and complex, beyond the capabilities of traditional cubic-scaling density functional theory. In this work, we present functionality to compute electron energy loss spectra within the onetep linear-scaling density functional theory code. We first demonstrate that simulated spectra agree with those computed using conventional plane wave pseudopotential methods to a high degree of precision. The ability of onetep to tackle large problems is then exploited to investigate convergence of spectra with respect to supercell size. Finally, we apply the novel functionality to a study of the electron energy loss spectra of defects on the (1 0 1) surface of an anatase slab and determine concentrations of defects which might be experimentally detectable. (paper)

Time dependent density functional theory of light absorption in dense plasmas: application to iron-plasma

International Nuclear Information System (INIS)

Grimaldi, F.; Grimaldi-Lecourt, A.; Dharma-Wardana, M.W.C.

1986-10-01

The objective of this paper is to present a simple time-dependent calculation of the light absorption cross section for a strongly coupled partially degenerate plasma so as to transcend the usual single-particle picture. This is achieved within the density functional theory (DFT) of plasmas by generalizing the method given by Zangwill and Soven for atomic calculations at zero temperature. The essential feature of the time dependent DFT is the correct treatment of the relaxation of the system under the external field. Exploratory calculations for a Fe-plasma at 100 eV show new features in the absorption cross section which are absent in the usual single particle theory. These arise from inter-shell correlations, channel mixing and self-energy effects. These many-body effects introduce significant modifications to the radiative properties of plasmas and are shown to be efficiently calculable by time dependent density functional theory (TD-DFT)

Time dependent density functional theory of light absorption in dense plasmas: application to iron-plasma

International Nuclear Information System (INIS)

Grimaldi, F.; Grimaldi-Lecourt, A.; Dharma-Wardana, M.W.C.

1985-02-01

The objective of this paper is to present a simple time-dependent calculation of the light absorption cross section for a strongly coupled partially degenerate plasma so as to transcend the usual single-particle picture. This is achieved within the density functional theory (DFT) of plasmas by generalizing the method given by Zangwill and Soven for atomic calculations at zero temperature. The essential feature of the time dependent DFT is the correct treatment of the relaxation of the system under the external field. Exploratory calculations for an Fe-plasma at 100 eV show new features in the absorption cross section which are absent in the usual single particle theory. These arise from inter-shell correlations, channel mixing and self-energy effects. These many-body effects introduce significant modifications to the radiative properties of plasma and are shown to be efficiently calculable by time dependent density functional theory (TD-DFT)

The oxygen reduction reaction mechanism on Pt(111) from density functional theory calculations

DEFF Research Database (Denmark)

Tripkovic, Vladimir; Skulason, Egill; Siahrostami, Samira

2010-01-01

We study the oxygen reduction reaction (ORR) mechanism on a Pt(1 1 1) surface using density functional theory calculations We find that at low overpotentials the surface is covered with a half dissociated water layer We estimate the barrier for proton transfer to this surface and the barrier...

Density functional theory for field emission from carbon nano-structures

Energy Technology Data Exchange (ETDEWEB)

Li, Zhibing, E-mail: stslzb@mail.sysu.edu.cn

2015-12-15

Electron field emission is understood as a quantum mechanical many-body problem in which an electronic quasi-particle of the emitter is converted into an electron in vacuum. Fundamental concepts of field emission, such as the field enhancement factor, work-function, edge barrier and emission current density, will be investigated, using carbon nanotubes and graphene as examples. A multi-scale algorithm basing on density functional theory is introduced. We will argue that such a first principle approach is necessary and appropriate for field emission of nano-structures, not only for a more accurate quantitative description, but, more importantly, for deeper insight into field emission. - Highlights: • Applications of DFT to electron field emission of nano-structures are reviewed. • Fundamental concepts of field emission are re-visited with emphasis on the many-body effects. • New insights to field emission of nano-structures are obtained by multi-scale DFT calculations. • It is shown that the exchange–correlation effect on the emission barrier is significant. • Spontaneous symmetry breaking in field emission of CNT has been predicted.

Covariant density functional theory beyond mean field and applications for nuclei far from stability

International Nuclear Information System (INIS)

Ring, P

2010-01-01

Density functional theory provides a very powerful tool for a unified microscopic description of nuclei all over the periodic table. It is not only successful in reproducing bulk properties of nuclear ground states such as binding energies, radii, or deformation parameters, but it also allows the investigation of collective phenomena, such as giant resonances and rotational excitations. However, it is based on the mean field concept and therefore it has its limits. We discuss here two methods based based on covariant density functional theory going beyond the mean field concept, (i) models with an energy dependent self energy allowing the coupling to complex configurations and a quantitative description of the width of giant resonances and (ii) methods of configuration mixing between Slater determinants with different deformation and orientation providing are very successful description of transitional nuclei and quantum phase transitions.

UV-Vis optoelectronic properties of α-SnWO4: A comparative experimental and density functional theory based study

KAUST Repository

Ziani, Ahmed

2015-09-03

We report a combined experimental and theoretical study on the optoelectronic properties of α-SnWO4 for UV-Vis excitation. The experimentally measured values for thin films were systematically compared with high-accuracy density functional theory and density functional perturbation theory using the HSE06 functional. The α-SnWO4 material shows an indirect bandgap of 1.52 eV with high absorption coefficient in the visible-light range (>2 × 105 cm−1). The results show relatively high dielectric constant (>30) and weak diffusion properties (large effective masses) of excited carriers.

Remarks on time-dependent [current]-density functional theory for open quantum systems.

Science.gov (United States)

Yuen-Zhou, Joel; Aspuru-Guzik, Alán

2013-08-14

Time-dependent [current]-density functional theory for open quantum systems (OQS) has emerged as a formalism that can incorporate dissipative effects in the dynamics of many-body quantum systems. Here, we review and clarify some formal aspects of these theories that have been recently questioned in the literature. In particular, we provide theoretical support for the following conclusions: (1) contrary to what we and others had stated before, within the master equation framework, there is in fact a one-to-one mapping between vector potentials and current densities for fixed initial state, particle-particle interaction, and memory kernel; (2) regardless of the first conclusion, all of our recently suggested Kohn-Sham (KS) schemes to reproduce the current and particle densities of the original OQS, and in particular, the use of a KS closed driven system, remains formally valid; (3) the Lindblad master equation maintains the positivity of the density matrix regardless of the time-dependence of the Hamiltonian or the dissipation operators; (4) within the stochastic Schrödinger equation picture, a one-to-one mapping from stochastic vector potential to stochastic current density for individual trajectories has not been proven so far, except in the case where the vector potential is the same for every member of the ensemble, in which case, it reduces to the Lindblad master equation picture; (5) master equations may violate certain desired properties of the density matrix, such as positivity, but they remain as one of the most useful constructs to study OQS when the environment is not easily incorporated explicitly in the calculation. The conclusions support our previous work as formally rigorous, offer new insights into it, and provide a common ground to discuss related theories.

Optimal control theory for quantum-classical systems: Ehrenfest molecular dynamics based on time-dependent density-functional theory

International Nuclear Information System (INIS)

Castro, A; Gross, E K U

2014-01-01

We derive the fundamental equations of an optimal control theory for systems containing both quantum electrons and classical ions. The system is modeled with Ehrenfest dynamics, a non-adiabatic variant of molecular dynamics. The general formulation, that needs the fully correlated many-electron wavefunction, can be simplified by making use of time-dependent density-functional theory. In this case, the optimal control equations require some modifications that we will provide. The abstract general formulation is complemented with the simple example of the H 2 + molecule in the presence of a laser field. (paper)

Time-dependent internal density functional theory formalism and Kohn-Sham scheme for self-bound systems

International Nuclear Information System (INIS)

Messud, Jeremie

2009-01-01

The stationary internal density functional theory (DFT) formalism and Kohn-Sham scheme are generalized to the time-dependent case. It is proven that, in the time-dependent case, the internal properties of a self-bound system (such as an atomic nuclei or a helium droplet) are all defined by the internal one-body density and the initial state. A time-dependent internal Kohn-Sham scheme is set up as a practical way to compute the internal density. The main difference from the traditional DFT formalism and Kohn-Sham scheme is the inclusion of the center-of-mass correlations in the functional.

Influence of van der Waals forces on the adsorption structure of benzene on silicon studied using density functional theory

DEFF Research Database (Denmark)

Johnston, Karen; Kleis, Jesper; Lundqvist, Bengt

2008-01-01

Two different adsorption configurations of benzene on the Si(001)-(2×1) surface, the tight-bridge and butterfly structures, were studied using density functional theory. Several exchange and correlation functionals were used, including the recently developed van der Waals density functional (vd...

Local-scaling density-functional method: Intraorbit and interorbit density optimizations

International Nuclear Information System (INIS)

Koga, T.; Yamamoto, Y.; Ludena, E.V.

1991-01-01

The recently proposed local-scaling density-functional theory provides us with a practical method for the direct variational determination of the electron density function ρ(r). The structure of ''orbits,'' which ensures the one-to-one correspondence between the electron density ρ(r) and the N-electron wave function Ψ({r k }), is studied in detail. For the realization of the local-scaling density-functional calculations, procedures for intraorbit and interorbit optimizations of the electron density function are proposed. These procedures are numerically illustrated for the helium atom in its ground state at the beyond-Hartree-Fock level

General atomistic approach for modeling metal-semiconductor interfaces using density functional theory and nonequilibrium Green's function

DEFF Research Database (Denmark)

Stradi, Daniele; Martinez, Umberto; Blom, Anders

2016-01-01

Metal-semiconductor contacts are a pillar of modern semiconductor technology. Historically, their microscopic understanding has been hampered by the inability of traditional analytical and numerical methods to fully capture the complex physics governing their operating principles. Here we introduce...... an atomistic approach based on density functional theory and nonequilibrium Green's function, which includes all the relevant ingredients required to model realistic metal-semiconductor interfaces and allows for a direct comparison between theory and experiments via I-Vbias curve simulations. We apply...... interfaces as it neglects electron tunneling, and that finite-size atomistic models have problems in describing these interfaces in the presence of doping due to a poor representation of space-charge effects. Conversely, the present method deals effectively with both issues, thus representing a valid...

Density functional theory for adsorption of gas mixtures in metal-organic frameworks.

Science.gov (United States)

Liu, Yu; Liu, Honglai; Hu, Ying; Jiang, Jianwen

2010-03-04

In this work, a recently developed density functional theory in three-dimensional space was extended to the adsorption of gas mixtures. Weighted density approximations to the excess free energy with different weighting functions were adopted for both repulsive and attractive contributions. An equation of state for hard-sphere mixtures and a modified Benedict-Webb-Rubin equation for Lennard-Jones mixtures were used to estimate the excess free energy of a uniform fluid. The theory was applied to the adsorption of CO(2)/CH(4) and CO(2)/N(2) mixtures in two metal-organic frameworks: ZIF-8 and Zn(2)(BDC)(2)(ted). To validate the theoretical predictions, grand canonical Monte Carlo simulations were also conducted. The predicted adsorption and selectivity from DFT were found to agree well with the simulation results. CO(2) has stronger adsorption than CH(4) and N(2), particularly in Zn(2)(BDC)(2)(ted). The selectivity of CO(2) over CH(4) or N(2) increases with increasing pressure as attributed to the cooperative interactions of adsorbed CO(2) molecules. The composition of the gas mixture exhibits a significant effect on adsorption but not on selectivity.

Density functional theory investigation of two-dimensional dipolar fermions in a harmonic trap

International Nuclear Information System (INIS)

Ustunel, Hande; Abedinpour, Saeed H; Tanatar, B

2014-01-01

We investigate the behavior of polarized dipolar fermions in a two-dimensional harmonic trap in the framework of the density functional theory (DFT) formalism using the local density approximation. We treat only a few particles interacting moderately. Important results were deduced concerning key characteristics of the system such as total energy and particle density. Our results indicate that, at variance with Coulombic systems, the exchange- correlation component was found to provide a large contribution to the total energy for a large range of interaction strengths and particle numbers. In addition, the density profiles of the dipoles are shown to display important features around the origin that is not possible to capture by earlier, simpler treatments of such systems

The structural and electronic properties of amine-functionalized boron nitride nanotubes via ammonia plasmas: a density functional theory study

International Nuclear Information System (INIS)

Cao Fenglei; Ji Yuemeng; Zhao Cunyuan; Ren Wei

2009-01-01

The reaction behavior of the chemical modification of boron nitride nanotubes (BNNTs) with ammonia plasmas has been investigated by density functional theory (DFT) calculations. Unlike previously studied functionalization with NH 3 and amino functional groups, we found that NH 2 * radicals involved in the ammonia plasmas can be covalently incorporated to BNNTs through a strong single B-N bond. Subsequently, the H * radicals also involved in the ammonia plasmas would prefer to combine with the N atoms neighboring the NH 2 -functionalized B atoms. Our study revealed that this reaction behavior can be elucidated using the frontier orbital theory. The calculated band structures and density of states (DOS) indicate that this modification is an effective method to modulate the electronic properties of BNNTs. We have discussed various defects on the surface of BNNTs generated by collisions of N 2 + ions. For most defects considered, the reactivity of the functionalization of BNNTs with NH 2 * are enhanced. Our conclusions are independent of the chirality, and the diameter dependence of the reaction energies is presented.

Corrections to the density-functional theory electronic spectrum: Copper phthalocyanine

DEFF Research Database (Denmark)

Vazquez, Hector; Jelinek, P.; Brandbyge, Mads

2009-01-01

A method for improving the electronic spectrum of standard Density-Functional Theory (DFT) calculations (i.e., LDA or GGA approximations) is presented, and its application is discussed for the case of the copper phthalocyanine (CuPc) molecule. The method is based on a treatment of exchange...... and correlation in a many-body Hamiltonian, and it leads to easy-to-evaluate corrections to the DFT eigenvalues. Self-interaction is largely corrected, so that the modified energy levels do not suffer from spurious crossings, as often encountered for CuPc in DFT, and they remedy the standard underestimation...... or semiempirical functionals for molecular levels, it can be easily applied to any local-orbital DFT approach, improving on several important limitations of standard DFT methods....

Dispersion- and Exchange-Corrected Density Functional Theory for Sodium Ion Hydration.

Science.gov (United States)

Soniat, Marielle; Rogers, David M; Rempe, Susan B

2015-07-14

A challenge in density functional theory is developing functionals that simultaneously describe intermolecular electron correlation and electron delocalization. Recent exchange-correlation functionals address those two issues by adding corrections important at long ranges: an atom-centered pairwise dispersion term to account for correlation and a modified long-range component of the electron exchange term to correct for delocalization. Here we investigate how those corrections influence the accuracy of binding free energy predictions for sodium-water clusters. We find that the dual-corrected ωB97X-D functional gives cluster binding energies closest to high-level ab initio methods (CCSD(T)). Binding energy decomposition shows that the ωB97X-D functional predicts the smallest ion-water (pairwise) interaction energy and larger multibody contributions for a four-water cluster than most other functionals - a trend consistent with CCSD(T) results. Also, ωB97X-D produces the smallest amounts of charge transfer and the least polarizable waters of the density functionals studied, which mimics the lower polarizability of CCSD. When compared with experimental binding free energies, however, the exchange-corrected CAM-B3LYP functional performs best (error <1 kcal/mol), possibly because of its parametrization to experimental formation enthalpies. For clusters containing more than four waters, "split-shell" coordination must be considered to obtain accurate free energies in comparison with experiment.

A Comparative Density Functional Theory and Density Functional Tight Binding Study of Phases of Nitrogen Including a High Energy Density Material N8

Directory of Open Access Journals (Sweden)

Nicholas Capel

2015-11-01

Full Text Available We present a comparative dispersion-corrected Density Functional Theory (DFT and Density Functional Tight Binding (DFTB-D study of several phases of nitrogen, including the well-known alpha, beta, and gamma phases as well as recently discovered highly energetic phases: covalently bound cubic gauche (cg nitrogen and molecular (vdW-bound N8 crystals. Among several tested parametrizations of N–N interactions for DFTB, we identify only one that is suitable for modeling of all these phases. This work therefore establishes the applicability of DFTB-D to studies of phases, including highly metastable phases, of nitrogen, which will be of great use for modelling of dynamics of reactions involving these phases, which may not be practical with DFT due to large required space and time scales. We also derive a dispersion-corrected DFT (DFT-D setup (atom-centered basis parameters and Grimme dispersion parameters tuned for accurate description simultaneously of several nitrogen allotropes including covalently and vdW-bound crystals and including high-energy phases.

Rationale for switching to nonlocal functionals in density functional theory.

Science.gov (United States)

Lazić, P; Atodiresei, N; Caciuc, V; Brako, R; Gumhalter, B; Blügel, S

2012-10-24

Density functional theory (DFT) has been steadily improving over the past few decades, becoming the standard tool for electronic structure calculations. The early local functionals (LDA) were eventually replaced by more accurate semilocal functionals (GGA) which are in use today. A major persisting drawback is the lack of the nonlocal correlation which is at the core of dispersive (van der Waals) forces, so that a large and important class of systems remains outside the scope of DFT. The vdW-DF correlation functional of Langreth and Lundqvist, published in 2004, was the first nonlocal functional which could be easily implemented. Beyond expectations, the nonlocal functional has brought significant improvement to systems that were believed not to be sensitive to nonlocal correlations. In this paper, we use the example of graphene nanodomes growing on the Ir(111) surface, where with an increase of the size of the graphene islands the character of the bonding changes from strong chemisorption towards almost pure physisorption. We demonstrate how the seamless character of the vdW-DF functionals makes it possible to treat all regimes self-consistently, proving to be a systematic and consistent improvement of DFT regardless of the nature of bonding. We also discuss the typical surface science example of CO adsorption on (111) surfaces of metals, which shows that the nonlocal correlation may also be crucial for strongly chemisorbed systems. We briefly discuss open questions, in particular the choice of the most appropriate exchange part of the functional. As the vdW-DF begins to appear implemented self-consistently in a number of popular DFT codes, with numerical costs close to the GGA calculations, we draw the attention of the DFT community to the advantages and benefits of the adoption of this new class of functionals.

Rationale for switching to nonlocal functionals in density functional theory

International Nuclear Information System (INIS)

Lazić, P; Atodiresei, N; Caciuc, V; Blügel, S; Brako, R; Gumhalter, B

2012-01-01

Density functional theory (DFT) has been steadily improving over the past few decades, becoming the standard tool for electronic structure calculations. The early local functionals (LDA) were eventually replaced by more accurate semilocal functionals (GGA) which are in use today. A major persisting drawback is the lack of the nonlocal correlation which is at the core of dispersive (van der Waals) forces, so that a large and important class of systems remains outside the scope of DFT. The vdW-DF correlation functional of Langreth and Lundqvist, published in 2004, was the first nonlocal functional which could be easily implemented. Beyond expectations, the nonlocal functional has brought significant improvement to systems that were believed not to be sensitive to nonlocal correlations. In this paper, we use the example of graphene nanodomes growing on the Ir(111) surface, where with an increase of the size of the graphene islands the character of the bonding changes from strong chemisorption towards almost pure physisorption. We demonstrate how the seamless character of the vdW-DF functionals makes it possible to treat all regimes self-consistently, proving to be a systematic and consistent improvement of DFT regardless of the nature of bonding. We also discuss the typical surface science example of CO adsorption on (111) surfaces of metals, which shows that the nonlocal correlation may also be crucial for strongly chemisorbed systems. We briefly discuss open questions, in particular the choice of the most appropriate exchange part of the functional. As the vdW-DF begins to appear implemented self-consistently in a number of popular DFT codes, with numerical costs close to the GGA calculations, we draw the attention of the DFT community to the advantages and benefits of the adoption of this new class of functionals.

Exchange-Correlation Effects for Noncovalent Interactions in Density Functional Theory.

Science.gov (United States)

Otero-de-la-Roza, A; DiLabio, Gino A; Johnson, Erin R

2016-07-12

In this article, we develop an understanding of how errors from exchange-correlation functionals affect the modeling of noncovalent interactions in dispersion-corrected density-functional theory. Computed CCSD(T) reference binding energies for a collection of small-molecule clusters are decomposed via a molecular many-body expansion and are used to benchmark density-functional approximations, including the effect of semilocal approximation, exact-exchange admixture, and range separation. Three sources of error are identified. Repulsion error arises from the choice of semilocal functional approximation. This error affects intermolecular repulsions and is present in all n-body exchange-repulsion energies with a sign that alternates with the order n of the interaction. Delocalization error is independent of the choice of semilocal functional but does depend on the exact exchange fraction. Delocalization error misrepresents the induction energies, leading to overbinding in all induction n-body terms, and underestimates the electrostatic contribution to the 2-body energies. Deformation error affects only monomer relaxation (deformation) energies and behaves similarly to bond-dissociation energy errors. Delocalization and deformation errors affect systems with significant intermolecular orbital interactions (e.g., hydrogen- and halogen-bonded systems), whereas repulsion error is ubiquitous. Many-body errors from the underlying exchange-correlation functional greatly exceed in general the magnitude of the many-body dispersion energy term. A functional built to accurately model noncovalent interactions must contain a dispersion correction, semilocal exchange, and correlation components that minimize the repulsion error independently and must also incorporate exact exchange in such a way that delocalization error is absent.

Towards a Density Functional Theory Exchange-Correlation Functional able to describe localization/delocalization

Science.gov (United States)

Mattsson, Ann E.; Wills, John M.

2013-03-01

The inability to computationally describe the physics governing the properties of actinides and their alloys is the poster child of failure of existing Density Functional Theory exchange-correlation functionals. The intricate competition between localization and delocalization of the electrons, present in these materials, exposes the limitations of functionals only designed to properly describe one or the other situation. We will discuss the manifestation of this competition in real materials and propositions on how to construct a functional able to accurately describe properties of these materials. I addition we will discuss both the importance of using the Dirac equation to describe the relativistic effects in these materials, and the connection to the physics of transition metal oxides. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Obtaining Hartree-Fock and density functional theory doubly excited states with Car-Parrinello density matrix search

Science.gov (United States)

Liang, Wenkel; Isborn, Christine M.; Li, Xiaosong

2009-11-01

The calculation of doubly excited states is one of the major problems plaguing the modern day excited state workhorse methodology of linear response time dependent Hartree-Fock (TDHF) and density function theory (TDDFT). We have previously shown that the use of a resonantly tuned field within real-time TDHF and TDDFT is able to simultaneously excite both the α and β electrons to achieve the two-electron excited states of minimal basis H2 and HeH+ [C. M. Isborn and X. Li, J. Chem. Phys. 129, 204107 (2008)]. We now extend this method to many electron systems with the use of our Car-Parrinello density matrix search (CP-DMS) with a first-principles fictitious mass method for wave function optimization [X. Li, C. L. Moss, W. Liang, and Y. Feng, J. Chem. Phys. 130, 234115 (2009)]. Real-time TDHF/TDDFT is used during the application of the laser field perturbation, driving the electron density toward the doubly excited state. The CP-DMS method then converges the density to the nearest stationary state. We present these stationary state doubly excited state energies and properties at the HF and DFT levels for H2, HeH+, lithium hydride, ethylene, and butadiene.

Total-energy Assisted Tight-binding Method Based on Local Density Approximation of Density Functional Theory

Science.gov (United States)

Fujiwara, Takeo; Nishino, Shinya; Yamamoto, Susumu; Suzuki, Takashi; Ikeda, Minoru; Ohtani, Yasuaki

2018-06-01

A novel tight-binding method is developed, based on the extended Hückel approximation and charge self-consistency, with referring the band structure and the total energy of the local density approximation of the density functional theory. The parameters are so adjusted by computer that the result reproduces the band structure and the total energy, and the algorithm for determining parameters is established. The set of determined parameters is applicable to a variety of crystalline compounds and change of lattice constants, and, in other words, it is transferable. Examples are demonstrated for Si crystals of several crystalline structures varying lattice constants. Since the set of parameters is transferable, the present tight-binding method may be applicable also to molecular dynamics simulations of large-scale systems and long-time dynamical processes.

Application of an excited state LDA exchange energy functional for the calculation of transition energy of atoms within time-independent density functional theory