Downloads & Free Reading Options - Results

A real-space formalism for density-functional perturbation theory (DFPT) is derived and applied for the computation of harmonic vibrational properties in molecules and solids. The practical implementation using numeric atom-centered orbitals as basis functions is demonstrated exemplarily for the all-electron Fritz Haber Institute ab initio molecular simulations (FHI-aims) package. The convergence of the calculations with respect to numerical parameters is carefully investigated and a systematic comparison with finite-difference approaches is performed both for finite (molecules) and extended (periodic) systems. Finally, the scaling tests and scalability tests on massively parallel computer systems demonstrate the computational efficiency.

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. 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--can be realized, for both all-electron and 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 pseudopotential calculations, and is comparable to the Gaussian basis for all-electron calculations. Further, we demonstrate the capability of the proposed approach to compute the electronic structure of materials systems containing a few thousand atoms using modest computational resources, and good scalability of the present implementation up to a few hundred processors.

We present calculations of ground state resonance structure associated with molecular clusters of the high explosives RDX, PETN, and TNT using density functional theory (DFT), which is for the construction of parameterized dielectric response functions for excitation by electromagnetic waves at compatible frequencies. These dielectric functions provide for different types of analyses concerning the dielectric response of explosives. In particular, these dielectric response functions provide quantitative initial estimates of spectral response features for subsequent adjustment with respect to additional information such as laboratory measurements and other types of theory based calculations. With respect to qualitative analysis, these spectra provide for the molecular level interpretation of response structure. The DFT software GAUSSIAN was used for the calculations of ground state resonance structure presented here.

Calculations are presented of vibrational and electronic excited-state absorption spectra for equilibrium structures of SixOy molecular clusters using density function theory (DFT) and time-dependent density functional theory (TD-DFT). The size of the clusters considered is relatively large compared to those considered in previous studies. DFT and TD-DFT can provide interpretation of absorption spectra with respect to molecular structure for excitation by electromagnetic waves at frequencies within the IR and UV-visible ranges. The absorption spectrum corresponding to excitation states of SixOy-nH2O molecular clusters consisting of relatively small numbers of atoms should be associated with response features that are intermediate between that of isolated molecules and that of a bulk system. DFT and TD-DFT calculated absorption spectra represent quantitative estimates that can be correlated with additional information obtained from laboratory measurements. The DFT software GAUSSIAN was used for the calculations of excitation states presented here.

The sensitivity of computed DFT (Density Functional Theory) molecular properties (including energetics, geometries, vibrational frequencies, and infrared intensities) to the radial and angular numerical integration grid meshes, as well as to the partitioning scheme, is discussed for a number of molecules using the Gaussian 98 program system. Problems with typical production grid sizes are particularly acute for third-row transition metal systems, but may still result in qualitatively incorrect results for a molecule as simple as CCH. Practical recommendations are made with respect to grid choices for the energy(+gradient) steps, as well as for the solution of the CPKS (Coupled Perturbed Kohn-Sham) equations.

We consider an atomistic model defined through an interaction field satisfying a variational principle, and can therefore be considered a toy model of (orbital free) density functional theory. We investigate atomistic-to-continuum coupling mechanisms for this atomistic model, paying special attention to the dependence of the atomistic subproblem on the atomistic region boundary and the boundary conditions. We rigorously prove first-order error estimates for two related coupling mechanisms.

We consider an atomistic model defined through an interaction field satisfying a variational principle, and can therefore be considered a toy model of (orbital free) density functional theory. We investigate atomistic-to-continuum coupling mechanisms for this atomistic model, paying special attention to the dependence of the atomistic subproblem on the atomistic region boundary and the boundary conditions. We rigorously prove first-order error estimates for two related coupling mechanisms.

