Downloads & Free Reading Options - Results

We present an implementation of the adiabatic spin-wave dynamics of Niu and Kleinman. This technique allows to decouple the spin and charge excitations of a many-electron system using a generalization of the adiabatic approximation. The only input for the spin-wave equations of motion are the energies and Berry curvatures of many-electron states describing frozen spin spirals. The latter are computed using a newly developed technique based on constrained density-functional theory, within the local spin density approximation and the pseudo-potential plane-wave method. Calculations for iron show an excellent agreement with experiments.

We present a novel method, which we call dual minima hopping method (DMHM), that allows us to find the global minimum of the potential energy surface (PES) within density functional theory for systems where a fast but less accurate calculation of the PES is possible. This method can rapidly find the ground state configuration of clusters and other complex systems with present day computer power by performing a systematic search. We apply the new method to silicon clusters. Even though these systems have already been extensively studied by other methods, we find new configurations that are lower in energy than the previously found.

Scanning tunneling spectroscopy measurements on the high temperature superconductor Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$ have reported an enhanced spectral gap in the neighborhood of O dopant atoms. We calculate, within density functional theory (DFT), the change in electronic structure due to such a dopant. We then construct and discuss the validity of several tight binding (TB) fits to the DFT bands with and without an O dopant. With the doping-modulated TB parameters, we finally evaluate the spin susceptibility and pairing interaction within spin fluctuation theory. The d-wave pairing eigenvalues are enhanced above the pure system without O dopant, supporting the picture of enhanced local pairing around such a defect.

We show that classical molecular density functional theory (MDFT), here in the homogeneous reference fluid approximation in which the functional is inferred from the properties of the bulk solvent, is a powerful new tool to study, at a fully molecular level, the solvation of complex surfaces and interfaces by polar solvents. This implicit solvent method allows for the determination of structural, orientational and energetic solvation properties that are on a par with all-atom molecular simulations performed for the same system, while reducing the computer time by two orders of magnitude. This is illustrated by the study of an atomistically-resolved clay surface composed of over a thousand atoms wetted by a molecular dipolar solvent. The high numerical efficiency of the method is exploited to carry a systematic analysis of the electrostatic and non-electrostatic components of the surface-solvent interaction within the popular CLAYFF force field. Solvent energetics and structure are found to depend weakly upon the atomic charges distribution of the clay surface, even for a rather polar solvent. We conclude on the consequences of such findings for force-field development.

We demonstrate the capabilities of time-dependent density functional theory (TDDFT) for strong-field, short wavelength (soft X-ray) physics, as compared to a formalism based on rate equations. We find that TDDFT provides a very good description of the total and individual ionization yields for Ne and Ar atoms exposed to strong laser pulses. We assess the reliability of different adiabatic density functionals and conclude that an accurate description of long-range interactions by the exchange and correlation potential is crucial for obtaining the correct ionization yield over a wide range of intensities ($10^{13}$ -- $5 \times 10^{15}$ W/cm$^2$). Our TDDFT calculations disentangle the contribution from each ionization channel based on the Kohn-Sham wavefunctions.

Approximate molecular calculations via standard Kohn-Sham Density Functional Theory are exactly reproduced by performing self-consistent calculations on isolated fragments via Partition Density Functional Theory [Phys. Rev. A 82, 024501 (2010)]. We illustrate this with the binding curves of small diatomic molecules. We ?nd that partition energies are in all cases qualitatively similar and numerically close to actual binding energies. We discuss qualitative features of the associated partition potentials.

A simple application of classical density functional theory is derived and applied to a system of polymers grafted to a plane. The system is assumed to have symmetry in directions parallel to the grafting plane hence it being a '1-dimensional' problem. A quick introduction to using propagators in the theory of chains in the presence of external fields and a numerical algorithm for solving the modified diffusion (or Fokker-Planck) equation are presented. Density profiles are found for hard-sphere chains which agree with results in the literature (see Murat,.., Milner,.., Muhukumar and others) and align with physical expectations of brush formation due to excluded volume. The linear scaling $h \approx (\sigma)^{1/3} N$, of brush height with polymerisation $N$, is reproduced here.

