No Arabic abstract
Finite temperature density functional theory provides, in principle, an exact description of the thermodynamical equilibrium of many-electron systems. In practical applications, however, the functionals must be approximated. Efficient and physically meaningful approximations can be developed if relevant properties of the exact functionals are known and taken into consideration as constraints. In this work, derivations of exact properties and scaling relations for the main quantities of finite temperature density functional theory are presented. In particular, a coordinate scaling transformation at finite temperature is introduced and its consequences are elucidated.
A previous analysis of scaling, bounds, and inequalities for the non-interacting functionals of thermal density functional theory is extended to the full interacting functionals. The results are obtained from analysis of the related functionals from the equilibrium statistical mechanics of thermodynamics for an inhomogeneous system. Their extension to the functionals of density functional theory is described.
Finite temperature density functional theory requires representations for the internal energy, entropy, and free energy as functionals of the local density field. A central formal difficulty for an orbital-free representation is construction of the corresponding functionals for non-interacting particles in an arbitrary external potential. That problem is posed here in the context of the equilibrium statistical mechanics of an inhomogeneous system. The density functionals are defined and shown to be equal to the extremal state for a functional of the reduced one-particle statistical operators. Convexity of the latter functionals implies a class of general inequalities. First, it is shown that the familiar von Weizsacker lower bound for zero temperature functionals applies at finite temperature as well. An upper bound is obtained in terms of a single-particle statistical operator corresponding to the Thomas-Fermi approximation. Next, the behavior of the density functionals under coordinate scaling is obtained. The inequalities are exploited to obtain a class of upper and lower bounds at constant temperature, and a complementary class at constant density. The utility of such constraints and their relationship to corresponding results at zero temperature are discussed.
Since the Time-Dependent Density Functional Theory is mathematically formulated through non-linear coupled time-dependent 3-dimensional partial differential equations it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrodinger equation for an interacting many-body system is however linear and (mathematically) the exact equations of the Density Functional Theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a Density Functional Theory framework. Whether for realistic implementations of the Time-Dependent Density Functional Theory the question of absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density functional theory is a unique tool allowing us the study of non-equilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Our analysis, for a number of quantum superfluid any-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom ${cal O}(10^{10})$ and larger, suggests that its maximum Lyapunov are negligible for all practical purposes.
We present an ab-initio approach for grand canonical ensembles in thermal equilibrium with local or nonlocal external potentials based on the one-reduced density matrix. We show that equilibrium properties of a grand canonical ensemble are determined uniquely by the eq-1RDM and establish a variational principle for the grand potential with respect to its one-reduced density matrix. We further prove the existence of a Kohn-Sham system capable of reproducing the one-reduced density matrix of an interacting system at finite temperature. Utilizing this Kohn-Sham system as an unperturbed system, we deduce a many-body approach to iteratively construct approximations to the correlation contribution of the grand potential.
The magnetic properties of the intermetallic compound FeAl are investigated using exact exchange density functional theory. This is implemented within a state of the art all-electron full potential method. We find that FeAl is magnetic with a moment of 0.70 $mu_B$, close to the LSDA result of 0.69 $mu_B$. A comparison with the non-magnetic density of states with experimental negative binding energy result shows a much better agreement than any previous calculations. We attribute this to the fine details of the exchange field, in particular its asymmetry, which is captured very well with the orbital dependent exchange potential.