Physical Review

Theory of a System Containing an Arbitrary Number of Different Interacting Bose Species

Richard David Etters, J. C. Raich

Published Jul 5, 1967• 2 citations

Abstract

A theory that describes the dynamic behavior of a system containing an arbitrary number of different interacting Bose species is developed. The excited eigenstates of the system are generated by the superposition of plane-wave creation- and annihilation-operator products ${{\ensuremath{\mu}}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}=\ensuremath{\Sigma}{\ensuremath{\gamma}}^{}\ensuremath{\Sigma}{p}^{}g(p,k,\ensuremath{\gamma}){{a}_{p+k}}^{\ifmmode\dagger\else\textdagger\fi{}}(\ensuremath{\gamma}){a}_{p}(\ensuremath{\gamma})$, where the $c$-number coefficients $g(p,k,\ensuremath{\gamma})$ are determined by solving the equation of motion for ${{\ensuremath{\mu}}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}$ in an approximate way. The sum over $\ensuremath{\gamma}$ extends over all different species of the system. The equation of motion is solved by a method similar to that originally derived by Pines and Bohm for the electron gas. Results are obtained for the energy spectrum of the elementary excitations in the system. In particular, explicit results are obtained for the energy spectrum of a two-species system---first for a weakly coupled gas, and secondly for the low-momentum collective modes of a strongly coupled fluid. The low-momentum excitations of the strongly coupled fluid are completely analogous to the "zero sound" modes commonly associated with Fermi liquids. The theory is then applied to a dilute mixture of ${\mathrm{He}}^{3}$ in ${\mathrm{He}}^{4}$ and the phonon velocity is determined as a function of ${\mathrm{He}}^{3}$ concentration. The probability per unit time $P(\mathbf{k},{\ensuremath{\omega}}_{k})$ that a particle transfers momentum k and energy ${\ensuremath{\omega}}_{k}$ to the system by scattering is also calculated in terms of the dynamic form factor ${S}_{\ensuremath{\gamma}}(\mathbf{k},\ensuremath{\omega})$.

DOI: 10.1103/physrev.159.146

Related papers

Dynamic many-body theory: Dynamics of strongly correlated Fermi fluids

Helga M. Böhm, Robert A. Holler, E. Krotscheck +1 more

arXiv:1010.1816Physical Review BDec 6, 201028 citations

Particle-number-conserving Bogoliubov method which demonstrates the validity of the time-dependent Gross-Pitaevskii equation for a highly condensed Bose gas

Crispin W. Gardiner

arXiv:quant-ph/9703005Physical Review AAug 1, 1997256 citations

Exactly Solvable Many-Boson Model

R. W. Richardson

Journal of Mathematical PhysicsSep 1, 196882 citations

Excitations in a Bose Gas at Finite Temperatures. II. Relation between Single-Particle and Density Fluctuations

Allan Griffin, T. H. Cheung

Physical Review AJun 1, 197336 citations

Collective Modes in Bose-Fermi Mixture Gases

Asano Yoji, Yoji Asano

Mar 17, 20200 citations

Excitation spectrum of a system of interacting bosons

Moorad Alexanian, R. A. Brito

Physical review. B, Condensed matterMay 1, 197811 citations

Physical Review

Theory of a System Containing an Arbitrary Number of Different Interacting Bose Species

Richard David Etters, J. C. Raich

Published Jul 5, 1967• 2 citations

Publisher

Abstract

{{\ensuremath{\mu}}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}=\ensuremath{\Sigma}{\ensuremath{\gamma}}^{}\ensuremath{\Sigma}{p}^{}g(p,k,\ensuremath{\gamma}){{a}_{p+k}}^{\ifmmode\dagger\else\textdagger\fi{}}(\ensuremath{\gamma}){a}_{p}(\ensuremath{\gamma})

, where the

c

-number coefficients

g(p,k,\ensuremath{\gamma})

are determined by solving the equation of motion for

{{\ensuremath{\mu}}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}

in an approximate way. The sum over

\ensuremath{\gamma}

extends over all different species of the system. The equation of motion is solved by a method similar to that originally derived by Pines and Bohm for the electron gas. Results are obtained for the energy spectrum of the elementary excitations in the system. In particular, explicit results are obtained for the energy spectrum of a two-species system---first for a weakly coupled gas, and secondly for the low-momentum collective modes of a strongly coupled fluid. The low-momentum excitations of the strongly coupled fluid are completely analogous to the "zero sound" modes commonly associated with Fermi liquids. The theory is then applied to a dilute mixture of