Starting from the observation that one of the most successful methods for solving the Kohn-Sham equations for periodic systems -- the plane-wave method -- is a spectral method based on eigenfunction expansion, we formulate a spectral method designed towards solving the Kohn-Sham equations for clusters. This allows for efficient calculation of the electronic structure of clusters (and molecules) with high accuracy and systematic convergence properties without the need for any artificial periodicity. The basis functions in this method form a complete orthonormal set and are expressible in terms of spherical harmonics and spherical Bessel functions. Computation of the occupied eigenstates of the discretized Kohn-Sham Hamiltonian is carried out using a combination of preconditioned block eigensolvers and Chebyshev polynomial filter accelerated subspace iterations. Several algorithmic and computational aspects of the method, including computation of the electrostatics terms and parallelization are discussed. We have implemented these methods and algorithms into an efficient and reliable package called ClusterES (Cluster Electronic Structure). A variety of benchmark calculations employing local and non-local pseudopotentials are carried out using our package and the results are compared to the literature. Convergence properties of the basis set are discussed through numerical examples. Computations involving large systems that contain thousands of electrons are demonstrated to highlight the efficacy of our methodology. The use of our method to study clusters with arbitrary point group symmetries is briefly discussed.

It is shown here that Kohn-Sham equations cannot be derived from Hohenberg-Kohn theory without an additional postulate. Assuming that a functional derivative with respect to total electron density exists leads in general to a theory inconsistent with the exclusion principle. A mathematically and physically correct variational theory of the Kohn-Sham model can be developed using functional derivatives with respect to orbital densities. These partial (G\^ateaux) derivatives can be constructed explicitly from general N-electron theory and are defined throughout the orbital Hilbert space. This theory is consistent with the local density approximation (LDA), but does not in general imply multiplicative local exchange-correlation potentials. Progress beyond the LDA in condensed-matter physics requires development of methodology for nonlocal exchange and correlation potentials.

Contents:Introduction; Basic formalism for electrons in time-dependent electric fields 2.1 One-to-one mapping between time-dependent potentials and time-dependent densities 2.2 Stationary-action principle 2.3 Time-dependent Kohn-Sham scheme 3 Motion of the nuclei 3.1 Quantum mechanical treatment of nuclear motion 3.2 Classical treatment of nuclear motion 4 Electrons in time-dependent electromagnetic fields 4.1 Coupling to spin 4.2 Coupling to orbital currents 5 Perturbative regime, basic equations 5.1 Time-dependent linear density response 5.2 Time-dependent higher-order response 6 The time-dependent exchange-correlation potential: Rigorous prop- erties and approximate functionals 6.1 Approximations based on the homogeneous electron gas 6.2 Time-dependent optimized effective potential 7 Applications within the perturbative regime 7.1 Photoresponse of finite and infinite Systems 7.2 Calculation of excitation energies 7.3 Van der Waals interactions 8 Applications beyond the perturbative regime: Atoms in strong femto-second laser pulses. Lecture Notes Collection FreeScience.info ID1948 Obtained from http://www.physik.fu-berlin.de/~ag-gross/articles/pdf/GDP96.pdf http://www.freescience.info/go.php?pagename=books&id=1948

MRS Advances © 2018 Materials Research Society . This is an Open Acces  article,  distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited. DOI: 10.1557/adv.2018.180 