A method for formulating the unit cell of arbitrarily stacked, two-dimensional (2D) transition metal dichalcogenides is presented. Geometrical considerations and genetic algorithms are used to minimize the number of unit cells utilized in the construction of a supercell which may accommodate the lattice constants of arbitrary 2D and other close-packed materials. Supercells for various combinations of layered 2D transition metal dichalcogenides are calculated and their electronic band structures are simulated using density functional theory. Results are compared with previously reported density functional theory simulations from the literature.

We formulate a time-dependent density-matrix functional theory (TDDMFT) approach for higher-order correlation effects like biexcitons in optical processes in solids based on the reduced two-particle density-matrix formalism within the normal orbital representation. A TDDMFT version of the Schr\"odinger equation for biexcitons in terms of one- and two-body reduced density matrices is derived, which leads to finite biexcitonic binding energies already with an adiabatic approximation. Biexcitonic binding energies for several bulk semiconductors are calculated using a contact biexciton model.

A fully symmetry unrestricted Time-Dependent Density Functional Theory extended to include pairing correlations is used to calculate properties of the isovector giant dipole resonances of the deformed open-shell nuclei 172Yb (axially deformed), 188Os (triaxially deformed), and 238U (axially deformed), and to demonstrate good agreement with experimental data on nuclear photo-absorption cross-sections for two different Skyrme force parametrizations of the energy density functional: SkP and SLy4.

This article is concerned with the mathematical analysis of the perturbation method for extended Kohn-Sham models, in which fractional occupation numbers are allowed. All our results are established in the framework of the reduced Hartree-Fock (rHF) model, but our approach can be used to study other kinds of extended Kohn-Sham models, under some assumptions on the mathematical structure of the exchange- correlation functional. The classical results of Density Functional Perturbation Theory in the non-degenerate case (that is when the Fermi level is not a degenerate eigenvalue of the mean-field Hamiltonian) are formalized, and a proof of Wigner's (2n + 1) rule is provided. We then focus on the situation when the Fermi level is a degenerate eigenvalue of the rHF Hamiltonian, which had not been considered so far.

We propose a novel scheme to bring reduced density matrix functional theory (RDMFT) into the realm of density functional theory (DFT) that preserves the accurate density functional description at equilibrium, while incorporating accurately static and left-right correlation effects in molecules and keeping the good computational performance of DFT-based schemes. The key ingredient is to relax the requirement that the local potential is the functional derivative of the energy with respect to the density. Instead, we propose to restrict the search for the approximate natural orbitals within a domain where these orbitals are eigenfunctions of a single-particle hamiltonian with a local effective potential. In this way, fractional natural occupation numbers are accommodated into Kohn-Sham equations allowing for the description of molecular dissociation without breaking spin symmetry. Additionally, our scheme provides a natural way to connect an energy eigenvalue spectrum to the approximate natural orbitals and this spectrum is found to represent accurately the ionization potentials of atoms and small molecules.

Two electrons at the threshold of ionization represent a severe test case for electronic structure theory. A pseudospectral method yields a very accurate density of the two-electron ion with nuclear charge close to the critical value. Highly accurate energy components and potentials of Kohn-Sham density functional theory are given, as well as a useful parametrization of the critical density. The challenges for density functional approximations and the strength of correlation are also discussed.

The van Leeuwen proof of linear-response time-dependent density functional theory (TDDFT) is generalized to thermal ensembles. This allows generalization to finite temperatures of the Gross-Kohn relation, the exchange-correlation kernel of TDDFT, and fluctuation dissipation theorem for DFT. This produces a natural method for generating new thermal exchange-correlation (XC) approximations.

