Lehtovaara, Lauri; Havu, Ville; Puska, Martti
2011-10-21
We present an all-electron method for time-dependent density functional theory which employs hierarchical nonuniform finite-element bases and the time-propagation approach. The method is capable of treating linear and nonlinear response of valence and core electrons to an external field. We also introduce (i) a preconditioner for the propagation equation, (ii) a stable way to implement absorbing boundary conditions, and (iii) a new kind of absorbing boundary condition inspired by perfectly matched layers. PMID:22029294
All-electron Kohn–Sham density functional theory on hierarchic finite element spaces
Schauer, Volker; Linder, Christian
2013-10-01
In this work, a real space formulation of the Kohn–Sham equations is developed, making use of the hierarchy of finite element spaces from different polynomial order. The focus is laid on all-electron calculations, having the highest requirement onto the basis set, which must be able to represent the orthogonal eigenfunctions as well as the electrostatic potential. A careful numerical analysis is performed, which points out the numerical intricacies originating from the singularity of the nuclei and the necessity for approximations in the numerical setting, with the ambition to enable solutions within a predefined accuracy. In this context the influence of counter-charges in the Poisson equation, the requirement of a finite domain size, numerical quadratures and the mesh refinement are examined as well as the representation of the electrostatic potential in a high order finite element space. The performance and accuracy of the method is demonstrated in computations on noble gases. In addition the finite element basis proves its flexibility in the calculation of the bond-length as well as the dipole moment of the carbon monoxide molecule.
Kuwahara, Riichi; Tadokoro, Yoichi; Ohno, Kaoru
2014-08-28
In this paper, we calculate kinetic and potential energy contributions to the electronic ground-state total energy of several isolated atoms (He, Be, Ne, Mg, Ar, and Ca) by using the local density approximation (LDA) in density functional theory, the Hartree-Fock approximation (HFA), and the self-consistent GW approximation (GWA). To this end, we have implemented self-consistent HFA and GWA routines in our all-electron mixed basis code, TOMBO. We confirm that virial theorem is fairly well satisfied in all of these approximations, although the resulting eigenvalue of the highest occupied molecular orbital level, i.e., the negative of the ionization potential, is in excellent agreement only in the case of the GWA. We find that the wave function of the lowest unoccupied molecular orbital level of noble gas atoms is a resonating virtual bound state, and that of the GWA spreads wider than that of the LDA and thinner than that of the HFA. PMID:25173006
NASA Astrophysics Data System (ADS)
Rivelino, Roberto; Malaspina, Thaciana; Fileti, Eudes E.
2009-01-01
We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C60(OH)n] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke’s three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C60(OH)36 . Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1 , 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C60(OH)24 isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4 , 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.
Gulans, Andris; Kontur, Stefan; Meisenbichler, Christian; Nabok, Dmitrii; Pavone, Pasquale; Rigamonti, Santiago; Sagmeister, Stephan; Werner, Ute; Draxl, Claudia
2014-09-10
Linearized augmented planewave methods are known as the most precise numerical schemes for solving the Kohn-Sham equations of density-functional theory (DFT). In this review, we describe how this method is realized in the all-electron full-potential computer package, exciting. We emphasize the variety of different related basis sets, subsumed as (linearized) augmented planewave plus local orbital methods, discussing their pros and cons and we show that extremely high accuracy (microhartrees) can be achieved if the basis is chosen carefully. As the name of the code suggests, exciting is not restricted to ground-state calculations, but has a major focus on excited-state properties. It includes time-dependent DFT in the linear-response regime with various static and dynamical exchange-correlation kernels. These are preferably used to compute optical and electron-loss spectra for metals, molecules and semiconductors with weak electron-hole interactions. exciting makes use of many-body perturbation theory for charged and neutral excitations. To obtain the quasi-particle band structure, the GW approach is implemented in the single-shot approximation, known as G(0)W(0). Optical absorption spectra for valence and core excitations are handled by the solution of the Bethe-Salpeter equation, which allows for the description of strongly bound excitons. Besides these aspects concerning methodology, we demonstrate the broad range of possible applications by prototypical examples, comprising elastic properties, phonons, thermal-expansion coefficients, dielectric tensors and loss functions, magneto-optical Kerr effect, core-level spectra and more. PMID:25135665
Mizutani, U; Inukai, M; Sato, H; Zijlstra, E S; Lin, Q
2014-05-16
There are three key electronic parameters in elucidating the physics behind the Hume–Rothery electron concentration rule: the square of the Fermi diameter (2kF)2, the square of the critical reciprocal lattice vector and the electron concentration parameter or the number of itinerant electrons per atom e/a. We have reliably determined these three parameters for 10 Rhombic Triacontahedron-type 2/1–2/1–2/1 (N = 680) and 1/1–1/1–1/1 (N = 160–162) approximants by making full use of the full-potential linearized augmented plane wave-Fourier band calculations based on all-electron density-functional theory. We revealed that the 2/1–2/1–2/1 approximants Al13Mg27Zn45 and Na27Au27Ga31 belong to two different sub-groups classified in terms of equal to 126 and 109 and could explain why they take different e/a values of 2.13 and 1.76, respectively. Among eight 1/1–1/1–1/1 approximants Al3Mg4Zn3, Al9Mg8Ag3, Al21Li13Cu6, Ga21Li13Cu6, Na26Au24Ga30, Na26Au37Ge18, Na26Au37Sn18 and Na26Cd40Pb6, the first two, the second two and the last four compounds were classified into three sub-groups with = 50, 46 and 42; and were claimed to obey the e/a = 2.30, 2.10–2.15 and 1.70–1.80 rules, respectively.
Mitin, Alexander V; van Wüllen, Christoph
2006-02-14
A two-component quasirelativistic Hamiltonian based on spin-dependent effective core potentials is used to calculate ionization energies and electron affinities of the heavy halogen atom bromine through the superheavy element 117 (eka-astatine) as well as spectroscopic constants of the homonuclear dimers of these atoms. We describe a two-component Hartree-Fock and density-functional program that treats spin-orbit coupling self-consistently within the orbital optimization procedure. A comparison with results from high-order Douglas-Kroll calculations--for the superheavy systems also with zeroth-order regular approximation and four-component Dirac results--demonstrates the validity of the pseudopotential approximation. The density-functional (but not the Hartree-Fock) results show very satisfactory agreement with theoretical coupled cluster as well as experimental data where available, such that the theoretical results can serve as an estimate for the hitherto unknown properties of astatine, element 117, and their dimers. PMID:16483205
Ching, W. Y.; Aryal, Sitram; Rulis, Paul; Schnick, Wolfgang
2011-04-15
Using density-functional-theory-based ab initio methods, the electronic structure and physical properties of the newly synthesized nitride BeP{sub 2}N{sub 4} with a phenakite-type structure and the predicted high-pressure spinel phase of BeP{sub 2}N{sub 4} are studied in detail. It is shown that both polymorphs are wide band-gap semiconductors with relatively small electron effective masses at the conduction-band minima. The spinel-type phase is more covalently bonded due to the increased number of P-N bonds for P at the octahedral sites. Calculations of mechanical properties indicate that the spinel-type polymorph is a promising superhard material with notably large bulk, shear, and Young's moduli. Also calculated are the Be K, P K, P L{sub 3}, and N K edges of the electron energy-loss near-edge structure for both phases. They show marked differences because of the different local environments of the atoms in the two crystalline polymorphs. These differences will be very useful for the experimental identification of the products of high-pressure syntheses targeting the predicted spinel-type phase of BeP{sub 2}N{sub 4}.
NASA Astrophysics Data System (ADS)
Betzinger, Markus; Friedrich, Christoph; Görling, Andreas; Blügel, Stefan
2015-12-01
We present a methodology to calculate frequency and momentum dependent all-electron response functions determined within Kohn-Sham density functional theory. It overcomes the main obstacle in calculating response functions in practice, which is the slow convergence with respect to the number of unoccupied states and the basis-set size. In this approach, the usual sum-over-states expression of perturbation theory is complemented by the response of the orbital basis functions, explicitly constructed by radial integrations of frequency-dependent Sternheimer equations. To an essential extent an infinite number of unoccupied states are included in this way. Furthermore, the response of the core electrons is treated virtually exactly, which is out of reach otherwise. The method is an extension of the recently introduced incomplete-basis-set correction (IBC) [Betzinger et al., Phys. Rev. B 85, 245124 (2012), 10.1103/PhysRevB.85.245124; Phys. Rev. B 88, 075130 (2013), 10.1103/PhysRevB.88.075130] to the frequency and momentum domain. We have implemented the generalized IBC within the all-electron full-potential linearized augmented-plane-wave method and demonstrate for rocksalt BaO the improved convergence of the dynamical Kohn-Sham polarizability. We apply this technique to compute (a) quasiparticle energies employing the COHSEX approximation for the self-energy of many-body perturbation theory and (b) all-electron RPA correlation energies. It is shown that the favorable convergence of the polarizability is passed over to the COHSEX and RPA calculation.
NASA Astrophysics Data System (ADS)
Betzinger, Markus; Friedrich, Christoph; Blügel, Stefan
2013-08-01
In a previous publication [Betzinger, Friedrich, Görling, and Blügel, Phys. Rev. BPRBMDO1098-012110.1103/PhysRevB.85.245124 85, 245124 (2012)] we presented a technique to compute accurate all-electron response functions, e.g., the density response function, within the full-potential linearized augmented-plane-wave (FLAPW) method. Response contributions that are not captured (completely) within the finite Hilbert space spanned by the LAPW basis are taken into account by an incomplete-basis-set correction (IBC). The latter is based on a formal response of the basis functions themselves, which is derived by exploiting their dependence on the effective potential. Its construction requires the solution of radial differential equations, having the form of Sternheimer equations, by numerical integration. The approach includes a formally exact treatment of the response contribution from the core states. While we restricted the formalism to spherical perturbations in the previous work, we here generalize the formalism to nonspherical perturbations. The improvements are demonstrated with exact-exchange optimized-effective-potential (EXX-OEP) calculations of antiferromagnetic NiO. It is shown that with the generalized IBC a basis-set convergence is realized that is as fast as in density-functional theory calculations using standard local or semilocal functionals. The EXX-OEP band gap, magnetic moment, and spectral function of NiO are in substantially better agreement with experiment than results obtained from calculations with local and semilocal functionals.
Hybrid functionals within the all-electron FLAPW method: Implementation and applications of PBE0
NASA Astrophysics Data System (ADS)
Betzinger, Markus; Friedrich, Christoph; Blügel, Stefan
2010-05-01
We present an efficient implementation of the Perdew-Burke-Ernzerhof hybrid functional PBE0 within the full-potential linearized augmented-plane-wave (FLAPW) method. The Hartree-Fock exchange term, which is a central ingredient of hybrid functionals, gives rise to a computationally expensive nonlocal potential in the one-particle Schrödinger equation. The matrix elements of this exchange potential are calculated with the help of an auxiliary basis that is constructed from products of FLAPW basis functions. By representing the Coulomb interaction in this basis the nonlocal exchange term becomes a Brillouin-zone sum over vector-matrix-vector products. The Coulomb matrix is calculated only once at the beginning of a self-consistent-field cycle. We show that it can be made sparse by a suitable unitary transformation of the auxiliary basis, which accelerates the computation of the vector-matrix-vector products considerably. Additionally, we exploit spatial and time-reversal symmetry to identify the nonvanishing exchange matrix elements in advance and to restrict the k summations for the nonlocal potential to an irreducible set of k points. Favorable convergence of the self-consistent-field cycle is achieved by a nested density-only and density-matrix iteration scheme. We discuss the convergence with respect to the parameters of our numerical scheme and show results for a variety of semiconductors and insulators, including the oxides ZnO, EuO, Al2O3 , and SrTiO3 , where the PBE0 hybrid functional improves the band gaps and the description of localized states in comparison with the PBE functional. Furthermore, we find that in contrast to conventional local exchange-correlation functionals ferromagnetic EuO is correctly predicted to be a semiconductor.
NASA Astrophysics Data System (ADS)
Matveev, Oleg; Shvaika, Andrij; Devereaux, Thomas; Freericks, James
The charge-density-wave phase of the Falicov-Kimball model displays a number of anomalous behavior including the appearance of subgap density of states as the temperature increases. These subgap states should have a significant impact on transport properties, particularly the nonlinear response of the system to a large dc electric field. Using the Kadanoff-Baym-Keldysh formalism, we employ nonequilibrium dynamical mean-field theory to exactly solve for this nonlinear response. We examine both the current and the order parameter of the conduction electrons as the ordered system is driven by a dc electric field. Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, Lviv, Ukraine.
NASA Astrophysics Data System (ADS)
Betzinger, Markus; Friedrich, Christoph; Görling, Andreas; Blügel, Stefan
2012-06-01
The optimized-effective-potential method is a special technique to construct local Kohn-Sham potentials from general orbital-dependent energy functionals. In a recent publication [M. Betzinger, C. Friedrich, S. Blügel, A. Görling, Phys. Rev. BPRBMDO1098-012110.1103/PhysRevB.83.045105 83, 045105 (2011)] we showed that uneconomically large basis sets were required to obtain a smooth local potential without spurious oscillations within the full-potential linearized augmented-plane-wave method. This could be attributed to the slow convergence behavior of the density response function. In this paper, we derive an incomplete-basis-set correction for the response, which consists of two terms: (1) a correction that is formally similar to the Pulay correction in atomic-force calculations and (2) a numerically more important basis response term originating from the potential dependence of the basis functions. The basis response term is constructed from the solutions of radial Sternheimer equations in the muffin-tin spheres. With these corrections the local potential converges at much smaller basis sets, at much fewer states, and its construction becomes numerically very stable. We analyze the improvements for rock-salt ScN and report results for BN, AlN, and GaN, as well as the perovskites CaTiO3, SrTiO3, and BaTiO3. The incomplete-basis-set correction can be applied to other electronic-structure methods with potential-dependent basis sets and opens the perspective to investigate a broad spectrum of problems in theoretical solid-state physics that involve response functions.
NASA Astrophysics Data System (ADS)
Levchenko, Sergey V.; Ren, Xinguo; Wieferink, Jürgen; Johanni, Rainer; Rinke, Patrick; Blum, Volker; Scheffler, Matthias
2015-07-01
We describe a framework to evaluate the Hartree-Fock exchange operator for periodic electronic-structure calculations based on general, localized atom-centered basis functions. The functionality is demonstrated by hybrid-functional calculations of properties for several semiconductors. In our implementation of the Fock operator, the Coulomb potential is treated either in reciprocal space or in real space, where the sparsity of the density matrix can be exploited for computational efficiency. Computational aspects, such as the rigorous avoidance of on-the-fly disk storage, and a load-balanced parallel implementation, are also discussed. We demonstrate linear scaling of our implementation with system size by calculating the electronic structure of a bulk semiconductor (GaAs) with up to 1,024 atoms per unit cell without compromising the accuracy.
Gedanken densities and exact constraints in density functional theory
Perdew, John P.; Ruzsinszky, Adrienn; Sun, Jianwei; Burke, Kieron
2014-05-14
Approximations to the exact density functional for the exchange-correlation energy of a many-electron ground state can be constructed by satisfying constraints that are universal, i.e., valid for all electron densities. Gedanken densities are designed for the purpose of this construction, but need not be realistic. The uniform electron gas is an old gedanken density. Here, we propose a spherical two-electron gedanken density in which the dimensionless density gradient can be an arbitrary positive constant wherever the density is non-zero. The Lieb-Oxford lower bound on the exchange energy can be satisfied within a generalized gradient approximation (GGA) by bounding its enhancement factor or simplest GGA exchange-energy density. This enhancement-factor bound is well known to be sufficient, but our gedanken density shows that it is also necessary. The conventional exact exchange-energy density satisfies no such local bound, but energy densities are not unique, and the simplest GGA exchange-energy density is not an approximation to it. We further derive a strongly and optimally tightened bound on the exchange enhancement factor of a two-electron density, which is satisfied by the local density approximation but is violated by all published GGA's or meta-GGA’s. Finally, some consequences of the non-uniform density-scaling behavior for the asymptotics of the exchange enhancement factor of a GGA or meta-GGA are given.
Sharkey, Keeper L.; Pavanello, Michele; Bubin, Sergiy; Adamowicz, Ludwik
2009-12-15
A new algorithm for calculating the Hamiltonian matrix elements with all-electron explicitly correlated Gaussian functions for quantum-mechanical calculations of atoms with two p electrons or a single d electron have been derived and implemented. The Hamiltonian used in the approach was obtained by rigorously separating the center-of-mass motion and it explicitly depends on the finite mass of the nucleus. The approach was employed to perform test calculations on the isotopes of the carbon atom in their ground electronic states and to determine the finite-nuclear-mass corrections for these states.
Sharkey, Keeper L; Kirnosov, Nikita; Adamowicz, Ludwik
2013-03-14
A new algorithm for quantum-mechanical nonrelativistic calculation of the Hamiltonian matrix elements with all-electron explicitly correlated Gaussian functions for atoms with an arbitrary number of s electrons and with three p electrons, or one p electron and one d electron, or one f electron is developed and implemented. In particular the implementation concerns atomic states with L = 3 and M = 0. The Hamiltonian used in the approach is obtained by rigorously separating the center-of-mass motion from the laboratory-frame all particle Hamiltonian, and thus it explicitly depends on the finite mass of the nucleus. The approach is employed to perform test calculations on the lowest (2)F state of the two main isotopes of the lithium atom, (7)Li and (6)Li. PMID:23514465
Partition density functional theory
NASA Astrophysics Data System (ADS)
Nafziger, Jonathan
Partition density functional theory (PDFT) is a method for dividing a molecular electronic structure calculation into fragment calculations. The molecular density and energy corresponding to Kohn Sham density-functional theory (KS-DFT) may be exactly recovered from these fragments. Each fragment acts as an isolated system except for the influence of a global one-body 'partition' potential which deforms the fragment densities. In this work, the developments of PDFT are put into the context of other fragment-based density functional methods. We developed three numerical implementations of PDFT: One within the NWChem computational chemistry package using basis sets, and the other two developed from scratch using real-space grids. It is shown that all three of these programs can exactly reproduce a KS-DFT calculation via fragment calculations. The first of our in-house codes handles non-interacting electrons in arbitrary one-dimensional potentials with any number of fragments. This code is used to explore how the exact partition potential changes for different partitionings of the same system and also to study features which determine which systems yield non-integer PDFT occupations and which systems are locked into integer PDFT occupations. The second in-house code, CADMium, performs real-space calculations of diatomic molecules. Features of the exact partition potential are studied for a variety of cases and an analytical formula determining singularities in the partition potential is derived. We introduce an approximation for the non-additive kinetic energy and show how this quantity can be computed exactly. Finally a PDFT functional is developed to address the issues of static correlation and delocalization errors in approximations within DFT. The functional is applied to the dissociation of H2 + and H2.
]: a density functional theory investigation
NASA Astrophysics Data System (ADS)
Ulian, Gianfranco; Tosoni, Sergio; Valdrè, Giovanni
2014-09-01
In this work, we modelled the structure, the compressional behaviour and the physical properties of talc over a wide range of pressure using a quantum mechanical approach based on periodic boundary conditions. We adopted the density functional theory using the B3LYP-D* functional, which includes a correction for the dispersive forces and all-electron Gaussian-type orbitals basis sets. An atomic level description of the athermal pressure-induced structural modification of talc is provided. From the compression results, we obtained the athermal ( T = 0 K) bulk modulus ( K T0), its first derivative ( K') and the athermal volume at zero pressure ( V 0) by a third-order Birch-Murnaghan equation with parameters K T0 = 56.25 GPa, K' = 5.66 and V 0 = 450.34 Å3. The mechanical behaviour is highly anisotropic, as observed by the axial compressibility. The presented data are in very good agreement with recent experimental results obtained by single-crystal neutron and X-ray diffraction experiments.
Evarestov, R A; Losev, M V
2009-12-01
For the first time the convergence of the phonon frequencies and dispersion curves in terms of the supercell size is studied in ab initio frozen phonon calculations on LiF crystal. Helmann-Feynman forces over atomic displacements are found in all-electron calculations with the localized atomic functions (LCAO) basis using CRYSTAL06 program. The Parlinski-Li-Kawazoe method and FROPHO program are used to calculate the dynamical matrix and phonon frequencies of the supercells. For fcc lattice, it is demonstrated that use of the full supercell space group (including the supercell inner translations) enables to reduce essentially the number of the displacements under consideration. For Hartree-Fock (HF), PBE and hybrid PBE0, B3LYP, and B3PW exchange-correlation functionals the atomic basis set optimization is performed. The supercells up to 216 atoms (3 x 3 x 3 conventional unit cells) are considered. The phonon frequencies using the supercells of different size and shape are compared. For the commensurate with supercell k-points the best agreement of the theoretical results with the experimental data is found for B3PW exchange-correlation functional calculations with the optimized basis set. The phonon frequencies at the most non-commensurate k-points converged for the supercell consisting of 4 x 4 x 4 primitive cells and ensures the accuracy 1-2% in the thermodynamic properties calculated (the Helmholtz free energy, entropy, and heat capacity at the room temperature). PMID:19382176
Density-dependent covariant energy density functionals
Lalazissis, G. A.
2012-10-20
Relativistic nuclear energy density functionals are applied to the description of a variety of nuclear structure phenomena at and away fromstability line. Isoscalar monopole, isovector dipole and isoscalar quadrupole giant resonances are calculated using fully self-consistent relativistic quasiparticle randomphase approximation, based on the relativistic Hartree-Bogoliubovmodel. The impact of pairing correlations on the fission barriers in heavy and superheavy nuclei is examined. The role of pion in constructing desnity functionals is also investigated.
Phenomenological Relativistic Energy Density Functionals
Lalazissis, G. A.; Kartzikos, S.; Niksic, T.; Paar, N.; Vretenar, D.; Ring, P.
2009-08-26
The framework of relativistic nuclear energy density functionals is applied to the description of a variety of nuclear structure phenomena, not only in spherical and deformed nuclei along the valley of beta-stability, but also in exotic systems with extreme isospin values and close to the particle drip-lines. Dynamical aspects of exotic nuclear structure is explored using the fully consistent quasiparticle random-phase approximation based on the relativistic Hartree-Bogoliubov model. Recent applications of energy density functionals with explicit density dependence of the meson-nucleon couplings are presented.
Density Functionals of Chemical Bonding
Putz, Mihai V.
2008-01-01
The behavior of electrons in general many-electronic systems throughout the density functionals of energy is reviewed. The basic physico-chemical concepts of density functional theory are employed to highlight the energy role in chemical structure while its extended influence in electronic localization function helps in chemical bonding understanding. In this context the energy functionals accompanied by electronic localization functions may provide a comprehensive description of the global-local levels electronic structures in general and of chemical bonds in special. Becke-Edgecombe and author’s Markovian electronic localization functions are discussed at atomic, molecular and solid state levels. Then, the analytical survey of the main workable kinetic, exchange, and correlation density functionals within local and gradient density approximations is undertaken. The hierarchy of various energy functionals is formulated by employing both the parabolic and statistical correlation degree of them with the electronegativity and chemical hardness indices by means of quantitative structure-property relationship (QSPR) analysis for basic atomic and molecular systems. PMID:19325846
Semiclassical origins of density functionals
NASA Astrophysics Data System (ADS)
Burke, Kieron
By careful numerical analysis of non-relativistic atomic correlation energies, we show that (a) the local density approximation becomes relatively exact for the correlation energy as the atomic number approaches infinity, (b) we find the leading correction, which is about 38.5 milliHartrees per atom, (c) show how this correction dominates for larger atoms and (d) how to construct a generalized gradient approximation that respects this limit (See KB, A. Cancio, T. Gould, S. Pittalis, arXiv:1409.4834). The relevance to density functional calculations will also be explained. Support provided by NSF CHE-1464795.
NASA Astrophysics Data System (ADS)
Sharkey, Keeper L.; Bubin, Sergiy; Adamowicz, Ludwik
2014-11-01
Accurate variational nonrelativistic quantum-mechanical calculations are performed for the five lowest 1D and four lowest 3D states of the 9Be isotope of the beryllium atom. All-electron explicitly correlated Gaussian (ECG) functions are used in the calculations and their nonlinear parameters are optimized with the aid of the analytical energy gradient determined with respect to these parameters. The effect of the finite nuclear mass is directly included in the Hamiltonian used in the calculations. The singlet-triplet energy gaps between the corresponding 1D and 3D states, are reported.
Density functionals from deep learning
NASA Astrophysics Data System (ADS)
McMahon, Jeffrey
Density-functional theory is a formally exact description of a many-body quantum system in terms of its density; in practice, however, approximations to the universal density functional (DF) are necessary. Machine learning has recently been proposed as a novel approach to discover such a DF (or components of it). Conventional machine learning algorithms, however, are limited in their ability to process data in their raw form, leading to invariance and/or sensitivity issues. In this presentation, an alternative approach based on deep learning will be demonstrated. Deep learning allows computational models that are capable of discovering intricate structure in large and/or high-dimensional data sets with multiple levels of abstraction, and do not suffer from the aforementioned issues. Results from the application of this approach to the prediction of the kinetic-energy DF of noninteracting electrons will be presented. Using theoretical results from computer science, a connection between the underlying model and the theorems of Hohenberg and Kohn will also be suggested.
Periodic Density Functional Theory Solver using Multiresolution Analysis with MADNESS
NASA Astrophysics Data System (ADS)
Harrison, Robert; Thornton, William
2011-03-01
We describe the first implementation of the all-electron Kohn-Sham density functional periodic solver (DFT) using multi-wavelets and fast integral equations using MADNESS (multiresolution adaptive numerical environment for scientific simulation; http://code.google.com/p/m-a-d-n-e-s-s). The multiresolution nature of a multi-wavelet basis allows for fast computation with guaranteed precision. By reformulating the Kohn-Sham eigenvalue equation into the Lippmann-Schwinger equation, we can avoid using the derivative operator which allows better control of overall precision for the all-electron problem. Other highlights include the development of periodic integral operators with low-rank separation, an adaptable model potential for nuclear potential, and an implementation for Hartree Fock exchange. This work was supported by NSF project OCI-0904972 and made use of resources at the Center for Computational Sciences at Oak Ridge National Laboratory under contract DE-AC05-00OR22725.
Density functional theory: Foundations reviewed
NASA Astrophysics Data System (ADS)
Kryachko, Eugene S.; Ludeña, Eduardo V.
2014-11-01
Guided by the above motto (quotation), we review a broad range of issues lying at the foundations of Density Functional Theory, DFT, a theory which is currently omnipresent in our everyday computational study of atoms and molecules, solids and nano-materials, and which lies at the heart of modern many-body computational technologies. The key goal is to demonstrate that there are definitely the ways to improve DFT. We start by considering DFT in the larger context provided by reduced density matrix theory (RDMT) and natural orbital functional theory (NOFT), and examine the implications that N-representability conditions on the second-order reduced density matrix (2-RDM) have not only on RDMT and NOFT but, also, by extension, on the functionals of DFT. This examination is timely in view of the fact that necessary and sufficient N-representability conditions on the 2-RDM have recently been attained. In the second place, we review some problems appearing in the original formulation of the first Hohenberg-Kohn theorem which is still a subject of some controversy. In this vein we recall Lieb's comment on this proof and the extension to this proof given by Pino et al. (2009), and in this context examine the conditions that must be met in order that the one-to-one correspondence between ground-state densities and external potentials remains valid for finite subspaces (namely, the subspaces where all Kohn-Sham solutions are obtained in practical applications). We also consider the issue of whether the Kohn-Sham equations can be derived from basic principles or whether they are postulated. We examine this problem in relation to ab initio DFT. The possibility of postulating arbitrary Kohn-Sham-type equations, where the effective potential is by definition some arbitrary mixture of local and non-local terms, is discussed. We also deal with the issue of whether there exists a universal functional, or whether one should advocate instead the construction of problem
Sharkey, Keeper L; Adamowicz, Ludwik
2014-05-01
An algorithm for quantum-mechanical nonrelativistic variational calculations of L = 0 and M = 0 states of atoms with an arbitrary number of s electrons and with three p electrons have been implemented and tested in the calculations of the ground (4)S state of the nitrogen atom. The spatial part of the wave function is expanded in terms of all-electrons explicitly correlated Gaussian functions with the appropriate pre-exponential Cartesian angular factors for states with the L = 0 and M = 0 symmetry. The algorithm includes formulas for calculating the Hamiltonian and overlap matrix elements, as well as formulas for calculating the analytic energy gradient determined with respect to the Gaussian exponential parameters. The gradient is used in the variational optimization of these parameters. The Hamiltonian used in the approach is obtained by rigorously separating the center-of-mass motion from the laboratory-frame all-particle Hamiltonian, and thus it explicitly depends on the finite mass of the nucleus. With that, the mass effect on the total ground-state energy is determined. PMID:24811630
Unexpected properties of a density functional
Karwowski, J.; Stanke, M.
2005-02-01
An observation on a pathological behavior of an exact density functional derived from either relativistic (Dirac) or nonrelativistic (Levy-Leblond) quantum-mechanical equation is reported. As expected, in the case of a one-electron atom the variational minimum of this functional is equal to the exact ground-state energy. However, apart from the correct density, this minimum is reached also by an infinite set of densities which do not correspond to the exact wave function. This paradoxical property of the functional is related to the multicomponent structure of both Dirac and Levy-Leblond wave functions. In particular, imposing the correct boundary conditions upon the trial densities removes only a part of the fake solutions. The results of this study demonstrate that in density-functional theories derived from models based on multicomponent wave functions, one should not expect any simple relation between the accuracy of the energy and the correctness of the corresponding density.
Orbital-free density functional theory implementation with the projector augmented-wave method
NASA Astrophysics Data System (ADS)
Lehtomäki, Jouko; Makkonen, Ilja; Caro, Miguel A.; Harju, Ari; Lopez-Acevedo, Olga
2014-12-01
We present a computational scheme for orbital-free density functional theory (OFDFT) that simultaneously provides access to all-electron values and preserves the OFDFT linear scaling as a function of the system size. Using the projector augmented-wave method (PAW) in combination with real-space methods, we overcome some obstacles faced by other available implementation schemes. Specifically, the advantages of using the PAW method are twofold. First, PAW reproduces all-electron values offering freedom in adjusting the convergence parameters and the atomic setups allow tuning the numerical accuracy per element. Second, PAW can provide a solution to some of the convergence problems exhibited in other OFDFT implementations based on Kohn-Sham (KS) codes. Using PAW and real-space methods, our orbital-free results agree with the reference all-electron values with a mean absolute error of 10 meV and the number of iterations required by the self-consistent cycle is comparable to the KS method. The comparison of all-electron and pseudopotential bulk modulus and lattice constant reveal an enormous difference, demonstrating that in order to assess the performance of OFDFT functionals it is necessary to use implementations that obtain all-electron values. The proposed combination of methods is the most promising route currently available. We finally show that a parametrized kinetic energy functional can give lattice constants and bulk moduli comparable in accuracy to those obtained by the KS PBE method, exemplified with the case of diamond.
Orbital-free density functional theory implementation with the projector augmented-wave method
Lehtomäki, Jouko; Makkonen, Ilja; Harju, Ari; Lopez-Acevedo, Olga; Caro, Miguel A.
2014-12-21
We present a computational scheme for orbital-free density functional theory (OFDFT) that simultaneously provides access to all-electron values and preserves the OFDFT linear scaling as a function of the system size. Using the projector augmented-wave method (PAW) in combination with real-space methods, we overcome some obstacles faced by other available implementation schemes. Specifically, the advantages of using the PAW method are twofold. First, PAW reproduces all-electron values offering freedom in adjusting the convergence parameters and the atomic setups allow tuning the numerical accuracy per element. Second, PAW can provide a solution to some of the convergence problems exhibited in other OFDFT implementations based on Kohn-Sham (KS) codes. Using PAW and real-space methods, our orbital-free results agree with the reference all-electron values with a mean absolute error of 10 meV and the number of iterations required by the self-consistent cycle is comparable to the KS method. The comparison of all-electron and pseudopotential bulk modulus and lattice constant reveal an enormous difference, demonstrating that in order to assess the performance of OFDFT functionals it is necessary to use implementations that obtain all-electron values. The proposed combination of methods is the most promising route currently available. We finally show that a parametrized kinetic energy functional can give lattice constants and bulk moduli comparable in accuracy to those obtained by the KS PBE method, exemplified with the case of diamond.
A Safari Through Density Functional Theory
NASA Astrophysics Data System (ADS)
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.
Magnetic fields and density functional theory
Salsbury Jr., Freddie
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.
Density functional theory for carbon dioxide crystal
Chang, Yiwen; Mi, Jianguo Zhong, Chongli
2014-05-28
We present a density functional approach to describe the solid−liquid phase transition, interfacial and crystal structure, and properties of polyatomic CO{sub 2}. Unlike previous phase field crystal model or density functional theory, which are derived from the second order direct correlation function, the present density functional approach is based on the fundamental measure theory for hard-sphere repulsion in solid. More importantly, the contributions of enthalpic interactions due to the dispersive attractions and of entropic interactions arising from the molecular architecture are integrated in the density functional model. Using the theoretical model, the predicted liquid and solid densities of CO{sub 2} at equilibrium triple point are in good agreement with the experimental values. Based on the structure of crystal-liquid interfaces in different planes, the corresponding interfacial tensions are predicted. Their respective accuracies need to be tested.
Accurate density functional thermochemistry for larger molecules.
Raghavachari, K.; Stefanov, B. B.; Curtiss, L. A.; Lucent Tech.
1997-06-20
Density functional methods are combined with isodesmic bond separation reaction energies to yield accurate thermochemistry for larger molecules. Seven different density functionals are assessed for the evaluation of heats of formation, Delta H 0 (298 K), for a test set of 40 molecules composed of H, C, O and N. The use of bond separation energies results in a dramatic improvement in the accuracy of all the density functionals. The B3-LYP functional has the smallest mean absolute deviation from experiment (1.5 kcal mol/f).
Locality of correlation in density functional theory
NASA Astrophysics Data System (ADS)
Burke, Kieron; Cancio, Antonio; Gould, Tim; Pittalis, Stefano
2016-08-01
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 EC → -AC ZlnZ + BCZ as Z → ∞, where Z is the atomic number, AC is known, and we estimate BC to be about 37 mhartree. The local density approximation yields AC exactly, but a very incorrect value for BC, 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 BC a functional of the TF density for the system. The implications for the construction of approximate density functionals are discussed.
Locality of correlation in density functional theory.
Burke, Kieron; Cancio, Antonio; Gould, Tim; Pittalis, Stefano
2016-08-01
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 EC → -AC ZlnZ + BCZ as Z → ∞, where Z is the atomic number, AC is known, and we estimate BC to be about 37 mhartree. The local density approximation yields AC exactly, but a very incorrect value for BC, 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 BC a functional of the TF density for the system. The implications for the construction of approximate density functionals are discussed. PMID:27497544
Particle conservation in dynamical density functional theory
NASA Astrophysics Data System (ADS)
de las Heras, Daniel; Brader, Joseph M.; Fortini, Andrea; Schmidt, Matthias
2016-06-01
We present the exact adiabatic theory for the dynamics of the inhomogeneous density distribution of a classical fluid. Erroneous particle number fluctuations of dynamical density functional theory are absent, both for canonical and grand canonical initial conditions. We obtain the canonical free energy functional, which yields the adiabatic interparticle forces of overdamped Brownian motion. Using an exact and one of the most advanced approximate hard core free energy functionals, we obtain excellent agreement with simulations. The theory applies to finite systems in and out of equilibrium.
Particle conservation in dynamical density functional theory.
de Las Heras, Daniel; Brader, Joseph M; Fortini, Andrea; Schmidt, Matthias
2016-06-22
We present the exact adiabatic theory for the dynamics of the inhomogeneous density distribution of a classical fluid. Erroneous particle number fluctuations of dynamical density functional theory are absent, both for canonical and grand canonical initial conditions. We obtain the canonical free energy functional, which yields the adiabatic interparticle forces of overdamped Brownian motion. Using an exact and one of the most advanced approximate hard core free energy functionals, we obtain excellent agreement with simulations. The theory applies to finite systems in and out of equilibrium. PMID:27115673
Nonlocal kinetic-energy-density functionals
Garcia-Gonzalez, P.; Alvarellos, J.E.; Chacon, E. |
1996-04-01
In this paper we present nonlocal kinetic-energy functionals {ital T}[{ital n}] within the average density approximation (ADA) framework, which do not require any extra input when applied to any electron system and recover the exact kinetic energy and the linear response function of a homogeneous system. In contrast with previous ADA functionals, these present good behavior of the long-range tail of the exact weight function. The averaging procedure for the kinetic functional (averaging the Fermi momentum of the electron gas, instead of averaging the electron density) leads to a functional without numerical difficulties in the calculation of extended systems, and it gives excellent results when applied to atoms and jellium surfaces. {copyright} {ital 1996 The American Physical Society.}
Density Functionals with Broad Applicability in Chemistry
Zhao, Yan; Truhlar, Donald G.
2008-02-01
The research described in this product was performed in part in the Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. Although density functional theory is widely used in the computational chemistry community, the most popular density functional, B3LYP, has some serious shortcomings: (i) it is better for main-group chemistry than for transition metals; (ii) it systematically underestimates reaction barrier heights; (iii) it is inaccurate for interactions dominated by mediumrange correlation energy, such as van der Waals attraction, aromatic-aromatic stacking, and alkane isomerization energies. We have developed a variety of databases for testing and designing new density functionals. We used these data to design new density functionals, called M06-class (and, earlier, M05-class) functionals, for which we enforced some fundamental exact constraints such as the uniform-electron-gas limit and the absence of self-correlation energy. Our M06-class functionals depend on spin-up and spin-down electron densities (i.e., spin densities), spin density gradients, spin kinetic energy densities, and, for nonlocal (also called hybrid) functionals, Hartree-Fock exchange. We have developed four new functionals that overcome the above-mentioned difficulties: (a) M06, a hybrid meta functional, is a functional with good accuracy “across-theboard” for transition metals, main group thermochemistry, medium-range correlation energy, and barrier heights; (b) M06- 2X, another hybrid meta functional, is not good for transition metals but has excellent performance for main group chemistry, predicts accurate valence and Rydberg electronic excitation energies, and is an excellent functional for aromatic-aromatic stacking interactions; (c) M06-L is not as accurate as M06 for barrier heights but is the most accurate
Periodic subsystem density-functional theory
Genova, Alessandro; Pavanello, Michele; Ceresoli, Davide
2014-11-07
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.
Connection formulas for thermal density functional theory
NASA Astrophysics Data System (ADS)
Pribram-Jones, A.; Burke, K.
2016-05-01
The adiabatic connection formula of ground-state density functional theory relates the correlation energy to a coupling-constant integral over a purely potential contribution, and is widely used to understand and improve approximations. The corresponding formula for thermal density functional theory is cast as an integral over temperatures instead, ranging upward from the system's physical temperature. We also show how to relate different correlation components to each other, either in terms of temperature or coupling-constant integrations. We illustrate our results on the uniform electron gas.
Koopmans' condition for density-functional theory
Dabo, Ismaila; Ferretti, Andrea; Poilvert, Nicolas; Marzari, Nicola; Li, Yanli; Cococcioni, Matteo
2010-09-15
In approximate Kohn-Sham density-functional theory, self-interaction manifests itself as the dependence of the energy of an orbital on its fractional occupation. This unphysical behavior translates into qualitative and quantitative errors that pervade many fundamental aspects of density-functional predictions. Here, we first examine self-interaction in terms of the discrepancy between total and partial electron removal energies, and then highlight the importance of imposing the generalized Koopmans' condition - that identifies orbital energies as opposite total electron removal energies - to resolve this discrepancy. In the process, we derive a correction to approximate functionals that, in the frozen-orbital approximation, eliminates the unphysical occupation dependence of orbital energies up to the third order in the single-particle densities. This non-Koopmans correction brings physical meaning to single-particle energies; when applied to common local or semilocal density functionals it provides results that are in excellent agreement with experimental data - with an accuracy comparable to that of GW many-body perturbation theory - while providing an explicit total energy functional that preserves or improves on the description of established structural properties.
Generalized local-spin-density-functional theory
NASA Astrophysics Data System (ADS)
Guo, Yufei; Whitehead, M. A.
1991-01-01
An alternative density-functional theory, the generalized local-spin-density-functional (G-LSD) theory, is proposed based on the boundary conditions and sum rule of the Fermi-correlation factor in the Hartree-Fock (HF) limit. It avoids the physical restrictions to the boundary conditions and the sum rule used in the generalized exchange local-spin-density-functional (GX-LSD) theory completely, the homogeneous electron-density approximation in the Hartree-Fock-Slater (HFS) theory and in the Gáspár-Kohn-Sham (GKS) theory partly, and the time-consuming step to search for the optimal exchange parameter for each atom or ion in the Xα and Ξa theories. The alternative G-LSD theory generates the GX-LSD, HFS, GKS, Xα, and Ξa theories, when additional approximations or certain Fermi-hole shapes or high electron-density limit in a system are used. Theoretically, the G-LSD theory is more rigorous than the GX-LSD, HFS, GKS, and Ξa theories. Numerically, the statistical total energies for atoms in the G-LSD theory are in excellent agreement with the HF results, when the Gopinathan, Whitehead, and Bogdanovic [Phys. Rev. A 14, 1 (1976)] Fermi-hole parameters are used.
Density-functional expansion methods: Grand challenges
Giese, Timothy J.; York, Darrin M.
2016-01-01
We discuss the source of errors in semiempirical density functional expansion (VE) methods. In particular, we show that VE methods are capable of well-reproducing their standard Kohn-Sham density functional method counterparts, but suffer from large errors upon using one or more of these approximations: the limited size of the atomic orbital basis, the Slater monopole auxiliary basis description of the response density, and the one- and two-body treatment of the core-Hamiltonian matrix elements. In the process of discussing these approximations and highlighting their symptoms, we introduce a new model that supplements the second-order density-functional tight-binding model with a self-consistent charge-dependent chemical potential equalization correction; we review our recently reported method for generalizing the auxiliary basis description of the atomic orbital response density; and we decompose the first-order potential into a summation of additive atomic components and many-body corrections, and from this examination, we provide new insights and preliminary results that motivate and inspire new approximate treatments of the core-Hamiltonian. PMID:27293378
Communication: Embedded fragment stochastic density functional theory
Neuhauser, Daniel; Baer, Roi; Rabani, Eran
2014-07-28
We develop a method in which the electronic densities of small fragments determined by Kohn-Sham density functional theory (DFT) are embedded using stochastic DFT to form the exact density of the full system. The new method preserves the scaling and the simplicity of the stochastic DFT but cures the slow convergence that occurs when weakly coupled subsystems are treated. It overcomes the spurious charge fluctuations that impair the applications of the original stochastic DFT approach. We demonstrate the new approach on a fullerene dimer and on clusters of water molecules and show that the density of states and the total energy can be accurately described with a relatively small number of stochastic orbitals.
Density functional theory studies of etoricoxib
NASA Astrophysics Data System (ADS)
Sachdeva, Ritika; Kaur, Prabhjot; Singh, V. P.; Saini, G. S. S.
2016-05-01
Etoricoxib is a COX-2 selective inhibitor drug with molecular formula C18H15ClN2O2S. It is primarily used for the treatment of arthritis(rheumatoid, psoriatic, osteoarthritis), ankylosing spondylitis, gout and chronic low back pain. Theoretical studies of the molecule including geometry optimization and vibrational frequency calculations were carried out with the help of density functional theory calculations using 6-311++ g (d, p) basis set and B3LYP functional.
Nuclear moments in covariant density functional theory
NASA Astrophysics Data System (ADS)
Meng, J.; Zhao, P. W.; Zhang, S. Q.; Hu, J. N.; Li, J.
2014-05-01
Recent progresses on microscopic and self-consistent description of the nuclear moments in covariant density functional theory based on a point-coupling interaction are briefly reviewed. In particular, the electric quadrupole moments of Cd isotopes and the magnetic moments of Pb isotopes are discussed.
Probability density function learning by unsupervised neurons.
Fiori, S
2001-10-01
In a recent work, we introduced the concept of pseudo-polynomial adaptive activation function neuron (FAN) and presented an unsupervised information-theoretic learning theory for such structure. The learning model is based on entropy optimization and provides a way of learning probability distributions from incomplete data. The aim of the present paper is to illustrate some theoretical features of the FAN neuron, to extend its learning theory to asymmetrical density function approximation, and to provide an analytical and numerical comparison with other known density function estimation methods, with special emphasis to the universal approximation ability. The paper also provides a survey of PDF learning from incomplete data, as well as results of several experiments performed on real-world problems and signals. PMID:11709808
Differentiability of Lieb functional in electronic density functional theory
NASA Astrophysics Data System (ADS)
Lammert, Paul E.
A solid understanding of the Lieb functional FL is important because of its centrality in the foundations of electronic density functional theory. A basic question is whether directional derivatives of FL at an ensemble-V-representable density are given by (minus) the potential. A widely accepted purported proof that FL is Gâteaux differentiable at EV-representable densities would say, ?yes.? But that proof is fallacious, as shown here. FL is not Gâteaux differentiable in the normal sense, nor is it continuous. By means of a constructive approach, however, we are able to show that the derivative of FL at an EV-representable density ?0 in the direction of ?1 is given by the potential if ?0 and ?1 are everywhere strictly greater than zero, and they and the ground state wave function have square integrable derivatives through second order.
Density gradient expansion of correlation functions
NASA Astrophysics Data System (ADS)
van Leeuwen, Robert
2013-04-01
We present a general scheme based on nonlinear response theory to calculate the expansion of correlation functions such as the pair-correlation function or the exchange-correlation hole of an inhomogeneous many-particle system in terms of density derivatives of arbitrary order. We further derive a consistency condition that is necessary for the existence of the gradient expansion. This condition is used to carry out an infinite summation of terms involving response functions up to infinite order from which it follows that the coefficient functions of the gradient expansion can be expressed in terms of the local density profile rather than the background density around which the expansion is carried out. We apply the method to the calculation of the gradient expansion of the one-particle density matrix to second order in the density gradients and recover in an alternative manner the result of Gross and Dreizler [Gross and Dreizler, Z. Phys. AZPAADB0340-219310.1007/BF01413038 302, 103 (1981)], which was derived using the Kirzhnits method. The nonlinear response method is more general and avoids the turning point problem of the Kirzhnits expansion. We further give a description of the exchange hole in momentum space and confirm the wave vector analysis of Langreth and Perdew [Langreth and Perdew, Phys. Rev. BPRBMDO1098-012110.1103/PhysRevB.21.5469 21, 5469 (1980)] for this case. This is used to derive that the second-order gradient expansion of the system averaged exchange hole satisfies the hole sum rule and to calculate the gradient coefficient of the exchange energy without the need to regularize divergent integrals.
Density-Functional Theory of Thermal Transport
NASA Astrophysics Data System (ADS)
Eich, F. G.; Principi, A.; di Ventra, M.; Vignale, G.
2014-03-01
We have recently introduced a non-equilibrium density-functional theory of local temperature and associated energy density that is suitable for the study of thermoelectric phenomena from first principles. This theory rests on a local temperature field coupled to the energy-density operator. Here we apply the theory to a simple two-terminal setup, in which the terminals are held at different temperatures. We show that our treatment becomes equivalent to the standard Landauer-Büttiker formulation of thermal transport in the non-interacting limit. We gratefully acknowledge support from DOE under Grant No. DE-FG02-05ER46203 (FGE, AP, GV) and DE-FG02-05ER46204 (MD).
NASA Astrophysics Data System (ADS)
Pribram-Jones, Aurora
Warm dense matter (WDM) is a high energy phase between solids and plasmas, with characteristics of both. It is present in the centers of giant planets, within the earth's core, and on the path to ignition of inertial confinement fusion. The high temperatures and pressures of warm dense matter lead to complications in its simulation, as both classical and quantum effects must be included. One of the most successful simulation methods is density functional theory-molecular dynamics (DFT-MD). Despite great success in a diverse array of applications, DFT-MD remains computationally expensive and it neglects the explicit temperature dependence of electron-electron interactions known to exist within exact DFT. Finite-temperature density functional theory (FT DFT) is an extension of the wildly successful ground-state DFT formalism via thermal ensembles, broadening its quantum mechanical treatment of electrons to include systems at non-zero temperatures. Exact mathematical conditions have been used to predict the behavior of approximations in limiting conditions and to connect FT DFT to the ground-state theory. An introduction to FT DFT is given within the context of ensemble DFT and the larger field of DFT is discussed for context. Ensemble DFT is used to describe ensembles of ground-state and excited systems. Exact conditions in ensemble DFT and the performance of approximations depend on ensemble weights. Using an inversion method, exact Kohn-Sham ensemble potentials are found and compared to approximations. The symmetry eigenstate Hartree-exchange approximation is in good agreement with exact calculations because of its inclusion of an ensemble derivative discontinuity. Since ensemble weights in FT DFT are temperature-dependent Fermi weights, this insight may help develop approximations well-suited to both ground-state and FT DFT. A novel, highly efficient approach to free energy calculations, finite-temperature potential functional theory, is derived, which has the
Jankowski, K; Nowakowski, K; Grabowski, I; Wasilewski, J
2009-04-28
The problem of linking the dynamic electron correlation effects defined in traditional ab initio methods [or wave function theories (WFTs)] with the structure of the individual density functional theory (DFT) exchange and correlation functionals has been analyzed for the Ne atom, for which nondynamic correlation effects play a negligible role. A density-based approach directly hinged on difference radial-density (DRD) distributions defined with respect the Hartree-Fock radial density has been employed for analyzing the impact of dynamic correlation effects on the density. Attention has been paid to the elimination of basis-set incompleteness errors. The DRD distributions calculated by several ab initio methods have been compared to their DFT counterparts generated for representatives of several generations of broadly used exchange-correlation functionals and for the recently developed orbital-dependent OEP2 exchange-correlation functional [Bartlett et al., J. Chem. Phys. 122, 034104 (2005)]. For the local, generalized-gradient, and hybrid functionals it has been found that the dynamic correlation effects are to a large extend accounted for by densities resulting from exchange-only calculations. Additional calculations with self-interaction corrected exchange potentials indicate that this finding cannot be explained as an artifact caused by the self-interaction error. It has been demonstrated that the VWN5 and LYP correlation functionals do not represent any substantial dynamical correlation effects on the electron density, whereas these effects are well represented by the orbital-dependent OEP2 correlation functional. Critical comparison of the present results with their counterparts reported in literature has been made. Some attention has been paid to demonstrating the differences between the energy- and density-based perspectives. They indicate the usefulness of density-based criteria for developing new exchange-correlation functionals. PMID:19405556
NASA Astrophysics Data System (ADS)
Jankowski, K.; Nowakowski, K.; Grabowski, I.; Wasilewski, J.
2009-04-01
The problem of linking the dynamic electron correlation effects defined in traditional ab initio methods [or wave function theories (WFTs)] with the structure of the individual density functional theory (DFT) exchange and correlation functionals has been analyzed for the Ne atom, for which nondynamic correlation effects play a negligible role. A density-based approach directly hinged on difference radial-density (DRD) distributions defined with respect the Hartree-Fock radial density has been employed for analyzing the impact of dynamic correlation effects on the density. Attention has been paid to the elimination of basis-set incompleteness errors. The DRD distributions calculated by several ab initio methods have been compared to their DFT counterparts generated for representatives of several generations of broadly used exchange-correlation functionals and for the recently developed orbital-dependent OEP2 exchange-correlation functional [Bartlett et al., J. Chem. Phys. 122, 034104 (2005)]. For the local, generalized-gradient, and hybrid functionals it has been found that the dynamic correlation effects are to a large extend accounted for by densities resulting from exchange-only calculations. Additional calculations with self-interaction corrected exchange potentials indicate that this finding cannot be explained as an artifact caused by the self-interaction error. It has been demonstrated that the VWN5 and LYP correlation functionals do not represent any substantial dynamical correlation effects on the electron density, whereas these effects are well represented by the orbital-dependent OEP2 correlation functional. Critical comparison of the present results with their counterparts reported in literature has been made. Some attention has been paid to demonstrating the differences between the energy- and density-based perspectives. They indicate the usefulness of density-based criteria for developing new exchange-correlation functionals.
Density functional calculations on model tyrosyl radicals.
Himo, F; Gräslund, A; Eriksson, L A
1997-01-01
A gradient-corrected density functional theory approach (PWP86) has been applied, together with large basis sets (IGLO-III), to investigate the structure and hyperfine properties of model tyrosyl free radicals. In nature, these radicals are observed in, e.g., the charge transfer pathways in photosystem II (PSII) and in ribonucleotide reductases (RNRs). By comparing spin density distributions and proton hyperfine couplings with experimental data, it is confirmed that the tyrosyl radicals present in the proteins are neutral. It is shown that hydrogen bonding to the phenoxyl oxygen atom, when present, causes a reduction in spin density on O and a corresponding increase on C4. Calculated proton hyperfine coupling constants for the beta-protons show that the alpha-carbon is rotated 75-80 degrees out of the plane of the ring in PSII and Salmonella typhimurium RNR, but only 20-30 degrees in, e.g., Escherichia coli, mouse, herpes simplex, and bacteriophage T4-induced RNRs. Furthermore, based on the present calculations, we have revised the empirical parameters used in the experimental determination of the oxygen spin density in the tyrosyl radical in E. coli RNR and of the ring carbon spin densities, from measured hyperfine coupling constants. Images FIGURE 1 FIGURE 5 PMID:9083661
Orbital-optimized density cumulant functional theory
Sokolov, Alexander Yu. Schaefer, Henry F.
2013-11-28
In density cumulant functional theory (DCFT) the electronic energy is evaluated from the one-particle density matrix and two-particle density cumulant, circumventing the computation of the wavefunction. To achieve this, the one-particle density matrix is decomposed exactly into the mean-field (idempotent) and correlation components. While the latter can be entirely derived from the density cumulant, the former must be obtained by choosing a specific set of orbitals. In the original DCFT formulation [W. Kutzelnigg, J. Chem. Phys. 125, 171101 (2006)] the orbitals were determined by diagonalizing the effective Fock operator, which introduces partial orbital relaxation. Here we present a new orbital-optimized formulation of DCFT where the energy is variationally minimized with respect to orbital rotations. This introduces important energy contributions and significantly improves the description of the dynamic correlation. In addition, it greatly simplifies the computation of analytic gradients, for which expressions are also presented. We offer a perturbative analysis of the new orbital stationarity conditions and benchmark their performance for a variety of chemical systems.
Density Functional Theory Models for Radiation Damage
NASA Astrophysics Data System (ADS)
Dudarev, S. L.
2013-07-01
Density functional theory models developed over the past decade provide unique information about the structure of nanoscale defects produced by irradiation and about the nature of short-range interaction between radiation defects, clustering of defects, and their migration pathways. These ab initio models, involving no experimental input parameters, appear to be as quantitatively accurate and informative as the most advanced experimental techniques developed for the observation of radiation damage phenomena. Density functional theory models have effectively created a new paradigm for the scientific investigation and assessment of radiation damage effects, offering new insight into the origin of temperature- and dose-dependent response of materials to irradiation, a problem of pivotal significance for applications.
Teaching Density Functional Theory Through Experiential Learning
NASA Astrophysics Data System (ADS)
Narasimhan, Shobhana
2015-09-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.
The problem of the universal density functional and the density matrix functional theory
Bobrov, V. B. Trigger, S. A.
2013-04-15
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.
Properties of the cosmological density distribution function
NASA Astrophysics Data System (ADS)
Bernardeau, Francis; Kofman, Lev
1995-04-01
The properties of the probability distribution function (PDF) of the cosmological continuous density field are studied. We focus our analysis on the quasi-linear regime where various calculations, based on dynamically motivated methods, have been presented: either by using the Zel'dovich approximation (ZA) or by using the perturbation theory to evaluate the behavior of the moments of the distribution function. We show how these two approaches are related to each other and that they can be used in a complementary way. For that respect, the one-dimensional dynamics, where the ZA is exact solution, has first been used as a testing ground. In particular, we show that, when the density PDF obtained with the ZA is regularized, its various moments exhibit the behavior expected by the perturbation theory applied to the ZA. We show that ZA approach can be used for arbitrary initial conditions (not only Gaussian) and that the nonlinear evolution of the moments can be obtained. The perturbation theory can be used for the exact dynamics. We take into account the final filtering of the density field both for ZA and perturbation theory. Applying these techniques, we got the generating function of the moments for the one-dimensional dynamics, the three-dimensional ZA, with and without smoothing effects. We also suggest methods to build PDFs. One is based on the Laplace inverse transform of the moment generating function. The other, the Edgeworth expansion, is obtained when the previous generating function is truncated at a given order and allows evaluation of the PDF out of limited number of moments. It provides insight on the relationship between the moments and the shape of the density PDF. In particular, it provides an alternative method to evaluate the skewness and kurtosis by measuring the PDF around its maximum. Eventually, results obtained from a numerical simulation with cold dark matter initial conditions have been used to validate the accuracy of the considered
Magnetic susceptibility of semiconductors by an all-electron first-principles approach
Ohno, K. |; Mauri, F.; Louie, S.G. |
1997-07-01
The magnetic susceptibility ({chi}) of the semiconductors (diamond, Si, GaAs, and GaP) and of the inert-gas solids (Ne, Ar, and Kr) are evaluated within density-functional theory in the local-density approximation, using a mixed-basis all-electron approach. In Si, GaAs, GaP, Ar, and Kr, the contribution of core electrons to {chi} is comparable to that of valence electrons. However, our results show that the contribution associated with the core states is independent of the chemical environment and can be computed from the isolated atoms. Moreover, our results indicate that the use of a {open_quotes}scissor operator{close_quotes} does not improve the agreement of the theoretical {chi} with experiments. {copyright} {ital 1997} {ital The American Physical Society}
Pseudospectral time-dependent density functional theory
NASA Astrophysics Data System (ADS)
Ko, Chaehyuk; Malick, David K.; Braden, Dale A.; Friesner, Richard A.; Martínez, Todd J.
2008-03-01
Time-dependent density functional theory (TDDFT) is implemented within the Tamm-Dancoff approximation (TDA) using a pseudospectral approach to evaluate two-electron repulsion integrals. The pseudospectral approximation uses a split representation with both spectral basis functions and a physical space grid to achieve a reduction in the scaling behavior of electronic structure methods. We demonstrate here that exceptionally sparse grids may be used in the excitation energy calculation, following earlier work employing the pseudospectral approximation for determining correlation energies in wavefunction-based methods with similar conclusions. The pseudospectral TDA-TDDFT method is shown to be up to ten times faster than a conventional algorithm for hybrid functionals without sacrificing chemical accuracy.
Optimization of constrained density functional theory
NASA Astrophysics Data System (ADS)
O'Regan, David D.; Teobaldi, Gilberto
2016-07-01
Constrained density functional theory (cDFT) is a versatile electronic structure method that enables ground-state calculations to be performed subject to physical constraints. It thereby broadens their applicability and utility. Automated Lagrange multiplier optimization is necessary for multiple constraints to be applied efficiently in cDFT, for it to be used in tandem with geometry optimization, or with molecular dynamics. In order to facilitate this, we comprehensively develop the connection between cDFT energy derivatives and response functions, providing a rigorous assessment of the uniqueness and character of cDFT stationary points while accounting for electronic interactions and screening. In particular, we provide a nonperturbative proof that stable stationary points of linear density constraints occur only at energy maxima with respect to their Lagrange multipliers. We show that multiple solutions, hysteresis, and energy discontinuities may occur in cDFT. Expressions are derived, in terms of convenient by-products of cDFT optimization, for quantities such as the dielectric function and a condition number quantifying ill definition in multiple constraint cDFT.
Recent progress in density functional theory
NASA Astrophysics Data System (ADS)
Truhlar, Donald
2014-03-01
Ongoing work involves several areas of density functional theory: new methods for computing electronic excitation energies, including a new way to remove spin contamination in the spin-flip Tamm-Dancoff approximation and a configuration-interaction-corrected Tamm-Dancoff Approximation for treating conical intersections; new ways to treat open-shell states, including a reinterpreted broken-symmetry method and multi-configuration Kohn-Sham theory; a new exchange-correlation functional; new tests of density functional theory against databases for electronic transition energies and molecules and solids containing metal atoms; and applications. A selection of results will be presented. I am grateful to the following collaborators for contributions to the ongoing work: Boris Averkiev, Rebecca Carlson, Laura Fernandez, Laura Gagliardi, Chad Hoyer, Francesc Illas, Miho Isegawa, Shaohong Li, Giovanni Li Manni, Sijie Luo, Dongxia Ma, Remi Maurice, Rubén Means-Pañeda, Roberto Peverati, Nora Planas, Prasenjit Seal, Pragya Verma, Bo Wang, Xuefei Xu, Ke R. Yang, Haoyu Yu, Wenjing Zhang, and Jingjing Zheng. Supported in part by the AFOSR and U.S. DOE.
Modulation Based on Probability Density Functions
NASA Technical Reports Server (NTRS)
Williams, Glenn L.
2009-01-01
A proposed method of modulating a sinusoidal carrier signal to convey digital information involves the use of histograms representing probability density functions (PDFs) that characterize samples of the signal waveform. The method is based partly on the observation that when a waveform is sampled (whether by analog or digital means) over a time interval at least as long as one half cycle of the waveform, the samples can be sorted by frequency of occurrence, thereby constructing a histogram representing a PDF of the waveform during that time interval.
Adiabatic corrections to density functional theory energies and wave functions.
Mohallem, José R; Coura, Thiago de O; Diniz, Leonardo G; de Castro, Gustavo; Assafrão, Denise; Heine, Thomas
2008-09-25
The adiabatic finite-nuclear-mass-correction (FNMC) to the electronic energies and wave functions of atoms and molecules is formulated for density-functional theory and implemented in the deMon code. The approach is tested for a series of local and gradient corrected density functionals, using MP2 results and diagonal-Born-Oppenheimer corrections from the literature for comparison. In the evaluation of absolute energy corrections of nonorganic molecules the LDA PZ81 functional works surprisingly better than the others. For organic molecules the GGA BLYP functional has the best performance. FNMC with GGA functionals, mainly BLYP, show a good performance in the evaluation of relative corrections, except for nonorganic molecules containing H atoms. The PW86 functional stands out with the best evaluation of the barrier of linearity of H2O and the isotopic dipole moment of HDO. In general, DFT functionals display an accuracy superior than the common belief and because the corrections are based on a change of the electronic kinetic energy they are here ranked in a new appropriate way. The approach is applied to obtain the adiabatic correction for full atomization of alcanes C(n)H(2n+2), n = 4-10. The barrier of 1 mHartree is approached for adiabatic corrections, justifying its insertion into DFT. PMID:18537228
NASA Astrophysics Data System (ADS)
Sicilia, E.; de Luca, G.; Chiodo, S.; Russo, N.; Calaminici, P.; Koster, A. M.; Jug, K.
Density functional calculations of the electric field gradient tensor at the nitrogen nucleus in 13 test molecules, containing 14 nitrogen sites, have been performed using the linear combination of Gaussian-type orbital Kohn-Sham density functional theory (LCGTO-KSDFT) approach. Local and gradient corrected functionals were used for all-electron calculations. All the molecular structures were optimized at their respective levels of theory with extended basis sets. Calibrated 14N nuclear quadrupole moments were obtained through a fitting procedure between calculated electric field gradients and experimental nuclear quadrupole coupling constants of the test set of molecules for each basis set and functional considered. With these calibrated 14N nuclear quadrupole moments, the nuclear quadrupole coupling constants of the following selected systems were determined: fluoromethylisonitrile, pyridine, pyrrole, imadazole, pyrazole, 1,8-bis(dimethyl-amino)naphthalene, cyclotetramethylenetetranitramine, cocaine and heroin.
Gaussian basis density functional theory for systems periodic in two or three dimensions
Hess, A.C.; Lin, Zijing; Clerc, D.G.
1996-12-31
We will present the formulation and results of a new, fully self-consistent all electron density functional theory that utilizes localized Gaussian basis functions to treat systems possessing 3- and 2- dimensional periodicity. Our method extends, in a fully rigorous fashion, many of the concepts commonly found in molecular density functional theories and can accurately be described as an {open_quotes}ab initio or first principles{close_quotes} approach. Results will be presented for several illustrative spin paired and unpaired insulating and semi-conducting bulk crystalline systems. The effect of employing local and non-local exchange and correlation functionals on the geometric, elastic, and electronic properties of diamond, silicon, silicon carbide, magnesium oxide and beryllium oxide will also be presented. Finally, the geometric and electronic structure of several oxide and semi-conducting surfaces will also be presented.
Density-functional theory of superconductivity
NASA Astrophysics Data System (ADS)
Gross, E. K. U.
2008-03-01
A prominent challenge of modern condensed-matter theory is to predict reliably material-specific properties of superconductors, such as the critical temperature. The traditional model of Bardeen, Cooper and Schrieffer (BCS) properly describes the universal features that all conventional superconductors have in common, but it is not able to make accurate predictions of material-specific properties. To tackle this problem, a density-functional formalism has been developed [1] which describes superconductors in thermal equilibrium in terms of three quantities: the ordinary density, the superconducting order parameter, and the nuclear N-body density. These three ``densities'' are determined self-consistently through a set of Kohn-Sham equations. Approximations of the universal exchange-correlation functional are derived on the basis of many-body perturbation theory. In this way, a true ab-initio description is achieved which does not contain any adjustable parameters such as the μ* of Eliashberg theory. Numerical results for the critical temperature, the isotope effect, the gap function and the jump of the specific heat will be presented for simple metals, for MgB2 [2] and CaBeSi, and for calcium intercalated graphite (CaC6) [3]. Furthermore, results for Li, Al, K, and H under pressure will be discussed. The calculations explain why Li and Al behave very differently, leading to a strong enhancement of superconductivity for Li and to a clear suppression for Al with increasing pressure [4]. For K we predict a behavior similar to Li, i.e. a strong increase of Tc with increasing pressure. Finally, hydrogen is found to be a three-gap superconductor whose critical temperature increases with increasing pressure until about 100K (at 500 GPa). [1] M. Lüders, M.A.L. Marques, N.N. Lathiotakis, A. Floris,G. Profeta, L. Fast, A.Continenza, S. Massidda, E.K.U. Gross, PRB 72, 024545 (2005). [2] A. Floris, G. Profeta, N.N. Lathiotakis, M. Lüders, M.A.L. Marques, C. Franchini, E
Local spin analyses using density functional theory
NASA Astrophysics Data System (ADS)
Abate, Bayileyegn; Peralta, Juan
Local spin analysis is a valuable technique in computational investigations magnetic interactions on mono- and polynuclear transition metal complexes, which play vital roles in catalysis, molecular magnetism, artificial photosynthesis, and several other commercially important materials. The relative size and complex electronic structure of transition metal complexes often prohibits the use of multi-determinant approaches, and hence, practical calculations are often limited to single-determinant methods. Density functional theory (DFT) has become one of the most successful and widely used computational tools for the electronic structure study of complex chemical systems; transition metal complexes in particular. Within the DFT formalism, a more flexible and complete theoretical modeling of transition metal complexes can be achieved by considering noncollinear spins, in which the spin density is 'allowed to' adopt noncollinear structures in stead of being constrained to align parallel/antiparallel to a universal axis of magnetization. In this meeting, I will present local spin analyses results obtained using different DFT functionals. Local projection operators are used to decompose the expectation value
Chemistry by Way of Density Functional Theory
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Ricca, Alessandra; Partridge, Harry; Langohff, Stephen R.; Arnold, James O. (Technical Monitor)
1996-01-01
In this work we demonstrate that density functional theory (DFT) methods make an important contribution to understanding chemical systems and are an important additional method for the computational chemist. We report calibration calculations obtained with different functionals for the 55 G2 molecules to justify our selection of the B3LYP functional. We show that accurate geometries and vibrational frequencies obtained at the B3LYP level can be combined with traditional methods to simplify the calculation of accurate heats of formation. We illustrate the application of the B3LYP approach to a variety of chemical problems from the vibrational frequencies of polycyclic aromatic hydrocarbons to transition metal systems. We show that the B3LYP method typically performs better than the MP2 method at a significantly lower computational cost. Thus the B3LYP method allows us to extend our studies to much larger systems while maintaining a high degree of accuracy. We show that for transition metal systems, the B3LYP bond energies are typically of sufficient accuracy that they can be used to explain experimental trends and even differentiate between different experimental values. We show that for boron clusters the B3LYP energetics are not as good as for many of the other systems presented, but even in this case the B3LYP approach is able to help understand the experimental trends.
Density functional studies of representative pericyclic reactions
NASA Astrophysics Data System (ADS)
Carpenter, John E.; Sosa, Carlos P.
1994-07-01
Density functional theory (DFT) has traditionally been shunned by computational chemists, but has long seen widespread use in the physics community. Recently, however, DFT has been adopted by the ab initio quantum chemistry community and much activity has been devoted to refining the methodology and exploring the range of its applicability. We investigate the use of DFT (both local (LDF) and non-local (NLDF) spin density approximations) to calculate transition and equilibrium structures for three representative pericyclic reactions: the electrocyclic ring opening reaction of cyclobutene, the [1,5] sigmatropic hydrogen shift reaction in (Z)-1,3-pentadiene, and the Diels-Alder cycloaddition reaction between ethylene and butadiene. LDF theory tends to overemphasize the stability of the ringed structures in each of these reactions. For example, LDF predicts a very low (6 kcal mol -1) barrier to reaction for the Deils-Alder reaction. NLDF theory substantially improves the calculated reaction barrier (20 kcal mol -1), but it is still low with respect to experiment.
Semiconductor Thermochemistry in Density Functional Calculations
Lany, S.
2008-01-01
The local-density and generalized gradient approximations (LDA and GGA) to density functional theory (DFT) exhibit incomplete error cancellation when energy differences are taken between chemically dissimilar systems. This energy inconsistency is manifested, e.g., in the tendency to underestimate the heat (enthalpy) of formation of semiconducting and insulating compounds in LDA and, even more so, in GGA. Considering a set of 61 compounds that can be formed from 14 elements (cations: Cu, Mg, Ca, Zn, Cd, Al, Ga, and In; anions: N, P, As, O, S, and Se), optimized elemental reference energies are determined by least-squares error minimization of an overdetermined set of linear equations. These elemental energies are 'optimally consistent' with the DFT energies of the semiconductor compounds and imply corrections of up to 1 eV compared to the respective LDA or GGA energies. While these 'corrections' are not to be understood to yield the correct absolute total energies of the elements, they are proposed to give appropriate bounds for the chemical potentials for thermodynamic processes in semiconductors and insulators, such as, e.g., defect formation, surface reconstruction, or catalytic processes. The present model allows to evaluate thermodynamic processes using DFT energy differences taken only between systems that are expected to show good error cancellation.
Carrier Modulation Via Waveform Probability Density Function
NASA Technical Reports Server (NTRS)
Williams, Glenn L.
2004-01-01
Beyond the classic modes of carrier modulation by varying amplitude (AM), phase (PM), or frequency (FM), we extend the modulation domain of an analog carrier signal to include a class of general modulations which are distinguished by their probability density function histogram. Separate waveform states are easily created by varying the pdf of the transmitted waveform. Individual waveform states are assignable as proxies for digital ONEs or ZEROs. At the receiver, these states are easily detected by accumulating sampled waveform statistics and performing periodic pattern matching, correlation, or statistical filtering. No fundamental natural laws are broken in the detection process. We show how a typical modulation scheme would work in the digital domain and suggest how to build an analog version. We propose that clever variations of the modulating waveform (and thus the histogram) can provide simple steganographic encoding.
Carrier Modulation Via Waveform Probability Density Function
NASA Technical Reports Server (NTRS)
Williams, Glenn L.
2006-01-01
Beyond the classic modes of carrier modulation by varying amplitude (AM), phase (PM), or frequency (FM), we extend the modulation domain of an analog carrier signal to include a class of general modulations which are distinguished by their probability density function histogram. Separate waveform states are easily created by varying the pdf of the transmitted waveform. Individual waveform states are assignable as proxies for digital one's or zero's. At the receiver, these states are easily detected by accumulating sampled waveform statistics and performing periodic pattern matching, correlation, or statistical filtering. No fundamental physical laws are broken in the detection process. We show how a typical modulation scheme would work in the digital domain and suggest how to build an analog version. We propose that clever variations of the modulating waveform (and thus the histogram) can provide simple steganographic encoding.
Phases of Polonium via Density Functional Theory
NASA Astrophysics Data System (ADS)
Verstraete, Matthieu J.
2010-01-01
The thermodynamical properties of the main phases of metallic polonium are examined using density functional theory. The exceptional nature of the solid-solid phase transition of α to β Po is underlined: it induces a lowering in symmetry, from cubic to rhombohedral, with increasing temperature. This is explained as the result of a delicate balance between bonding and entropic effects. Overall agreement with existing experimental data is good by state-of-the-art standards. The phonons of Po present Kohn anomalies, and it is shown that the effect of spin-orbit interactions is the inverse of that in normal metals: due to the nonspherical nature of the Fermi Surface, spin-orbit effects reduce nesting and harden most phonon frequencies.
Effective potential in density matrix functional theory.
Nagy, A; Amovilli, C
2004-10-01
In the previous paper it was shown that in the ground state the diagonal of the spin independent second-order density matrix n can be determined by solving a single auxiliary equation of a two-particle problem. Thus the problem of an arbitrary system with even electrons can be reduced to a two-particle problem. The effective potential of the two-particle equation contains a term v(p) of completely kinetic origin. Virial theorem and hierarchy of equations are derived for v(p) and simple approximations are proposed. A relationship between the effective potential u(p) of the shape function equation and the potential v(p) is established. PMID:15473719
Improving Density Functionals with Quantum Harmonic Oscillators
NASA Astrophysics Data System (ADS)
Tkatchenko, Alexandre
2013-03-01
Density functional theory (DFT) is the most widely used and successful approach for electronic structure calculations. However, one of the pressing challenges for DFT is developing efficient functionals that can accurately capture the omnipresent long-range electron correlations, which determine the structure and stability of many molecules and materials. Here we show that, under certain conditions, the problem of computing the long-range correlation energy of interacting electrons can be mapped to a system of coupled quantum harmonic oscillators (QHOs). The proposed model allows us to synergistically combine concepts from DFT, quantum chemistry, and the widely discussed random-phase approximation for the correlation energy. In the dipole limit, the interaction energy for a system of coupled QHOs can be calculated exactly, thereby leading to an efficient and accurate model for the many-body dispersion energy of complex molecules and materials. The studied examples include intermolecular binding energies, the conformational hierarchy of DNA structures, the geometry and stability of molecular crystals, and supramolecular host-guest complexes (A. Tkatchenko, R. A. DiStasio Jr., R. Car, M. Scheffler, Phys. Rev. Lett. 108, 236402 (2012); R. A. DiStasio Jr., A. von Lilienfeld, A. Tkatchenko, PNAS 109, 14791 (2012); A. Tkatchenko, D. Alfe, K. S. Kim, J. Chem. Theory and Comp. (2012), doi: 10.1021/ct300711r; A. Tkatchenko, A. Ambrosetti, R. A. DiStasio Jr., arXiv:1210.8343v1).
Building a Universal Nuclear Energy Density Functional
Carlson, Joe A.; Furnstahl, Dick; Horoi, Mihai; Lust, Rusty; Nazaewicc, Witek; Ng, Esmond; Thompson, Ian; Vary, James
2012-12-30
During the period of Dec. 1 2006 – Jun. 30, 2012, the UNEDF collaboration carried out a comprehensive study of all nuclei, based on the most accurate knowledge of the strong nuclear interaction, the most reliable theoretical approaches, the most advanced algorithms, and extensive computational resources, with a view towards scaling to the petaflop platforms and beyond. The long-term vision initiated with UNEDF is to arrive at a comprehensive, quantitative, and unified description of nuclei and their reactions, grounded in the fundamental interactions between the constituent nucleons. We seek to replace current phenomenological models of nuclear structure and reactions with a well-founded microscopic theory that delivers maximum predictive power with well-quantified uncertainties. Specifically, the mission of this project has been three-fold: First, to find an optimal energy density functional (EDF) using all our knowledge of the nucleonic Hamiltonian and basic nuclear properties; Second, to apply the EDF theory and its extensions to validate the functional using all the available relevant nuclear structure and reaction data; Third, to apply the validated theory to properties of interest that cannot be measured, in particular the properties needed for reaction theory.
A density functional for sparse matter
NASA Astrophysics Data System (ADS)
Langreth, D. C.; Lundqvist, B. I.; Chakarova-Käck, S. D.; Cooper, V. R.; Dion, M.; Hyldgaard, P.; Kelkkanen, A.; Kleis, J.; Kong, Lingzhu; Li, Shen; Moses, P. G.; Murray, E.; Puzder, A.; Rydberg, H.; Schröder, E.; Thonhauser, T.
2009-02-01
Sparse matter is abundant and has both strong local bonds and weak nonbonding forces, in particular nonlocal van der Waals (vdW) forces between atoms separated by empty space. It encompasses a broad spectrum of systems, like soft matter, adsorption systems and biostructures. Density-functional theory (DFT), long since proven successful for dense matter, seems now to have come to a point, where useful extensions to sparse matter are available. In particular, a functional form, vdW-DF (Dion et al 2004 Phys. Rev. Lett. 92 246401; Thonhauser et al 2007 Phys. Rev. B 76 125112), has been proposed for the nonlocal correlations between electrons and applied to various relevant molecules and materials, including to those layered systems like graphite, boron nitride and molybdenum sulfide, to dimers of benzene, polycyclic aromatic hydrocarbons (PAHs), doped benzene, cytosine and DNA base pairs, to nonbonding forces in molecules, to adsorbed molecules, like benzene, naphthalene, phenol and adenine on graphite, alumina and metals, to polymer and carbon nanotube (CNT) crystals, and hydrogen storage in graphite and metal-organic frameworks (MOFs), and to the structure of DNA and of DNA with intercalators. Comparison with results from wavefunction calculations for the smaller systems and with experimental data for the extended ones show the vdW-DF path to be promising. This could have great ramifications.
Semiclassical origins of density functional theory
NASA Astrophysics Data System (ADS)
Burke, Kieron
2009-03-01
Until the seminal work of Hohenberg, Kohn, and Sham of the mid 60's, most density functional theory (DFT) was derived from semiclassical approximations. This non-empirical approach shows an intrinsic difference between solids (for which DFT was originally developed) and molecules, and explains many of its more mysterious manifestations. For example, the success of DFT for molecules has nothing to do with the uniform gas. Results include [1] a derivation of the empirical parameter in the B88 exchange functional, [2] PBEsol, a new GGA that restores the exchange gradient expansion and improves lattice constants in solids, [3] a novel approach to ``orbital-free'' DFT that, in preliminary tests, is 40 times more accurate than its DFT counterpart. The talk is aimed at a general theoretical audience. Detailed technical knowledge of DFT is neither needed, nor desirable. [4pt] [1] J.P. Perdew, L.A. Constantin, E. Sagvolden, and KB, Phys. Rev. Lett. 97, 223002 (2006). [0pt] [2] J.P. Perdew, A. Ruzsinszky, G.I. Csonka, O.A. Vydrov, G.E. Scuseria, L.A. Constantin, X. Zhou, and KB, Phys. Rev. Lett. 100, 136406 (2008). [0pt] [3] Peter Elliott, Donghyung Lee, Attila Cangi, KB, Phys. Rev. Lett. 100, 256406 (2008).
Rationale for switching to nonlocal functionals in density functional theory
NASA Astrophysics Data System (ADS)
Lazić, P.; Atodiresei, N.; Caciuc, V.; Brako, R.; Gumhalter, B.; Blügel, S.
2012-10-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.
Density Functional Theory study of the equilibrium density of water at normal conditions
NASA Astrophysics Data System (ADS)
Wang, Jue; Roman-Perez, Guillermo; Soler, Jose M.; Artacho, Emilio; Fernandez-Serra, Marivi
2010-03-01
Ab initio molecular dynamics of liquid water with the use of density functional theory (DFT) currently underperform experimental equilibrium density 1g/cm^3 under room temperature. At constant density, not much is known about the equilibrium density of commonly used GGA functionals in liquid water simulations. We present a DFT-based AIMD study of liquid water at different densities and analyze the structure and diffusivity of water with different exchange and correlation functionals. We show that all current GGA functionals fail to reproduce experimental density, however, the explicit description of long range correlations through a Van der Walls density functional (DRSLL)footnotetextM. Dion, H. Rydberg, E. Schr"oder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004) can potentially transform our current understanding of the structure of liquid water. Our results shows that this new functional improves density, with only 2% error to experiment. But it underperforms GGA functionals in terms of structure.
Embedded-cluster calculations in a numeric atomic orbital density-functional theory framework
Berger, Daniel Oberhofer, Harald; Reuter, Karsten; Logsdail, Andrew J. Farrow, Matthew R.; Catlow, C. Richard A.; Sokol, Alexey A.; Sherwood, Paul; Blum, Volker
2014-07-14
We integrate the all-electron electronic structure code FHI-aims into the general ChemShell package for solid-state embedding quantum and molecular mechanical (QM/MM) calculations. A major undertaking in this integration is the implementation of pseudopotential functionality into FHI-aims to describe cations at the QM/MM boundary through effective core potentials and therewith prevent spurious overpolarization of the electronic density. Based on numeric atomic orbital basis sets, FHI-aims offers particularly efficient access to exact exchange and second order perturbation theory, rendering the established QM/MM setup an ideal tool for hybrid and double-hybrid level density functional theory calculations of solid systems. We illustrate this capability by calculating the reduction potential of Fe in the Fe-substituted ZSM-5 zeolitic framework and the reaction energy profile for (photo-)catalytic water oxidation at TiO{sub 2}(110)
Flexoelectricity from density-functional perturbation theory
NASA Astrophysics Data System (ADS)
Stengel, Massimiliano
2013-11-01
We derive the complete flexoelectric tensor, including electronic and lattice-mediated effects, of an arbitrary insulator in terms of the microscopic linear response of the crystal to atomic displacements. The basic ingredient, which can be readily calculated from first principles in the framework of density-functional perturbation theory, is the quantum-mechanical probability current response to a long-wavelength acoustic phonon. Its second-order Taylor expansion in the wave vector q around the Γ (q=0) point in the Brillouin zone naturally yields the flexoelectric tensor. At order one in q we recover Martin's theory of piezoelectricity [Martin, Phys. Rev. B 5, 1607 (1972)], thus providing an alternative derivation thereof. To put our derivations on firm theoretical grounds, we perform a thorough analysis of the nonanalytic behavior of the dynamical matrix and other response functions in a vicinity of Γ. Based on this analysis, we find that there is an ambiguity in the specification of the “zero macroscopic field” condition in the flexoelectric case; such arbitrariness can be related to an analytic band-structure term, in close analogy to the theory of deformation potentials. As a by-product, we derive a rigorous generalization of the Cochran-Cowley formula [Cochran and Cowley, J. Phys. Chem. Solids 23, 447 (1962)] to higher orders in q. This can be of great utility in building reliable atomistic models of electromechanical phenomena, as well as for improving the accuracy of the calculation of phonon dispersion curves. Finally, we discuss the physical interpretation of the various contributions to the flexoelectric response, either in the static or dynamic regime, and we relate our findings to earlier theoretical works on the subject.
Bone mineral density, adiposity, and cognitive functions
Sohrabi, Hamid R.; Bates, Kristyn A.; Weinborn, Michael; Bucks, Romola S.; Rainey-Smith, Stephanie R.; Rodrigues, Mark A.; Bird, Sabine M.; Brown, Belinda M.; Beilby, John; Howard, Matthew; Criddle, Arthur; Wraith, Megan; Taddei, Kevin; Martins, Georgia; Paton, Athena; Shah, Tejal; Dhaliwal, Satvinder S.; Mehta, Pankaj D.; Foster, Jonathan K.; Martins, Ian J.; Lautenschlager, Nicola T.; Mastaglia, Francis; Laws, Simon M.; Martins, Ralph N.
2015-01-01
Cognitive decline and dementia due to Alzheimer's disease (AD) have been associated with genetic, lifestyle, and environmental factors. A number of potentially modifiable risk factors should be taken into account when preventive or ameliorative interventions targeting dementia and its preclinical stages are investigated. Bone mineral density (BMD) and body composition are two such potentially modifiable risk factors, and their association with cognitive decline was investigated in this study. 164 participants, aged 34–87 years old (62.78 ± 9.27), were recruited for this longitudinal study and underwent cognitive and clinical examinations at baseline and after 3 years. Blood samples were collected for apolipoprotein E (APOE) genotyping and dual energy x-ray absorptiometry (DXA) was conducted at the same day as cognitive assessment. Using hierarchical regression analysis, we found that BMD and lean body mass, as measured using DXA were significant predictors of episodic memory. Age, gender, APOE status, and premorbid IQ were controlled for. Specifically, the List A learning from California Verbal Learning Test was significantly associated with BMD and lean mass both at baseline and at follow up assessment. Our findings indicate that there is a significant association between BMD and lean body mass and episodic verbal learning. While the involvement of modifiable lifestyle factors in human cognitive function has been examined in different studies, there is a need for further research to understand the potential underlying mechanisms. PMID:25741279
Bone mineral density, adiposity, and cognitive functions.
Sohrabi, Hamid R; Bates, Kristyn A; Weinborn, Michael; Bucks, Romola S; Rainey-Smith, Stephanie R; Rodrigues, Mark A; Bird, Sabine M; Brown, Belinda M; Beilby, John; Howard, Matthew; Criddle, Arthur; Wraith, Megan; Taddei, Kevin; Martins, Georgia; Paton, Athena; Shah, Tejal; Dhaliwal, Satvinder S; Mehta, Pankaj D; Foster, Jonathan K; Martins, Ian J; Lautenschlager, Nicola T; Mastaglia, Francis; Laws, Simon M; Martins, Ralph N
2015-01-01
Cognitive decline and dementia due to Alzheimer's disease (AD) have been associated with genetic, lifestyle, and environmental factors. A number of potentially modifiable risk factors should be taken into account when preventive or ameliorative interventions targeting dementia and its preclinical stages are investigated. Bone mineral density (BMD) and body composition are two such potentially modifiable risk factors, and their association with cognitive decline was investigated in this study. 164 participants, aged 34-87 years old (62.78 ± 9.27), were recruited for this longitudinal study and underwent cognitive and clinical examinations at baseline and after 3 years. Blood samples were collected for apolipoprotein E (APOE) genotyping and dual energy x-ray absorptiometry (DXA) was conducted at the same day as cognitive assessment. Using hierarchical regression analysis, we found that BMD and lean body mass, as measured using DXA were significant predictors of episodic memory. Age, gender, APOE status, and premorbid IQ were controlled for. Specifically, the List A learning from California Verbal Learning Test was significantly associated with BMD and lean mass both at baseline and at follow up assessment. Our findings indicate that there is a significant association between BMD and lean body mass and episodic verbal learning. While the involvement of modifiable lifestyle factors in human cognitive function has been examined in different studies, there is a need for further research to understand the potential underlying mechanisms. PMID:25741279
Density functional theory in the solid state
Hasnip, Philip J.; Refson, Keith; Probert, Matt I. J.; Yates, Jonathan R.; Clark, Stewart J.; Pickard, Chris J.
2014-01-01
Density functional theory (DFT) has been used in many fields of the physical sciences, but none so successfully as in the solid state. From its origins in condensed matter physics, it has expanded into materials science, high-pressure physics and mineralogy, solid-state chemistry and more, powering entire computational subdisciplines. Modern DFT simulation codes can calculate a vast range of structural, chemical, optical, spectroscopic, elastic, vibrational and thermodynamic phenomena. The ability to predict structure–property relationships has revolutionized experimental fields, such as vibrational and solid-state NMR spectroscopy, where it is the primary method to analyse and interpret experimental spectra. In semiconductor physics, great progress has been made in the electronic structure of bulk and defect states despite the severe challenges presented by the description of excited states. Studies are no longer restricted to known crystallographic structures. DFT is increasingly used as an exploratory tool for materials discovery and computational experiments, culminating in ex nihilo crystal structure prediction, which addresses the long-standing difficult problem of how to predict crystal structure polymorphs from nothing but a specified chemical composition. We present an overview of the capabilities of solid-state DFT simulations in all of these topics, illustrated with recent examples using the CASTEP computer program. PMID:24516184
Photostriction in Ferroelectrics from Density Functional Theory
NASA Astrophysics Data System (ADS)
Paillard, Charles; Xu, Bin; Dkhil, Brahim; Geneste, Grégory; Bellaiche, L.
2016-06-01
An ab initio procedure allowing the computation of the deformation of ferroelectric-based materials under light is presented. This numerical scheme consists in structurally relaxing the system under the constraint of a fixed ne concentration of electrons photoexcited into a specific conduction band edge state from a chosen valence band state, via the use of a constrained density functional theory method. The resulting change in lattice constant along a selected crystallographic direction is then calculated for a reasonable estimate of ne. This method is applied to bulk multiferroic BiFeO3 and predicts a photostriction effect of the same order of magnitude than the ones recently observed. A strong dependence of photostrictive response on both the reached conduction state and the crystallographic direction (along which this effect is determined) is also revealed. Furthermore, analysis of the results demonstrates that the photostriction mechanism mostly originates from the screening of the spontaneous polarization by the photoexcited electrons in combination with the inverse piezoelectric effect.
Dynamical density functional theory for microswimmers
NASA Astrophysics Data System (ADS)
Menzel, Andreas M.; Saha, Arnab; Hoell, Christian; Löwen, Hartmut
2016-01-01
Dynamical density functional theory (DDFT) has been successfully derived and applied to describe on one hand passive colloidal suspensions, including hydrodynamic interactions between individual particles. On the other hand, active "dry" crowds of self-propelled particles have been characterized using DDFT. Here, we go one essential step further and combine these two approaches. We establish a DDFT for active microswimmer suspensions. For this purpose, simple minimal model microswimmers are introduced. These microswimmers self-propel by setting the surrounding fluid into motion. They hydrodynamically interact with each other through their actively self-induced fluid flows and via the common "passive" hydrodynamic interactions. An effective soft steric repulsion is also taken into account. We derive the DDFT starting from common statistical approaches. Our DDFT is then tested and applied by characterizing a suspension of microswimmers, the motion of which is restricted to a plane within a three-dimensional bulk fluid. Moreover, the swimmers are confined by a radially symmetric trapping potential. In certain parameter ranges, we find rotational symmetry breaking in combination with the formation of a "hydrodynamic pumping state," which has previously been observed in the literature as a result of particle-based simulations. An additional instability of this pumping state is revealed.
Photostriction in Ferroelectrics from Density Functional Theory.
Paillard, Charles; Xu, Bin; Dkhil, Brahim; Geneste, Grégory; Bellaiche, L
2016-06-17
An ab initio procedure allowing the computation of the deformation of ferroelectric-based materials under light is presented. This numerical scheme consists in structurally relaxing the system under the constraint of a fixed n_{e} concentration of electrons photoexcited into a specific conduction band edge state from a chosen valence band state, via the use of a constrained density functional theory method. The resulting change in lattice constant along a selected crystallographic direction is then calculated for a reasonable estimate of n_{e}. This method is applied to bulk multiferroic BiFeO_{3} and predicts a photostriction effect of the same order of magnitude than the ones recently observed. A strong dependence of photostrictive response on both the reached conduction state and the crystallographic direction (along which this effect is determined) is also revealed. Furthermore, analysis of the results demonstrates that the photostriction mechanism mostly originates from the screening of the spontaneous polarization by the photoexcited electrons in combination with the inverse piezoelectric effect. PMID:27367406
The Nuclear Energy Density Functional Formalism
NASA Astrophysics Data System (ADS)
Duguet, T.
The present document focuses on the theoretical foundations of the nuclear energy density functional (EDF) method. As such, it does not aim at reviewing the status of the field, at covering all possible ramifications of the approach or at presenting recent achievements and applications. The objective is to provide a modern account of the nuclear EDF formalism that is at variance with traditional presentations that rely, at one point or another, on a Hamiltonian-based picture. The latter is not general enough to encompass what the nuclear EDF method represents as of today. Specifically, the traditional Hamiltonian-based picture does not allow one to grasp the difficulties associated with the fact that currently available parametrizations of the energy kernel E[g',g] at play in the method do not derive from a genuine Hamilton operator, would the latter be effective. The method is formulated from the outset through the most general multi-reference, i.e. beyond mean-field, implementation such that the single-reference, i.e. "mean-field", derives as a particular case. As such, a key point of the presentation provided here is to demonstrate that the multi-reference EDF method can indeed be formulated in a mathematically meaningful fashion even if E[g',g] does not derive from a genuine Hamilton operator. In particular, the restoration of symmetries can be entirely formulated without making any reference to a projected state, i.e. within a genuine EDF framework. However, and as is illustrated in the present document, a mathematically meaningful formulation does not guarantee that the formalism is sound from a physical standpoint. The price at which the latter can be enforced as well in the future is eventually alluded to.
A half century of density functional theory
Zangwill, Andrew
2015-07-15
Today’s most popular method for calculating the electronic structure of atoms, molecules, liquids, solids, and plasmas began as a bold hypothesis: The electron density distribution completely characterizes the ground state of a many-electron system.
Ions in solution: Density corrected density functional theory (DC-DFT)
Kim, Min-Cheol; Sim, Eunji; Burke, Kieron
2014-05-14
Standard density functional approximations often give questionable results for odd-electron radical complexes, with the error typically attributed to self-interaction. In density corrected density functional theory (DC-DFT), certain classes of density functional theory calculations are significantly improved by using densities more accurate than the self-consistent densities. We discuss how to identify such cases, and how DC-DFT applies more generally. To illustrate, we calculate potential energy surfaces of HO·Cl{sup −} and HO·H{sub 2}O complexes using various common approximate functionals, with and without this density correction. Commonly used approximations yield wrongly shaped surfaces and/or incorrect minima when calculated self consistently, while yielding almost identical shapes and minima when density corrected. This improvement is retained even in the presence of implicit solvent.
Optimizing a parametrized Thomas-Fermi-Dirac-Weizsäcker density functional for atoms.
Espinosa Leal, L A; Karpenko, A; Caro, M A; Lopez-Acevedo, O
2015-12-21
Because of issues with accuracy and transferability of existing orbital-free (OF) density functionals, OF functional development remains an active research area. However, due to numerical difficulties, all-electron self-consistent assessment of OF functionals is limited. Using an all-electron radial OFDFT code, we evaluate the performance of a parametrized OF functional for a wide range in parameter space. Specifically, we combine the parametrized Thomas-Fermi-Weizsäcker kinetic model (λ and γ for the fractions of Weizsäcker and Thomas-Fermi functionals, respectively) with a local density approximation (LDA) for the exchange-correlation functional. In order to obtain the converged results for λ values other than λ = 1, we use the potential scaling introduced in previous work. Because we work within a wide region in parameter space, this strategy provides an effective route towards better understanding of the parameter interplay that allows us to achieve good agreement with the Kohn-Sham (KS) model. Here, our interest lies in total energy, Euler equation eigenvalue, and electronic densities when the parameters are varied between 0.2 and 1.5. We observe that a one-to-one relation between λ and γ defines a region in parameter space that allows the atomic energies to be approximated with a very small average error (less than 3% percent for all the atoms studied) with respect to the KS reference energies. For each atom, the reference KS HOMO eigenvalue can also be reproduced with a similar error, but the one-to-one correspondence between λ and γ belongs to a different region of the same parameter space. Contrary to both properties, the atomic density behaves more smoothly and the error in reproducing the KS reference densities appears more insensitive to variation of the parameters (with mostly an average integrated difference of 0.15-0.20 |e| per electron). These results pave the way towards testing of parameter transferability and further systematic
a Renormalization Group Calculation of the Velocity - and Density-Density Correlation Functions.
NASA Astrophysics Data System (ADS)
Cowan, Mark Timothy
The velocity-velocity correlation function of a free field theory is obtained. The renormalization group, along with a 4-varepsilon expansion, is then used to find the leading order behavior of the velocity-velocity correlation function for an interacting field theory in the high temperature phase near the critical point. The details of the calculation of the density-density correlation function for Hedgehogs, in the context of a free field theory, is presented next. Finally the renormalization group, along with a 4-varepsilon expansion, is used to find the leading order behavior of the density-density correlation function for Hedgehogs in an interacting field theory near the critical point.
Fermionic density functional at a Feshbach resonance
Seidl, Michael; Bhaduri, Rajat K.
2007-05-15
We consider a dilute gas of neutral unpolarized fermionic atoms at zero temperature. The atoms interact via a short-range (tunable) attractive interaction. We demonstrate analytically a curious property of the gas at unitarity. Namely, the correlation energy of the gas, evaluated by second-order perturbation theory, has the same density dependence as the first-order exchange energy, and the two almost exactly cancel each other at a Feshbach resonance irrespective of the shape of the potential, provided ({mu}r{sub s})>>1. Here ({mu}){sup -1} is the range of the two-body potential, and r{sub s} is defined through the number density, n=3/(4{pi}r{sub s}{sup 3}). The implications of this result for universality are discussed.
Density functional theory for Yukawa fluids
NASA Astrophysics Data System (ADS)
Hatlo, Marius M.; Banerjee, Priyanka; Forsman, Jan; Lue, Leo
2012-08-01
We develop an approximate field theory for particles interacting with a generalized Yukawa potential. This theory improves and extends a previous splitting field theory, originally developed for counterions around a fixed charge distribution. The resulting theory bridges between the second virial approximation, which is accurate at low particle densities, and the mean-field approximation, accurate at high densities. We apply this theory to charged, screened ions in bulk solution, modeled to interact with a Yukawa potential; the theory is able to accurately reproduce the thermodynamic properties of the system over a broad range of conditions. The theory is also applied to "dressed counterions," interacting with a screened electrostatic potential, contained between charged plates. It is found to work well from the weak coupling to the strong coupling limits. The theory is able to reproduce the counterion profiles and force curves for closed and open systems obtained from Monte Carlo simulations.
Density-functional Monte-Carlo simulation of CuZn order-disorder transition
Khan, Suffian N.; Eisenbach, Markus
2016-01-25
We perform a Wang-Landau Monte Carlo simulation of a Cu0.5Zn0.5 order-disorder transition using 250 atoms and pairwise atom swaps inside a 5 x 5 x 5 BCC supercell. Each time step uses energies calculated from density functional theory (DFT) via the all-electron Korringa-Kohn- Rostoker method and self-consistent potentials. Here we find CuZn undergoes a transition from a disordered A2 to an ordered B2 structure, as observed in experiment. Our calculated transition temperature is near 870 K, comparing favorably to the known experimental peak at 750 K. We also plot the entropy, temperature, specific-heat, and short-range order as a function ofmore » internal energy.« less
Label-free all-electronic biosensing in microfluidic systems
NASA Astrophysics Data System (ADS)
Stanton, Michael A.
Label-free, all-electronic detection techniques offer great promise for advancements in medical and biological analysis. Electrical sensing can be used to measure both interfacial and bulk impedance changes in conducting solutions. Electronic sensors produced using standard microfabrication processes are easily integrated into microfluidic systems. Combined with the sensitivity of radiofrequency electrical measurements, this approach offers significant advantages over competing biological sensing methods. Scalable fabrication methods also provide a means of bypassing the prohibitive costs and infrastructure associated with current technologies. We describe the design, development and use of a radiofrequency reflectometer integrated into a microfluidic system towards the specific detection of biologically relevant materials. We developed a detection protocol based on impedimetric changes caused by the binding of antibody/antigen pairs to the sensing region. Here we report the surface chemistry that forms the necessary capture mechanism. Gold-thiol binding was utilized to create an ordered alkane monolayer on the sensor surface. Exposed functional groups target the N-terminus, affixing a protein to the monolayer. The general applicability of this method lends itself to a wide variety of proteins. To demonstrate specificity, commercially available mouse anti- Streptococcus Pneumoniae monoclonal antibody was used to target the full-length recombinant pneumococcal surface protein A, type 2 strain D39 expressed by Streptococcus Pneumoniae. We demonstrate the RF response of the sensor to both the presence of the surface decoration and bound SPn cells in a 1x phosphate buffered saline solution. The combined microfluidic sensor represents a powerful platform for the analysis and detection of cells and biomolecules.
Molecular dynamics simulation of liquid water: Hybrid density functionals
Todorova, T; Seitsonen, A; Hutter, J; Kuo, W; Mundy, C
2005-09-12
The structure, dynamical and electronic properties of liquid water utilizing different hybrid density functionals were tested within the plane wave framework of first principles molecular dynamics simulations. The computational approach, which employs modified functionals with short-ranged Hartree-Fock exchange, was first tested in calculations of the structural and bonding properties of the water dimer and cyclic water trimer. Liquid water simulations were performed at the state point of 350 K at the experimental density. Simulations included three different hybrid functionals, a meta functional, four gradient corrected functionals, the local density and Hartree-Fock approximation. It is found that hybrid functionals are superior in reproducing the experimental structure and dynamical properties as measured by the radial distribution function and self diffusion constant when compared to the pure density functionals. The local density and Hartree-Fock approximations show strongly over- and under-structured liquids, respectively. Hydrogen bond analysis shows that the hybrid functionals give slightly smaller averaged numbers of hydrogen bonds and similar hydrogen bond populations as pure density functionals. The average molecular dipole moments in the liquid from the three hybrid functionals are lower than from the corresponding pure density functionals.
Phonons in nonlocal van der Waals density functional theory
NASA Astrophysics Data System (ADS)
Sabatini, Riccardo; Küçükbenli, Emine; Pham, Cong Huy; de Gironcoli, Stefano
2016-06-01
We extend the formulation of density functional perturbation theory to treat nonlocal density functionals, accounting for van der Waals interactions, in a rigorous and efficient way. We provide a general formalism, suitable for any functional in this family, and give specific equations for the most widely used ones. We then study the lattice dynamics of graphite, comparing several nonlocal functionals and the local density approximation, showing that our recent revision of the VV10 functional [R. Sabatini et al., Phys. Rev. B 87, 041108(R) (2013), 10.1103/PhysRevB.87.041108] gives the best comparison with experiments.
Exact conditions on the temperature dependence of density functionals
NASA Astrophysics Data System (ADS)
Burke, K.; Smith, J. C.; Grabowski, P. E.; Pribram-Jones, A.
2016-05-01
Universal exact conditions guided the construction of most ground-state density functional approximations in use today. We derive the relation between the entropy and Mermin free energy density functionals for thermal density functional theory. Both the entropy and sum of kinetic and electron-electron repulsion functionals are shown to be monotonically increasing with temperature, while the Mermin functional is concave downwards. Analogous relations are found for both exchange and correlation. The importance of these conditions is illustrated in two extremes: the Hubbard dimer and the uniform gas.
Density-potential mapping in time-dependent density-functional theory
Maitra, N. T.; Todorov, T. N.; Woodward, C.; Burke, K.
2010-04-15
The key questions of uniqueness and existence in time-dependent density-functional theory are usually formulated only for potentials and densities that are analytic in time. Simple examples, standard in quantum mechanics, lead, however, to nonanalyticities. We reformulate these questions in terms of a nonlinear Schroedinger equation with a potential that depends nonlocally on the wave function.
Spin constraints on nuclear energy density functionals
NASA Astrophysics Data System (ADS)
Robledo, L. M.; Bernard, R. N.; Bertsch, G. F.
2014-02-01
The Gallagher-Moszkowski rule in the spectroscopy of odd-odd nuclei imposes a new spin constraint on the energy functionals for self-consistent mean field theory. The commonly used parametrization of the effective three-body interaction in the Gogny and Skyrme families of energy functionals is ill suited to satisfy the spin constraint. In particular, the Gogny parametrization of the three-body interaction has the spin dependence opposite to that required by the observed spectra. The two-body part has a correct sign, but in combination the rule is violated as often as not. We conclude that a new functional form is needed for the effective three-body interaction that can take into better account the different spin-isospin channels of the interaction.
Introduction to Classical Density Functional Theory by a Computational Experiment
ERIC Educational Resources Information Center
Jeanmairet, Guillaume; Levy, Nicolas; Levesque, Maximilien; Borgis, Daniel
2014-01-01
We propose an in silico experiment to introduce the classical density functional theory (cDFT). Density functional theories, whether quantum or classical, rely on abstract concepts that are nonintuitive; however, they are at the heart of powerful tools and active fields of research in both physics and chemistry. They led to the 1998 Nobel Prize in…
Nuclear collective excitations: A relativistic density functional approach
NASA Astrophysics Data System (ADS)
Piekarewicz, J.
2015-08-01
Density functional theory provides the most promising, and likely unique, microscopic framework to describe nuclear systems ranging from finite nuclei to neutron stars. Properly optimized energy density functionals define a new paradigm in nuclear theory where predictive capability is possible and uncertainty quantification is demanded. Moreover, density functional theory offers a consistent approach to the linear response of the nuclear ground state. In this paper, we review the fundamental role played by nuclear collective modes in uncovering novel excitations and in guiding the optimization of the density functional. Indeed, without collective excitations the determination of the density functional remains incomplete. Without collective excitations, the equation of state of neutron-rich matter continues to be poorly constrained. We conclude with a discussion of some of the remaining challenges in this field and propose a path forward to address these challenges.
Petrenko, Taras; Kossmann, Simone; Neese, Frank
2011-02-01
In this paper, we present the implementation of efficient approximations to time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation (TDA) for hybrid density functionals. For the calculation of the TDDFT/TDA excitation energies and analytical gradients, we combine the resolution of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calculation of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approximation leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 molecules with the configuration interaction singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equilibrium geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calculated excited state normal coordinate displacements is of the order of 0.01, and relative error in the calculated excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approximation are, thus, insignificant compared to the errors related to the approximate nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calculations on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ~26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ~27-29 for 30 processors). The solution of the Z-vector equations shows a speedup of ~24 on 30 processors. The
NASA Astrophysics Data System (ADS)
Petrenko, Taras; Kossmann, Simone; Neese, Frank
2011-02-01
In this paper, we present the implementation of efficient approximations to time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation (TDA) for hybrid density functionals. For the calculation of the TDDFT/TDA excitation energies and analytical gradients, we combine the resolution of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calculation of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approximation leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 molecules with the configuration interaction singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equilibrium geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calculated excited state normal coordinate displacements is of the order of 0.01, and relative error in the calculated excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approximation are, thus, insignificant compared to the errors related to the approximate nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calculations on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ˜26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ˜27-29 for 30 processors). The solution of the Z-vector equations shows a speedup of ˜24 on 30 processors. The
NASA Astrophysics Data System (ADS)
Bochevarov, Arteum D.; Friesner, Richard A.
2008-01-01
We investigate one of the fundamental observables, electronic charge density, as produced by a number of popular functionals of the density functional theory (DFT): SVWN5, B3LYP, B3LYP, OLYP, O3LYP, BP86, B3P86, O3P86, and PBE using restricted and unrestricted orbitals. Measuring and comparing the quality of the densities could tell us more about the physical soundness of the functional models. The study is performed on the small molecules He, H2, LiH, H4 in an extensive range of correlation-consistent basis sets. We compare DFT densities to those of full configuration interaction (FCI) under the assumption that the FCI density in the largest employed basis set is sufficiently close to the exact one. For LiH and H4, we also compare the DFT densities to those of CCSD. The SVWN5 functional consistently shows the worst performance. The OPTX exchange functional regularly beats the Becke exchange. Among the best performers are all the hybrid functionals, the novel O3P86 being the most accurate in most cases. The popular functional B3LYP was consistently outmatched by O3LYP, and produced, in fact, some of the poorest densities among the hybrids. CCSD was found to produce much more accurate densities than any DFT functional in the case of LiH in equilibrium geometry, but was sometimes outperformed by DFT in the case of slightly stretched H4, where CCSD theory itself starts to break down. Surprisingly, as one stretches the H2 molecule, BP86 and PBE improve the description of density although such behavior is not observed in other systems. We conclude by reasoning how functionals such as B3LYP, despite being quite average for density, could still be very successful in predicting thermodynamic properties.
Molecular density functional theory of water including density-polarization coupling.
Jeanmairet, Guillaume; Levy, Nicolas; Levesque, Maximilien; Borgis, Daniel
2016-06-22
We present a three-dimensional molecular density functional theory of water derived from first-principles that relies on the particle's density and multipolar polarization density and includes the density-polarization coupling. This brings two main benefits: (i) scalar density and vectorial multipolar polarization density fields are much more tractable and give more physical insight than the full position and orientation densities, and (ii) it includes the full density-polarization coupling of water, that is known to be non-vanishing but has never been taken into account. Furthermore, the theory requires only the partial charge distribution of a water molecule and three measurable bulk properties, namely the structure factor and the Fourier components of the longitudinal and transverse dielectric susceptibilities. PMID:27116250
The benchmark of gutzwiller density functional theory in hydrogen systems
Yao, Y.; Wang, Cai-Zhuang; Ho, Kai-Ming
2012-02-23
We propose an approximate form of the exchange-correlation energy functional for the Gutzwiller density functional theory. It satisfies certain physical constraints in both weak and strong electron correlation limits. We benchmark the Gutzwiller density functional approximation in the hydrogen systems, where the static correlation error is shown to be negligible. The good transferability is demonstrated by applications to the hydrogen molecule and some crystal structures.
The Benchmark of Gutzwiller Density Functional Theory in Hydrogen Systems
Yao, Yongxin; Wang, Cai-Zhuang; Ho, Kai-Ming
2011-01-13
We propose an approximate form of the exchange-correlation energy functional for the Gutzwiller density functional theory. It satisfies certain physical constraints in both weak and strong electron correlation limits. We benchmark the Gutzwiller density functional approximation in the hydrogen systems, where the static correlation error is shown to be negligible. The good transferability is demonstrated by applications to the hydrogen molecule and some crystal structures. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012
Impact ionization in GaAs: A screened exchange density-functional approach
Picozzi, S.; Asahi, R.; Geller, C. B.; Continenza, A.; Freeman, A. J.
2001-08-13
Results are presented of a fully ab initio calculation of impact ionization rates in GaAs within the density functional theory framework, using a screened-exchange formalism and the highly precise all-electron full-potential linearized augmented plane wave method. The calculated impact ionization rates show a marked orientation dependence in k space, indicating the strong restrictions imposed by the conservation of energy and momentum. This anisotropy diminishes as the impacting electron energy increases. A Keldysh type fit performed on the energy-dependent rate shows a rather soft edge and a threshold energy greater than the direct band gap. The consistency with available Monte Carlo and empirical pseudopotential calculations shows the reliability of our approach and paves the way to ab initio calculations of pair production rates in new and more complex materials.
Generalized van der Waals density functional theory for nonuniform polymers
Patra, Chandra N.; Yethiraj, Arun
2000-01-15
A density functional theory is presented for the effect of attractions on the structure of polymers at surfaces. The theory treats the ideal gas functional exactly, and uses a weighted density approximation for the hard chain contribution to the excess free energy functional. The attractive interactions are treated using a van der Waals approximation. The theory is in good agreement with computer simulations for the density profiles at surfaces for a wide range of densities and temperatures, except for low polymer densities at low temperatures where it overestimates the depletion of chains from the surface. This deficiency is attributed to the neglect of liquid state correlations in the van der Waals term of the free energy functional. (c) 2000 American Institute of Physics.
Exact maps in density functional theory for lattice models
NASA Astrophysics Data System (ADS)
Dimitrov, Tanja; Appel, Heiko; Fuks, Johanna I.; Rubio, Angel
2016-08-01
In the present work, we employ exact diagonalization for model systems on a real-space lattice to explicitly construct the exact density-to-potential and graphically illustrate the complete exact density-to-wavefunction map that underly the Hohenberg–Kohn theorem in density functional theory. Having the explicit wavefunction-to-density map at hand, we are able to construct arbitrary observables as functionals of the ground-state density. We analyze the density-to-potential map as the distance between the fragments of a system increases and the correlation in the system grows. We observe a feature that gradually develops in the density-to-potential map as well as in the density-to-wavefunction map. This feature is inherited by arbitrary expectation values as functional of the ground-state density. We explicitly show the excited-state energies, the excited-state densities, and the correlation entropy as functionals of the ground-state density. All of them show this exact feature that sharpens as the coupling of the fragments decreases and the correlation grows. We denominate this feature as intra-system steepening and discuss how it relates to the well-known inter-system derivative discontinuity. The inter-system derivative discontinuity is an exact concept for coupled subsystems with degenerate ground state. However, the coupling between subsystems as in charge transfer processes can lift the degeneracy. An important conclusion is that for such systems with a near-degenerate ground state, the corresponding cut along the particle number N of the exact density functionals is differentiable with a well-defined gradient near integer particle number.
NASA Astrophysics Data System (ADS)
Painter, G. S.; Averill, F. W.
1983-11-01
The principle of augmentation, used to introduce inner-atom core structure into slowly varying basis functions, is applied to Gaussian orbitals to define a new basis set for highly accurate total-energy calculations for atomic clusters within the density-functional formalism. Diffuse Gaussian-orbital tails are matched continuously and differentiably to inner-atom numeric radial functions at the atomic-sphere radius. Major advantages of Gaussian-orbital basis sets are acquired without the need for numerous Gaussians of large exponent for the core region. The numeric functions used inside the atom permit essentially exact solutions for that region. Procedures are described which recover use of the efficient integral algorithms for the Gaussian-orbital-tail matrix elements. The interactions over the structured inner-atom region are treated by efficient integrand smoothing and integration procedures for the sphere. The new augmented Gaussian basis removes the primary limitations on the use of Gaussian orbitals for heavy atoms. As an illustration the method is applied to the copper dimer in an all-electron framework within the local-spin-density approximation (LSDA). The calculated binding energy, equilibrium separation, and first ionization potential of Cu2 are within 2% of experiment within the Xα model. Excitation energies are better described within more recent refined exchange-correlation functionals. These all-electron results show the LSDA model predicts a slightly contracted bond length for Cu2, consistent with bulk LSDA calculations for the 3d transition-metal series.
Postfragmentation density function for bacterial aggregates in laminar flow
NASA Astrophysics Data System (ADS)
Byrne, Erin; Bortz, David M.; Dzul, Steve; Solomon, Michael; Younger, John
2011-04-01
The postfragmentation probability density of daughter flocs is one of the least well-understood aspects of modeling flocculation. We use three-dimensional positional data of Klebsiella pneumoniae bacterial flocs in suspension and the knowledge of hydrodynamic properties of a laminar flow field to construct a probability density function of floc volumes after a fragmentation event. We provide computational results which predict that the primary fragmentation mechanism for large flocs is erosion. The postfragmentation probability density function has a strong dependence on the size of the original floc and indicates that most fragmentation events result in clumps of one to three bacteria eroding from the original floc. We also provide numerical evidence that exhaustive fragmentation yields a limiting density inconsistent with the log-normal density predicted in the literature, most likely due to the heterogeneous nature of K. pneumoniae flocs. To support our conclusions, artificial flocs were generated and display similar postfragmentation density and exhaustive fragmentation.
Predicting Pt-195 NMR chemical shift using new relativistic all-electron basis set.
Paschoal, D; Guerra, C Fonseca; de Oliveira, M A L; Ramalho, T C; Dos Santos, H F
2016-10-01
Predicting NMR properties is a valuable tool to assist the experimentalists in the characterization of molecular structure. For heavy metals, such as Pt-195, only a few computational protocols are available. In the present contribution, all-electron Gaussian basis sets, suitable to calculate the Pt-195 NMR chemical shift, are presented for Pt and all elements commonly found as Pt-ligands. The new basis sets identified as NMR-DKH were partially contracted as a triple-zeta doubly polarized scheme with all coefficients obtained from a Douglas-Kroll-Hess (DKH) second-order scalar relativistic calculation. The Pt-195 chemical shift was predicted through empirical models fitted to reproduce experimental data for a set of 183 Pt(II) complexes which NMR sign ranges from -1000 to -6000 ppm. Furthermore, the models were validated using a new set of 75 Pt(II) complexes, not included in the descriptive set. The models were constructed using non-relativistic Hamiltonian at density functional theory (DFT-PBEPBE) level with NMR-DKH basis set for all atoms. For the best model, the mean absolute deviation (MAD) and the mean relative deviation (MRD) were 150 ppm and 6%, respectively, for the validation set (75 Pt-complexes) and 168 ppm (MAD) and 5% (MRD) for all 258 Pt(II) complexes. These results were comparable with relativistic DFT calculation, 200 ppm (MAD) and 6% (MRD). © 2016 Wiley Periodicals, Inc. PMID:27510431
Integration, Continuity and a Connection with Probability Density Functions
ERIC Educational Resources Information Center
Samuels, M.
2006-01-01
This note considers functions of two variables which are continuous on a possibly unbounded closed region in [vertical bar]R[squared], and the functions of one variable obtained by integrating out the other variable over this region. The question of continuity of these functions is investigated, as are connections with joint density and marginal…
NASA Astrophysics Data System (ADS)
Knuth, Franz; Carbogno, Christian; Atalla, Viktor; Blum, Volker; Scheffler, Matthias
2015-05-01
We derive and implement the strain derivatives of the total energy of solids, i.e., the analytic stress tensor components, in an all-electron, numeric atom-centered orbital based density-functional formalism. We account for contributions that arise in the semi-local approximation (LDA/GGA) as well as in the generalized Kohn-Sham case, in which a fraction of exact exchange (hybrid functionals) is included. In this work, we discuss the details of the implementation including the numerical corrections for sparse integrations grids which allow to produce accurate results. We validate the implementation for a variety of test cases by comparing to strain derivatives performed via finite differences. Additionally, we include the detailed definition of the overlapping atom-centered integration formalism used in this work to obtain total energies and their derivatives.
Optical properties of alkali halide crystals from all-electron hybrid TD-DFT calculations
Webster, R. Harrison, N. M.; Bernasconi, L.
2015-06-07
We present a study of the electronic and optical properties of a series of alkali halide crystals AX, with A = Li, Na, K, Rb and X = F, Cl, Br based on a recent implementation of hybrid-exchange time-dependent density functional theory (TD-DFT) (TD-B3LYP) in the all-electron Gaussian basis set code CRYSTAL. We examine, in particular, the impact of basis set size and quality on the prediction of the optical gap and exciton binding energy. The formation of bound excitons by photoexcitation is observed in all the studied systems and this is shown to be correlated to specific features of the Hartree-Fock exchange component of the TD-DFT response kernel. All computed optical gaps and exciton binding energies are however markedly below estimated experimental and, where available, 2-particle Green’s function (GW-Bethe-Salpeter equation, GW-BSE) values. We attribute this reduced exciton binding to the incorrect asymptotics of the B3LYP exchange correlation ground state functional and of the TD-B3LYP response kernel, which lead to a large underestimation of the Coulomb interaction between the excited electron and hole wavefunctions. Considering LiF as an example, we correlate the asymptotic behaviour of the TD-B3LYP kernel to the fraction of Fock exchange admixed in the ground state functional c{sub HF} and show that there exists one value of c{sub HF} (∼0.32) that reproduces at least semi-quantitatively the optical gap of this material.
Testing for parity violation in nuclei using spin density matrices for nuclear density functionals
NASA Astrophysics Data System (ADS)
Barrett, B. R.; Giraud, B. G.
2015-06-01
The spin density matrix (SDM) used in atomic and molecular physics is revisited for nuclear physics, in the context of the radial density functional theory. The vector part of the SDM defines a ‘hedgehog’ situation, which exists only if nuclear states contain some amount of parity violation. A toy model is given as an illustrative example.
Fluctuation-dissipation theorem density-functional theory
NASA Astrophysics Data System (ADS)
Furche, Filipp; Van Voorhis, Troy
2005-04-01
Using the fluctuation-dissipation theorem (FDT) in the context of density-functional theory (DFT), one can derive an exact expression for the ground-state correlation energy in terms of the frequency-dependent density response function. When combined with time-dependent density-functional theory, a new class of density functionals results that use approximations to the exchange-correlation kernel fxc as input. This FDT-DFT scheme holds promise to solve two of the most distressing problems of conventional Kohn-Sham DFT: (i) It leads to correlation energy functionals compatible with exact exchange, and (ii) it naturally includes dispersion. The price is a moderately expensive O(N6) scaling of computational cost and a slower basis set convergence. These general features of FDT-DFT have all been recognized previously. In this paper, we present the first benchmark results for a set of molecules using FDT-DFT beyond the random-phase approximation (RPA)—that is, the first such results with fxc≠0. We show that kernels derived from the adiabatic local-density approximation and other semilocal functionals suffer from an "ultraviolet catastrophe," producing a pair density that diverges at small interparticle distance. Nevertheless, dispersion interactions can be treated accurately if hybrid functionals are employed, as is demonstrated for He2 and HeNe. We outline constraints that future approximations to fxc should satisfy and discuss the prospects of FDT-DFT.
Charge and spin fluctuations in the density functional theory
Gyoerffy, B.L.; Barbieri, A. . H.H. Wills Physics Lab.); Staunton, J.B. . Dept. of Physics); Shelton, W.A.; Stocks, G.M. )
1990-01-01
We introduce a conceptual framework which allow us to treat charge and spin fluctuations about the Local density Approximation (LDA) to the Density Functional Theory (DFT). We illustrate the approach by explicit study of the Disordered Local Moment (DLM) state in Fe above the Curie Temperature {Tc} and the Mott insulating state in MnO. 27 refs., 6 figs.
Density Functional Model for Nondynamic and Strong Correlation.
Kong, Jing; Proynov, Emil
2016-01-12
A single-term density functional model for the left-right nondynamic/strong electron correlation is presented based on single-determinant Kohn-Sham density functional theory. It is derived from modeling the adiabatic connection for kinetic correlation energy based on physical arguments, with the correlation potential energy based on the Becke'13 model ( Becke, A.D. J. Chem. Phys . 2013 , 138 , 074109 ). This functional satisfies some known scaling relationships for correlation functionals. The fractional spin error is further reduced substantially with a new density-functional correction. Preliminary tests with self-consistent-field implementation show that the model, with only three empirical parameters, recovers the majority of left-right nondynamic/strong correlation upon bond dissociation and performs reasonably well for atomization energies and singlet-triplet energy splittings. This study also demonstrates the feasibility of developing DFT functionals for nondynamic and strong correlation within the single-determinant KS scheme. PMID:26636190
Density and Spin Response Functions in Ultracold Fermionic Atom Gases
Mihaila, Bogdan; Blagoev, Krastan B.; Balatsky, Alexander V.; Smith, Darryl L.; Gaudio, Sergio; Littlewood, Peter B.
2005-08-26
We propose a new method of detecting the onset of superfluidity in a two-component ultracold fermionic gas of atoms governed by an attractive short-range interaction. By studying the two-body correlation functions we find that a measurement of the momentum distribution of the density and spin-response functions allows one to access separately the normal and anomalous densities. The change in sign at low momentum transfer of the normal-ordered part of the density response function signals the transition between a BEC and a BCS regime, characterized by small and large pairs, respectively. This change in sign of the density response function represents an unambiguous signature of the BEC-to-BCS crossover. Spin rotational symmetry breaking due to the magnetic field, if observed, can be used to validate the one-channel model.
Density Functional Theory with Dissipation: Transport through Single Molecules
Kieron Burke
2012-04-30
A huge amount of fundamental research was performed on this grant. Most of it focussed on fundamental issues of electronic structure calculations of transport through single molecules, using density functional theory. Achievements were: (1) First density functional theory with dissipation; (2) Pseudopotential plane wave calculations with master equation; (3) Weak bias limit; (4) Long-chain conductance; and (5) Self-interaction effects in tunneling.
Gutzwiller density functional theory for correlated electron systems
Ho, K. M.; Schmalian, J.; Wang, C. Z.
2008-02-04
We develop a density functional theory (DFT) and formalism for correlated electron systems by taking as reference an interacting electron system that has a ground state wave function which exactly obeys the Gutzwiller approximation for all one-particle operators. The solution of the many-electron problem is mapped onto the self-consistent solution of a set of single-particle Schroedinger equations, analogously to standard DFT-local density approximation calculations.
Preface: Special Topic on Advances in Density Functional Theory
Yang, Weitao
2014-05-14
This Special Topic Issue on the Advances in Density Functional Theory, published as a celebration of the fifty years of density functional theory, contains a retrospective article, a perspective article, and a collection of original research articles that showcase recent theoretical advances in the field. It provides a timely discussion reflecting a cross section of our understanding, and the theoretical and computational developments, which have significant implications in broad areas of sciences and engineering.
Preface: Special topic on advances in density functional theory.
Yang, Weitao
2014-05-14
This Special Topic Issue on the Advances in Density Functional Theory, published as a celebration of the fifty years of density functional theory, contains a retrospective article, a perspective article, and a collection of original research articles that showcase recent theoretical advances in the field. It provides a timely discussion reflecting a cross section of our understanding, and the theoretical and computational developments, which have significant implications in broad areas of sciences and engineering. PMID:24832306
Reflection-Asymmetric Nuclear Deformations within the Density Functional Theory
Olsen, E; Erler, J; Nazarewicz, W.; Stoitsov, M
2012-01-01
Within the nuclear density functional theory (DFT) we study the effect of reflection- asymmetric shapes on ground-state binding energies and binding energy differences. To this end, we developed the new DFT solver axialhfb that uses an approximate second-order gradient to solve the Hartree-Fock-Bogoliubov equations of superconducting DFT with the quasi-local Skyrme energy density functionals. Illustrative calculations are carried out for even- even isotopes of radium and thorium.
Density functional calculations for prediction of ultra-thin film structure and properties
Trickey, S.B.; Boettger, J.C.
1995-09-01
Most of the published theoretical and calculational effort on unsupported, ordered, ultra-thin films (``UTF``) in vacuo has focused on the thickest computationally feasible systems as models for surface properties of semi-infinite slabs. Crystalline periodic length scales in two cartesian dimensions combined with molecular- scale thickness in the third, however, make UTF`s strong candidates for the occurrence of quantum interference effects. Many UTF properties were predicted first from jellium slab models. A noteworthy prediction was that there would be large oscillations in the work function as a function of layer number. Extensive calculations on a variety of N-layers (N = 1,2,3... atomic planes) using all-electron, full-potential, local-spin-density approximation techniques show that the work function oscillation is weaker than expected but that there are significant layer-number dependences in the equilibrium lattice parameters, inter-planar spacings, electronic structure, density of states, and electronic stopping power. This paper reviews our own calculations as well as some others. Our objectives include the discernment of systematics within UTF`s, systematics in relationship to their counterpart crystals, relationship with surface properties, and appraisal of challenges to current models and methods.
Energy Densities in the Strong-Interaction Limit of Density Functional Theory.
Mirtschink, André; Seidl, Michael; Gori-Giorgi, Paola
2012-09-11
We discuss energy densities in the strong-interaction limit of density functional theory, deriving an exact expression within the definition (gauge) of the electrostatic potential of the exchange-correlation hole. Exact results for small atoms and small model quantum dots (Hooke's atoms) are compared with available approximations defined in the same gauge. The idea of a local interpolation along the adiabatic connection is discussed, comparing the energy densities of the Kohn-Sham, the physical, and the strong-interacting systems. We also use our results to analyze the local version of the Lieb-Oxford bound, widely used in the construction of approximate exchange-correlation functionals. PMID:26605721
Comments on the locality in density-functional theory
Lindgren, Ingvar; Salomonson, Sten
2003-05-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.
Density functional theory for polymeric systems in 2D
NASA Astrophysics Data System (ADS)
Słyk, Edyta; Roth, Roland; Bryk, Paweł
2016-06-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.
Density functional theory for polymeric systems in 2D.
Słyk, Edyta; Roth, Roland; Bryk, Paweł
2016-06-22
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. PMID:27115343
Franca, Vivian V.; D'Amico, Irene
2011-04-15
We derive an analytical density functional for the single-site entanglement of the one-dimensional homogeneous Hubbard model by means of an approximation to the linear entropy. We show that this very simple density functional reproduces quantitatively the exact results. We then use this functional as input for a local-density approximation to the single-site entanglement of inhomogeneous systems. We illustrate the power of this approach in a harmonically confined system, which could simulate recent experiments with ultracold atoms in optical lattices as well as in a superlattice and in an impurity system. The impressive quantitative agreement with numerical calculations--which includes reproducing subtle signatures of the particle density stages--shows that our density functional can provide entanglement calculations for actual experiments via density measurements. Next we use our functional to calculate the entanglement in disordered systems. We find that, in contrast with the expectation that disorder destroys the entanglement, there exist regimes for which the entanglement remains almost unaffected by the presence of disordered impurities.
Reactivity of Graphene Investigated by Density-Functional Theory
NASA Astrophysics Data System (ADS)
Soni, Himadri; Gebhardt, Julian; Görling, Andreas; Chair of Theoretical Chemistry Team
Using spin-polarized density-functional theory, we study the adsorption and reaction of hydrogen and fluorine with graphene. Graphene has a bipartite lattice with two different sublattices and hence, due to Lieb's theorem, the inequality between two sublattices should lead to a net magnetic moment upon adsorption of hydrogen or fluorine. Our calculations using density-functional theory with the generalized gradient approximation predict a magnetic moment of 1 µB for a single hydrogen adsorbed on graphene but not for a single fluorine atom adsorbed on graphene. Switching to hybrid density-functional theory with the HSE functional, we obtain a magnetic moment of 1 µB for of a single fluorine atom adsorption on graphene. This is in line with work of Kim et al., who also found in density-functional theory calculations with the HSE exchange-correlation functional spin-polarization for a fluorine adatom on graphene. Here, we present a systematic study of the reactivity and relevant adsorption mechanism for single-sided graphene, i.e., a graphene sheet which is accessible by an adsorbate from only one side with hydrogen and fluorine using hybrid density-functional theory. German Research Council (DFG) by the Collaborative Research Center 953.
Range Separation and Local Hybridization in Density Functional Theory†
Henderson, Thomas M.; Janesko, Benjamin G.; Scuseria, Gustavo E.
2016-01-01
Kohn–Sham density functional theory has become a standard method for modeling energetic, spectroscopic, and chemical reactivity properties of large molecules and solids. Density functional theory provides a rigorous theoretical framework for modeling the many-body exchange-correlation effects that dominate the computational cost of traditional wave function approaches. The advent of hybrid exchange-correlation functionals which incorporate a fraction of nonlocal exact exchange has solidified the prominence of density functional theory within computational chemistry. Hybrids provide accurate treatments of properties such as thermochemistry and molecular geometry. But they also exhibit some rather spectacular failures, and often contain multiple empirical parameters. This article reviews our work on developing novel exchange-correlation functionals that build upon the successes of global hybrids. We focus on more flexible functional forms, including local and range-separated hybrid functionals, constructed to obey known exact constraints and (ideally) to incorporate a minimum of empirical parametrization. The article places our work within the context of some other new approximate density functionals and discusses prospects for future work. PMID:19006280
Optimal-transport formulation of electronic density-functional theory
NASA Astrophysics Data System (ADS)
Buttazzo, Giuseppe; De Pascale, Luigi; Gori-Giorgi, Paola
2012-06-01
The most challenging scenario for Kohn-Sham density-functional theory, that is, when the electrons move relatively slowly trying to avoid each other as much as possible because of their repulsion (strong-interaction limit), is reformulated here as an optimal transport (or mass transportation theory) problem, a well-established field of mathematics and economics. In practice, we show that to solve the problem of finding the minimum possible internal repulsion energy for N electrons in a given density ρ(r) is equivalent to find the optimal way of transporting N-1 times the density ρ into itself, with the cost function given by the Coulomb repulsion. We use this link to set the strong-interaction limit of density-functional theory on firm ground and to discuss the potential practical aspects of this reformulation.
Relaxation of Actinide Surfaces: An All Electron Study
NASA Astrophysics Data System (ADS)
Atta-Fynn, Raymond; Dholabhai, Pratik; Ray, Asok
2006-10-01
Fully relativistic full potential density functional calculations with a linearized augmented plane wave plus local orbitals basis (LAPW + lo) have been performed to investigate the relaxations of heavy actinide surfaces, namely the (111) surface of fcc δ-Pu and the (0001) surface of dhcp Am using WIEN2k. This code uses the LAPW + lo method with the unit cell divided into non-overlapping atom-centered spheres and an interstitial region. The APW+lo basis is used to describe all s, p, d, and f states and LAPW basis to describe all higher angular momentum states. Each surface was modeled by a three-layer periodic slab separated by 60 Bohr vacuum with four atoms per surface unit cell. In general, we have found a contraction of the interlayer separations for both Pu and Am. We will report, in detail, the electronic and geometric structures of the relaxed surfaces and comparisons with the respective non-relaxed surfaces.
Optimization of an exchange-correlation density functional for water.
Fritz, Michelle; Fernández-Serra, Marivi; Soler, José M
2016-06-14
We describe a method, that we call data projection onto parameter space (DPPS), to optimize an energy functional of the electron density, so that it reproduces a dataset of experimental magnitudes. Our scheme, based on Bayes theorem, constrains the optimized functional not to depart unphysically from existing ab initio functionals. The resulting functional maximizes the probability of being the "correct" parameterization of a given functional form, in the sense of Bayes theory. The application of DPPS to water sheds new light on why density functional theory has performed rather poorly for liquid water, on what improvements are needed, and on the intrinsic limitations of the generalized gradient approximation to electron exchange and correlation. Finally, we present tests of our water-optimized functional, that we call vdW-DF-w, showing that it performs very well for a variety of condensed water systems. PMID:27305990
Optimization of an exchange-correlation density functional for water
NASA Astrophysics Data System (ADS)
Fritz, Michelle; Fernández-Serra, Marivi; Soler, José M.
2016-06-01
We describe a method, that we call data projection onto parameter space (DPPS), to optimize an energy functional of the electron density, so that it reproduces a dataset of experimental magnitudes. Our scheme, based on Bayes theorem, constrains the optimized functional not to depart unphysically from existing ab initio functionals. The resulting functional maximizes the probability of being the "correct" parameterization of a given functional form, in the sense of Bayes theory. The application of DPPS to water sheds new light on why density functional theory has performed rather poorly for liquid water, on what improvements are needed, and on the intrinsic limitations of the generalized gradient approximation to electron exchange and correlation. Finally, we present tests of our water-optimized functional, that we call vdW-DF-w, showing that it performs very well for a variety of condensed water systems.
Density functional for ternary non-additive hard sphere mixtures.
Schmidt, Matthias
2011-10-19
Based on fundamental measure theory, a Helmholtz free energy density functional for three-component mixtures of hard spheres with general, non-additive interaction distances is constructed. The functional constitutes a generalization of the previously given theory for binary non-additive mixtures. The diagrammatic structure of the spatial integrals in both functionals is of star-like (or tree-like) topology. The ternary diagrams possess a higher degree of complexity than the binary diagrams. Results for partial pair correlation functions, obtained via the Ornstein-Zernike route from the second functional derivatives of the excess free energy functional, agree well with Monte Carlo simulation data. PMID:21946780
Modal density function and number of propagating modes in ducts
NASA Technical Reports Server (NTRS)
Rice, E. J.
1976-01-01
Often raised questions in duct sound propagation studies involve the total number of propagating modes, the number of propagating radial modes for a particular spinning lobe number, and the number of modes possible between two given values of cutoff ratio or eigenvalue. These questions can be answered approximately by using the modal distribution function which is the integral of the modal density function for ducts in a manner similar to that previously published for architectural acoustics. The modal density functions are derived for rectangular and circular ducts with a uniform steady flow. Results from this continuous theory are compared to the actual (discrete) modal distributions.
Van der Waals density functional applied to adsorption systems
NASA Astrophysics Data System (ADS)
Hamada, Ikutaro
2013-03-01
The van der Waals density functional (vdW-DF) is a promising density functional to describe the van der Waals forces within density functional theory. However, despite the recent efforts, there is still room for further improvement, especially for describing molecular adsorption on metal surfaces. I will show that by choosing appropriate exchange and nonlocal correlation functionals, it is possible to calculate geometries and electronic structures for adsorption systems accurately within the framework of vdW-DF. Applicability of the present approach will be illustrated with its applications to graphene/metal, fullerene/metal, and water/graphene interfaces. This work is partly supported by a Grant-in-Aid for Scientific Research on Innovative Area (No. 23104501). AIMR was established by the World Premier International Research Center Initiative (WPI), MEXT, Japan.
A relativistic time-dependent density functional study of the excited states of the mercury dimer.
Kullie, Ossama
2014-01-14
In previous works on Zn2 and Cd2 dimers we found that the long-range corrected CAMB3LYP gives better results than other density functional approximations for the excited states, especially in the asymptotic region. In this paper, we use it to present a time-dependent density functional (TDDFT) study for the ground-state as well as the excited states corresponding to the (6s(2) + 6s6p), (6s(2) + 6s7s), and (6s(2) + 6s7p) atomic asymptotes for the mercury dimer Hg2. We analyze its spectrum obtained from all-electron calculations performed with the relativistic Dirac-Coulomb and relativistic spinfree Hamiltonian as implemented in DIRAC-PACKAGE. A comparison with the literature is given as far as available. Our result is excellent for the most of the lower excited states and very encouraging for the higher excited states, it shows generally good agreements with experimental results and outperforms other theoretical results. This enables us to give a detailed analysis of the spectrum of the Hg2 including a comparative analysis with the lighter dimers of the group 12, Cd2, and Zn2, especially for the relativistic effects, the spin-orbit interaction, and the performance of CAMB3LYP and is enlightened for similar systems. The result shows, as expected, that spinfree Hamiltonian is less efficient than Dirac-Coulomb Hamiltonian for systems containing heavy elements such as Hg2. PMID:24437874
A relativistic time-dependent density functional study of the excited states of the mercury dimer
NASA Astrophysics Data System (ADS)
Kullie, Ossama
2014-01-01
In previous works on Zn2 and Cd2 dimers we found that the long-range corrected CAMB3LYP gives better results than other density functional approximations for the excited states, especially in the asymptotic region. In this paper, we use it to present a time-dependent density functional (TDDFT) study for the ground-state as well as the excited states corresponding to the (6s2 + 6s6p), (6s2 + 6s7s), and (6s2 + 6s7p) atomic asymptotes for the mercury dimer Hg2. We analyze its spectrum obtained from all-electron calculations performed with the relativistic Dirac-Coulomb and relativistic spinfree Hamiltonian as implemented in DIRAC-PACKAGE. A comparison with the literature is given as far as available. Our result is excellent for the most of the lower excited states and very encouraging for the higher excited states, it shows generally good agreements with experimental results and outperforms other theoretical results. This enables us to give a detailed analysis of the spectrum of the Hg2 including a comparative analysis with the lighter dimers of the group 12, Cd2, and Zn2, especially for the relativistic effects, the spin-orbit interaction, and the performance of CAMB3LYP and is enlightened for similar systems. The result shows, as expected, that spinfree Hamiltonian is less efficient than Dirac-Coulomb Hamiltonian for systems containing heavy elements such as Hg2.
Dal Corso, A.; Pasquarello, A.; Baldereschi, A.; Car, R.
1996-01-01
Pseudopotential calculations within density-functional theory for a few selected solids (Si, GaAs, and Al) are used to assess the validity of two generalized-gradient approximations (GGA{close_quote}s), the one proposed by Becke and Perdew (BP) and the more recent one proposed by Perdew and Wang (PW) in comparison with the currently used local-density approximation (LDA). The GGA{close_quote}s give total energies of atoms and cohesive energies of solids that are closer to experiment than the LDA results. Lattice constants are reproduced with the same accuracy as in LDA, while bulk moduli and zone-center phonon frequencies are underestimated with respect to both the LDA and the experimental values. Comparison to all-electron results shows that in both GGA schemes the validity of the pseudopotential approach is as good as in LDA. The predictions of the two GGA{close_quote}s are similar, the PW functional yielding results marginally, but systematically, closer to experiment. The calculated values of the transition pressure in Si between the diamond and the {beta}-tin structure are 72 kbar (LDA), 164 kbar (BP), and 135 kbar (PW), to be compared with the two available experimental values of 103 and 125 kbar. {copyright} {ital 1996 The American Physical Society.}
Force Density Function Relationships in 2-D Granular Media
NASA Technical Reports Server (NTRS)
Youngquist, Robert C.; Metzger, Philip T.; Kilts, Kelly N.
2004-01-01
An integral transform relationship is developed to convert between two important probability density functions (distributions) used in the study of contact forces in granular physics. Developing this transform has now made it possible to compare and relate various theoretical approaches with one another and with the experimental data despite the fact that one may predict the Cartesian probability density and another the force magnitude probability density. Also, the transforms identify which functional forms are relevant to describe the probability density observed in nature, and so the modified Bessel function of the second kind has been identified as the relevant form for the Cartesian probability density corresponding to exponential forms in the force magnitude distribution. Furthermore, it is shown that this transform pair supplies a sufficient mathematical framework to describe the evolution of the force magnitude distribution under shearing. Apart from the choice of several coefficients, whose evolution of values must be explained in the physics, this framework successfully reproduces the features of the distribution that are taken to be an indicator of jamming and unjamming in a granular packing. Key words. Granular Physics, Probability Density Functions, Fourier Transforms
Improved association in a classical density functional theory for water.
Krebs, Eric J; Schulte, Jeff B; Roundy, David
2014-03-28
We present a modification to our recently published statistical associating fluid theory-based classical density functional theory for water. We have recently developed and tested a functional for the averaged radial distribution function at contact of the hard-sphere fluid that is dramatically more accurate at interfaces than earlier approximations. We now incorporate this improved functional into the association term of our free energy functional for water, improving its description of hydrogen bonding. We examine the effect of this improvement by studying two hard solutes (a hard hydrophobic rod and a hard sphere) and a Lennard-Jones approximation of a krypton atom solute. The improved functional leads to a moderate change in the density profile and a large decrease in the number of hydrogen bonds broken in the vicinity of the hard solutes. We find an improvement of the partial radial distribution for a krypton atom in water when compared with experiment. PMID:24697459
Improved association in a classical density functional theory for water
Krebs, Eric J.; Schulte, Jeff B.; Roundy, David
2014-03-28
We present a modification to our recently published statistical associating fluid theory-based classical density functional theory for water. We have recently developed and tested a functional for the averaged radial distribution function at contact of the hard-sphere fluid that is dramatically more accurate at interfaces than earlier approximations. We now incorporate this improved functional into the association term of our free energy functional for water, improving its description of hydrogen bonding. We examine the effect of this improvement by studying two hard solutes (a hard hydrophobic rod and a hard sphere) and a Lennard-Jones approximation of a krypton atom solute. The improved functional leads to a moderate change in the density profile and a large decrease in the number of hydrogen bonds broken in the vicinity of the hard solutes. We find an improvement of the partial radial distribution for a krypton atom in water when compared with experiment.
Rapid Bacterial Detection via an All-Electronic CMOS Biosensor.
Nikkhoo, Nasim; Cumby, Nichole; Gulak, P Glenn; Maxwell, Karen L
2016-01-01
The timely and accurate diagnosis of infectious diseases is one of the greatest challenges currently facing modern medicine. The development of innovative techniques for the rapid and accurate identification of bacterial pathogens in point-of-care facilities using low-cost, portable instruments is essential. We have developed a novel all-electronic biosensor that is able to identify bacteria in less than ten minutes. This technology exploits bacteriocins, protein toxins naturally produced by bacteria, as the selective biological detection element. The bacteriocins are integrated with an array of potassium-selective sensors in Complementary Metal Oxide Semiconductor technology to provide an inexpensive bacterial biosensor. An electronic platform connects the CMOS sensor to a computer for processing and real-time visualization. We have used this technology to successfully identify both Gram-positive and Gram-negative bacteria commonly found in human infections. PMID:27618185
Addressing spectroscopic quality of covariant density functional theory
NASA Astrophysics Data System (ADS)
Afanasjev, A. V.
2015-03-01
The spectroscopic quality of covariant density functional theory has been accessed by analyzing the accuracy and theoretical uncertainties in the description of spectroscopic observables. Such analysis is first presented for the energies of the single-particle states in spherical and deformed nuclei. It is also shown that the inclusion of particle-vibration coupling improves the description of the energies of predominantly single-particle states in medium and heavy-mass spherical nuclei. However, the remaining differences between theory and experiment clearly indicate missing physics and missing terms in covariant energy density functionals. The uncertainties in the predictions of the position of two-neutron drip line sensitively depend on the uncertainties in the prediction of the energies of the single-particle states. On the other hand, many spectroscopic observables in well deformed nuclei at ground state and finite spin only weakly depend on the choice of covariant energy density functional.
Kinetic-energy density functional: Atoms and shell structure
Garcia-Gonzalez, P.; Alvarellos, J.E.; Chacon, E. |
1996-09-01
We present a nonlocal kinetic-energy functional which includes an anisotropic average of the density through a symmetrization procedure. This functional allows a better description of the nonlocal effects of the electron system. The main consequence of the symmetrization is the appearance of a clear shell structure in the atomic density profiles, obtained after the minimization of the total energy. Although previous results with some of the nonlocal kinetic functionals have given incipient structures for heavy atoms, only our functional shows a clear shell structure for most of the atoms. The atomic total energies have a good agreement with the exact calculations. Discussion of the chemical potential and the first ionization potential in atoms is included. The functional is also extended to spin-polarized systems. {copyright} {ital 1996 The American Physical Society.}
Exploration of a modified density dependence in the Skyrme functional
Erler, J.; Reinhard, P.-G.; Kluepfel, P.
2010-10-15
A variant of the basic Skyrme-Hartree-Fock functional is considered dealing with a new form of density dependence. It employs only integer powers and thus will allow a more sound basis for projection schemes (particle number, angular momentum). We optimize the new functional with exactly the same adjustment strategy as used in an earlier study with a standard Skyrme functional. This allows direct comparisons of the performance of the new functional relative to the standard one. We discuss various observables: bulk properties of finite nuclei, nuclear matter, giant resonances, superheavy elements, and energy systematics. The new functional performs at least as well as the standard one, but offers a wider range of applicability (e.g., for projection) and more flexibility in the regime of high densities.
An Accurate Density Functional from Exchange-Correlation Hole
NASA Astrophysics Data System (ADS)
Tao, Jianmin; Mo, Yuxiang
The exchange-correlation hole is most fundamentally important in the development and understanding of density functional theory (DFT). However, due to the nonlocal nature of the exchange-correlation hole, development of DFT from the underlying hole presents a great challenge, and the works along this direction are limited. Here I will discuss a novel nonempirical DFT based on a semilocal hole, which is obtained from the density matrix expansion. Extensive tests on molecules and solids show that this functional can achieve remarkable accuracy for wide-ranging properties in condensed matter physics and quantum chemistry. This work was supported by NSF under Grant No. CHE-1261918.
Nonlinear eigenvalue problems in Density Functional Theory calculations
Fattebert, J
2009-08-28
Developed in the 1960's by W. Kohn and coauthors, Density Functional Theory (DFT) is a very popular quantum model for First-Principles simulations in chemistry and material sciences. It allows calculations of systems made of hundreds of atoms. Indeed DFT reduces the 3N-dimensional Schroedinger electronic structure problem to the search for a ground state electronic density in 3D. In practice it leads to the search for N electronic wave functions solutions of an energy minimization problem in 3D, or equivalently the solution of an eigenvalue problem with a non-linear operator.
Constrained Density Functional Calculations of alpha and delta Pu
NASA Astrophysics Data System (ADS)
Eriksson, Olle
2002-03-01
The electronic structure of α and δ Pu are described using a modified density functional theory that incorporates localization effects of the 5f shell. It is argued that a Russel-Saunders coupled state involving a 5f^4 multiplet, together with one itinerant 5f electron explains most of the observed ground state properties of δ Pu (equilibrium volume, elastic constants, near degeneracy with the alpha phase). This 5f electrons in the α phase are argued to form itinerant states, that are well described in density functional theory. The two distinctly different electronic ground states give rise to different excitation spectra and a comparison with experimental data is made.
Direct propagation of probability density functions in hydrological equations
NASA Astrophysics Data System (ADS)
Kunstmann, Harald; Kastens, Marko
2006-06-01
Sustainable decisions in hydrological risk management require detailed information on the probability density function ( pdf) of the model output. Only then probabilities for the failure of a specific management option or the exceedance of critical thresholds (e.g. of pollutants) can be derived. A new approach of uncertainty propagation in hydrological equations is developed that directly propagates the probability density functions of uncertain model input parameters into the corresponding probability density functions of model output. The basics of the methodology are presented and central applications to different disciplines in hydrology are shown. This work focuses on the following basic hydrological equations: (1) pumping test analysis (Theis-equation, propagation of uncertainties in recharge and transmissivity), (2) 1-dim groundwater contaminant transport equation (Gauss-equation, propagation of uncertainties in decay constant and dispersivity), (3) evapotranspiration estimation (Penman-Monteith-equation, propagation of uncertainty in roughness length). The direct propagation of probability densities is restricted to functions that are monotonically increasing or decreasing or that can be separated in corresponding monotonic branches so that inverse functions can be derived. In case no analytic solutions for inverse functions could be derived, semi-analytical approximations were used. It is shown that the results of direct probability density function propagation are in perfect agreement with results obtained from corresponding Monte Carlo derived frequency distributions. Direct pdf propagation, however, has the advantage that is yields exact solutions for the resulting hydrological pdfs rather than approximating discontinuous frequency distributions. It is additionally shown that the type of the resulting pdf depends on the specific values (order of magnitude, respectively) of the standard deviation of the input pdf. The dependency of skewness and kurtosis
Higher-order adaptive finite-element methods for Kohn–Sham density functional theory
Motamarri, P.; Nowak, M.R.; Leiter, K.; Knap, J.; Gavini, V.
2013-11-15
We present an efficient computational approach to perform real-space electronic structure calculations using an adaptive higher-order finite-element discretization of Kohn–Sham density-functional theory (DFT). To this end, we develop an a priori mesh-adaption technique to construct a close to optimal finite-element discretization of the problem. We further propose an efficient solution strategy for solving the discrete eigenvalue problem by using spectral finite-elements in conjunction with Gauss–Lobatto quadrature, and a Chebyshev acceleration technique for computing the occupied eigenspace. The proposed approach has been observed to provide a staggering 100–200-fold computational advantage over the solution of a generalized eigenvalue problem. Using the proposed solution procedure, we investigate the computational efficiency afforded by higher-order finite-element discretizations of the Kohn–Sham DFT problem. Our studies suggest that staggering computational savings—of the order of 1000-fold—relative to linear finite-elements can be realized, for both all-electron and local pseudopotential calculations, by using higher-order finite-element discretizations. On all the benchmark systems studied, we observe diminishing returns in computational savings beyond the sixth-order for accuracies commensurate with chemical accuracy, suggesting that the hexic spectral-element may be an optimal choice for the finite-element discretization of the Kohn–Sham DFT problem. A comparative study of the computational efficiency of the proposed higher-order finite-element discretizations suggests that the performance of finite-element basis is competing with the plane-wave discretization for non-periodic local pseudopotential calculations, and compares to the Gaussian basis for all-electron calculations to within an order of magnitude. Further, we demonstrate the capability of the proposed approach to compute the electronic structure of a metallic system containing 1688
Higher-order adaptive finite-element methods for Kohn-Sham density functional theory
NASA Astrophysics Data System (ADS)
Motamarri, P.; Nowak, M. R.; Leiter, K.; Knap, J.; Gavini, V.
2013-11-01
We present an efficient computational approach to perform real-space electronic structure calculations using an adaptive higher-order finite-element discretization of Kohn-Sham density-functional theory (DFT). To this end, we develop an a priori mesh-adaption technique to construct a close to optimal finite-element discretization of the problem. We further propose an efficient solution strategy for solving the discrete eigenvalue problem by using spectral finite-elements in conjunction with Gauss-Lobatto quadrature, and a Chebyshev acceleration technique for computing the occupied eigenspace. The proposed approach has been observed to provide a staggering 100-200-fold computational advantage over the solution of a generalized eigenvalue problem. Using the proposed solution procedure, we investigate the computational efficiency afforded by higher-order finite-element discretizations of the Kohn-Sham DFT problem. Our studies suggest that staggering computational savings-of the order of 1000-fold-relative to linear finite-elements can be realized, for both all-electron and local pseudopotential calculations, by using higher-order finite-element discretizations. On all the benchmark systems studied, we observe diminishing returns in computational savings beyond the sixth-order for accuracies commensurate with chemical accuracy, suggesting that the hexic spectral-element may be an optimal choice for the finite-element discretization of the Kohn-Sham DFT problem. A comparative study of the computational efficiency of the proposed higher-order finite-element discretizations suggests that the performance of finite-element basis is competing with the plane-wave discretization for non-periodic local pseudopotential calculations, and compares to the Gaussian basis for all-electron calculations to within an order of magnitude. Further, we demonstrate the capability of the proposed approach to compute the electronic structure of a metallic system containing 1688 atoms using
The exact density functional for two electrons in one dimension
NASA Astrophysics Data System (ADS)
Cohen, Aron; Mori-Sanchez, Paula
The exact universal density functional F [ ρ ] is calculated for real space two-electron densities in one dimension ρ (x) with a soft-Coulomb interaction. It is calculated by the Levy constrained search F [ ρ ] =minΨ-->ρ < Ψ | \\Tcirc +\\Vcircee | Ψ > over wavefunctions of a two-dimensional Hilbert space Ψ (x1 ,x2) --> ρ (x1) and can be directly visualized. We do an approximate constrained search via density matrices and a direct approximation to natural orbitals. This allows us to make an accurate approximation to the exact functional that is calculated using a search over potentials. We investigate the exact functional and the performance of many approximations on some of the most challenging electronic structure in two-electron systems, from strongly-correlated electron transfer to the description of a localized-delocalized transition. The exact Kohn-Sham potential, vs (x) , and exact Kohn-Sham eigenvalues, ɛi, are calculated and this allows us to discuss the band-gap problem versus the perspective of the exact density functional F [ ρ ] for all numbers of electrons. We calculate the derivative discontinuity of the exact functional in an example of a Mott-Insulator, one-dimensional stretched H2.
Bakosi, Jozsef; Ristorcelli, Raymond J
2010-01-01
Probability density function (PDF) methods are extended to variable-density pressure-gradient-driven turbulence. We apply the new method to compute the joint PDF of density and velocity in a non-premixed binary mixture of different-density molecularly mixing fluids under gravity. The full time-evolution of the joint PDF is captured in the highly non-equilibrium flow: starting from a quiescent state, transitioning to fully developed turbulence and finally dissipated by molecular diffusion. High-Atwood-number effects (as distinguished from the Boussinesq case) are accounted for: both hydrodynamic turbulence and material mixing are treated at arbitrary density ratios, with the specific volume, mass flux and all their correlations in closed form. An extension of the generalized Langevin model, originally developed for the Lagrangian fluid particle velocity in constant-density shear-driven turbulence, is constructed for variable-density pressure-gradient-driven flows. The persistent small-scale anisotropy, a fundamentally 'non-Kolmogorovian' feature of flows under external acceleration forces, is captured by a tensorial diffusion term based on the external body force. The material mixing model for the fluid density, an active scalar, is developed based on the beta distribution. The beta-PDF is shown to be capable of capturing the mixing asymmetry and that it can accurately represent the density through transition, in fully developed turbulence and in the decay process. The joint model for hydrodynamics and active material mixing yields a time-accurate evolution of the turbulent kinetic energy and Reynolds stress anisotropy without resorting to gradient diffusion hypotheses, and represents the mixing state by the density PDF itself, eliminating the need for dubious mixing measures. Direct numerical simulations of the homogeneous Rayleigh-Taylor instability are used for model validation.
On the Evolution of the Density Probability Density Function in Strongly Self-gravitating Systems
NASA Astrophysics Data System (ADS)
Girichidis, Philipp; Konstandin, Lukas; Whitworth, Anthony P.; Klessen, Ralf S.
2014-02-01
The time evolution of the probability density function (PDF) of the mass density is formulated and solved for systems in free-fall using a simple approximate function for the collapse of a sphere. We demonstrate that a pressure-free collapse results in a power-law tail on the high-density side of the PDF. The slope quickly asymptotes to the functional form PV (ρ)vpropρ-1.54 for the (volume-weighted) PDF and PM (ρ)vpropρ-0.54 for the corresponding mass-weighted distribution. From the simple approximation of the PDF we derive analytic descriptions for mass accretion, finding that dynamically quiet systems with narrow density PDFs lead to retarded star formation and low star formation rates (SFRs). Conversely, strong turbulent motions that broaden the PDF accelerate the collapse causing a bursting mode of star formation. Finally, we compare our theoretical work with observations. The measured SFRs are consistent with our model during the early phases of the collapse. Comparison of observed column density PDFs with those derived from our model suggests that observed star-forming cores are roughly in free-fall.
On the evolution of the density probability density function in strongly self-gravitating systems
Girichidis, Philipp; Konstandin, Lukas; Klessen, Ralf S.; Whitworth, Anthony P.
2014-02-01
The time evolution of the probability density function (PDF) of the mass density is formulated and solved for systems in free-fall using a simple approximate function for the collapse of a sphere. We demonstrate that a pressure-free collapse results in a power-law tail on the high-density side of the PDF. The slope quickly asymptotes to the functional form P{sub V} (ρ)∝ρ{sup –1.54} for the (volume-weighted) PDF and P{sub M} (ρ)∝ρ{sup –0.54} for the corresponding mass-weighted distribution. From the simple approximation of the PDF we derive analytic descriptions for mass accretion, finding that dynamically quiet systems with narrow density PDFs lead to retarded star formation and low star formation rates (SFRs). Conversely, strong turbulent motions that broaden the PDF accelerate the collapse causing a bursting mode of star formation. Finally, we compare our theoretical work with observations. The measured SFRs are consistent with our model during the early phases of the collapse. Comparison of observed column density PDFs with those derived from our model suggests that observed star-forming cores are roughly in free-fall.
All-electron scalar relativistic calculation of water molecule adsorption onto small gold clusters.
Kuang, Xiang-Jun; Wang, Xin-Qiang; Liu, Gao-Bin
2011-08-01
An all-electron scalar relativistic calculation was performed on Au( n )H(2)O (n = 1-13) clusters using density functional theory (DFT) with the generalized gradient approximation at PW91 level. The calculation results reveal that, after adsorption, the small gold cluster would like to bond with oxygen and the H(2)O molecule prefers to occupy the single fold coordination site. Reflecting the strong scalar relativistic effect, Au( n ) geometries are distorted slightly but still maintain a planar structure. The Au-Au bond is strengthened and the H-O bond is weakened, as manifested by the shortening of the Au-Au bond-length and the lengthening of the H-O bond-length. The H-O-H bond angle becomes slightly larger. The enhancement of reactivity of the H(2)O molecule is obvious. The Au-O bond-lengths, adsorption energies, VIPs, HLGs, HOMO (LUMO) energy levels, charge transfers and the highest vibrational frequencies of the Au-O mode for Au( n )H(2)O clusters exhibit an obvious odd-even oscillation. The most favorable adsorption between small gold clusters and the H(2)O molecule takes place when the H(2)O molecule is adsorbed onto an even-numbered Au( n ) cluster and becomes an Au( n )H(2)O cluster with an even number of valence electrons. The odd-even alteration of magnetic moments is observed in Au( n )H(2)O clusters and may serve as material with a tunable code capacity of "0" and "1" by adsorbing a H(2)O molecule onto an odd or even-numbered small gold cluster. PMID:21140279
Subsystem real-time time dependent density functional theory
NASA Astrophysics Data System (ADS)
Krishtal, Alisa; Ceresoli, Davide; Pavanello, Michele
2015-04-01
We present the extension of Frozen Density Embedding (FDE) formulation of subsystem Density Functional Theory (DFT) to real-time Time Dependent Density Functional Theory (rt-TDDFT). FDE is a DFT-in-DFT embedding method that allows to partition a larger Kohn-Sham system into a set of smaller, coupled Kohn-Sham systems. Additional to the computational advantage, FDE provides physical insight into the properties of embedded systems and the coupling interactions between them. The extension to rt-TDDFT is done straightforwardly by evolving the Kohn-Sham subsystems in time simultaneously, while updating the embedding potential between the systems at every time step. Two main applications are presented: the explicit excitation energy transfer in real time between subsystems is demonstrated for the case of the Na4 cluster and the effect of the embedding on optical spectra of coupled chromophores. In particular, the importance of including the full dynamic response in the embedding potential is demonstrated.
External orthogonality in subsystem time-dependent density functional theory.
Chulhai, Dhabih V; Jensen, Lasse
2016-08-01
Subsystem density functional theory (subsystem DFT) is a DFT partitioning method that is exact in principle, but depends on approximations to the kinetic energy density functional (KEDF). One may avoid the use of approximate KEDFs by ensuring that the inter-subsystem molecular orbitals are orthogonal, termed external orthogonality (EO). We present a method that extends a subsystem DFT method, that includes EO, into the time-dependent DFT (TDDFT) regime. This method therefore removes the need for approximations to the kinetic energy potential and kernel, and we show that it can accurately reproduce the supermolecular TDDFT results for weakly and strongly coupled subsystems, and for systems with strongly overlapping densities (where KEDF approximations traditionally fail). PMID:26932176
Applying Density Functional Theory for Atomic Vacancies in Solids
NASA Astrophysics Data System (ADS)
Zhou, Xiaolan; Perdew, John P.
2008-03-01
We use a new generalized gradient approximation of density functional theory -- PBEsol, a revised Perdew-Burke-Ernzerhof GGA, to calculate the vacancy formation energies and other properties of metals and semiconductors: Al, Pt, Pd, GaN [1], etc. By restoring the gradient expansion over a wide range of density gradients, PBEsol [2] yields excellent jellium exchange and correlation surface energies. We expect that this new functional will improve the description of vacancies in real materials, since the vacancy formation energy is essentially the work needed to create an interior surface. [1] Thomas R.Mattsson and Ann E. Mattson. Phys. Rev. B 66, 214410 (2002). [2] John P. Perdew, Adrienn Ruzsinszky, Gabor I. Csonka, Oleg A. Vydrov, Gustavo E. Scuseria, Lucian A. Constantin, Xiaolan Zhou, and Kieron Burke, Restoring the density-gradient expansion for exchange in solids and surfaces, http://arxiv.org/abs/0711.0156
Catalytic activities of platinum nanotubes: a density functional study
NASA Astrophysics Data System (ADS)
Mukherjee, Prajna; Gupta, Bikash C.; Jena, Puru
2015-10-01
In this work we investigate the catalytic properties of platinum nanotubes using density functional theory based calculations. In particular, we study the dissociation of hydrogen and oxygen molecules as well as oxidation of CO molecules. The results indicate that platinum nanotubes have good catalytic properties and can be effectively used in converting CO molecule to CO2.
NASA Astrophysics Data System (ADS)
Manzhos, Sergei
2016-01-01
A comparative DFTB (density functional tight binding)-DFT (density functional theory) study of several adsorption modes of 2-anthroic acid on titania is presented. Two parameterizations of DFTB previously used for dye-TiO2 interfaces are tested. DFTB predicts adsorption energies which differ from those computed by DFT not only in magnitude (by up to 0.5 eV) but also in the order among different configurations. The band alignment computed with DFTB is not consistent with DFT results and with experimental data. The strategy of geometry optimization with DFTB followed by single-point DFT calculations also does not necessarily result in plausible adsorption energies.
NASA Astrophysics Data System (ADS)
Vidya, R.; Ravindran, P.; Fjellvåg, H.; Svensson, B. G.; Monakhov, E.; Ganchenkova, M.; Nieminen, R. M.
2011-01-01
Formation energies of various intrinsic defects and defect complexes in ZnO have been calculated using a density-functional-theory-based pseudopotential all-electron method. The various defects considered are oxygen vacancy (VO), zinc vacancy (VZn), oxygen at an interstitial site (Oi), Zn at an interstitial site (Zni), Zn at VO (ZnO), O at VZn(OZn), and an antisite pair (combination of the preceding two defects). In addition, defect complexes like (VO+Zni) and Zn-vacancy clusters are studied. The Schokkty pair (VO+VZn) and Frenkel pairs [(VO+Oi) and (VZn+Zni)] are considered theoretically for the first time. Upon comparing the formation energies of these defects, we find that VO would be the dominant intrinsic defect under both Zn-rich and O-rich conditions and it is a deep double donor. Both ZnO and Zni are found to be shallow donors. The low formation energy of donor-type intrinsic defects could lead to difficulty in achieving p-type conductivity in ZnO. Defect complexes have charge transitions deep inside the band gap. The red, yellow, and green photoluminescence peaks of undoped samples can be assigned to some of the defect complexes considered. It is believed that the red luminescence originates from an electronic transition in VO, but we find that it can originate from the antisite ZnO defect. Charge density and electron-localization function analyses have been used to understand the effect of these defects on the ZnO lattice. The electronic structure of ZnO with intrinsic defects has been studied using density-of-states and electronic band structure plots. The acceptor levels introduced by VZn are relatively localized, making it difficult to achieve p-type conductivity with sufficient hole mobility.
Reduced density-matrix functionals applied to the Hubbard dimer
NASA Astrophysics Data System (ADS)
Kamil, Ebad; Schade, Robert; Pruschke, Thomas; Blöchl, Peter E.
2016-02-01
Common density-matrix functionals, the Müller and the power functional, have been benchmarked for the half-filled Hubbard dimer, which allows us to model the bond dissociation problem and the transition from the weakly to the strongly correlated limit. Unbiased numerical calculations are combined with analytical results. Despite the well known successes of the Müller functional, the ground state is degenerate with a one-dimensional manifold of ferromagnetic solutions. The resulting infinite magnetic susceptibility indicates another qualitative flaw of the Müller functional. The derivative discontinuity with respect to particle number is not present indicating an incorrect metal-like behavior. The power functional actually favors the ferromagnetic state for weak interaction. Analogous to the Hartree-Fock approximation, the power functional undergoes a transition beyond a critical interaction strength, in this case, however, to a noncollinear antiferromagnetic state.
Density Functional Study of Perovskite Superconductor MgCNi3
NASA Astrophysics Data System (ADS)
Kumar, Jagdish; Sharma, Devina; Kumar, Ranjan; Awana, V. P. S.; Ahluwalia, P. K.
2011-12-01
We here report the first principle density functional study of MgCNi3 which crystallize in cubic perovskite structure having critical transition temperature of 8 K. The interesting aspect of this compound is that in normal state it is non magnetic in nature despite conduction electrons in it are derived from partially filled Ni d states, which typically lead to ferromagnetism in metallic Ni and many Ni-based binary alloys. To investigate the detailed microscopic origin of the non magnetic nature we have done density functional based calculations on this compound. The lattice constant is calculated using minimum energy criteria from total energy versus lattice constant plot. By taking the calculated values of lattice constant we have done the precise calculations on the compound using Full Potential Linear Augmented Plane Wave (FP-LAPW) method implemented in ELK code. The electronic density of states is found spin degenerate that corresponds to a non-magnetic ground state. The density of states (DOS) at Fermi level, N(EF) is dominated by Ni-d states. The sharp peak observed just below Fermi level corresponds to van Hove singularity (vHs). The projected density of states (PDOS) suggests a strong hybridization of Ni-3d and C-2p states which is responsible for the observed non magnetic nature of MgCNi3.
Density functional calculations of point defects in InAs
NASA Astrophysics Data System (ADS)
Moussa, Jonathan; Schultz, Peter
2013-03-01
Standard semilocal density functionals do not generate a gap in the Kohn-Sham eigenvalues for InAs, a semiconductor with an experimental gap of 0.4 eV. Without a theoretical band gap, it becomes difficult to identify, specify, and characterize pure localized states of point defects with energy levels within the experimental band gap. The bulk band gap problem can be alleviated with screened hybrid density functionals, such as the Heyd-Scuseria-Ernzerhof (HSE) functional, that open the generalized Kohn-Sham eigenvalue gap of InAs to near the experimental value. However, even without a Kohn-Sham gap, the local moment countercharge (LMCC) method [Phys. Rev. Lett. 96, 246401 (2006)] is able to predict charge transition energy levels of localized defect states, using standard semi-local functionals. We present an LMCC-based study of standard point defects in InAs using semilocal density functionals and compare the results to HSE-based calculations to assess the validity of LMCC calculations in this situation. 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. DOE's National Nuclear Security Administration under contract DE-AC04-94AL85000.
Density Functional Theory of Biologically Relevant Metal Centers
NASA Astrophysics Data System (ADS)
Siegbahn, Per E. M.; Blomberg, Margareta R. A.
1999-10-01
Recent applications of density functional theory to biologically relevant metal centers are reviewed. The emphasis is on reaction mechanisms, structures, and modeling. The accuracy of different functionals is discussed for standard benchmark tests of first- and second-row molecules and for transition metal systems. Modeling aspects of the protein metal complexes are discussed regarding both the size of the model being treated quantum mechanically and the treatment of the protein surrounding it. To illustrate the effects, structures computed without the effects of the protein are compared with experimental structures from enzymes, and results from simple dielectric models of the protein for electron transfer processes are described. The choice of spin state is discussed for multimetal complexes. Examples of mechanisms studied recently by density functional theory are described, such as O2 and methane activation in methane monooxygenase and O2 formation in photosystem II.
Ab initio derivation of model energy density functionals
NASA Astrophysics Data System (ADS)
Dobaczewski, Jacek
2016-08-01
I propose a simple and manageable method that allows for deriving coupling constants of model energy density functionals (EDFs) directly from ab initio calculations performed for finite fermion systems. A proof-of-principle application allows for linking properties of finite nuclei, determined by using the nuclear nonlocal Gogny functional, to the coupling constants of the quasilocal Skyrme functional. The method does not rely on properties of infinite fermion systems but on the ab initio calculations in finite systems. It also allows for quantifying merits of different model EDFs in describing the ab initio results.
Perspective: Fundamental aspects of time-dependent density functional theory
NASA Astrophysics Data System (ADS)
Maitra, Neepa T.
2016-06-01
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.
Partial ionization in dense plasmas: comparisons among average-atom density functional models.
Murillo, Michael S; Weisheit, Jon; Hansen, Stephanie B; Dharma-wardana, M W C
2013-06-01
Nuclei interacting with electrons in dense plasmas acquire electronic bound states, modify continuum states, generate resonances and hopping electron states, and generate short-range ionic order. The mean ionization state (MIS), i.e, the mean charge Z of an average ion in such plasmas, is a valuable concept: Pseudopotentials, pair-distribution functions, equations of state, transport properties, energy-relaxation rates, opacity, radiative processes, etc., can all be formulated using the MIS of the plasma more concisely than with an all-electron description. However, the MIS does not have a unique definition and is used and defined differently in different statistical models of plasmas. Here, using the MIS formulations of several average-atom models based on density functional theory, we compare numerical results for Be, Al, and Cu plasmas for conditions inclusive of incomplete atomic ionization and partial electron degeneracy. By contrasting modern orbital-based models with orbital-free Thomas-Fermi models, we quantify the effects of shell structure, continuum resonances, the role of exchange and correlation, and the effects of different choices of the fundamental cell and boundary conditions. Finally, the role of the MIS in plasma applications is illustrated in the context of x-ray Thomson scattering in warm dense matter. PMID:23848795
Gonis, A.; Zhang, X. G.; Stocks, G. M.; Nicholson, D. M.
2015-10-23
Density functional theory for the case of general, N-representable densities is reformulated in terms of density functional derivatives of expectation values of operators evaluated with wave functions leading to a density, making no reference to the concept of potential. The developments provide a complete solution of the v-representability problem by establishing a mathematical procedure that determines whether a density is v-representable and in the case of an affirmative answer determines the potential (within an additive constant) as a derivative with respect to the density of a constrained search functional. It also establishes the existence of an energy functional of the density that, for v-representable densities, assumes its minimum value at the density describing the ground state of an interacting many-particle system. The theorems of Hohenberg and Kohn emerge as special cases of the formalism.
Gonis, A.; Zhang, X. G.; Stocks, G. M.; Nicholson, D. M.
2015-10-23
Density functional theory for the case of general, N-representable densities is reformulated in terms of density functional derivatives of expectation values of operators evaluated with wave functions leading to a density, making no reference to the concept of potential. The developments provide a complete solution of the v-representability problem by establishing a mathematical procedure that determines whether a density is v-representable and in the case of an affirmative answer determines the potential (within an additive constant) as a derivative with respect to the density of a constrained search functional. It also establishes the existence of an energy functional of themore » density that, for v-representable densities, assumes its minimum value at the density describing the ground state of an interacting many-particle system. The theorems of Hohenberg and Kohn emerge as special cases of the formalism.« less
Modal density function and number of propagating modes in ducts
NASA Technical Reports Server (NTRS)
Rice, E. J.
1976-01-01
The question of the number of propagating modes within a small range of mode cut off ratio was raised. The population density of modes were shown to be greatest near cut off and least for the well propagating modes. It was shown that modes of nearly the same cut off ratio behave nearly the same in a sound absorbing duct as well as in the way they propagate to the far. Handling all of the propagating modes individually, they can be grouped into several cut off ratio ranges. It is important to know the modal density function to estimate acoustic power distribution.
Numerical density-to-potential inversions in time-dependent density functional theory.
Jensen, Daniel S; Wasserman, Adam
2016-08-01
We treat the density-to-potential inverse problem of time-dependent density functional theory as an optimization problem with a partial differential equation constraint. The unknown potential is recovered from a target density by applying a multilevel optimization method controlled by error estimates. We employ a classical optimization routine using gradients efficiently computed by the discrete adjoint method. The inverted potential has both a real and imaginary part to reduce reflections at the boundaries and other numerical artifacts. We demonstrate this method on model one-dimensional systems. The method can be straightforwardly extended to a variety of numerical solvers of the time-dependent Kohn-Sham equations and to systems in higher dimensions. PMID:26984427
SURFACE SYMMETRY ENERGY OF NUCLEAR ENERGY DENSITY FUNCTIONALS
Nikolov, N; Schunck, N; Nazarewicz, W; Bender, M; Pei, J
2010-12-20
We study the bulk deformation properties of the Skyrme nuclear energy density functionals. Following simple arguments based on the leptodermous expansion and liquid drop model, we apply the nuclear density functional theory to assess the role of the surface symmetry energy in nuclei. To this end, we validate the commonly used functional parametrizations against the data on excitation energies of superdeformed band-heads in Hg and Pb isotopes, and fission isomers in actinide nuclei. After subtracting shell effects, the results of our self-consistent calculations are consistent with macroscopic arguments and indicate that experimental data on strongly deformed configurations in neutron-rich nuclei are essential for optimizing future nuclear energy density functionals. The resulting survey provides a useful benchmark for further theoretical improvements. Unlike in nuclei close to the stability valley, whose macroscopic deformability hangs on the balance of surface and Coulomb terms, the deformability of neutron-rich nuclei strongly depends on the surface-symmetry energy; hence, its proper determination is crucial for the stability of deformed phases of the neutron-rich matter and description of fission rates for r-process nucleosynthesis.
Steady-State Density Functional Theory for Finite Bias Conductances.
Stefanucci, G; Kurth, S
2015-12-01
In the framework of density functional theory, a formalism to describe electronic transport in the steady state is proposed which uses the density on the junction and the steady current as basic variables. We prove that, in a finite window around zero bias, there is a one-to-one map between the basic variables and both local potential on as well as bias across the junction. The resulting Kohn-Sham system features two exchange-correlation (xc) potentials, a local xc potential, and an xc contribution to the bias. For weakly coupled junctions the xc potentials exhibit steps in the density-current plane which are shown to be crucial to describe the Coulomb blockade diamonds. At small currents these steps emerge as the equilibrium xc discontinuity bifurcates. The formalism is applied to a model benzene junction, finding perfect agreement with the orthodox theory of Coulomb blockade. PMID:26571349
Computational complexity of time-dependent density functional theory
NASA Astrophysics Data System (ADS)
Whitfield, J. D.; Yung, M.-H.; Tempel, D. G.; Boixo, S.; Aspuru-Guzik, A.
2014-08-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.
Scrutinizing "Invisible" astatine: A challenge for modern density functionals.
Sergentu, Dumitru-Claudiu; David, Grégoire; Montavon, Gilles; Maurice, Rémi; Galland, Nicolas
2016-06-01
The main-group 6p elements did not receive much attention in the development of recent density functionals. In many cases it is still difficult to choose among the modern ones a relevant functional for various applications. Here, we illustrate the case of astatine species (At, Z = 85) and we report the first, and quite complete, benchmark study on several properties concerning such species. Insights on geometries, transition energies and thermodynamic properties of a set of 19 astatine species, for which reference experimental or theoretical data has been reported, are obtained with relativistic (two-component) density functional theory calculations. An extensive set of widely used functionals is employed. The hybrid meta-generalized gradient approximation (meta-GGA) PW6B95 functional is overall the best choice. It is worth noting that the range-separated HSE06 functional as well as the old and very popular B3LYP and PBE0 hybrid-GGAs appear to perform quite well too. Moreover, we found that astatine chemistry in solution can accurately be predicted using implicit solvent models, provided that specific parameters are used to build At cavities. © 2016 Wiley Periodicals, Inc. PMID:27059181
What Density Functional Theory could do for Quantum Information
NASA Astrophysics Data System (ADS)
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.
Dynamical density functional theory for colloidal dispersions including hydrodynamic interactions
NASA Astrophysics Data System (ADS)
Rex, M.; Löwen, H.
2009-02-01
A dynamical density functional theory (DDFT) for translational Brownian dynamics is derived which includes hydrodynamic interactions. The theory reduces to the simple Brownian DDFT proposed by Marconi and Tarazona (U. Marini Bettolo Marconi and P. Tarazona, J. Chem. Phys. 110, 8032 (1999); J. Phys.: Condens. Matter 12, A413 (2000)) when hydrodynamic interactions are neglected. The derivation is based on Smoluchowski’s equation for the time evolution of the probability density with pairwise hydrodynamic interactions. The theory is applied to hard-sphere colloids in an oscillating spherical optical trap which switches periodically in time from a stable confining to an unstable potential. Rosenfeld’s fundamental measure theory for the equilibrium density functional is used and hydrodynamics are incorporated on the Rotne-Prager level. The results for the time-dependent density profiles are compared to extensive Brownian dynamics simulations which are performed on the same Rotne-Prager level and excellent agreement is obtained. It is further found that hydrodynamic interactions damp and slow the dynamics of the confined colloid cluster in comparison to the same situation with neglected hydrodynamic interactions.
Precise effective masses from density functional perturbation theory
NASA Astrophysics Data System (ADS)
Laflamme Janssen, J.; Gillet, Y.; Poncé, S.; Martin, A.; Torrent, M.; Gonze, X.
2016-05-01
The knowledge of effective masses is a key ingredient to analyze numerous properties of semiconductors, like carrier mobilities, (magneto)transport properties, or band extrema characteristics yielding carrier densities and density of states. Currently, these masses are usually calculated using finite-difference estimation of density functional theory (DFT) electronic band curvatures. However, finite differences require an additional convergence study and are prone to numerical noise. Moreover, the concept of effective mass breaks down at degenerate band extrema. We assess the former limitation by developing a method that allows to obtain the Hessian of DFT bands directly, using density functional perturbation theory. Then, we solve the latter issue by adapting the concept of "transport equivalent effective mass" to the k .p ̂ framework. The numerical noise inherent to finite-difference methods is thus eliminated, along with the associated convergence study. The resulting method is therefore more general, more robust, and simpler to use, which makes it especially appropriate for high-throughput computing. After validating the developed techniques, we apply them to the study of silicon, graphane, and arsenic. The formalism is implemented into the abinit software and supports the norm-conserving pseudopotential approach, the projector augmented-wave method, and the inclusion of spin-orbit coupling. The derived expressions also apply to the ultrasoft pseudopotential method.
Density Functional Theory for Steady-State Nonequilibrium Molecular Junctions
NASA Astrophysics Data System (ADS)
Liu, Shuanglong; Nurbawono, Argo; Zhang, Chun
2015-10-01
We present a density functional theory (DFT) for steady-state nonequilibrium quantum systems such as molecular junctions under a finite bias. Based on the steady-state nonequilibrium statistics that maps nonequilibrium to an effective equilibrium, we show that ground-state DFT (GS-DFT) is not applicable in this case and two densities, the total electron density and the density of current-carrying electrons, are needed to uniquely determine the properties of the corresponding nonequilibrium system. A self-consistent mean-field approach based on two densities is then derived. The theory is implemented into SIESTA computational package and applied to study nonequilibrium electronic/transport properties of a realistic carbon-nanotube (CNT)/Benzene junction. Results obtained from our steady-state DFT (SS-DFT) are compared with those of conventional GS-DFT based transport calculations. We show that SS-DFT yields energetically more stable nonequilibrium steady state, predicts significantly lower electric current, and is able to produce correct electronic structures in local equilibrium under a limiting case.
Semilocal density functional theory with correct surface asymptotics
NASA Astrophysics Data System (ADS)
Constantin, Lucian A.; Fabiano, Eduardo; Pitarke, J. M.; Della Sala, Fabio
2016-03-01
Semilocal density functional theory is the most used computational method for electronic structure calculations in theoretical solid-state physics and quantum chemistry of large systems, providing good accuracy with a very attractive computational cost. Nevertheless, because of the nonlocality of the exchange-correlation hole outside a metal surface, it was always considered inappropriate to describe the correct surface asymptotics. Here, we derive, within the semilocal density functional theory formalism, an exact condition for the imagelike surface asymptotics of both the exchange-correlation energy per particle and potential. We show that this condition can be easily incorporated into a practical computational tool, at the simple meta-generalized-gradient approximation level of theory. Using this tool, we also show that the Airy-gas model exhibits asymptotic properties that are closely related to those at metal surfaces. This result highlights the relevance of the linear effective potential model to the metal surface asymptotics.
Excitations and benchmark ensemble density functional theory for two electrons
Pribram-Jones, Aurora; Burke, Kieron; Yang, Zeng-hui; Ullrich, Carsten A.; Trail, John R.; Needs, Richard J.
2014-05-14
A new method for extracting ensemble Kohn-Sham potentials from accurate excited state densities is applied to a variety of two-electron systems, exploring the behavior of exact ensemble density functional theory. The issue of separating the Hartree energy and the choice of degenerate eigenstates is explored. A new approximation, spin eigenstate Hartree-exchange, is derived. Exact conditions that are proven include the signs of the correlation energy components and the asymptotic behavior of the potential for small weights of the excited states. Many energy components are given as a function of the weights for two electrons in a one-dimensional flat box, in a box with a large barrier to create charge transfer excitations, in a three-dimensional harmonic well (Hooke's atom), and for the He atom singlet-triplet ensemble, singlet-triplet-singlet ensemble, and triplet bi-ensemble.
Analyzing the financial crisis using the entropy density function
NASA Astrophysics Data System (ADS)
Oh, Gabjin; Kim, Ho-yong; Ahn, Seok-Won; Kwak, Wooseop
2015-02-01
The risk that is created by nonlinear interactions among subjects in economic systems is assumed to increase during an abnormal state of a financial market. Nevertheless, investigating the systemic risk in financial markets following the global financial crisis is not sufficient. In this paper, we analyze the entropy density function in the return time series for several financial markets, such as the S&P500, KOSPI, and DAX indices, from October 2002 to December 2011 and analyze the variability in the entropy value over time. We find that the entropy density function of the S&P500 index during the subprime crisis exhibits a significant decrease compared to that in other periods, whereas the other markets, such as those in Germany and Korea, exhibit no significant decrease during the market crisis. These findings demonstrate that the S&P500 index generated a regular pattern in the return time series during the financial crisis.
Covariant density functional theory: The role of the pion
Lalazissis, G. A.; Karatzikos, S.; Serra, M.; Otsuka, T.; Ring, P.
2009-10-15
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 {sigma} 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.
Neutrinoless double-beta decay in covariant density functional theory
Ring, P.; Yao, J. M.; Song, L. S.; Hagino, K.; Meng, J.
2015-10-15
We use covariant density functional theory beyond mean field in order to describe neutrinoless double-beta decay in a fully relativistic way. The dynamic effects of particle-number and angular-momentum conservations as well as shape fluctuations of quadrupole character are taken into account within the generator coordinate method for both initial and final nuclei. The calculations are based on the full relativistic transition operator. The nuclear matrix elements (NME’s) for a large number of possible transitions are investigated. The results are compared with various non-relativistic calculations, in particular also with the density functional theory based on the Gogny force. We find that the non-relativistic approximation is justified and that the total NME’s can be well approximated by the pure axial-vector coupling term. This corresponds to a considerable reduction of the computational effort.
Density-functional theory of inhomogeneous systems of hard spherocylinders
Velasco; Mederos; Sullivan
2000-09-01
The smectic-A phase boundaries of a hard-spherocylinder fluid are calculated using a density-functional theory based on one proposed earlier by Somoza and Tarazona [Phys. Rev. A 41, 965 (1990)]. Our calculations do not employ the translation-rotation decoupling approximation used in previous density-functional theories. The calculated phase boundaries agree well with computer simulation results up to aspect ratios L/D approximately 5 and are in better agreement with the simulations than are previous theories. We generalize the model fluid by including long-range interactions with quadrupolar orientational symmetry, which are taken into account by mean-field approximation. For sufficiently large strength, these interactions produce a smectic-C phase, which undergoes either a continuous or weakly first-order transition to the smectic-A phase. The theory and numerical methods discussed here can be applied to the analysis of interfacial phenomena. PMID:11088887
Neutrinoless double-beta decay in covariant density functional theory
NASA Astrophysics Data System (ADS)
Ring, P.; Yao, J. M.; Song, L. S.; Hagino, K.; Meng, J.
2015-10-01
We use covariant density functional theory beyond mean field in order to describe neutrinoless double-beta decay in a fully relativistic way. The dynamic effects of particle-number and angular-momentum conservations as well as shape fluctuations of quadrupole character are taken into account within the generator coordinate method for both initial and final nuclei. The calculations are based on the full relativistic transition operator. The nuclear matrix elements (NME's) for a large number of possible transitions are investigated. The results are compared with various non-relativistic calculations, in particular also with the density functional theory based on the Gogny force. We find that the non-relativistic approximation is justified and that the total NME's can be well approximated by the pure axial-vector coupling term. This corresponds to a considerable reduction of the computational effort.
Density versus spin-density functional in DFT+U and DFT+DMFT
NASA Astrophysics Data System (ADS)
Park, Hyowon; Millis, Andrew; Marianetti, Chris
2015-03-01
The construction of multi-variable effective action theories such as DFT+U and DFT+DMFT requires the choice of a local subspace of correlated orbitals and an additional variable being either the charge density or spin density. This talk examines the differences between using charge-only and spin-dependent exchange-correlation functionals with the aim of providing guidance for constructing more sophisticated beyond-density functional theories. The widely used spin-dependent approximations to the exchange-correlation functional are found to lead to a large and in some cases unphysical effective exchange coupling within the correlated subspace. Additionally, the differences between Wannier and Projector based definitions of the correlated orbitals are examined, and only small differences are found provided that the orbitals are orthonormal and strongly localized. These results are documented in the context of the rare earth nickelates. This work is supported under the Grant DOE-ER-046169 and under the FAME grant, one of six centers of STARnet, a Semiconductor Research Corporation program sponsored by MARCO and DARPA.
Ś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. PMID:25903880
Autoionization in time-dependent density-functional theory
NASA Astrophysics Data System (ADS)
Kapoor, V.
2016-06-01
We compute the exact exchange-correlation potential of the time-dependent density-functional theory (TDDFT) for the correlated process of autoionization. The potential develops barriers which regulate the autoionization rate. TDDFT employing known and practicable exchange-correlation potentials does not capture any autoionization dynamics. Approximate exchange-correlation potentials capturing such dynamics would necessarily require memory effects and are unlikely to be developed, as will be illustrated.
Density functional theory across chemistry, physics and biology
van Mourik, Tanja; Bühl, Michael; Gaigeot, Marie-Pierre
2014-01-01
The past decades have seen density functional theory (DFT) evolve from a rising star in computational quantum chemistry to one of its major players. This Theme Issue, which comes half a century after the publication of the Hohenberg–Kohn theorems that laid the foundations of modern DFT, reviews progress and challenges in present-day DFT research. Rather than trying to be comprehensive, this Theme Issue attempts to give a flavour of selected aspects of DFT. PMID:24516181
Association between lung function and airway wall density
NASA Astrophysics Data System (ADS)
Leader, J. Ken; Zheng, Bin; Fuhrman, Carl R.; Tedrow, John; Park, Sang C.; Tan, Jun; Pu, Jiantao; Drescher, John M.; Gur, David; Sciurba, Frank C.
2009-02-01
Computed tomography (CT) examination is often used to quantify the relation between lung function and airway remodeling in chronic obstructive pulmonary disease (COPD). In this preliminary study, we examined the association between lung function and airway wall computed attenuation ("density") in 200 COPD screening subjects. Percent predicted FVC (FVC%), percent predicted FEV1 (FEV1%), and the ratio of FEV1 to FVC as a percentage (FEV1/FVC%) were measured post-bronchodilator. The apical bronchus of the right upper lobe was manually selected from CT examinations for evaluation. Total airway area, lumen area, wall area, lumen perimeter and wall area as fraction of the total airway area were computed. Mean HU (meanHU) and maximum HU (maxHU) values were computed across pixels assigned membership in the wall and with a HU value greater than -550. The Pearson correlation coefficients (PCC) between FVC%, FEV1%, and FEV1/FVC% and meanHU were -0.221 (p = 0.002), -0.175 (p = 0.014), and -0.110 (p = 0.123), respectively. The PCCs for maxHU were only significant for FVC%. The correlations between lung function and the airway morphometry parameters were slightly stronger compared to airway wall density. MeanHU was significantly correlated with wall area (PCC = 0.720), airway area (0.498) and wall area percent (0.611). This preliminary work demonstrates that airway wall density is associated with lung function. Although the correlations in our study were weaker than a recent study, airway wall density initially appears to be an important parameter in quantitative CT analysis of COPD.
Probability Density Function Method for Langevin Equations with Colored Noise
Wang, Peng; Tartakovsky, Alexandre M.; Tartakovsky, Daniel M.
2013-04-05
We present a novel method to derive closed-form, computable PDF equations for Langevin systems with colored noise. The derived equations govern the dynamics of joint or marginal probability density functions (PDFs) of state variables, and rely on a so-called Large-Eddy-Diffusivity (LED) closure. We demonstrate the accuracy of the proposed PDF method for linear and nonlinear Langevin equations, describing the classical Brownian displacement and dispersion in porous media.
Density functional theory across chemistry, physics and biology.
van Mourik, Tanja; Bühl, Michael; Gaigeot, Marie-Pierre
2014-03-13
The past decades have seen density functional theory (DFT) evolve from a rising star in computational quantum chemistry to one of its major players. This Theme Issue, which comes half a century after the publication of the Hohenberg-Kohn theorems that laid the foundations of modern DFT, reviews progress and challenges in present-day DFT research. Rather than trying to be comprehensive, this Theme Issue attempts to give a flavour of selected aspects of DFT. PMID:24516181
Probability density function modeling for sub-powered interconnects
NASA Astrophysics Data System (ADS)
Pater, Flavius; Amaricǎi, Alexandru
2016-06-01
This paper proposes three mathematical models for reliability probability density function modeling the interconnect supplied at sub-threshold voltages: spline curve approximations, Gaussian models,and sine interpolation. The proposed analysis aims at determining the most appropriate fitting for the switching delay - probability of correct switching for sub-powered interconnects. We compare the three mathematical models with the Monte-Carlo simulations of interconnects for 45 nm CMOS technology supplied at 0.25V.
Next generation high density self assembling functional protein arrays
Ramachandran, Niroshan; Raphael, Jacob V.; Hainsworth, Eugenie; Demirkan, Gokhan; Fuentes, Manuel G.; Rolfs, Andreas; Hu, Yanhui; LaBaer, Joshua
2009-01-01
We report a high-density self assembling protein microarray that displays thousands of proteins, produced and captured in situ from immobilized cDNA templates. Over 1500 unique cDNAs were tested with > 90% success with nearly all proteins displaying yields within 2 fold of the mean, minimal sample variation and good day to day reproducibility. The displayed proteins revealed selective protein interactions. This method will enable various experimental approaches to study protein function in high throughput. PMID:18469824
Plato: A localised orbital based density functional theory code
NASA Astrophysics Data System (ADS)
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
Advances in time-dependent current-density functional theory
NASA Astrophysics Data System (ADS)
Berger, Arjan
In this work we solve the problem of the gauge dependence of molecular magnetic properties (magnetizabilities, circular dichroism) using time-dependent current-density functional theory [1]. We also present a new functional that accurately describes the optical absorption spectra of insulators, semiconductors and metals [2] N. Raimbault, P.L. de Boeij, P. Romaniello, and J.A. Berger Phys. Rev. Lett. 114, 066404 (2015) J.A. Berger, Phys. Rev. Lett. 115, 137402 (2015) This study has been partially supported through the Grant NEXT No. ANR-10-LABX-0037 in the framework of the Programme des Investissements d'Avenir.
Thermal Corrections to Density Functional Simulations of Warm Dense Matter
NASA Astrophysics Data System (ADS)
Smith, Justin; Pribram-Jones, Aurora; Burke, Kieron
Present density functional calculations of warm dense matter often use the Mermin-Kohn-Sham (MKS) scheme at finite temperature, but employ ground-state approximations to the exchange-correlation (XC) free energy. In the simplest solvable non-trivial model, an asymmetric Hubbard dimer, we calculate the exact many-body energies, the exact Mermin-Kohn-Sham functionals for this system, and extract the exact XC free energy. For moderate temperatures and weak correlation, we show this approximation is excellent, but fails for stronger correlations. Additionally, we use this system to test various conditions that must be satisfied.
Linear Scaling Density Functional Calculations with Gaussian Orbitals
NASA Technical Reports Server (NTRS)
Scuseria, Gustavo E.
1999-01-01
Recent advances in linear scaling algorithms that circumvent the computational bottlenecks of large-scale electronic structure simulations make it possible to carry out density functional calculations with Gaussian orbitals on molecules containing more than 1000 atoms and 15000 basis functions using current workstations and personal computers. This paper discusses the recent theoretical developments that have led to these advances and demonstrates in a series of benchmark calculations the present capabilities of state-of-the-art computational quantum chemistry programs for the prediction of molecular structure and properties.
Conformational Analysis of Thioether Musks Using Density Functional Theory
Setzer, William N.
2009-01-01
A conformational analysis of nine macrocyclic thioether musks has been carried out using molecular mechanics (MMFF), density functional theory (DFT) using both B3LYP and M06 functionals, as well as Hartree-Fock and post-Hartree-Fock (MP2) ab initio methods. 6-Thia-, 10-thia- and 4-methyl-5-thia-14-tetradecananolide, 4-thia-, 7-thia-, 11-thia- and 12-thia-15-pentadecanolide and 6-thia- and 12-thia-16-hexadecanolide were modeled. Unfortunately, there was little agreement between the computational methods at the levels of theory used in this study. PMID:20111690
Atomistic force field for alumina fit to density functional theory
Sarsam, Joanne; Finnis, Michael W.; Tangney, Paul
2013-11-28
We present a force field for bulk alumina (Al{sub 2}O{sub 3}), which has been parametrized by fitting the energies, forces, and stresses of a large database of reference configurations to those calculated with density functional theory (DFT). We use a functional form that is simpler and computationally more efficient than some existing models of alumina parametrized by a similar technique. Nevertheless, we demonstrate an accuracy of our potential that is comparable to those existing models and to DFT. We present calculations of crystal structures and energies, elastic constants, phonon spectra, thermal expansion, and point defect formation energies.
Franco-Pérez, Marco; Ayers, Paul W; Gázquez, José L; Vela, Alberto
2015-12-28
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. PMID:26723661
Franco-Pérez, Marco E-mail: ayers@mcmaster.ca E-mail: avela@cinvestav.mx; Ayers, Paul W. E-mail: ayers@mcmaster.ca E-mail: avela@cinvestav.mx; Gázquez, José L. E-mail: ayers@mcmaster.ca E-mail: avela@cinvestav.mx; Vela, Alberto E-mail: ayers@mcmaster.ca E-mail: avela@cinvestav.mx
2015-12-28
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.
A Wigner Monte Carlo approach to density functional theory
NASA Astrophysics Data System (ADS)
Sellier, J. M.; Dimov, I.
2014-08-01
In order to simulate quantum N-body systems, stationary and time-dependent density functional theories rely on the capacity of calculating the single-electron wave-functions of a system from which one obtains the total electron density (Kohn-Sham systems). In this paper, we introduce the use of the Wigner Monte Carlo method in ab-initio calculations. This approach allows time-dependent simulations of chemical systems in the presence of reflective and absorbing boundary conditions. It also enables an intuitive comprehension of chemical systems in terms of the Wigner formalism based on the concept of phase-space. Finally, being based on a Monte Carlo method, it scales very well on parallel machines paving the way towards the time-dependent simulation of very complex molecules. A validation is performed by studying the electron distribution of three different systems, a Lithium atom, a Boron atom and a hydrogenic molecule. For the sake of simplicity, we start from initial conditions not too far from equilibrium and show that the systems reach a stationary regime, as expected (despite no restriction is imposed in the choice of the initial conditions). We also show a good agreement with the standard density functional theory for the hydrogenic molecule. These results demonstrate that the combination of the Wigner Monte Carlo method and Kohn-Sham systems provides a reliable computational tool which could, eventually, be applied to more sophisticated problems.
Density functional theory and chromium: Insights from the dimers
Würdemann, Rolf; Kristoffersen, Henrik H.; Moseler, Michael; Walter, Michael
2015-03-28
The binding in small Cr clusters is re-investigated, where the correct description of the dimer in three charge states is used as criterion to assign the most suitable density functional theory approximation. The difficulty in chromium arises from the subtle interplay between energy gain from hybridization and energetic cost due to exchange between s and d based molecular orbitals. Variations in published bond lengths and binding energies are shown to arise from insufficient numerical representation of electron density and Kohn-Sham wave-functions. The best functional performance is found for gradient corrected (GGA) functionals and meta-GGAs, where we find severe differences between functionals from the same family due to the importance of exchange. Only the “best fit” from Bayesian error estimation is able to predict the correct energetics for all three charge states unambiguously. With this knowledge, we predict small bond-lengths to be exclusively present in Cr{sub 2} and Cr{sub 2}{sup −}. Already for the dimer cation, solely long bond-lengths appear, similar to what is found in the trimer and in chromium bulk.
Curvature and Frontier Orbital Energies in Density Functional Theory
Stein, Tamar; Autschbach, Jochen; Govind, Niranjan; Kronik, Leeor; Baer, Roi
2012-12-20
Perdew et al. [Phys. Rev. Lett 49, 1691 (1982)] discovered and proved two different properties that exact Kohn-Sham density functional theory (DFT) must obey: (i) The exact total energy versus particle number must be a series of linear segments between integer electron points; (ii) Across an integer number of electrons, the exchange-correlation potential may ``jump’’ by a constant, known as the derivative discontinuity (DD). Here, we show analytically that in both the original and the generalized Kohn-Sham formulation of DFT, the two are in fact two sides of the same coin. Absence of a derivative discontinuity necessitates deviation from piecewise linearity, and the latter can be used to correct for the former, thereby restoring the physical meaning of the orbital energies. Using selected small molecules, we show that this results in a simple correction scheme for any underlying functional, including semi-local and hybrid functionals as well as Hartree-Fock theory, suggesting a practical correction for the infamous gap problem of density functional theory. Moreover, we show that optimally-tuned range-separated hybrid functionals can inherently minimize both DD and curvature, thus requiring no correction, and show that this can be used as a sound theoretical basis for novel tuning strategies.
Van der Waals Density Functional Theory with Applications
NASA Astrophysics Data System (ADS)
Langreth, David C.
2004-03-01
We discuss the development of electronic density functionals that are applicable for weakly bound systems where the van der Waals interaction and its ramifications become important. Our current functionals approach the correct asymptotic dependence at large distances and are seamless at small distances. The first form of the functional, appropriate for layered systems, has been recently applied to graphite, boron nitride, and molybdenum sulfide [H. Rydberg et al., Phys. Rev. Lett. 91, 126402 (2003) and D. C. Langreth, Int. J. Quant. Chem. (submitted), see http//:www.physics.rutgers.edu/ ˜langreth/preprints/dft2003.pdf]. The second form of the functional [M. Dion it et al. (to be published)] is appropriate for arbitrary geometries. Recent results on rare gas dimers and the benzene dimer suggest promise for this method as well.
A relativistic time-dependent density functional study of the excited states of the mercury dimer
Kullie, Ossama E-mail: ossama.kullie@unistra.fr
2014-01-14
In previous works on Zn{sub 2} and Cd{sub 2} dimers we found that the long-range corrected CAMB3LYP gives better results than other density functional approximations for the excited states, especially in the asymptotic region. In this paper, we use it to present a time-dependent density functional (TDDFT) study for the ground-state as well as the excited states corresponding to the (6s{sup 2} + 6s6p), (6s{sup 2} + 6s7s), and (6s{sup 2} + 6s7p) atomic asymptotes for the mercury dimer Hg{sub 2}. We analyze its spectrum obtained from all-electron calculations performed with the relativistic Dirac-Coulomb and relativistic spinfree Hamiltonian as implemented in DIRAC-PACKAGE. A comparison with the literature is given as far as available. Our result is excellent for the most of the lower excited states and very encouraging for the higher excited states, it shows generally good agreements with experimental results and outperforms other theoretical results. This enables us to give a detailed analysis of the spectrum of the Hg{sub 2} including a comparative analysis with the lighter dimers of the group 12, Cd{sub 2}, and Zn{sub 2}, especially for the relativistic effects, the spin-orbit interaction, and the performance of CAMB3LYP and is enlightened for similar systems. The result shows, as expected, that spinfree Hamiltonian is less efficient than Dirac-Coulomb Hamiltonian for systems containing heavy elements such as Hg{sub 2}.
Generalized Pauli constraints in reduced density matrix functional theory
Theophilou, Iris; Helbig, Nicole; Lathiotakis, Nektarios N.; Marques, Miguel A. L.
2015-04-21
Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman’s ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.
Density-dependence of functional spiking networks in vitro
Ham, Michael I; Gintautuas, Vadas; Rodriguez, Marko A; Bettencourt, Luis M A; Bennett, Ryan; Santa Maria, Cara L
2008-01-01
During development, the mammalian brain differentiates into specialized regions with unique functional abilities. While many factors contribute to this functional specialization, we explore the effect neuronal density can have on neuronal interactions. Two types of networks, dense (50,000 neurons and glia support cells) and sparse (12,000 neurons and glia support cells), are studied. A competitive first response model is applied to construct activation graphs that represent pairwise neuronal interactions. By observing the evolution of these graphs during development in vitro we observe that dense networks form activation connections earlier than sparse networks, and that link-!llltropy analysis of the resulting dense activation graphs reveals that balanced directional connections dominate. Information theoretic measures reveal in addition that early functional information interactions (of order 3) are synergetic in both dense and sparse networks. However, during development in vitro, such interactions become redundant in dense, but not sparse networks. Large values of activation graph link-entropy correlate strongly with redundant ensembles observed in the dense networks. Results demonstrate differences between dense and sparse networks in terms of informational groups, pairwise relationships, and activation graphs. These differences suggest that variations in cell density may result in different functional specialization of nervous system tissue also in vivo.
Downlink Probability Density Functions for EOS-McMurdo Sound
NASA Technical Reports Server (NTRS)
Christopher, P.; Jackson, A. H.
1996-01-01
The visibility times and communication link dynamics for the Earth Observations Satellite (EOS)-McMurdo Sound direct downlinks have been studied. The 16 day EOS periodicity may be shown with the Goddard Trajectory Determination System (GTDS) and the entire 16 day period should be simulated for representative link statistics. We desire many attributes of the downlink, however, and a faster orbital determination method is desirable. We use the method of osculating elements for speed and accuracy in simulating the EOS orbit. The accuracy of the method of osculating elements is demonstrated by closely reproducing the observed 16 day Landsat periodicity. An autocorrelation function method is used to show the correlation spike at 16 days. The entire 16 day record of passes over McMurdo Sound is then used to generate statistics for innage time, outage time, elevation angle, antenna angle rates, and propagation loss. The levation angle probability density function is compared with 1967 analytic approximation which has been used for medium to high altitude satellites. One practical result of this comparison is seen to be the rare occurrence of zenith passes. The new result is functionally different than the earlier result, with a heavy emphasis on low elevation angles. EOS is one of a large class of sun synchronous satellites which may be downlinked to McMurdo Sound. We examine delay statistics for an entire group of sun synchronous satellites ranging from 400 km to 1000 km altitude. Outage probability density function results are presented three dimensionally.
Revealing the density of encoded functions in a viral RNA
Patel, Nikesh; Dykeman, Eric C.; Coutts, Robert H. A.; Lomonossoff, George P.; Rowlands, David J.; Phillips, Simon E. V.; Ranson, Neil; Twarock, Reidun; Tuma, Roman; Stockley, Peter G.
2015-01-01
We present direct experimental evidence that assembly of a single-stranded RNA virus occurs via a packaging signal-mediated mechanism. We show that the sequences of coat protein recognition motifs within multiple, dispersed, putative RNA packaging signals, as well as their relative spacing within a genomic fragment, act collectively to influence the fidelity and yield of capsid self-assembly in vitro. These experiments confirm that the selective advantages for viral yield and encapsidation specificity, predicted from previous modeling of packaging signal-mediated assembly, are found in Nature. Regions of the genome that act as packaging signals also function in translational and transcriptional enhancement, as well as directly coding for the coat protein, highlighting the density of encoded functions within the viral RNA. Assembly and gene expression are therefore direct molecular competitors for different functional folds of the same RNA sequence. The strongest packaging signal in the test fragment, encodes a region of the coat protein that undergoes a conformational change upon contact with packaging signals. A similar phenomenon occurs in other RNA viruses for which packaging signals are known. These contacts hint at an even deeper density of encoded functions in viral RNA, which if confirmed, would have profound consequences for the evolution of this class of pathogens. PMID:25646435
Current density partitioning in time-dependent current density functional theory
Mosquera, Martín A.; Wasserman, Adam; Department of Physics, Purdue University, West Lafayette, Indiana 47907
2014-05-14
We adapt time-dependent current density functional theory to allow for a fragment-based solution of the many-electron problem of molecules in the presence of time-dependent electric and magnetic fields. Regarding a molecule as a set of non-interacting subsystems that individually evolve under the influence of an auxiliary external electromagnetic vector-scalar potential pair, the partition 4-potential, we show that there are one-to-one mappings between this auxiliary potential, a sharply-defined set of fragment current densities, and the total current density of the system. The partition electromagnetic (EM) 4-potential is expressed in terms of the real EM 4-potential of the system and a gluing EM 4-potential that accounts for exchange-correlation effects and mutual interaction forces between fragments that are required to yield the correct electron dynamics. We prove the zero-force theorem for the fragmented system, establish a variational formulation in terms of action functionals, and provide a simple illustration for a charged particle in a ring.
Storchi, Loriano; Rampino, Sergio; Belpassi, Leonardo; Tarantelli, Francesco; Quiney, Harry M
2013-12-10
We propose a new complete memory-distributed algorithm, which significantly improves the parallel implementation of the all-electron four-component Dirac-Kohn-Sham (DKS) module of BERTHA (J. Chem. Theory Comput. 2010, 6, 384). We devised an original procedure for mapping the DKS matrix between an efficient integral-driven distribution, guided by the structure of specific G-spinor basis sets and by density fitting algorithms, and the two-dimensional block-cyclic distribution scheme required by the ScaLAPACK library employed for the linear algebra operations. This implementation, because of the efficiency in the memory distribution, represents a leap forward in the applicability of the DKS procedure to arbitrarily large molecular systems and its porting on last-generation massively parallel systems. The performance of the code is illustrated by some test calculations on several gold clusters of increasing size. The DKS self-consistent procedure has been explicitly converged for two representative clusters, namely Au20 and Au34, for which the density of electronic states is reported and discussed. The largest gold cluster uses more than 39k basis functions and DKS matrices of the order of 23 GB. PMID:26592273
Stochastic Optimally Tuned Range-Separated Hybrid Density Functional Theory.
Neuhauser, Daniel; Rabani, Eran; Cytter, Yael; Baer, Roi
2016-05-19
We develop a stochastic formulation of the optimally tuned range-separated hybrid density functional theory that enables significant reduction of the computational effort and scaling of the nonlocal exchange operator at the price of introducing a controllable statistical error. Our method is based on stochastic representations of the Coulomb convolution integral and of the generalized Kohn-Sham density matrix. The computational cost of the approach is similar to that of usual Kohn-Sham density functional theory, yet it provides a much more accurate description of the quasiparticle energies for the frontier orbitals. This is illustrated for a series of silicon nanocrystals up to sizes exceeding 3000 electrons. Comparison with the stochastic GW many-body perturbation technique indicates excellent agreement for the fundamental band gap energies, good agreement for the band edge quasiparticle excitations, and very low statistical errors in the total energy for large systems. The present approach has a major advantage over one-shot GW by providing a self-consistent Hamiltonian that is central for additional postprocessing, for example, in the stochastic Bethe-Salpeter approach. PMID:26651840
Density functional theory: Its origins, rise to prominence, and future
NASA Astrophysics Data System (ADS)
Jones, R. O.
2015-07-01
In little more than 20 years, the number of applications of the density functional (DF) formalism in chemistry and materials science has grown in an astonishing fashion. The number of publications alone shows that DF calculations make up a huge success story, and many younger colleagues are surprised to learn that the widespread application of density functional methods, particularly in chemistry, began only after 1990. This is indeed unexpected, because the origins are usually traced to the papers of Hohenberg, Kohn, and Sham more than a quarter of a century earlier. The DF formalism, its applications, and prospects were reviewed for this journal in 1989. About the same time, the combination of DF calculations with molecular dynamics promised to provide an efficient way to study structures and reactions in molecules and extended systems. This paper reviews the development of density-related methods back to the early years of quantum mechanics and follows the breakthrough in their application after 1990. The two examples from biochemistry and materials science are among the many current applications that were simply far beyond expectations in 1990. The reasons why—50 years after its modern formulation and after two decades of rapid expansion—some of the most cited practitioners in the field are concerned about its future are discussed.
Subsystem real-time time dependent density functional theory.
Krishtal, Alisa; Ceresoli, Davide; Pavanello, Michele
2015-04-21
We present the extension of Frozen Density Embedding (FDE) formulation of subsystem Density Functional Theory (DFT) to real-time Time Dependent Density Functional Theory (rt-TDDFT). FDE is a DFT-in-DFT embedding method that allows to partition a larger Kohn-Sham system into a set of smaller, coupled Kohn-Sham systems. Additional to the computational advantage, FDE provides physical insight into the properties of embedded systems and the coupling interactions between them. The extension to rt-TDDFT is done straightforwardly by evolving the Kohn-Sham subsystems in time simultaneously, while updating the embedding potential between the systems at every time step. Two main applications are presented: the explicit excitation energy transfer in real time between subsystems is demonstrated for the case of the Na4 cluster and the effect of the embedding on optical spectra of coupled chromophores. In particular, the importance of including the full dynamic response in the embedding potential is demonstrated. PMID:25903875
Applications of large-scale density functional theory in biology
NASA Astrophysics Data System (ADS)
Cole, Daniel J.; Hine, Nicholas D. M.
2016-10-01
Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality.
Applications of large-scale density functional theory in biology.
Cole, Daniel J; Hine, Nicholas D M
2016-10-01
Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality. PMID:27494095
Van der Waals Interactions in Density Functional Theory: Intermolecular Complexes
NASA Astrophysics Data System (ADS)
Kannemann, Felix; Becke, Axel
2010-03-01
Conventional density functional theory (GGA and hybrid functionals) fails to account for dispersion interactions and is therefore not applicable to systems where van der Waals interactions play a dominant role, such as intermolecular complexes and biomolecules. The exchange-hole dipole moment (XDM) dispersion model of Becke and Johnson [A. D. Becke and E. R. Johnson, J. Chem. Phys. 127, 154108 (2007)] corrects for this deficiency. We have previously shown that the XDM dispersion model can be combined with standard GGA functionals (PW86 for exchange and PBE for correlation) to give accurate binding energy curves for rare-gas diatomics [F. O. Kannemann and A. D. Becke, J. Chem. Theory Comput. 5, 719 (2009)]. Here we present further tests of the GGA-XDM method using benchmark sets including hydrogen bonding, electrostatic, dispersion and stacking interactions, and systems ranging from rare-gas diatomics to biomolecular complexes.
Dipole polarizability of 120Sn and nuclear energy density functionals
NASA Astrophysics Data System (ADS)
Hashimoto, T.; Krumbholz, A. M.; Reinhard, P.-G.; Tamii, A.; von Neumann-Cosel, P.; Adachi, T.; Aoi, N.; Bertulani, C. A.; Fujita, H.; Fujita, Y.; GanioÇ§lu, E.; Hatanaka, K.; Ideguchi, E.; Iwamoto, C.; Kawabata, T.; Khai, N. T.; Krugmann, A.; Martin, D.; Matsubara, H.; Miki, K.; Neveling, R.; Okamura, H.; Ong, H. J.; Poltoratska, I.; Ponomarev, V. Yu.; Richter, A.; Sakaguchi, H.; Shimbara, Y.; Shimizu, Y.; Simonis, J.; Smit, F. D.; Süsoy, G.; Suzuki, T.; Thies, J. H.; Yosoi, M.; Zenihiro, J.
2015-09-01
The electric dipole strength distribution in 120Sn between 5 and 22 MeV has been determined at the Research Center for Nuclear Physics, Osaka, from polarization transfer observables measured in proton inelastic scattering at E0=295 MeV and forward angles including 0∘. Combined with photoabsorption data, a highly precise electric dipole polarizability αD(120Sn) =8.93 (36 ) fm3 is extracted. The dipole polarizability as isovector observable par excellence carries direct information on nuclear symmetry energy and its density dependence. The correlation of the new value with the well-established αD(208Pb) serves as a test of its prediction by nuclear energy density functionals. Models based on modern Skyrme interactions describe the data fairly well while most calculations based on relativistic Hamiltonians cannot.
Dynamic density functional theory of solid tumor growth: Preliminary models.
Chauviere, Arnaud; Hatzikirou, Haralambos; Kevrekidis, Ioannis G; Lowengrub, John S; Cristini, Vittorio
2012-03-01
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. PMID:22489279
Sublinear scaling for time-dependent stochastic density functional theory
Gao, Yi; Neuhauser, Daniel; Baer, Roi; Rabani, Eran
2015-01-21
A stochastic approach to time-dependent density functional theory is developed for computing the absorption cross section and the random phase approximation (RPA) correlation energy. The core idea of the approach involves time-propagation of a small set of stochastic orbitals which are first projected on the occupied space and then propagated in time according to the time-dependent Kohn-Sham equations. The evolving electron density is exactly represented when the number of random orbitals is infinite, but even a small number (≈16) of such orbitals is enough to obtain meaningful results for absorption spectrum and the RPA correlation energy per electron. We implement the approach for silicon nanocrystals using real-space grids and find that the overall scaling of the algorithm is sublinear with computational time and memory.
Hubbard operator density functional theory for Fermionic lattice models
NASA Astrophysics Data System (ADS)
Cheng, Zhengqian; Marianetti, Chris
We formulate an effective action as a functional of Hubbard operator densities whose stationary point delivers all local static information of the interacting lattice model. Using the variational principle, we get a self-consistent equation for Hubbard operator densities. The computational cost of our approach is set by diagonalizing the local Fock space. We apply our method to the one and two band Hubbard model (including crystal field and on-site exchange) in infinite dimensions where the exact solution is known. Excellent agreement is obtained for the one-band model. In the two-band model, good agreement is obtained in the metallic region of the phase diagram in addition to the metal-insulator transition. While our approach does not address frequency dependent observables, it has a negligible computational cost as compared to dynamical mean field theory and could be highly applicable in the context total energies of strongly correlated materials and molecules.
Density Functional Calculations for Atoms, Molecules and Clusters
NASA Astrophysics Data System (ADS)
Gunnarsson, O.; Jones, R. O.
1980-01-01
The density functional formalism provides a framework for including exchange and correlation effects in the calculation of ground state properties of many-electron systems. The reduction of the problem to the solution of single-particle equations leads to important numerical advantages over other ab initio methods of incorporating correlation effects. The essential features of the scheme are outlined and results obtained for atomic and molecular systems are surveyed. The local spin density (LSD) approximation gives generally good results for systems where the bonding involves s and p electrons, but results are less satisfactory for d-bonded systems. Non-local modifications to the LSD approximation have been tested on atomic systems yielding much improved total energies.
Relativistic density functional theory for finite nuclei and neutron stars
NASA Astrophysics Data System (ADS)
Piekarewicz, Jorge
In 1939 Oppenheimer and Volkoff demonstrated using Einstein's theory of general relativity that a neutron star supported exclusively by neutron degeneracy pressure will collapse into a black hole if its mass exceeds seven tenths of a solar mass. Seventy five years after such a pioneering prediction the existence of neutron stars with masses as large as two solar masses has been firmly established. This fact alone highlights the critical role that nuclear interactions play in explaining the structure of neutron stars. Indeed, a neutron star is a gold mine for the study of nuclear phenomena that span an enormous range of densities and neutron-proton asymmetries. Physical phenomena over such diverse scales are best described by a formalism based on Relativistic Density Functional Theory. In this contribution I focus on the synergy between theory, experiment, and observation that is needed to elucidate the myriad of exotic states of matter that are believed to exist in a neutron star.
Nuclear chiral and magnetic rotation in covariant density functional theory
NASA Astrophysics Data System (ADS)
Meng, Jie; Zhao, Pengwei
2016-05-01
Excitations of chiral rotation observed in triaxial nuclei and magnetic and/or antimagnetic rotations (AMR) seen in near-spherical nuclei have attracted a lot of attention. Unlike conventional rotation in well-deformed or superdeformed nuclei, here the rotational axis is not necessary coinciding with any principal axis of the nuclear density distribution. Thus, tilted axis cranking (TAC) is mandatory to describe these excitations self-consistently in the framework of covariant density functional theory (CDFT). We will briefly introduce the formalism of TAC–CDFT and its application for magnetic and AMR phenomena. Configuration-fixed CDFT and its predictions for nuclear chiral configurations and for favorable triaxial deformation parameters are also presented, and the discoveries of the multiple chiral doublets in 133Ce and 103Rh are discussed.
Sissay, Adonay; Abanador, Paul; Mauger, François; Gaarde, Mette; Schafer, Kenneth J; Lopata, Kenneth
2016-09-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. PMID:27608987
Progress at the interface of wave-function and density-functional theories
Gidopoulos, Nikitas I.
2011-04-15
The Kohn-Sham (KS) potential of density-functional theory (DFT) emerges as the minimizing effective potential in a variational scheme that does not involve fixing the unknown single-electron density. Using Rayleigh Schroedinger (RS) perturbation theory (PT), we construct ab initio approximations for the energy difference, the minimization of which determines the KS potential directly - thereby bypassing DFT's traditional algorithm to search for the density that minimizes the total energy. From second-order RS PT, we obtain variationally stable energy differences to be minimized, solving the severe problem of variational collapse of orbital-dependent exchange-correlation functionals based on second-order RS PT.
Characterizing the Spatial Density Functions of Neural Arbors
NASA Astrophysics Data System (ADS)
Teeter, Corinne Michelle
Recently, it has been proposed that a universal function describes the way in which all arbors (axons and dendrites) spread their branches over space. Data from fish retinal ganglion cells as well as cortical and hippocampal arbors from mouse, rat, cat, monkey and human provide evidence that all arbor density functions (adf) can be described by a Gaussian function truncated at approximately two standard deviations. A Gaussian density function implies that there is a minimal set of parameters needed to describe an adf: two or three standard deviations (depending on the dimensionality of the arbor) and an amplitude. However, the parameters needed to completely describe an adf could be further constrained by a scaling law found between the product of the standard deviations and the amplitude of the function. In the following document, I examine the scaling law relationship in order to determine the minimal set of parameters needed to describe an adf. First, I find that the at, two-dimensional arbors of fish retinal ganglion cells require only two out of the three fundamental parameters to completely describe their density functions. Second, the three-dimensional, volume filling, cortical arbors require four fundamental parameters: three standard deviations and the total length of an arbor (which corresponds to the amplitude of the function). Next, I characterize the shape of arbors in the context of the fundamental parameters. I show that the parameter distributions of the fish retinal ganglion cells are largely homogenous. In general, axons are bigger and less dense than dendrites; however, they are similarly shaped. The parameter distributions of these two arbor types overlap and, therefore, can only be differentiated from one another probabilistically based on their adfs. Despite artifacts in the cortical arbor data, different types of arbors (apical dendrites, non-apical dendrites, and axons) can generally be differentiated based on their adfs. In addition, within
SUMO1 Affects Synaptic Function, Spine Density and Memory
Matsuzaki, Shinsuke; Lee, Linda; Knock, Erin; Srikumar, Tharan; Sakurai, Mikako; Hazrati, Lili-Naz; Katayama, Taiichi; Staniszewski, Agnieszka; Raught, Brian; Arancio, Ottavio; Fraser, Paul E.
2015-01-01
Small ubiquitin-like modifier-1 (SUMO1) plays a number of roles in cellular events and recent evidence has given momentum for its contributions to neuronal development and function. Here, we have generated a SUMO1 transgenic mouse model with exclusive overexpression in neurons in an effort to identify in vivo conjugation targets and the functional consequences of their SUMOylation. A high-expressing line was examined which displayed elevated levels of mono-SUMO1 and increased high molecular weight conjugates in all brain regions. Immunoprecipitation of SUMOylated proteins from total brain extract and proteomic analysis revealed ~95 candidate proteins from a variety of functional classes, including a number of synaptic and cytoskeletal proteins. SUMO1 modification of synaptotagmin-1 was found to be elevated as compared to non-transgenic mice. This observation was associated with an age-dependent reduction in basal synaptic transmission and impaired presynaptic function as shown by altered paired pulse facilitation, as well as a decrease in spine density. The changes in neuronal function and morphology were also associated with a specific impairment in learning and memory while other behavioral features remained unchanged. These findings point to a significant contribution of SUMO1 modification on neuronal function which may have implications for mechanisms involved in mental retardation and neurodegeneration. PMID:26022678
SUMO1 Affects Synaptic Function, Spine Density and Memory.
Matsuzaki, Shinsuke; Lee, Linda; Knock, Erin; Srikumar, Tharan; Sakurai, Mikako; Hazrati, Lili-Naz; Katayama, Taiichi; Staniszewski, Agnieszka; Raught, Brian; Arancio, Ottavio; Fraser, Paul E
2015-01-01
Small ubiquitin-like modifier-1 (SUMO1) plays a number of roles in cellular events and recent evidence has given momentum for its contributions to neuronal development and function. Here, we have generated a SUMO1 transgenic mouse model with exclusive overexpression in neurons in an effort to identify in vivo conjugation targets and the functional consequences of their SUMOylation. A high-expressing line was examined which displayed elevated levels of mono-SUMO1 and increased high molecular weight conjugates in all brain regions. Immunoprecipitation of SUMOylated proteins from total brain extract and proteomic analysis revealed ~95 candidate proteins from a variety of functional classes, including a number of synaptic and cytoskeletal proteins. SUMO1 modification of synaptotagmin-1 was found to be elevated as compared to non-transgenic mice. This observation was associated with an age-dependent reduction in basal synaptic transmission and impaired presynaptic function as shown by altered paired pulse facilitation, as well as a decrease in spine density. The changes in neuronal function and morphology were also associated with a specific impairment in learning and memory while other behavioral features remained unchanged. These findings point to a significant contribution of SUMO1 modification on neuronal function which may have implications for mechanisms involved in mental retardation and neurodegeneration. PMID:26022678
The performance of semilocal and hybrid density functionals in 3d transition-metal chemistry
NASA Astrophysics Data System (ADS)
Furche, Filipp; Perdew, John P.
2006-01-01
We investigate the performance of contemporary semilocal and hybrid density functionals for bond energetics, structures, dipole moments, and harmonic frequencies of 3d transition-metal (TM) compounds by comparison with gas-phase experiments. Special attention is given to the nonempirical metageneralized gradient approximation (meta-GGA) of Tao, Perdew, Staroverov, and Scuseria (TPSS) [Phys. Rev. Lett. 91, 146401 (2003)], which has been implemented in TURBOMOLE for the present work. Trends and error patterns for classes of homologous compounds are analyzed, including dimers, monohydrides, mononitrides, monoxides, monofluorides, polyatomic oxides and halogenides, carbonyls, and complexes with organic π ligands such as benzene and cyclopentadienyl. Weakly bound systems such as Ca2, Mn2, and Zn2 are discussed. We propose a reference set of reaction energies for benchmark purposes. Our all-electron results with quadruple zeta valence basis sets validate semilocal density-functional theory as the workhorse of computational TM chemistry. Typical errors in bond energies are substantially larger than in (organic) main group chemistry, however. The Becke-Perdew'86 [Phys. Rev. A 38, 3098 (1988); Phys. Rev. B 33, 8822 (1986)] GGA and the TPSS meta-GGA have the best price/performance ratio, while the TPSS hybrid functional achieves a slightly lower mean absolute error in bond energies. The popular Becke three-parameter hybrid B3LYP underbinds significantly and tends to overestimate bond distances; we give a possible explanation for this. We further show that hybrid mixing does not reduce the width of the error distribution on our reference set. The error of a functional for the s-d transfer energy of a TM atom does not predict its error for TM bond energies and bond lengths. For semilocal functionals, self-interaction error in one- and three-electron bonds appears to be a major source of error in TM reaction energies. Nevertheless, TPSS predicts the correct ground