Characterization of Phase Transition in Heisenberg Fluids from Density Functional Theory 作者： Li Liang-Sheng 1 Li Li 1 Chen Xiao-Song 1 作者单位： 1. Institute of Theoretical Physics, Chinese Academy of sciences 提交时间： 2016-05-08 摘要: The phase transition of Heisenberg fluid has been investigated with the density functional theory in mean-field approximation (MF). The matrix of the second derivatives of the grand canonical potential Omega with respect to the particle density fluctuations and the magnetization fluctuations has been investigated and diagonalized. The smallest eigenvalue being 0 signalizes the phase instability and the related eigenvector characterizes this phase transition. We find a Curie line where the order parameter is pure magnetization and a spinodal where the order parameter is a mixture of particle density and magnetization. Along the spinodal, the character of phase in stability changes continuously from predominant condensation to predominant ferromagnetic phase transition with the decrease of total density. The spinodal meets the Curie line at the critical end point with the reduced density rho* = rho sigma(3) = 0.224 and the reduced temperature T* = kT/epsilon = 1.87 (sigma is the diameter of Heisenberg hard sphere and epsilon is the coupling constant). critical phenomena order-disorder magnetic liquid 期刊： COMMUNICATIONS IN THEORETICAL PHYSICS 分类： 物理学 >> 基本粒子与场物理学 引用： ChinaXiv:201605.01061 (或此版本 ChinaXiv:201605.01061V1 ) doi:10.12074/201605.01061 CSTR:32003.36.ChinaXiv.201605.01061.V1 推荐引用方式： Li Liang-Sheng,Li Li,Chen Xiao-Song.(2016).Characterization of Phase Transition in Heisenberg Fluids from Density Functional Theory.COMMUNICATIONS IN THEORETICAL PHYSICS.[ChinaXiv:201605.01061] 版本历史 [V1] 2016-05-08 22:39:08 ChinaXiv:201605.01061V1 下载全文

The exact interaction energy of a many-electron system is determined by the electron pair density, which is not well-approximated in standard Kohn-Sham density functional models. Here we study the (complicated but well-defined) exact universal map from density to pair density. We survey how many common functionals, including the most basic version of the LDA (Dirac exchange with no correlation contribution), arise from particular approximations of this map. We develop an algorithm to compute the map numerically, and apply it to one-parameter families {a*rho(a*x)} of one-dimensional homogeneous and inhomogeneous single-particle densities. We observe that the pair density develops remarkable multiscale patterns which strongly depend on both the particle number and the "width" 1/a of the single-particle density. The simulation results are confirmed by rigorous asymptotic results in the limiting regimes a>>1 and a <

In this work we give a comprehensive derivation of an exact and numerically feasible method to perform ab-initio calculations of quantum particles interacting with a quantized electromagnetic field. We present a hierachy of density-functional-type theories that describe the interaction of charged particles with photons and introduce the appropriate Kohn-Sham schemes. We show how the evolution of a system described by quantum electrodynamics in Coulomb gauge is uniquely determined by its initial state and two reduced quantities. These two fundamental observables, the polarization of the Dirac field and the vector potential of the photon field, can be calculated by solving two coupled, non-linear evolution equations without the need to explicitly determine the (numerically infeasible) many-body wave function of the coupled quantum system. To find reliable approximations to the implicit functionals we present the according Kohn-Sham construction. In the non-relativistic limit this density-functional-type theory of quantum electrodynamics reduces to the density-functional reformulation of the Pauli-Fierz Hamiltonian, which is based on the current density of the electrons and the vector potential of the photon field. By making further approximations, e.g. restricting the allowed modes of the photon field, we derive further density functional-type theories of coupled matter-photon systems for the corresponding approximate Hamiltonians. In the limit of only two sites and one mode we deduce the according effective theory for the two-site Hubbard model coupled to one photonic mode. This model system is used to illustrate the basic ideas of a density-functional reformulation in great detail and we present the exact Kohn-Sham potentials for our coupled matter-photon model system.

The Hohenberg-Kohn density functional was long ago shown to reduce to the Thomas-Fermi 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 supports the conjecture that this is also true for correlation, but much less relevant to atoms. We illustrate how expansions around 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 the correlation energy tends to $-A_c Z ln Z + B_c Z$ as $Z$ tends to infinity, where $Z$ is the atomic number, $A_c$ is known, and we estimate $B_c$ to be about 37 millihartrees. The local density approximation yields $A_c$ exactly, but a very incorrect value for $B_c$, showing that the local approximation is less relevant for correlation alone. This limit is a benchmark for the non-empirical construction of density functional approximations. We conjecture that, beyond atoms, the leading correction to the local density approximation in the large-$Z$ limit generally takes this form, but with $B_c$ a functional of the TF density for the system. The implications for construction of approximate density functionals are discussed.