Numerous technological applications of Ni-based Zr and Hf intermetallic alloys promoted comprehensive studies in Zr$_8$Ni$_{21}$ and Hf$_8$Ni$_{21}$ by perturbed angular correlation (PAC) spectroscopy, which were not studied earlier until this report. The different phases produced in the samples have been identified by PAC and X-ray diffraction (XRD) measurements. Using $^{181}$Hf probe, two non-equivalent Zr/Hf sites have been observed in both Zr$_8$Ni$_{21}$ and Hf$_8$Ni$_{21}$ compounds. From present PAC measurements in Zr$_8$Ni$_{21}$, a component due to the production of Zr$_7$Ni$_{10}$ by eutectic reaction from the liquid metals is also observed. The phase Zr$_7$Ni$_{10}$, however, is not found from the XRD measurement. In Zr$_8$Ni$_{21}$, while the results do not change appreciably up to 973 K exhibit drastic changes at 1073 K. In Hf$_8$Ni$_{21}$, similar results for the two non-equivalent sites have been found but site fractions are in reverse order. In this alloy, a different contaminating phase, possibly due to HfNi$_3$, has been found from PAC measurements but is not found from XRD measurement. Density functional theory (DFT) based calculations of electric field gradient (EFG) and asymmetry parameter ($\eta$) at the sites of $^{181}$Ta probe nucleus allowed us to assign the observed EFG fractions to the various lattice sites in (Zr/Hf)$_8$Ni$_{21}$ compounds.

A number of diverse bulk properties of the zincblende and wurtzite III-V nitrides AlN, GaN, and InN, are predicted from first principles within density functional theory using the plane-wave ultrasoft pseudopotential method, within both the LDA (local density) and GGA (generalized gradient) approximations to the exchange-correlation functional. Besides structure and cohesion, we study formation enthalpies (a key ingredient in predicting defect solubilities and surface stability), spontaneous polarizations and piezoelectric constants (central parameters for nanostructure modeling), and elastic constants. Our study bears out the relative merits of the two density functional approaches in describing diverse properties of the III-V nitrides (and of the parent species N$_2$, Al, Ga, and In), and leads us to conclude that the GGA approximation, associated with high-accuracy techniques such as multiprojector ultrasoft pseudopotentials or modern all-electron methods, is to be preferred in the study of III-V nitrides.

Using density-functional theory within the generalized gradient approximation, we show that Si-based heterostructures with 1/4 layer $\delta$-doping of {\em interstitial} Mn (Mn$_{\mathrm int}$) are half-metallic. For Mn$_{\mathrm int}$ concentrations of 1/2 or 1 layer, the states induced in the band gap of $\delta$-doped heterostructures still display high spin polarization, about 85% and 60%, respectively. The proposed heterostructures are more stable than previously assumed $\delta$-layers of {\em substitutional} Mn. Contrary to wide-spread belief, the present study demonstrates that {\em interstitial} Mn can be utilized to tune the magnetic properties of Si, and thus provides a new clue for Si-based spintronics materials.

The effective action for the charge density and the photon field is proposed as a generalization of the density functional. A simple definition is given for the density functional, as the functional Legendre transform of the generator functional of connected Green functions for the density and the photon field, offering systematic approximation schemes. The leading order of the perturbation expansion reproduces the Hartree-Fock equation. A renormalization group motivated method is introduced to turn on the Coulomb interaction gradually and to find corrections to the Hartree-Fock and the Kohn-Sham schemes.

We show that the time-dependent particle density $n(\vec r,t)$ and the current density ${\vec j}(\vec r,t)$ of a many-particle system that evolves under the action of external scalar and vector potentials $V(\vec r,t)$ and $\vec A(\vec r,t)$ and is initially in the quantum state $|\psi (0)>$, can always be reproduced (under mild assumptions) in another many-particle system, with different two-particle interaction, subjected to external potentials $V'(\vec r,t)$ and $\vec A'(\vec r,t)$, starting from an initial state $|\psi' (0)>$, which yields the same density and current as $|\psi (0)>$. Given the initial state of this other many-particle system, the potentials $V'(\vec r,t)$ and $\vec A'(\vec r,t)$ are uniquely determined up to gauge transformations that do not alter the initial state. As a special case, we obtain a new and simpler proof of the Runge-Gross theorem for time-dependent current density functional theory. This theorem provides a formal basis for the application of time-dependent current density functional theory to transport problems.