{\mathrm{He}}^{3}

{\mathrm{He}}^{4}

and the phonon velocity is determined as a function of

{\mathrm{He}}^{3}

concentration. The probability per unit time

P(\mathbf{k},{\ensuremath{\omega}}_{k})

that a particle transfers momentum k and energy

{\ensuremath{\omega}}_{k}

to the system by scattering is also calculated in terms of the dynamic form factor

{S}_{\ensuremath{\gamma}}(\mathbf{k},\ensuremath{\omega})

DOI: 10.1103/physrev.159.146

Related papers

Dynamic many-body theory: Dynamics of strongly correlated Fermi fluids

Helga M. Böhm, Robert A. Holler, E. Krotscheck +1 more

We develop a systematic theory of multiparticle excitations in strongly interacting Fermi systems. Our work is the generalization of the time-honored work by Jackson, Feenberg, and Campbell for bosons, that provides, in its most advanced implementation, quantitative predictions for the dynamic structure function in the whole experimentally accessible energy/momentum regime. Our view is that the same physical effects, namely, fluctuations of the wave function at an atomic length scale are responsible for the correct energetics of the excitations in both Bose and Fermi fluids. Besides a comprehensive derivation of the fermion version of the theory and discussion of the approximations made, we present results for homogeneous $^{3}\text{H}\text{e}$ and electrons in three dimensions. We find indeed a significant lowering of the zero-sound mode in $^{3}\text{H}\text{e}$ and a broadening of the collective mode due to the coupling to two-particle--two-hole excitations in good agreement with experiments. The most visible effect in electronic systems is the appearance of a ``double-plasmon'' excitation.We develop a systematic theory of multiparticle excitations in strongly interacting Fermi systems. Our work is the generalization of the time-honored work by Jackson, Feenberg, and Campbell for bosons, that provides, in its most advanced implementation, quantitative predictions for the dynamic structure function in the whole experimentally accessible energy/momentum regime. Our view is that the same physical effects, namely, fluctuations of the wave function at an atomic length scale are responsible for the correct energetics of the excitations in both Bose and Fermi fluids. Besides a comprehensive derivation of the fermion version of the theory and discussion of the approximations made, we present results for homogeneous $^{3}\text{H}\text{e}$ and electrons in three dimensions. We find indeed a significant lowering of the zero-sound mode in $^{3}\text{H}\text{e}$ and a broadening of the collective mode due to the coupling to two-particle--two-hole excitations in good agreement with experiments. The most visible effect in electronic systems is the appearance of a ``double-plasmon'' excitation.

arXiv:1010.1816Physical Review BDec 6, 201028 citations

Particle-number-conserving Bogoliubov method which demonstrates the validity of the time-dependent Gross-Pitaevskii equation for a highly condensed Bose gas

Crispin W. Gardiner