This paper overviews various phenomena related to the concept of isospin symmetry. The focus is on N~Z nuclei, which are excellent laboratories of isospin physics. The theoretical framework applied is nuclear Density Functional Theory and its isospin- and angular-momentum projected extensions, as well as symmetry-projected multi-reference models. The topics covered include: isospin impurities, superallowed beta decays, beta-transitions in mirror nuclei, isospin breaking hadronic interactions, mirror and triplet binding energy differences, and isoscalar pairing.

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

A new method is presented for performing first-principles molecular-dynamics simulations of systems with variable occupancies. We adopt a matrix representation for the one-particle statistical operator Gamma, to introduce a ``projected'' free energy functional G that depends on the Kohn-Sham orbitals only and that is invariant under their unitary transformations. The Liouville equation [ Gamma , H ] = 0 is always satisfied, guaranteeing a very efficient and stable variational minimization algorithm that can be extended to non-conventional entropic formulations or fictitious thermal distributions.

A new method is presented for performing first-principles molecular-dynamics simulations of systems with variable occupancies. We adopt a matrix representation for the one-particle statistical operator Gamma, to introduce a ``projected'' free energy functional G that depends on the Kohn-Sham orbitals only and that is invariant under their unitary transformations. The Liouville equation [ Gamma , H ] = 0 is always satisfied, guaranteeing a very efficient and stable variational minimization algorithm that can be extended to non-conventional entropic formulations or fictitious thermal distributions.

For the calculation of neutral excitations, time-dependent density functional theory (TDDFT) is an exact reformulation of the many-body time-dependent Schr ̈ dinger equation, based o on knowledge of the density instead of the many-body wavefunction. The density can be determined in an efficient scheme by solving one-particle non-interacting Schr ̈ dinger o equations�the Kohn�Sham equations. The complication of the problem is hidden in the� unknown�time-dependent exchange and correlation potential that appears in the Kohn�Sham equations and for which it is essential to find good approximations. Many approximations have been suggested and tested for finite systems, where even the very simple adiabatic local-density approximation (ALDA) has often proved to be successful. In the case of solids, ALDA fails to reproduce optical absorption spectra, which are instead well described by solving the Bethe� Salpeter equation of many-body perturbation theory (MBPT). On the other hand, ALDA can lead to excellent results for loss functions (at vanishing and finite momentum transfer). In view of this and thanks to recent successful developments of improved linear-response kernels derived from MBPT, TDDFT is today considered a promising alternative to MBPT for the calculation of electronic spectra, even for solids. After reviewing the fundamentals of TDDFT within linear response, we discuss different approaches and a variety of applications to extended systems. Lecture Notes Collection FreeScience.info ID2183 Obtained from http://www.fhi-berlin.mpg.de/th/publications/RPP-70-357-2007.pdf http://www.freescience.info/go.php?pagename=books&id=2183

Using first-principles calculations, we systematically study the adsorption behavior of a single molecular H$_{2}$O on the Be(0001) surface. We find that the favored molecular adsorption site is the top site with an adsorption energy of about 0.3 eV, together with the detailed electronic structure analysis, suggesting a weak binding strength of the H$_{2}$O/Be(0001) surface. The adsorption interaction is mainly contributed by the overlapping between the $s$ and $p_{z}$ states of the top-layer Be atom and the molecular orbitals 1$b_{1}$ and 3$a_{1}$ of H$_{2}$O. The activation energy for H$_{2}$O diffusion on the surface is about 0.3 eV. Meanwhile, our study indicates that no dissociation state exists for the H$_{2}$O/Be(0001) surface.

