An analytical expression for the lower bound of the average number of hyperfine component
μ=0 particles in the ground state of spin-1 condensates (denoted as
ρ0__) under a magnetic field has been derived. In the derivation the total magnetization
M is kept rigorously conserved. Numerical examples are given to show the applicability of the analytical expression. It was found that, in a broad domain of parameters specified in the paper, the lower bound is very close to the actual
ρ0__. Thereby, in this domain, $\overset{\_\_}{% ρ_{0}}$ can be directly evaluated simply by using the analytical expression.
The variation of the hyperfine populations $N{\ensuremath{\rho}}_{\ensuremath{\mu}}$ of the ground state (g.s.) of spin-2 condensates against a magnetic field
B is studied. A
N-body theory which is rigorous in dealing with the spin degrees of freedom has been adopted. The conservation of the total magnetization
M is assumed. A combined hyperfine population $N{\ensuremath{\rho}}_{\text{comb}}=N({\ensuremath{\rho}}_{0}\ensuremath{-}3{\ensuremath{\rho}}_{2}\ensuremath{-}3{\ensuremath{\rho}}_{\ensuremath{-}2})$, which is an observable, is defined. The lower bound for this population has been derived analytically. It turns out that the lower bound is identical to the actual population when
B=0 and $\ensuremath{\infty}$, and the population is monotonically uprising in between. Thus, via the analytical form of the lower bound, the variation of the population against
B can be roughly known in advance. Numerical examples are given to demonstrate the applicability of the bound. When
M=0, the g.s. in the polar phase was found to be highly unstable against
B because the neighboring level density is extremely high. When
M=N, the g.s. in the cyclic phase was also found to be highly unstable. This is because the appearance of
B, even if it is very weak, will break the high degeneracy of the g.s. and result in having numerous levels emerging in the neighborhood of the g.s.
Based on perturbation theory up to second order, analytical asymptotic expressions for the variation of the population of hyperfine component $\ensuremath{\mu}=0$ particles in the ground state of spin-1 condensates against a magnetic field
B have been derived. The ranges of
B in which the asymptotic expressions are applicable have been clarified via a comparison of the numerical results from the analytical expressions and from a diagonalization of the Hamiltonian in a complete spin space. It was found that, for
87Rb, the two analytical expressions, one for a weak field and the other one for a strong field, together cover the whole range of
B from zero to infinity. For Na, the analytical expressions are valid only if
B is very weak or sufficiently strong.
The Feynman inequality
EG\ensuremath{\le}
Etrial+${\mathrm{lim}}_{\mathrm{\ensuremath{\beta}}\ensuremath{\rightarrow}\mathrm{\ensuremath{\infty}}}$〈S-
Strial〉/\ensuremath{\beta} for path integrals provides a powerful upper bound on the ground-state energy
EG of a large variety of systems.
Etrial is the ground-state energy of some trial system with action
Strial for imaginary values of the time variable, and S is the action (also expressed in imaginary time variables) of the system under study. \ensuremath{\beta}=1/
kBT, where
kB is the Boltzmann constant and T the temperature. However, the Feynman inequality is not a priori justified for a system in a magnetic field, because imaginary terms subsist in the action also after transforming to imaginary time variables. Replacing or extending this inequality when magnetic fields are present has therefore been a long-standing problem. In the present paper we solve this problem. We first derive an inequality, providing an upper bound for the ground-state energy, that is valid even in the case of a nonzero magnetic field,
EG\ensuremath{\le}
Etrial+〈\ensuremath{\infty}\ensuremath{\Vert}scrT{
markUtrial (\ensuremath{\infty},-\ensuremath{\infty})[V(0)-
Vtrial(0)]}-\ensuremath{\infty}〉,for a system with Hamiltonian
H0+V. T time-ordering operator, and
Utrial is the time evolution operator of a trial system with Hamiltonian
H0+
Vtrial in the interaction representation, with the interactions V(t) and
Vtrial(t) switched on adiabatically. Because of the time ordering, retardation effects are also properly taken into account. The contribution of the magnetic field is included in the unperturbed Hamiltonian
H0. If the time-dependent integrands occurring in the matrix element in the right-hand side of our generalized inequality satisfy certain analyticity conditions in the complex-time plane, this inequality reduces to the Feynman inequality for path integrals. If these analyticity conditions are not satisfied, our generalized inequality may introduce supplementary terms
EDB in the right-hand side of the Feynman upper bound,
EG\ensuremath{\le}
Etial+${\mathrm{lim}}_{\mathrm{\ensuremath{\beta}}\ensuremath{\rightarrow}\mathrm{\ensuremath{\infty}}}$〈S-
Strial〉/\ensuremath{\beta} +
EDB,because different branch lines or singularities have to be taken into account in the transformation to imaginary time variables. As an important illustration, our generalized inequality is applied to the problem of the Fr\"ohlich polaron in a magnetic field. From the generalization of the Feynman an inequality derived in the present paper, we determine the conditions to be imposed on the variational parameters in the trial action, such that the Feynman upper bound in its original form remains valid for a polaron in a magnetic field. Some limiting cases are studied analytically to illustrate the relevance of our additional constraints on the variational parameters of the trial system. In the free-particle limit and for a particular value of one of the variational parameters, we explicitly derive the contributions from the branch lines in the complex-time plane which arise if these additional constraints are not satisfied.