The Bogoliubov method for the excitation spectrum of a Bose-condensed gas is generalized to apply to a gas with an exact large number $N$ of particles. This generalization yields a description of the Schr\"odinger picture field operators as the product of an annihilation operator $A$ for the total number of particles and the sum of a ``condensate wave function'' $\ensuremath{\xi}(\mathbf{x})$ and a phonon field operator $\ensuremath{\chi}(\mathbf{x})$ in the form $\ensuremath{\psi}(\mathbf{x})\ensuremath{\approx}A{\ensuremath{\xi}(\mathbf{x})+\ensuremath{\chi}(\mathbf{x})/\sqrt{N}}$ when the field operator acts on the $N$ particle subspace. It is then possible to expand the Hamiltonian in decreasing powers of $\sqrt{N}$ , and thus obtain solutions for eigenvalues and eigenstates as an asymptotic expansion of the same kind. It is also possible to compute all matrix elements of field operators between states of different $N$ . The excitation spectrum can be obtained by essentially the same method as Bogoliubov only if $\ensuremath{\xi}(\mathbf{x})$ is a solution of the time-independent Gross-Pitaevskii equation for $N$ particles and any chemical potential $\ensuremath{\mu}$, which yields a valid and stable solution of the Gross-Pitaevskii equation. The treatment within a subspace of fixed $N$ is identical in form to that usually used, but the interpretation of the operators is slightly different. A time-dependent generalization is then made, yielding an asymptotic expansion in decreasing powers of $\sqrt{N}$ for the equations of motion. In this expansion the condensate wave function has the time-dependent form $\ensuremath{\xi}(\mathbf{x},t)$, and the condition for the validity of the expansion is that $\ensuremath{\xi}(\mathbf{x},t)$ satisfies the time-dependent Gross-Pitaevskii equation $\ensuremath{\partial}\ensuremath{\xi}/\ensuremath{\partial}t=\ensuremath{-}({\ensuremath{\Elzxh}}^{2}/2m){\ensuremath{\nabla}}^{2}\ensuremath{\xi}+V\ensuremath{\xi}+Nu|\ensuremath{\xi}{|}^{2}\ensuremath{\xi}$. The physics is then described in a kind of interaction picture, called the condensate picture, in which the phonon operator can be expressed as $\ensuremath{\chi}(\mathbf{x},t)={\ensuremath{\sum}}_{k}{\ensuremath{\xi}}_{k}(\mathbf{x},t){\ensuremath{\alpha}}_{k},$ where the operators ${\ensuremath{\alpha}}_{k}$ are time-independent annihilation operators, and the state vector has a time evolution described by a Schr\"odinger equation in which the Hamiltonian is a time-dependent quadratic form in the phonon creation and annihilation operators, whose coefficients are explicitly determined in terms of the time-dependent condensate wave function $\ensuremath{\xi}(\mathbf{x},t)$.The Bogoliubov method for the excitation spectrum of a Bose-condensed gas is generalized to apply to a gas with an exact large number $N$ of particles. This generalization yields a description of the Schr\"odinger picture field operators as the product of an annihilation operator $A$ for the total number of particles and the sum of a ``condensate wave function'' $\ensuremath{\xi}(\mathbf{x})$ and a phonon field operator $\ensuremath{\chi}(\mathbf{x})$ in the form $\ensuremath{\psi}(\mathbf{x})\ensuremath{\approx}A{\ensuremath{\xi}(\mathbf{x})+\ensuremath{\chi}(\mathbf{x})/\sqrt{N}}$ when the field operator acts on the $N$ particle subspace. It is then possible to expand the Hamiltonian in decreasing powers of $\sqrt{N}$, and thus obtain solutions for eigenvalues and eigenstates as an asymptotic expansion of the same kind. It is also possible to compute all matrix elements of field operators between states of different $N$. The excitation spectrum can be obtained by essentially the same method as Bogoliubov only if $\ensuremath{\xi}(\mathbf{x})$ is a solution of the time-independent Gross-Pitaevskii equation for $N$ particles and any chemical potential $\ensuremath{\mu}$, which yields a valid and stable solution of the Gross-Pitaevskii equation. The treatment within a subspace of fixed $N$ is identical in form to that usually used, but the interpretation of the operators is slightly different. A time-dependent generalization is then made, yielding an asymptotic expansion in decreasing powers of $\sqrt{N}$ for the equations of motion. In this expansion the condensate wave function has the time-dependent form $\ensuremath{\xi}(\mathbf{x},t)$, and the condition for the validity of the expansion is that $\ensuremath{\xi}(\mathbf{x},t)$ satisfies the time-dependent Gross-Pitaevskii equation $\ensuremath{\partial}\ensuremath{\xi}/\ensuremath{\partial}t=\ensuremath{-}({\ensuremath{\Elzxh}}^{2}/2m){\ensuremath{\nabla}}^{2}\ensuremath{\xi}+V\ensuremath{\xi}+Nu|\ensuremath{\xi}{|}^{2}\ensuremath{\xi}$. The physics is then described in a kind of interaction picture, called the condensate picture, in which the phonon operator can be expressed as $\ensuremath{\chi}(\mathbf{x},t)={\ensuremath{\sum}}_{k}{\ensuremath{\xi}}_{k}(\mathbf{x},t){\ensuremath{\alpha}}_{k},$ where the operators ${\ensuremath{\alpha}}_{k}$ are time-independent annihilation operators, and the state vector has a time evolution described by a Schr\"odinger equation in which the Hamiltonian is a time-dependent quadratic form in the phonon creation and annihilation operators, whose coefficients are explicitly determined in terms of the time-dependent condensate wave function $\ensuremath{\xi}(\mathbf{x},t)$.

arXiv:quant-ph/9703005Physical Review AAug 1, 1997256 citations

Exactly Solvable Many-Boson Model

R. W. Richardson

