We present an ab initio approach for grand-canonical ensembles in thermal equilibrium (eq) with local or nonlocal external potentials based on the one-reduced density matrix (1RDM). 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 1RDM. We further prove the existence of a Kohn-Sham system capable of reproducing the 1RDM 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.
We present a method to compute the Fermi function of the Hamiltonian for a system of independent fermions based on an exact decomposition of the grand-canonical potential. This scheme does not rely on the localization of the orbitals and is insensitive to ill-conditioned Hamiltonians. It lends itself naturally to linear scaling as soon as the sparsity of the system's density matrix is exploited. By using a combination of polynomial expansion and Newton-like iterative techniques, an arbitrarily large number of terms can be employed in the expansion, overcoming some of the difficulties encountered in previous papers. Moreover, this hybrid approach allows us to obtain a very favorable scaling of the computational cost with increasing inverse temperature, which makes the method competitive with other Fermi operator expansion techniques. After performing an in-depth theoretical analysis of computational cost and accuracy, we test our approach on the density functional theory Hamiltonian for the metallic phase of the LiAl alloy.
We study how the grand-canonical density matrix arises in macroscopic quantum systems. ``Canonical typicality'' is the known statement that for a typical wave function
Ψ from a micro-canonical energy shell of a quantum system
S weakly coupled to a large but finite quantum system
B, the reduced density matrix
ρ^ΨS=trB∣Ψ⟩⟨Ψ∣ is approximately equal to the canonical density matrix
ρ^can=Zcan−1exp(−βH^S). Here, we discuss the analogous statement and related questions for the \emph{grand-canonical} density matrix
ρ^gc=Zgc−1exp(−β(H^S−μ1N^1S−…−μrN^rS)) with
N^iS the number operator for molecules of type
i in the system
S. This includes (i) the case of chemical reactions (which requires some novel considerations) and (ii) that of systems
S defined by a spatial region which particles may enter or leave. It includes statements about how
ρ^gc arises from the density matrix of the appropriate (generalized micro-canonical) Hilbert subspace
Hgmc⊂HS⊗HB (defined by a micro-canonical interval of total energy and suitable particle number sectors) or from typical
Ψ in
Hgmc, as well as statements about the distribution of the (conditional) wave function
ψS of
S, which turns out to be a so-called GAP or Scrooge measure. That is, we discuss the foundation and justification of both the density matrix and the distribution of the wave function in the grand-canonical case. To this end (particularly for the chemical reactions), we also need to extend these considerations to the so-called generalized Gibbs ensembles, which apply to systems for which some macroscopic observables are conserved.
In a previous study, the author formulated the density functional theory in terms of functional integrals. It was valid at zero and finite temperature. It was possible to derive the Hohenberg and Kohn formulation at zero temperature and the Mermin formulation at finite temperature of the density functional theory, which states that the energy or the grand potential are functionals of the true density of the system considered. In particular, the Kohn and Sham equations are proven to appear naturally by performing a saddle-point evaluation of a specific functional integral. This result is valid at zero or finite temperature. Unfortunately, the expression of the grand potential given in our previous work differs from the usual expression found in the literature. In this short paper, we derive the common expression of the grand potential in the framework of the density functional theory by starting from the expression given in this previous work. This completes the formulation of the density functional theory using functional integrals. This work could be of interest to people working in the field of quantum Monte Carlo methods at finite temperature.
We introduce a quantum Monte Carlo method at finite temperature for interacting fermionic models in the canonical ensemble, where the conservation of the particle number is enforced. Although general thermodynamic arguments ensure the equivalence of the canonical and the grand-canonical ensembles in the thermodynamic limit, their approach to the infinite-volume limit is distinctively different. Observables computed in the canonical ensemble generically display a finite-size correction proportional to the inverse volume, whereas in the grand-canonical ensemble the approach is exponential in the ratio of the linear size over the correlation length. We verify these predictions by quantum Monte Carlo simulations of the Hubbard model in one and two dimensions in the grand-canonical and the canonical ensemble. We prove an exact formula for the finite-size part of the free energy density, energy density and other observables in the canonical ensemble and relate this correction to a susceptibility computed in the corresponding grand-canonical ensemble. This result is confirmed by an exact computation of the one-dimensional classical Ising model in the canonical ensemble, which for classical models corresponds to the so-called fixed-magnetization ensemble. Our method is useful for simulating finite systems which are not coupled to a particle bath, such as in nuclear or cold atom physics.