We report the formation of a non-magnetic band insulator at the isopolar interface between the antiferromagnetic Mott-Hubbard insulator LaTiO3 and the antiferromagnetic charge transfer insulator LaFeO3. By density functional theory calculations, we find that the formation of this interface state is driven by the combination of O band alignment and crystal field splitting energy of the t2g and eg bands. As a result of these two driving forces, the Fe 3d bands rearrange and electrons are transferred from Ti to Fe. This picture is supported by x-ray photoelectron spectroscopy, which confirms the rearrangement of the Fe 3d bands and reveals an unprecedented charge transfer up to 1.2+/-0.2 e-/interface unit cell in our LaTiO3/LaFeO3 heterostructures.

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

We present a method capable of calculating elastic scanning tunneling microscopy (STM) currents from localized atomic orbital density functional theory (DFT). To overcome the poor accuracy of the localized orbital description of the wave functions far away from the atoms, we propagate the wave functions, using the total DFT potential. From the propagated wave functions, the Bardeen's perturbative approach provides the tunneling current. To illustrate the method we investigate carbon monoxide adsorbed on a Cu(111) surface and recover the depression/protrusion observed experimentally with normal/CO-functionalized STM tips. The theory furthermore allows us to discuss the significance of s- and p-wave tips.

We examine the performance of a recently developed nonlocal density functional in predicting a model noncovalent interaction, the weak bond between an aromatic $\pi$ system and an aliphatic C-H group. The new functional is a significant improvement over traditional density functionals, providing results which compare favorably to high-level quantum-chemistry techniques but at considerably lower computational cost. Interaction energies in several model C-H/$\pi$ systems are in generally good agreement with coupled-cluster calculations, though equilibrium distances are consistently overpredicted when using the revPBE functional for exchange. The new functional correctly predicts changes in energy upon addition of halogen substituents.

We describe a procedure for mapping a self-consistent mean-field theory (also known as density functional theory) into a shell model Hamiltonian that includes quadrupole-quadrupole and monopole pairing interactions in a truncated space. We test our method in the deformed N=Z sd-shell nuclei Ne-20, Mg-24 and Ar-36, starting from the Hartree-Fock plus BCS approximation of the USD shell model interaction. A similar procedure is then followed using the SLy4 Skyrme energy density functional in the particle-hole channel plus a zero-range density-dependent force in the pairing channel. Using the ground-state solution of this density functional theory at the Hartree-Fock plus BCS level, an effective shell model Hamiltonian is constructed. We use this mapped Hamiltonian to extract quadrupolar and pairing correlation energies beyond the mean field approximation. The rescaling of the mass quadrupole operator in the truncated shell model space is found to be almost independent of the coupling strength used in the pairing channel of the underlying mean-field theory.

Density-functional theory is applied to compute the ground-state energies of quantum hard-sphere solids. The modified weighted-density approximation is used to map both the Bose and the Fermi solid onto a corresponding uniform Bose liquid, assuming negligible exchange for the Fermi solid. The required liquid-state input data are obtained from a paired phonon analysis and the Feynman approximation, connecting the static structure factor and the linear response function. The Fermi liquid is treated by the Wu-Feenberg cluster expansion, which approximately accounts for the effects of antisymmetry. Liquid-solid transitions for both systems are obtained with no adjustment of input data. Limited quantitative agreement with simulation indicates a need for further improvement of the liquid-state input through practical alternatives to the Feynman approximation.

Time-dependent density functional theory is thought to work well for the test cases of He and Be atoms. We perform a quantum defect analysis of the s, p, and d Rydberg states of Be with accurate ground state Kohn-Sham potentials. The s- and p-quantum defects are well described by the ALDA, but fails badly for the d-quantum defect. The same failure is observed in case of He. This provides a new challenge for functional development in time-dependent density functional theory.

The onset of condensation of hard spheres in a gravitational field is studied using density functional theory. In particular, we find that the local density approximation yields results identical to those obtained previously using the kinetic theory [Physica A 271, 192, (1999)], and a weighted density functional theory gives qualitatively similar results, namely, that the temperature at which condensation begins at the bottom scales linearly with weight, diameter, and number of layers of particles.