A many-boson model is formulated and expressions for its exact eigenstates and energies are obtained for both an arbitrary finite and an infinite number of bosons. The Hamiltonian of the model contains interactions between bosons whose momenta have equal magnitudes but opposite directions. The matrix elements of this interaction are taken to be a constant over a range of momenta surrounding k = 0. The ground state of the 2N-particle system is shown to be a product of N pair-creation operators acting on the vacuum state. Each of these pair-creation operators depends upon one of N parameters which are called pair energies. The N pair energies are shown to satisfy a coupled system of nonlinear algebraic equations. The energy of the state is the sum of the pair energies and the occupation probabilities of the single-particle levels are given as simple functions of the pair energies. Similar results are derived for the excited states of the system and for the states of an odd number of particles. These results are valid for both a repulsive and an attractive interaction, since they only depend upon the form of the interaction. The equations are solved algebraically for two model systems. The first of these is one whose single-particle kinetic energy takes on only one value. The equations for this system are solved for an arbitrary interaction strength and it is shown that the pair energies are proportional to the zeros of certain Laguerre polynomials. The second system is one in which the single-particle kinetic energy can take on two values. The equations for this system are solved in the strong repulsive-interaction limit and it is shown that the pair energies are proportional to the zeros of certain Jacobi polynomials. The excitation energies of this second system are shown to be proportional to 1/n and the occupations of the two single-particle levels in the ground state are shown to be proportional to n, where n is the total number of particles. For a repulsive interaction and an arbitrary single-particle spectrum, the algebraic equations for the pair energies are converted into an approximate integral equation for the density of roots which is accurate to order 1/n. This integral equation is solved for a strong interaction which, in the context of this model, means an interaction whose strength is greater than a constant times 1/V⅔ in the limit of a large volume. From this solution, the following results are obtained: (1) the lowest two single-particle levels have occupations of order n; (2) the excitation spectrum is that of a set of noninteracting quasiparticles; and (3) the quasiparticle spectrum has two zeros corresponding to the lowest two single-particle levels. Apart from the presence of two zeros, the quasiparticle spectrum does not differ significantly from that of the noninteracting particles.A many-boson model is formulated and expressions for its exact eigenstates and energies are obtained for both an arbitrary finite and an infinite number of bosons. The Hamiltonian of the model contains interactions between bosons whose momenta have equal magnitudes but opposite directions. The matrix elements of this interaction are taken to be a constant over a range of momenta surrounding k = 0. The ground state of the 2N-particle system is shown to be a product of N pair-creation operators acting on the vacuum state. Each of these pair-creation operators depends upon one of N parameters which are called pair energies. The N pair energies are shown to satisfy a coupled system of nonlinear algebraic equations. The energy of the state is the sum of the pair energies and the occupation probabilities of the single-particle levels are given as simple functions of the pair energies. Similar results are derived for the excited states of the system and for the states of an odd number of particles. These results are valid for both a repulsive and an attractive interaction, since they only depend upon the form of the interaction. The equations are solved algebraically for two model systems. The first of these is one whose single-particle kinetic energy takes on only one value. The equations for this system are solved for an arbitrary interaction strength and it is shown that the pair energies are proportional to the zeros of certain Laguerre polynomials. The second system is one in which the single-particle kinetic energy can take on two values. The equations for this system are solved in the strong repulsive-interaction limit and it is shown that the pair energies are proportional to the zeros of certain Jacobi polynomials. The excitation energies of this second system are shown to be proportional to 1/n and the occupations of the two single-particle levels in the ground state are shown to be proportional to n, where n is the total number of particles. For a repulsive interaction and an arbitrary single-particle spectrum, the algebraic equations for the pair energies are converted into an approximate integral equation for the density of roots which is accurate to order 1/n. This integral equation is solved for a strong interaction which, in the context of this model, means an interaction whose strength is greater than a constant times 1/V⅔ in the limit of a large volume. From this solution, the following results are obtained: (1) the lowest two single-particle levels have occupations of order n; (2) the excitation spectrum is that of a set of noninteracting quasiparticles; and (3) the quasiparticle spectrum has two zeros corresponding to the lowest two single-particle levels. Apart from the presence of two zeros, the quasiparticle spectrum does not differ significantly from that of the noninteracting particles.

Journal of Mathematical PhysicsSep 1, 196882 citations

Excitations in a Bose Gas at Finite Temperatures. II. Relation between Single-Particle and Density Fluctuations