The effect of epitaxial strain on the cation distribution in spinel ferrites CoFe2O4 and NiFe2O4 is investigated by GGA+U total energy calculations. We obtain a very strong (moderate) tendency for cation inversion in NiFe2O4 (CoFe2O4), in agreement with experimental bulk studies. This preference for the inverse spinel structure is reduced by tensile epitaxial strain, which can lead to strong sensitivity of the cation distribution on specific growth conditions in thin films. Furthermore, we obtain significant energy differences between different cation arrangements with the same degree of inversion, providing further evidence for recently proposed short range B site order in NiFe2O4.

In this tutorial presentation, we give a comprehensive introduction into the Gutzwiller variational approach and its merger with the density functional theory. The merits of this method are illustrated by a discussion of results for two-band Hubbard models and realistic multi-band Hamiltonians for iron-pnictide systems.

In this tutorial presentation, we give a comprehensive introduction into the Gutzwiller variational approach and its merger with the density functional theory. The merits of this method are illustrated by a discussion of results for two-band Hubbard models and realistic multi-band Hamiltonians for iron-pnictide systems.

In this article we obtain a fundamental measure functional for the model of aligned hard hexagons in the plane. Our aim is not just to provide a functional for a new, admittedly academic, model, but to investigate the structure of fundamental measure theory. A model of aligned hard hexagons has similarities with the hard disk model. Both share "lost cases", i.e. admit configurations of three particles in which there is pairwise overlap but not triple overlap. These configurations are known to be problematic for fundamental measure functionals, which are not able to capture their contribution correctly. This failure lies in the inability of these functionals to yield a correct low density limit of the third order direct correlation function. Here we derive the functional by projecting aligned hard cubes on the plane x+y+z=0. The correct dimensional crossover behavior of these functionals permits us to follow this strategy. The functional of aligned hard cubes, however, does not have lost cases, so neither had the resulting functional for aligned hard hexagons. The latter exhibits, in fact, a peculiar structure as compared to the one for hard disks. It depends on a uniparametric family of weighted densities through a new term not appearing in the functional for hard disks. Apart from studying the freezing of this system, we discuss the implications of the functional structure for new developments of fundamental measure theory.

A new method for calculating optical absorption spectra within linear-scaling density-functional theory (LS-DFT) is presented, incorporating a scheme for optimizing a set of localized orbitals to accurately represent unoccupied Kohn-Sham states. Three different schemes are compared and the most promising of these, based on the use of a projection operator, has been implemented in a fully-functional LS-DFT code. The method has been applied to the calculation of optical absorption spectra for the metal-free phthalocyanine molecule and the conjugated polymer poly(para-phenylene). Excellent agreement with results from a traditional DFT code is obtained.

The study of the spin excitations in antiferromagnetic (AFM) and ferromagnetic (FM) phases of FeRh is reported. We demonstrate that although the Fe atomic moments are well defined there is a number of important phenomena absent in the Heisenberg description: Landau damping of spin waves, large Rh moments induced by the AFM magnons, the formation of the optical magnons terminated by Stoner excitations. We relate the properties of the spin-wave damping to the features of the Stoner continuum and compare the chirality of the spin excitations in AFM, FM and paramagnetic (PM) systems.

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

Reduced density-matrix functional theory (RDMFT) has become an appealing alternative to density-functional theory to describe electronic properties of highly-correlated systems. Here we derive exact conditions for the suitability of RDMFT to describe the two-dimensional homogeneous electron gas, which is the base system for, e.g., semiconductor quantum dots and quantum Hall devices. Following the method of Cioslowski and Pernal [J. Chem. Phys. {\bf 111}, 3396 (1999)] we focus on the properties of power functionals of the form $f(n,n')=(n n')^\alpha$ for the scaling function in the exchange-correlation energy. We show that in order to have stable and analytic solutions, and for $f$ to satisfy the homogeneous scaling constraint, the power is restricted to $1/4 \leq \alpha \leq 3/4$. Applying a reasonable ansatz for the momentum distribution and the lower bound for the exchange-correlation energy tightens the physical regime further to $0.64 \lesssim \alpha \leq 0.75$.