Sensitivity impacts owing to the variations in the type of zero-range pairing forces on the fission properties using the density functional theory 作者： Yang Su 1 Ze-Yu Li 1 Li-Le Liu 1 Guo-xiang Dong 2 Xiaobao Wang 2 Yong-Jing Chen 1 作者单位： 1. China Nuclear Data Center, China Institute of Atomic Energy, Beijing 102413, China 2. School of Science, Huzhou University, Huzhou 313000, China 通讯作者： Xiaobao Wang Email:[email protected] 提交时间： 2024-02-21 13:08:44 摘要: By the Skyrme density functional theory (DFT), potential energy surfaces (PES) of $^{240}$Pu with constraints on the axial quadrupole and octupole deformations($q_{20}$ and $q_{30}$) are calculated. The volume-like, and the surface-like pairing forces, and the mixing between these two forces are used, within the Hartree-Fock-Bogoliubov (HFB) approximation. The variations of the least-energy fission path, fission barrier, paring energy, total kinetic energy, scission line, and mass distribution of fission fragments by the different forms of paring force are analyzed and discussed. The fission dynamics is studied based on the time-dependent generator coordinate method (TDGCM) plus Gaussian overlap approximation (GOA). The results show a sensitivity of the mass and charge distributions of fission fragments to the form of the pairing force. Based on the investigation of the neutron-induced fission of $^{239}$Pu, among the volume, mixed, and surface pairing forces, the mixed-type of pairing force can give a good reproduction of experimental data. Nuclear fission density functional theory paring force potential energy surfaces fission fragment distribution 来自： 小保王 分类： 物理学 >> 核物理学 投稿状态： 已投稿期刊 引用： ChinaXiv:202402.00221 (或此版本 ChinaXiv:202402.00221V1 ) DOI:10.12074/202402.00221V1 CSTR:32003.36.ChinaXiv.202402.00221.V1 推荐引用方式： Yang Su,Ze-Yu Li,Li-Le Liu,Guo-xiang Dong,Xiaobao Wang,Yong-Jing Chen.(2024).Sensitivity impacts owing to the variations in the type of zero-range pairing forces on the fission properties using the density functional theory.中国科学院科技论文预发布平台.doi:10.12074/T202402.00268V1 版本历史 [V1] 2024-02-21 13:08:44 ChinaXiv:202402.00221V1 下载全文

Density functional theory for a simple model of dendrimers is proposed. The theory is based on fundamental measure theory which accounts for the hard-sphere repulsion of the segments and on the Wertheim first-order perturbation theory for the correlations due to connectivity. Set of the recurrence formulae for the ideal chain contribution involving simple integrals is derived. By using perturbation theory dispersion forces can be easily included.

Calculations are presented of vibrational and electronic excited-state absorption spectra for equilibrium structures of SixOy-nH2O molecular clusters using density function theory (DFT) and time-dependent density functional theory (TD-DFT). The size of the clusters considered is relatively large compared to those considered in previous studies. DFT and TD-DFT can provide interpretation of absorption spectra with respect to molecular structure for excitation by electromagnetic waves at frequencies within the IR and UV-visible ranges. The absorption spectrum corresponding to excitation states of SixOy-nH2O molecular clusters consisting of relatively small numbers of atoms should be associated with response features that are intermediate between that of isolated molecules and that of a bulk system. DFT and TD-DFT calculated absorption spectra represent quantitative estimates that can be correlated with additional information obtained from laboratory measurements. The DFT software GAUSSIAN was used for the calculations of excitation states presented here.

Reliable predictions of nuclear properties are needed as much to answer fundamental science questions as in applications such as reactor physics or data evaluation. Nuclear density functional theory is currently the only microscopic, global approach to nuclear structure that is applicable throughout the nuclear chart. In the past few years, a lot of effort has been devoted to setting up a general methodology to assess theoretical uncertainties in nuclear DFT calculations. In this paper, we summarize some of the recent progress in this direction. Most of the new material discussed here will be be published in separate articles.