Allan Griffin, T. H. Cheung

By making a systematic analysis of correlation functions in terms of irreducible and reducible parts, Ma and others have developed a dielectric formulation of the dynamics of a condensed Bose system which fully includes the background (or excited) atoms. We show that this formulation quite generally leads to the density and single-particle-correlation functions both having two resonances ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ and ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$ although with different weights. If the condensate ${n}_{0}=0$ , the ${\ensuremath{\omega}}_{2}$ mode only appears in the single-particle spectrum while ${\ensuremath{\omega}}_{1}$ only appears in the density-fluctuation spectrum. For finite values of ${n}_{0}$ , these modes are coupled and renormalized to ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ and ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$. As a specific illustration, we use the shielded-potential approximation (SPA) for the reducible self-energies. In this model, we find the free-particle excitations (${\ensuremath{\omega}}_{2}$) are coupled to the zero-sound density fluctuations (${\ensuremath{\omega}}_{1}$) through the action of the condensate. In the SPA, the reducible self-energies have a pole at ${\ensuremath{\omega}}_{1}$. The usual Bogoliubov, Hartree-Fock, and ordinary $t$ -matrix approximations involve nonsingular self-energies and hence do not exhibit the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ mode. Experimentally, it appears that the high-frequency ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$ mode has never been detected, but this is probably owing to the fact that it is strongly damped as a result of its coupling to the zero-sound mode. Using the two-fluid hydrodynamic equations, one can argue that at low frequencies the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ excitation corresponds to first sound and the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$ excitation corresponds to second sound. Finally, we briefly discuss the possibility of singularities in the irreducible self-energies at high frequencies and the resulting two-roton bound states.By making a systematic analysis of correlation functions in terms of irreducible and reducible parts, Ma and others have developed a dielectric formulation of the dynamics of a condensed Bose system which fully includes the background (or excited) atoms. We show that this formulation quite generally leads to the density and single-particle-correlation functions both having two resonances ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ and ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$ although with different weights. If the condensate ${n}_{0}=0$, the ${\ensuremath{\omega}}_{2}$ mode only appears in the single-particle spectrum while ${\ensuremath{\omega}}_{1}$ only appears in the density-fluctuation spectrum. For finite values of ${n}_{0}$, these modes are coupled and renormalized to ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ and ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$. As a specific illustration, we use the shielded-potential approximation (SPA) for the reducible self-energies. In this model, we find the free-particle excitations (${\ensuremath{\omega}}_{2}$) are coupled to the zero-sound density fluctuations (${\ensuremath{\omega}}_{1}$) through the action of the condensate. In the SPA, the reducible self-energies have a pole at ${\ensuremath{\omega}}_{1}$. The usual Bogoliubov, Hartree-Fock, and ordinary $t$-matrix approximations involve nonsingular self-energies and hence do not exhibit the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ mode. Experimentally, it appears that the high-frequency ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$ mode has never been detected, but this is probably owing to the fact that it is strongly damped as a result of its coupling to the zero-sound mode. Using the two-fluid hydrodynamic equations, one can argue that at low frequencies the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ excitation corresponds to first sound and the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$ excitation corresponds to second sound. Finally, we briefly discuss the possibility of singularities in the irreducible self-energies at high frequencies and the resulting two-roton bound states.

Physical Review AJun 1, 197336 citations

Collective Modes in Bose-Fermi Mixture Gases

Asano Yoji, Yoji Asano