Density Functional Resonance Theory (DFRT) is a complex-scaled version of ground-state Density Functional Theory (DFT) that allows one to calculate the resonance energies and lifetimes of metastable anions. In this formalism, the exact energy and lifetime of the lowest-energy resonance of unbound systems is encoded into a complex "density" that can be obtained via complex-coordinate scaling. This complex density is used as the primary variable in a DFRT calculation just as the ground-state density would be used as the primary variable in DFT. As in DFT, there exists a mapping of the N-electron interacting system to a Kohn-Sham system of N non-interacting particles in DFRT. This mapping facilitates self consistent calculations with an initial guess for the complex density, as illustrated with an exactly-solvable model system. Whereas DFRT yields in principle the exact resonance energy and lifetime of the interacting system, we find that neglecting the complex-correlation contribution leads to errors of similar magnitude to those of standard scattering close-coupling calculations under the bound-state approximation.

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

A new "on the fly" method to perform Born-Oppenheimer ab initio molecular dynamics (AIMD) is presented. Inspired by Ehrenfest dynamics in time-dependent density functional theory, the electronic orbitals are evolved by a Schroedinger-like equation, where the orbital time derivative is multiplied by a parameter. This parameter controls the time scale of the fictitious electronic motion and speeds up the calculations with respect to standard Ehrenfest dynamics. In contrast to other methods, wave function orthogonality needs not be imposed as it is automatically preserved, which is of paramount relevance for large scale AIMD simulations.

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.

A coarse grained model for flexible polymers end-grafted to repulsive spherical nanoparticles is studied for various chain lengths and grafting densities under good solvent conditions, by Molecular Dynamics methods and density functional theory. With increasing chain length the monomer density profile exhibits a crossover to the star polymer limit. The distribution of polymer ends and the linear dimensions of individual polymer chains are obtained, while the inhomogeneous stretching of the chains is characterized by the local persistence lengths. The results on the structure factor of both single chain and full spherical brush as well as the range of applicability of the different theoretical tools are presented. Eventually an outlook on experiments is given.

We have studied solid hydrogen up to pressures of 300 GPa and temperatures of 350 K using density functional theory methods and have found "mixed structures" that are more stable than those predicted earlier. Mixed structures consist of alternate layers of strongly bonded molecules and weakly bonded graphene-like sheets. Quasi-harmonic vibrational calculations show that mixed structures are the most stable at room temperature over the pressure range 250-295 GPa. These structures are stabilized with respect to strongly-bonded molecular phases at room temperature by the presence of lower-frequency vibrational modes arising from the graphene-like sheets. Our results for the mixed structures are consistent with the experimental Raman data [M. I. Eremets and I. A. Troyan, Nature Mater. 10 927 (2011) and R. T. Howie et al. Phys. Rev. Lett. 108 125501 (2012)]. We find that mixed phases are reasonable structural models for phase IV of hydrogen.

We consider two perturbative schemes to calculate excitation energies, each employing the Kohn-Sham Hamiltonian as the unperturbed system. Using accurate exchange-correlation potentials generated from essentially exact densities and their exchange components determined by a recently proposed method, we evaluate energy differences between the ground state and excited states in first-order perturbation theory for the Helium, ionized Lithium and Beryllium atoms. It was recently observed that the zeroth-order excitations energies, simply given by the difference of the Kohn-Sham eigenvalues, almost always lie between the singlet and triplet experimental excitations energies, corrected for relativistic and finite nuclear mass effects. The first-order corrections provide about a factor of two improvement in one of the perturbative schemes but not in the other. The excitation energies within perturbation theory are compared to the excitations obtained within $\Delta$SCF and time-dependent density functional theory. We also calculate the excitation energies in perturbation theory using approximate functionals such as the local density approximation and the optimized effective potential method with and without the Colle-Salvetti correlation contribution.

Recent experiments have shown that the striking structure formation in dewetting films of evaporating colloidal nanoparticle suspensions occurs in an ultrathin `postcursor' layer that is left behind by a mesoscopic dewetting front. Various phase change and transport processes occur in the postcursor layer, that may lead to nanoparticle deposits in the form of labyrinthine, network or strongly branched `finger' structures. We develop a versatile dynamical density functional theory to model this system which captures all these structures and may be employed to investigate the influence of evaporation/condensation, nanoparticle transport and solute transport in a differentiated way. We highlight, in particular, the influence of the subtle interplay of decomposition in the layer and contact line motion on the observed particle-induced transverse instability of the dewetting front.

