A variational criterion is used to find approximate eigenfunctions and eigenvalues of the Liouville operator in classical many-body systems. The trial functions are taken to be sums over molecules of functions depending on the position and momentum of a single molecule. In a harmonic lattice, this approach leads to exact eigenfunctions and eigenvalues. In a fluid, the eigenvalue spectrum is continuous, and the eigenfunctions are related to those found by Van Kampen in his study of the linearized Vlasov equation for a plasma. The time dependence of the fluid current density is found by means of these eigenfunctions and eigenvalues. The results show persistent free-particle propagation and damped sound-wave propagation, with relative importance depending on the magnitude of the sound velocity.
A variational criterion is used to find approximate eigenfunctions and eigenvalues of the Liouville operator in classical many-body systems. The trial functions are taken to be sums over molecules of functions depending on the position and momentum of a single molecule. In a harmonic lattice, this approach leads to exact eigenfunctions and eigenvalues. In a fluid, the eigenvalue spectrum is continuous, and the eigenfunctions are related to those found by Van Kampen in his study of the linearized Vlasov equation for a plasma. The time dependence of the fluid current density is found by means of these eigenfunctions and eigenvalues. The results show persistent free-particle propagation and damped sound-wave propagation, with relative importance depending on the magnitude of the sound velocity.
A variational solution of the Liouville equation is used to obtain collective coordinates for simple classical fluids. Comparison is made with macroscopic analogs determined from relevant hydrodynamic equations. Approximate eigenfunctions of the Liouville operator, when constructed from linear combinations of the spatial Fourier transforms (fluctuations) of the mass, momentum, and energy-density operators for the fluid, are found appropriate to low-frequency disturbances. When time derivatives of fluctuations are included among trial functions, the variational procedure provides fluid-state counterparts of usual solidstate phonons.
By means of a quantum mechanical phase space distribution function
introduced by von Roos, the Schroedinger equation for a non-relativistic
system of N identical particles with scalar interactions is transformed
into a quantum mechanical generalization of the Liouville equation,
thereby formulating the problem in terms of a generalized density in
phase space, a quantity of primary interest in most treatments of the
corresponding classical system (or plasma). This transformation permits
a parallel development of the theories of classical and quantum plasmas
and thus allows the quantum many-body problem to be discussed virtually
completely in classical terms. In particular, a kinetic theory of
quantum plasmas is obtained by deriving the quantum analogue of the BBGKY
hierarchy, and applying thereto approximation techniques similar to those
of Rostocker and Rosenbluth, and Bogoliubov. The point of departure from
similar previous studies based on the Wigner distribution function is that
the proper exchange symmetry can be tractably introduced into the formalism. Attention is first focused on the Hartree and Hartree-Fock approximations,
in which case the quantum BBGKY system reduces to a simple quantum
generalization of the Vlasov equation. This equation is used to study
the response of spatially homogeneous systems to weak external forces,
and the associated problems of plasmon and spin-wave excitations. It is
also used to derive the quantum and exchange corrected equations of
inviscid hydrodynamical transport which are then applied to the problem
of sound propagation in the degenerate electron gas. The second part of the study is concerned with the theory of the
many-electron atom in the Hartree and Hartree-Fock approximations.
The relevant quantum Vlasov equations lead naturally to a statistical
theory of the atom which reduces to the Thomas-Fermi-Amaldi and Thomas Fermi
models (respectively) as ħ → 0. For ħ ≠ 0, the quantum and
exchange corrections to these models are simultaneously generated. The
quantum hydrodynamical theory developed earlier is used to determine
the influence of these corrections on the boundary conditions of the
model, and a theory of the compressed atom is consequently obtained.
Considered in somewhat less detail are the effects of non-zero temperature,
net orbital angular momentum, relativity and correlations, as well as
time dependent processes. The final part deals with the problem of the degenerate electron
gas with a uniform neutralizing background. Going beyond the Hartree-Fock approximation, the pair functions for particles with
parallel and anti-parallel spin are obtained by neglecting three particle
correlations. From these functions, a quantum-mechanical collision
integral is derived which differs from that obtained by Silin
and Guernsey and conjectured by Wyld and Pines in that dynamical exchange
effects are included. Also obtained from the pair function
is an expression for the correlation which reduces in the high
density limit to the result of Gell-Mann and Brueckner. At intermediate
densities an additional term appears in the energy due to the screening
of the exchange interaction by the dielectric properties of the medium.
It is evaluated in the high density limit and found to be -0.151 r_s ln r_s
Ryd/electron in marked disagreement with the corresponding value obtained
by DuBois.
We derive a closed equation of motion for the current density of an inhomogeneous quantum many-body system under the assumption that the time-dependent wave function can be described as a geometric deformation of the ground-state wave function. By describing the many-body system in terms of a single collective field we provide an alternative to traditional approaches, which emphasize one-particle orbitals. We refer to our approach as continuum mechanics for quantum many-body systems. In the linear response regime, the equation of motion for the displacement field becomes a linear fourth-order integrodifferential equation, whose only inputs are the one-particle density matrix and the pair-correlation function of the ground state. The complexity of this equation remains essentially unchanged as the number of particles increases. We show that our equation of motion is a Hermitian eigenvalue problem, which admits a complete set of orthonormal eigenfunctions under a scalar product that involves the ground-state density. Further, we show that the excitation energies derived from this approach satisfy a sum rule which guarantees the exactness of the integrated spectral strength. Our formulation becomes exact for systems consisting of a single particle and for any many-body system in the high-frequency limit. The theory is illustrated by explicit calculations for simple one- and two-particle systems.
Variational calculations play an important role in quantum mechanics, particularly in determining the ground-state energy of physical systems. Didactic examples in which one can see how the exact energy and wave function are approached by using a family of trial functions are not very numerous. An example was given in a recent book [M. Moshinsky, The Harmonic Oscillator in Modern Physics: From Atoms to Quarks (Gordon and Breach, New York, 1969)] which discusses the ground state of the hydrogen atom using as trial wave function a linear combination of harmonic-oscillator states. In the present paper we carry out a similar analysis for the ground state of a particle of mass m in a three-dimensional square-well potential. We discuss not only the approach to the exact energy when we vary the frequency and the number of oscillator states but also analyze the overlap of the exact and variational wave functions and compare the exact and approximate form factors.
Exact solutions to the Schrödinger equation for realistic nanoscale systems are beyond reach, hence, different strategies for approximating the solutions are necessary. Perturbation theory relies on a “zero-order Hamiltonian” to express the desired eigenvalues and eigenvectors of “the full Hamiltonian.” We derive working equations for the Rayleigh–Schrödinger perturbation theory, and the validity of the approach is analyzed for a generic two-level system. Applications are given to atoms perturbed by point charges or static fields, and for electrons in quantum wells. An alternative strategy is the variation method, which replaces exact solutions by their projection on a reduced space of “trial functions,” varied to minimize the associated error. Particularly important is the method of linear variation, which can potentially converge to the exact Hamiltonian eigenstates. The mean-field approximation, commonly used for many-particle systems, is derived by optimizing a trial function in the form of product of single-particle functions.
Abstract We present a kinetic description of superfluid currents in ring-shaped Bose-Einstein condensates based on the Wigner phase-space formalism. Starting from the Gross-Pitaevskii equation in a toroidal geometry, we derive a Vlasov-type equation for the angular Wigner function, in which the mean-field interaction generates an effective force proportional to the density gradient. Within this framework, we obtain the dispersion relation of collective modes and recover the Bogoliubov spectrum in the long-wavelength limit. We show that the Landau criterion for superfluidity can be interpreted as the absence of resonant phase-space trajectories satisfying the condition ω=qvℓ ω = q v ℓ . In a ring geometry, the quantization of angular momentum leads to a discrete set of velocities, which suppresses the availability of resonant states and strongly inhibits Landau damping. In contrast, in the continuous limit R→∞ R → ∞ , the spectrum becomes quasi-continuous and the standard Landau damping mechanism is recovered, establishing a direct connection between kinetic resonances and the energetic criterion for superfluidity. We further analyze the role of Bogoliubov depletion by considering a finite-width angular momentum distribution. Although resonant states formally exist in this case, we show that, for flow velocities below the sound velocity, the phase-space distribution does not provide the gradients required for energy transfer, and the superfluid current remains dynamically stable. Our results provide a unified phase-space interpretation of superfluidity, highlighting the role of angular momentum quantization and the structure of the distribution function in determining the stability of persistent currents. Graphic abstract Phase-space interpretation of superfluid stability in a ring Bose-Einstein condensate. For flow velocities below the sound velocity, resonant phase-space trajectories are unavailable and the superflow remains stable, whereas above the sound velocity resonant coupling becomes possible, leading to kinetic instability