This thesis theoretically deals with collective modes in the normal phase of a Bose-Fermi mixture composed of single-species bosons and single-species fermions.The mixture is experimentally established as an ultracold atomic gas.While the community studies the mixture including a superfluid phase, relevant to Bose-Einstein condensation of bosons and Cooper paring of fermions, dynamical behaviors in the normal phase have been not fully understood so that we are interested in.We consider the mixture in both spatially uniform and harmonically trapped systems with the Boltzmann equation at finite temperatures.We perform the analysis under the assumption that the physical quantities are near to equilibrium.After the linerization, we apply the moment method to the equation, obtaining a closed set of coupled equations to describe sound or oscillation.The collisional property of the system characterizes the collective modes in the sense of the hydrodynamic regime, collisionless regime, and transition between these regimes.In the uniform system, we investigate sound mode and the dynamic density structure factor.With decreasing temperature, we find that the sound velocity deviates from that of the hydrodynamic mode.Even at a temperature much lower than the Fermi temperature, the sound velocity does not reach the collisionless limit analogous to the zero sound in a Fermi system.However, we show the existence of a long-lived sound at a low temperature as distinct from the first and zero sounds.Nature of the long-lived sound seems to be originated by Pauli blocking of fermions; this statistical effect makes the interspecies collision suppress in the Bose-Fermi mixture.In the trapped system, we mainly study dipole oscillation.We consistently define the relaxation time originated by the interspecies scattering.We find that two oscillating modes in the collisionless regime have distinct fates in the hydrodynamic regime; one collisionless mode shows a crossover to a hydrodynamic in-phase mode, and the other collisionless mode shows a transition to two purely-damped modes.We show the time evolution of the three types of hydrodynamic modes with suitable initial conditions.We discuss temperature dependence of the dipole modes.We observe that the relaxation time dominates the dependence compared to the mean-field contribution.We also observe the nonmonotonic behavior of the relaxation time when Pauli blocking of fermions is not negligible.i List of Figures 1.1 Two schematic images of collisions with a period of oscillation . . . . . . . . . . . . . . . . 1 1.2 Measurement of the zero sound in pure liquid 3 He . . . . . . . . . . . . . . . . . . . . . . . 2 1.3 Transition between the first sound and the zero sound in a two-component Fermi gas . . . 3 1.4 Dipole oscillations of the two-component gas of 40 K . . . . . . . . . . . . . . . . . . . . . . 4 1.5 Frequencies of dipole mode in the two-component gas of 40 K . . . . . . . . . . . . . . . . . 4 1.6 Damping times of dipole mode in the two-component gas of 40 K . . . . . . . . . . . . . . . 5 1.7 Dipole oscillations of the mixture of 87 Rb and 40 K in the hydrodynamic regime . . . . . . .6 1.8 Dipole oscillations of the mixture of 87 Rb and 40 K in the deeper hydrodynamic regime . . .6 1.9 Dipole oscillations of the mixture of 87 Rb and 40 K in the collisionless regime . . . . . . . .7 2.1 Results obtained by solving the coupled moment equations for sound mode . . . . . . . . .26 3.1 Time evolution of the center-of-mass positions and velocity fields . . . . . . . . . . . . . .39 3.2 Dipole modes as a function of the relaxation time with N B = N F . . . . . . . . . . . . . . .40 3.3 Dipole modes as a function of the relaxation time with N B /N F = 0.2 . . . . . . . . . . . . .41 3.4 Ratio of the displacement of the center of mass to the velocity field as a function of the relaxation time . . . . . . . . . . . . . . . . . . .This thesis theoretically deals with collective modes in the normal phase of a Bose-Fermi mixture composed of single-species bosons and single-species fermions.The mixture is experimentally established as an ultracold atomic gas.While the community studies the mixture including a superfluid phase, relevant to Bose-Einstein condensation of bosons and Cooper paring of fermions, dynamical behaviors in the normal phase have been not fully understood so that we are interested in.We consider the mixture in both spatially uniform and harmonically trapped systems with the Boltzmann equation at finite temperatures.We perform the analysis under the assumption that the physical quantities are near to equilibrium.After the linerization, we apply the moment method to the equation, obtaining a closed set of coupled equations to describe sound or oscillation.The collisional property of the system characterizes the collective modes in the sense of the hydrodynamic regime, collisionless regime, and transition between these regimes.In the uniform system, we investigate sound mode and the dynamic density structure factor.With decreasing temperature, we find that the sound velocity deviates from that of the hydrodynamic mode.Even at a temperature much lower than the Fermi temperature, the sound velocity does not reach the collisionless limit analogous to the zero sound in a Fermi system.However, we show the existence of a long-lived sound at a low temperature as distinct from the first and zero sounds.Nature of the long-lived sound seems to be originated by Pauli blocking of fermions; this statistical effect makes the interspecies collision suppress in the Bose-Fermi mixture.In the trapped system, we mainly study dipole oscillation.We consistently define the relaxation time originated by the interspecies scattering.We find that two oscillating modes in the collisionless regime have distinct fates in the hydrodynamic regime; one collisionless mode shows a crossover to a hydrodynamic in-phase mode, and the other collisionless mode shows a transition to two purely-damped modes.We show the time evolution of the three types of hydrodynamic modes with suitable initial conditions.We discuss temperature dependence of the dipole modes.We observe that the relaxation time dominates the dependence compared to the mean-field contribution.We also observe the nonmonotonic behavior of the relaxation time when Pauli blocking of fermions is not negligible.i List of Figures 1.1 Two schematic images of collisions with a period of oscillation . . . . . . . . . . . . . . . . 1 1.2 Measurement of the zero sound in pure liquid 3 He . . . . . . . . . . . . . . . . . . . . . . . 2 1.3 Transition between the first sound and the zero sound in a two-component Fermi gas . . . 3 1.4 Dipole oscillations of the two-component gas of 40 K . . . . . . . . . . . . . . . . . . . . . . 4 1.5 Frequencies of dipole mode in the two-component gas of 40 K . . . . . . . . . . . . . . . . . 4 1.6 Damping times of dipole mode in the two-component gas of 40 K . . . . . . . . . . . . . . . 5 1.7 Dipole oscillations of the mixture of 87 Rb and 40 K in the hydrodynamic regime . . . . . . .6 1.8 Dipole oscillations of the mixture of 87 Rb and 40 K in the deeper hydrodynamic regime . . .6 1.9 Dipole oscillations of the mixture of 87 Rb and 40 K in the collisionless regime . . . . . . . .7 2.1 Results obtained by solving the coupled moment equations for sound mode . . . . . . . . .26 3.1 Time evolution of the center-of-mass positions and velocity fields . . . . . . . . . . . . . .39 3.2 Dipole modes as a function of the relaxation time with N B = N F . . . . . . . . . . . . . . .40 3.3 Dipole modes as a function of the relaxation time with N B /N F = 0.2 . . . . . . . . . . . . .41 3.4 Ratio of the displacement of the center of mass to the velocity field as a function of the relaxation time . . . . . . . . . . . . . . . . . . .