We present a new time-dependent Density Functional approach to study the relaxational dynamics of an assembly of interacting particles subject to thermal noise. Starting from the Langevin stochastic equations of motion for the velocities of the particles we are able by means of an approximated closure to derive a self-consistent deterministic equation for the temporal evolution of the average particle density. The closure is equivalent to assuming that the equal-time two-point correlation function out of equilibrium has the same properties as its equilibrium version. The changes in time of the density depend on the functional derivatives of the grand canonical free energy functional $F[\rho]$ of the system. In particular the static solutions of the equation for the density correspond to the exact equilibrium profiles provided one is able to determine the exact form of $F[\rho]$. In order to assess the validity of our approach we performed a comparison between the Langevin dynamics and the dynamic density functional method for a one-dimensional hard-rod system in three relevant cases and found remarkable agreement, with some interesting exceptions, which are discussed and explained. In addition, we consider the case where one is forced to use an approximate form of $F[\rho]$. Finally we compare the present method with the stochastic equation for the density proposed by other authors [Kawasaki,Kirkpatrick etc.] and discuss the role of the thermal fluctuations.

The density-functional (DF) theory provides a simple method for calculating the properties of an interacting system under an external potential by associating it with a corresponding non-interacting system. Here, we find some relations in this non-interacting system, which can enable us to set up the thermodynamical relations for a neutral electron-nucleus mixture in terms of quantities of the non-interacting system and the exchange-correlation effect. In this way, Andersen's force theorems are clearly described and easily proved in conjunction with other types of force theorems, Janak's theorem extended to a finite-temperature system, pressure formulas and some thermodynamical relations. For this purpose, the thermodynamics in the DF theory is presented in a systematic and explicit way.

Based on the work of Gorling and that of Levy and Nagy, density-functional formalism for many Fermionic excited-states is explored through a careful and rigorous analysis of the excited-state density to external potential mapping. It is shown that the knowledge of the ground-state density is a must to fix the mapping from an excited-state density to the external potential. This is the excited-state counterpart of the Hohenberg-Kohn theorem, where instead of the ground-state density the density of the excited-state gives the true many-body wavefunctions of the system. Further, the excited-state Kohn-Sham system is defined by comparing it's non-interacting kinetic energy with the true kinetic energy. The theory is demonstrated by studying a large number of atomic systems.