Using both dynamical density functional theory and particle-resolved Brownian dynamics simulations, we explore the flow of two-dimensional colloidal solids and fluids driven through a linear channel with a geometric constriction. The flow is generated by a constant external force acting on all colloids. The initial configuration is equilibrated in the absence of flow and then the external force is switched on instantaneously. Upon starting the flow, we observe four different scenarios: a complete blockade, a monotonic decay to a constant particle flux (typical for a fluid), a damped oscillatory behaviour in the particle flux, and a long-lived stop-and-go behaviour in the flow (typical for a solid). The dynamical density functional theory describes all four situations but predicts infinitely long undamped oscillations in the flow which are always damped in the simulations. We attribute the mechanisms of the underlying stop-and-go flow to symmetry conditions on the flowing solid. Our predictions are verifiable in real-space experiments on magnetic colloidal monolayers which are driven through structured microchannels and can be exploited to steer the flow throughput in microfluidics.

This article is from Scientific Reports , volume 4 . Abstract The detection of partial discharge by analyzing the components of SF6 gas in gas-insulated switchgears is important to the diagnosis and assessment of the operational state of power equipment. A gas sensor based on anatase TiO2 is used to detect decomposed gases in SF6. In this paper, first-principle density functional theory calculations are adopted to analyze the adsorption of SO2, SOF2, and SO2F2, the primary decomposition by-products of SF6 under partial discharge, on anatase (101) and (001) surfaces. Simulation results show that the perfect anatase (001) surface has a stronger interaction with the three gases than that of anatase (101), and both surfaces are more sensitive and selective to SO2 than to SOF2 and SO2F2. The selection of a defect surface to SO2, SOF2, and SO2F2 differs from that of a perfect surface. This theoretical result is corroborated by the sensing experiment using a TiO2 nanotube array (TNTA) gas sensor. The calculated values are analyzed to explain the results of the Pt-doped TNTA gas sensor sensing experiment. The results imply that the deposited Pt nanoparticles on the surface increase the active sites of the surface and the gas molecules may decompose upon adsorption on the active sites.

We present calculations of ground state resonance structure associated with molecular clusters of HMX using density functional theory (DFT), which is for the construction of parameterized dielectric response functions for excitation by electromagnetic waves at compatible frequencies. These dielectric functions provide for different types of analyses concerning the dielectric response of explosives. In particular, these dielectric response functions provide quantitative initial estimates of spectral response features for subsequent adjustment with respect to additional information such as laboratory measurements and other types of theory-based calculations. With respect to qualitative analysis, these spectra provide for the molecular-level interpretation of response structure. The DFT software GAUSSIAN was used for the calculations of ground state resonance structure presented here.

We discuss the relationship between modern time-dependent density functional theory and earlier time-periodic versions, and why the criticisms in a recent paper (Chem. Phys. Lett. {\bf 433} (2006) 204) of our earlier analysis (Chem. Phys. Lett. {\bf 359} (2002) 237) are incorrect.

Time-dependent density functional theory is widely used to describe excitations of many-fermion systems. In its many applications, 3D coordinate-space representation is used, and infinite-domain calculations are limited to a finite volume represented by a box. For finite quantum systems (atoms, molecules, nuclei), the commonly used periodic or reflecting boundary conditions introduce spurious quantization of the continuum states and artificial reflections from boundary; hence, an incorrect treatment of evaporated particles. These artifacts can be practically cured by introducing absorbing boundary conditions (ABC) through an absorbing potential in a certain boundary region sufficiently far from the described system. But also the calculations of infinite matter (crystal electrons, quantum fluids, neutron star crust) suffer artifacts from a finite computational box. In this regime, twist- averaged boundary conditions (TABC) have been used successfully to diminish the finite-volume effects. In this work, we extend TABC to time-dependent framework and apply it to resolve the box artifacts for finite quantum systems using as test case small- and large-amplitude nuclear vibrations. We demonstrate that by using such a method, one can reduce finite volume effects drastically without adding any additional parameters. While they are almost equivalent in the linear regime, TABC and ABC differ in the nonlinear regime in their treatment of evaporated particles.