Mar 17, 20200 citations

Excitation spectrum of a system of interacting bosons

Moorad Alexanian, R. A. Brito

The Bogoliubov treatment of a weakly interacting Bose gas is extended by allowing macroscopic occupation of many single-particle quantum states. This generalization of having more than a single condensed mode gives rise to a nonuniform condensate and is being considered as a possible mechanism for the description of a dense strongly interacting Bose gas at low temperatures. The model with a repulsive $\ensuremath{\delta}$-function interparticle potential gives an energy spectrum of elementary excitations or quasiparticles with both phonon and roton features. The values obtained for the roton minimum energy $\ensuremath{\Delta}$ and the speed of (first) sound agree reasonably well with inelastic-neutron-scattering data.The Bogoliubov treatment of a weakly interacting Bose gas is extended by allowing macroscopic occupation of many single-particle quantum states. This generalization of having more than a single condensed mode gives rise to a nonuniform condensate and is being considered as a possible mechanism for the description of a dense strongly interacting Bose gas at low temperatures. The model with a repulsive $\ensuremath{\delta}$-function interparticle potential gives an energy spectrum of elementary excitations or quasiparticles with both phonon and roton features. The values obtained for the roton minimum energy $\ensuremath{\Delta}$ and the speed of (first) sound agree reasonably well with inelastic-neutron-scattering data.

Physical review. B, Condensed matterMay 1, 197811 citations

{{\ensuremath{\mu}}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}=\ensuremath{\Sigma}{\ensuremath{\gamma}}^{}\ensuremath{\Sigma}{p}^{}g(p,k,\ensuremath{\gamma}){{a}_{p+k}}^{\ifmmode\dagger\else\textdagger\fi{}}(\ensuremath{\gamma}){a}_{p}(\ensuremath{\gamma})

, where the

c

-number coefficients

g(p,k,\ensuremath{\gamma})

are determined by solving the equation of motion for

{{\ensuremath{\mu}}_{k}}^{\ifmmode\dagger\else\textdagger\fi{}}

in an approximate way. The sum over

\ensuremath{\gamma}

{\mathrm{He}}^{3}

{\mathrm{He}}^{4}

and the phonon velocity is determined as a function of

{\mathrm{He}}^{3}

concentration. The probability per unit time

P(\mathbf{k},{\ensuremath{\omega}}_{k})

that a particle transfers momentum k and energy

{\ensuremath{\omega}}_{k}

to the system by scattering is also calculated in terms of the dynamic form factor

{S}_{\ensuremath{\gamma}}(\mathbf{k},\ensuremath{\omega})

^{3}\text{H}\text{e}

and electrons in three dimensions. We find indeed a significant lowering of the zero-sound mode in