We study the interactions between two negatively charged macroscopic surfaces confining positive counterions. A density-functional approach is introduced which, besides the usual mean-field interactions, takes into account the correlations in the positions of counterions. The excess free energy is derived in the framework of the Debye-H{\"u}ckel theory of the one-component plasma, with the homogeneous density replaced by a weighted density. The minimization of the total free energy yields the density profile of the microions. The pressure is calculated and compared with the simulations and the results derived from integral equations theories. We find that the interaction between the two plates becomes attractive when their separation distance is sufficiently small and the surface charge density is larger than a threshold value.

Even though chiral nematic phases were the first liquid crystals experimentally observed more than a century ago, the origin of the thermodynamic stability of cholesteric states is still unclear. In this Letter we address the problem by means of a novel density functional theory for the equilibrium pitch of chiral particles. When applied to right-handed hard helices, our theory predicts an entropy-driven cholesteric phase, which can be either right- or left-handed, depending not only on the particle shape but also on the thermodynamic state. We explain the origin of the chiral ordering as an interplay between local nematic alignment and excluded-volume differences between left- and right-handed particle pairs.

Due to the large value of the scattering length in nuclear systems, standard density--functional theories based on effective interactions usually fail to reproduce the nuclear Fermi liquid behavior both at very low densities and close to equilibrium. Guided on one side by the success of the Skyrme density functional and, on the other side, by resummation techniques used in Effective Field Theories for systems with large scattering lengths, a new energy--density functional is proposed. This functional, adjusted on microscopic calculations, reproduces the nuclear equations of state of neutron and symmetric matter at various densities. Furthermore, it provides reasonable saturation properties as well as an appropriate density dependence for the symmetry energy.

We introduce a new class of exchange-correlation potentials for a static and time-dependent Density Functional Theory of strongly correlated systems in 3D. The potentials are obtained via Dynamical Mean Field Theory and, for strong enough interactions, exhibit a discontinuity at half filling density, a signature of the Mott transition. For time-dependent perturbations, the dynamics is described in the adiabatic local density approximation. Results from the new scheme compare very favorably to exact ones in clusters. As an application, we study Bloch oscillations in the 3D Hubbard model.

This article is from International Journal of Molecular Sciences , volume 13 . Abstract Density functional theory has been employed to investigate the deformation behaviors of glassy Fe-Si-B model systems prepared by ab initio molecular dynamics. The atomistic deformation defects which are closely related to the local dilation volumes or excess volumes and unstable bonding have been systematically analyzed. It has been found that the icosahedral structures are relatively stable under shear deformation until fracture occurs. Plastic flow is indicated by interruption of percolating icosahedral structures, caused by unstable Fe-Si bonding of p-s hybridization in nature.

The current status of the application of covariant density functional theory to the description of fission barriers in actinides and superheavy nuclei is reviewed. The achievements and open problems are discussed.

Electronic structure calculations within screened-exchange hybrid density functional theory show that Cu(In,Ga) antisites in both CuInSe2 and CuGaSe2 are localized hole traps, which can be attributed to the experimentally observed N2 level. In contrast, GaCu antisites and their defect complexes with copper vacancies exhibit an electron trap level, which can limit the open-circuit voltage and efficiency in Ga-rich Cu(In,Ga)Se2 alloys. Low-temperature photoluminescence measurements in CuGaSe2 thin-film solar cells show a free-to-bound transition at an energy of 1.48 eV, in very good agreement with the calculated transition energy for the GaCu antisite. Since the intrinsic DX center does not exhibit a pinning level within the band gap of CuInSe2, metastable DX behaviour can only be expected for GaCu antisites.

In search of the potential realization of novel normal-state phases on the surface of Sr2RuO4 - those stemming from either topological bulk properties or the interplay between spin-orbit coupling (SO) and the broken symmetry of the surface - we revisit the electronic structure of the top-most layers by ARPES with improved data quality as well as ab-initio LDA slab calculations. We find that the current model of a single surface layer (\surd2x\surd2)R45{\deg} reconstruction does not explain all detected features. The observed depth-dependent signal degradation, together with the close quantitative agreement with LDA+SO slab calculations based on the LEED-determined surface crystal structure, reveal that (at a minimum) the sub-surface layer also undergoes a similar although weaker reconstruction. This points to a surface-to-bulk progression of the electronic states driven by structural instabilities, with no evidence for Dirac and Rashba-type states or surface magnetism.

Ensemble density-functional theory (eDFT) suffers from the so-called "ghost interaction" error when approximate exchange-correlation functionals are used. In this work, we present a rigorous ghost interaction correction (GIC) scheme in the context of range-separated eDFT. The method relies on an exact decomposition of the ensemble short-range exchange-correlation energy into a multideterminantal exact exchange term, which involves the long-range interacting ensemble density matrix instead of the Kohn--Sham (KS) one, and a complementary density-functional correlation energy. A generalized adiabatic connection formula is derived for the latter. In order to perform practical calculations, the complementary correlation functional has been simply modeled by its ground-state local density approximation (LDA) while long-range interacting ground- and excited-state wavefunctions have been obtained self-consistently by combining a long-range configuration interaction calculation with a short-range LDA potential. We show that GIC reduces the curvature of approximate range-separated ensemble energies drastically while providing considerably more accurate excitation energies, even for charge-transfer and double excitations. Interestingly, the method performs well also in the context of standard KS-eDFT, which is recovered when the range-separation parameter is set to zero.