Time-dependent density-functional response theory for molecules. Lecture Notes Collection FreeScience.info ID2122 Obtained from http://tddft.org/TDDFT2004/PracticalSessions/papers/casida_1996.ps http://www.freescience.info/go.php?pagename=books&id=2122

Improving the accuracy and thus broadening the applicability of electronic density functional theory (DFT) is crucial to many research areas, from material science, to theoretical chemistry, biophysics and biochemistry. In the last three years, the mathematical structure of the strong-interaction limit of density functional theory has been uncovered, and exact information on this limit has started to become available. The aim of this paper is to give a perspective on how this new piece of exact information can be used to treat situations that are problematic for standard Kohn-Sham DFT. One way to use the strong-interaction limit, more relevant for solid-state physical devices, is to define a new framework to do practical, non-conventional, DFT calculations in which a strong-interacting reference system is used instead of the traditional non-interacting one of Kohn and Sham. Another way to proceed, more related to chemical applications, is to include the exact treatment of the strong-interaction limit into approximate exchange-correlation energy density functionals in order to describe difficult situations such as the breaking of the chemical bond.

We outline a Kohn-Sham-Dirac density-functional-theory (DFT) scheme for graphene sheets that treats slowly-varying inhomogeneous external potentials and electron-electron interactions on an equal footing. The theory is able to account for the the unusual property that the exchange-correlation contribution to chemical potential increases with carrier density in graphene. Consequences of this property, and advantages and disadvantages of using the DFT approach to describe it, are discussed. The approach is illustrated by solving the Kohn-Sham-Dirac equations self-consistently for a model random potential describing charged point-like impurities located close to the graphene plane. The influence of electron-electron interactions on these non-linear screening calculations is discussed at length, in the light of recent experiments reporting evidence for the presence of electron-hole puddles in nearly-neutral graphene sheets.

Forty-five years after the point de d\'epart [1] of density functional theory, its applications in chemistry and the study of electronic structures keep steadily growing. However, the precise form of the "divine" energy functional in terms of the electron density [2] still eludes us --and possibly will do so forever [3]. In what follows we examine a formulation in the same spirit with phase space variables. The validity of Hohenberg-Kohn-Levy-type theorems on phase space is recalled. We study the representability problem for reduced Wigner functions, and proceed to analyze properties of the new functional. Along the way, new results on states in the phase-space formalism of quantum mechanics are established. Natural Wigner orbital theory is developed in depth, with the final aim of constructing accurate correlation-exchange functionals on phase space. A new proof of the overbinding property of the Mueller functional is given. This exact theory supplies its home at long last to that illustrious ancestor, the Thomas-Fermi model.

Kohn-Sham density functional theory is one of the most widely used electronic structure theories. In the pseudopotential framework, uniform discretization of the Kohn-Sham Hamiltonian generally results in a large number of basis functions per atom in order to resolve the rapid oscillations of the Kohn-Sham orbitals around the nuclei. Previous attempts to reduce the number of basis functions per atom include the usage of atomic orbitals and similar objects, but the atomic orbitals generally require fine tuning in order to reach high accuracy. We present a novel discretization scheme that adaptively and systematically builds the rapid oscillations of the Kohn-Sham orbitals around the nuclei as well as environmental effects into the basis functions. The resulting basis functions are localized in the real space, and are discontinuous in the global domain. The continuous Kohn-Sham orbitals and the electron density are evaluated from the discontinuous basis functions using the discontinuous Galerkin (DG) framework. Our method is implemented in parallel and the current implementation is able to handle systems with at least thousands of atoms. Numerical examples indicate that our method can reach very high accuracy (less than 1meV) with a very small number ($4\sim 40$) of basis functions per atom.