^{3}\text{H}\text{e}

and a broadening of the collective mode due to the coupling to two-particle--two-hole excitations in good agreement with experiments. The most visible effect in electronic systems is the appearance of a ``double-plasmon'' excitation.

The Bogoliubov method for the excitation spectrum of a Bose-condensed gas is generalized to apply to a gas with an exact large number

N

of particles. This generalization yields a description of the Schr\"odinger picture field operators as the product of an annihilation operator

A

for the total number of particles and the sum of a ``condensate wave function'' $\ensuremath{\xi}(\mathbf{x})$ and a phonon field operator $\ensuremath{\chi}(\mathbf{x})$ in the form $\ensuremath{\psi}(\mathbf{x})\ensuremath{\approx}A{\ensuremath{\xi}(\mathbf{x})+\ensuremath{\chi}(\mathbf{x})/\sqrt{N}}$ when the field operator acts on the

N

particle subspace. It is then possible to expand the Hamiltonian in decreasing powers of

\sqrt{N}

, and thus obtain solutions for eigenvalues and eigenstates as an asymptotic expansion of the same kind. It is also possible to compute all matrix elements of field operators between states of different

N

. The excitation spectrum can be obtained by essentially the same method as Bogoliubov only if $\ensuremath{\xi}(\mathbf{x})$ is a solution of the time-independent Gross-Pitaevskii equation for

N

particles and any chemical potential $\ensuremath{\mu}$, which yields a valid and stable solution of the Gross-Pitaevskii equation. The treatment within a subspace of fixed

N

is identical in form to that usually used, but the interpretation of the operators is slightly different. A time-dependent generalization is then made, yielding an asymptotic expansion in decreasing powers of

\sqrt{N}

for the equations of motion. In this expansion the condensate wave function has the time-dependent form $\ensuremath{\xi}(\mathbf{x},t)$, and the condition for the validity of the expansion is that $\ensuremath{\xi}(\mathbf{x},t)$ satisfies the time-dependent Gross-Pitaevskii equation $\ensuremath{\partial}\ensuremath{\xi}/\ensuremath{\partial}t=\ensuremath{-}({\ensuremath{\Elzxh}}^{2}/2m){\ensuremath{\nabla}}^{2}\ensuremath{\xi}+V\ensuremath{\xi}+Nu|\ensuremath{\xi}{|}^{2}\ensuremath{\xi}$. The physics is then described in a kind of interaction picture, called the condensate picture, in which the phonon operator can be expressed as $\ensuremath{\chi}(\mathbf{x},t)={\ensuremath{\sum}}_{k}{\ensuremath{\xi}}_{k}(\mathbf{x},t){\ensuremath{\alpha}}_{k},$ where the operators ${\ensuremath{\alpha}}_{k}$ are time-independent annihilation operators, and the state vector has a time evolution described by a Schr\"odinger equation in which the Hamiltonian is a time-dependent quadratic form in the phonon creation and annihilation operators, whose coefficients are explicitly determined in terms of the time-dependent condensate wave function $\ensuremath{\xi}(\mathbf{x},t)$.

{n}_{0}=0

, the ${\ensuremath{\omega}}_{2}$ mode only appears in the single-particle spectrum while ${\ensuremath{\omega}}_{1}$ only appears in the density-fluctuation spectrum. For finite values of

{n}_{0}

, these modes are coupled and renormalized to ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ and ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$. As a specific illustration, we use the shielded-potential approximation (SPA) for the reducible self-energies. In this model, we find the free-particle excitations (${\ensuremath{\omega}}_{2}$) are coupled to the zero-sound density fluctuations (${\ensuremath{\omega}}_{1}$) through the action of the condensate. In the SPA, the reducible self-energies have a pole at ${\ensuremath{\omega}}_{1}$. The usual Bogoliubov, Hartree-Fock, and ordinary

t

-matrix approximations involve nonsingular self-energies and hence do not exhibit the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ mode. Experimentally, it appears that the high-frequency ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$ mode has never been detected, but this is probably owing to the fact that it is strongly damped as a result of its coupling to the zero-sound mode. Using the two-fluid hydrodynamic equations, one can argue that at low frequencies the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{1}$ excitation corresponds to first sound and the ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\omega}}}_{2}$ excitation corresponds to second sound. Finally, we briefly discuss the possibility of singularities in the irreducible self-energies at high frequencies and the resulting two-roton bound states.