We reformulate and generalize the uniqueness and existence proofs of time-dependent density-functional theory. The central idea is to restate the fundamental one-to-one correspondence between densities and potentials as a global fixed point question for potentials on a given time-interval. We show that the unique fixed point, i.e. the unique potential generating a given density, is reached as the limiting point of an iterative procedure. The one-to-one correspondence between densities and potentials is a straightforward result provided that the response function of the divergence of the internal forces is bounded. The existence, i.e. the v-representability of a density, can be proven as well provided that the operator norms of the response functions of the members of the iterative sequence of potentials have an upper bound. The densities under consideration have second time-derivatives that are required to satisfy a condition slightly weaker than being square-integrable. This approach avoids the usual restrictions of Taylor-expandability in time of the uniqueness theorem by Runge and Gross [Phys.Rev.Lett.52, 997 (1984)] and of the existence theorem by van Leeuwen [Phys.Rev.Lett. 82, 3863 (1999)]. Owing to its generality, the proof not only answers basic questions in density-functional theory but also has potential implications in other fields of physics.

Local and semilocal density-functional approximations for the exchange-correlation energy fail badly in the zero-thickness limit of a quasi-two-dimensional electron gas, where the density variation is rapid almost everywhere. Here we show that a fully nonlocal fifth-rung functional, the inhomogeneous Singwi-Tosi-Land-Sj\"olander (STLS) approach, which employs both occupied and unoccupied Kohn-Sham orbitals, recovers the true two-dimensional STLS limit and appears to be remarkably accurate for any thickness of the slab (and thus for the dimensional crossover). We also show that this good behavior is only partly due to the use of the full exact exchange energy.

We introduce an orbital-optimized double-hybrid (DH) scheme using the optimized-effective-potential (OEP) method. The orbitals are optimized using a local potential corresponding to the complete exchange-correlation energy expression including the second-order M{{\o}}ller-Plesset (MP2) correlation contribution. We have implemented a one-parameter version of this OEP-based self-consistent DH scheme using the BLYP density-functional approximation and compared it to the corresponding non-self-consistent DH scheme for calculations on a few closed-shell atoms and molecules. While the OEP-based self-consistency does not provide any improvement for the calculations of ground-state total energies and ionization potentials, it does improve the accuracy of electron affinities and restores the meaning of the LUMO orbital energy as being connected to a neutral excitation energy. Moreover, the OEP-based self-consistent DH scheme provides reasonably accurate exchange-correlation potentials and correlated densities.

Wavelets

A spin angular gradient approximation for the exchange correlation magnetic field in the density functional formalism is proposed. The usage of such corrections leads to a consistent spin dynamical approach beyond the local approximation. The proposed technique does not contain any approximations for the form of potential and can be used in modern full potential band structure methods. The obtained results indicate that the direct 'potential' exchange in 3d magnets is rather small compared to the indirect 'kinetic' exchange, thus justifies the dynamical aspect of the local density approximation in 3d metals.

The group theory framework developed by Fukutome for a systematic analysis of the various broken symmetry types of Hartree-Fock solutions exhibiting spin structures is here extended to the general many body context using spinor-Green function formalism for describing magnetic systems. Consequences of this theory are discussed for examining the magnetism of itinerant electrons in nanometric systems of current interest as well as bulk systems where a vector spin-density form is required, by specializing our work to spin-density-functional formalism. We also formulate the linear response theory for such a system and compare and contrast them with the recent results obtained for localized electron systems. The various phenomenological treatments of itinerant magnetic systems are here unified in this group-theoretical description.

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

The multi-band Gutzwiller method, combined with calculations based on density functional theory, is employed to study total energy curves of the ferromagnetic ground state of Ni. A new method is presented which allows flow of charge between d and s, p type orbitals in an approximate way. Further it is emphasized that the missing repulsive contribution to the total energy at large magnetic moments can be estimated from an analysis of specific DFT calculations.

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