Bound States at Threshold Resulting from Coulomb Repulsion
The eigenvalue absorption for a many–particle Hamiltonian depending on a parameter is analyzed in the framework of non–relativistic quantum mechanics. The long–range part of pair potentials is assumed to be pure Coulomb and no restriction on the particle statistics is imposed. It is proved that if the lowest dissociation threshold corresponds to the decay into two likewise non–zero charged clusters then the bound state, which approaches the threshold, does not spread and eventually becomes the bound state at threshold. The obtained results have applications in atomic and nuclear physics. In particular, we prove that an atomic ion with the critical charge and electrons has a bound state at threshold given that , whereby the electrons are treated as fermions and the mass of the nucleus is finite.
On leave from: ] Institute of Physics, St. Petersburg State University, Ulyanovskaya 1, 198504 Russia
In Refs. 1, ; 2, it was proved that a critically bound N–body system, where none of the subsystems has bound states with and particle pairs have no zero energy resonances, has a square integrable state at zero energy. The condition on the absence of 2–body zero energy resonances was shown to be essential in the three–body case 1 (). Here we consider the –particle system, where particles can be charged and apart from short–range pair–interactions may also interact via Coulomb attraction/repulsion. The formation of bound states at threshold in the two–particle case when the particles Coulomb repel each other is well–studied 3 (); gest (). In the three–particle case there is a well–known proof ostenhof () that a two–electron ion with an infinitely heavy nucleus has a bound state at threshold, when the nuclear charge becomes critical.
Our aim here is to investigate the general many–particle case. Here we generalize the result in Ref. ostenhof, to the case of many electron ions with Fermi statistics and finite nuclear mass. In the proofs we shall use the bounds on Green’s functions from Ref. 3, as well as the technique of spreading sequences from Ref. 1, , that is we prove the eigenvalue absorption by demonstrating that the wave functions corresponding to bound states do not spread, c.f. Theorem 1 in Ref. 1, . A different approach based on the calculus of variations was recently developed in Ref. frank, , where the authors give an alternative proof to the result in Ref. ostenhof, . The authors in Ref. frank, indicate that their approach could be generalized to the many–particle case. In the present paper as well as in Refs. ostenhof, ; frank, one uses essentially the same idea, namely, one uses the fact that the weak limit of ground state wave functions is a solution to the Schrödinger equation at the threshold. The hardest part is to prove that the weak limit is not identically zero. Our approach differs from the ones in Refs. ostenhof, ; frank, in that we use the upper bounds on the two–particle Green’s functions3 ().
The paper is organized as follows. In Sec. II we introduce notations, formulate the main theorem and prove a number of technical lemmas. In Sec. III we derive an upper bound on the Green’s function, which is used in Sec. IV for the proof of Theorem 1. In Sec. V we discuss two main applications of Theorem 1 concerning the stability diagram of three Coulomb charges (Theorem 2 in Sec. V.1) and negative atomic ions (Theorem 3 in Sec. V.2). In Appendix A we derive various criteria for non–spreading sequences.
Let us mention physical applications. The effect when a size of a bound system increases near the threshold and by far exceeds the scales set by attractive parts of potentials was discovered in neutron halos, helium dimer, Efimov states, for discussion see Refs. fedorov, ; zhukov, ; efimov, ; hansen0, . Here we demonstrate that in a many–particle system similarly to the two–body case 3 (); gest () a Coulomb repulsion between possible decay products blocks the spreading of bound states and forces an bound state at threshold. In nuclear physics, this, in particular, explains why contrary to neutron halos no proton halos are found hansen ().
Ii Formulation of the Main Theorem
We consider the –particle Hamiltonian ()
where is a parameter, is the kinetic energy operator with the center of mass removed, are particles’ position vectors and denote the particles’ charges depending on . We shall assume that for each given . Here denotes the space of bounded Borel functions vanishing at infinity. We shall also take particle spins into account, though we shall consider only spin–independent Hamiltonians. The Hamiltonian acts in , where the direct sum has summands and denotes the spin of particle . Similar notation for the Hilbert space can be found in Refs. quotesimon, ; quotethaller, . By Kato’s theorem reed (); teschl () is self–adjoint on , where and denotes the corresponding Sobolev space teschl (); liebloss (). A function depends explicitly on the arguments as , where and are the spin variables.
We treat the particles with integer spins as bosons and particles with half–integer spin as fermions. denotes the orthogonal projection operator on the subspace of functions, which are symmetric with respect to the interchange of bosons and antisymmetric with respect to the interchange of fermions. We denote the bottom of the continuous spectrum by
We shall use the function , which determines the asymptotic behavior at infinity
where and always denotes the characteristic function of the set . Note that is continuous and . We make the following assumptions
is defined for an infinite sequence of parameter values and , where . For all there is such that , where , and . Besides, .
and , where is such that and , are fixed constants. Additionally, for all .
Let label all the distinct ways ims () of partitioning particles into two non–empty clusters and . We define the Jacobi intercluster coordinates for the clusters as and respectively, where and (the symbol denotes the number of particles in the corresponding cluster). By we denote the full set of intercluster coordinates and we set
points from the center of mass of to the center of mass of . The full set of Jacobi coordinates is .
We denote the sum of interaction cross terms between the clusters by
The product of net charges of the clusters is defined as
The projection operators on the proper symmetry subspace for the particles within clusters and are and respectively. Namely, projects on a subspace of functions, which are antisymmetric with respect to the interchange of fermions in and symmetric with respect to the interchange of bosons in (). Naturally, and . We also define . The Hamiltonian (1) can be decomposed in the following way
where is the Hamiltonian of the clusters’ intrinsic motion and denotes the reduced mass derived from clusters’ total masses. From now on without loss of generality we set .
It is convenient to treat the Hilbert space as the tensor product , where the first term in the product corresponds to the space associated with coordinates and spin variables, while the second one refers to the space associated with the coordinate. In such case the operator has the form , where is the restriction of to . The coordinate is unaffected by permutations of particles within the clusters or . Therefore, , where denotes the restriction of to the space associated with coordinates and spin variables.
The set of assumptions is continued as follows.
For and one has . There is such that the following inequalities hold for
The requirement R3 says that the bottom of the continuous spectrum of is set by the decomposition into those two clusters that correspond to any of the decompositions . Inequality (9) introduces a gap between the ground state energy of the two clusters and other excited states. For and we define the projection operator acting on
where are spectral projections of . Note that by R3 the projection operators , have a finite dimensional range.
The last assumption introduces the uniform control over the fall off of clusters’ wave functions
There are constants such that
for and .
Due to R3 there must exist orthonormal for such that
Suppose that the orthonormal set of function is such that and for , where are constants. If then
where is the Lieb’s constant in the Cwikel–Lieb–Rosenbljum bound.
From it follows that
where we set . Hence,
By the min–max principle does not exceed the number of negative energy bound states of the operator in square brackets in (16). This number, in turn, is equal to the number of negative energy bound states of the operator in square brackets considered in times due to the spin degeneracy. Now (14) follows from the Cwikel–Lieb–Rosebljum bound reed (); cwikel (). ∎
Now we can formulate the main theorem.
Suppose that satisfies and
Then the sequence does not spread and there exists such that , where and .
Let us remark that the term spreading was defined in Ref. 1, for sequences in . We shall say that a sequence spreads if spreads for all possible fixed values of the spin variables. We postpone the proof of Theorem 1 to Sec. IV.
Together with the upper bound on the Green’s function derived in the next section the following lemma is the key ingredient in the proof of Theorem 1.
There is independent of such that
for defined in R1, defined in R2 and defined in R4.
The statement of the lemma is based on the following inequality, which can be checked directly. For all
For fixed the term on the lhs of (19) falls off like . We write
For any cluster decomposition and
where , are numerical coefficients depending on masses. It is easy to see that the coefficient in front of is always by fixing and taking . Therefore, by (19) we have
Using R2 it is easy to see that . Because we have . ∎
Iii Upper Bound on the Two Particle Green’s Function
Consider the following integral operator on
for , whose integral kernel we denote as (the superscript “c” refers to “Coulomb”). Note that away from if either or , c. f. Lemma 1 in Ref. 3, . The following Lemma uses the upper bound on a two particle Green’s function from Ref. 3, .
There is such that for all ,
where the norm on the lhs is the operator norm.
where the last inequality follows from being an integral operator with a positive kernel 3 (); lpestim (). We shall derive the following estimates and for , from which the the statement of the Lemma follows. The first term on the rhs of (27) is the norm of the self–adjoint operator, which can be estimated as follows
where we have used the inequality for non–negative self–adjoint operators and (see, for example Ref. glimmjaffe, , Proposition A.2.5 on page 131). Thus as claimed.
Lets us set equal to the positive root of the equation . Then we get
where have to be chosen to satisfy the following inequalities
From the inequality (31) we obtain the bound
The integral in (36) can be calculated explicitly and we obtain as claimed. ∎
We shall need the following corollary of Lemma 3
For fixed the following inequality holds
For an arbitrary we have
Now using that and applying the Cauchy-Schwartz inequality we get from Eq. (39)
For the series on the rhs of Eq. (40) obviously converge. ∎
Iv Proof of the Main Theorem
We shall need an analogue of the IMS localization formula, see Ref. ims, . The functions form the partition of unity and are homogeneous of degree zero in the exterior of the unit sphere, i.e. for , (this makes fall off at infinity). Additionally, there exists a constant such that
The functions of the IMS decomposition can be chosen sigalsays1 (); sigalsays2 () invariant under permutations of particle coordinates both in and in , hence . Note also that for all one has , c. f. Lemma 7.4 in Ref. liebloss, (the proof in Ref. liebloss, easily extends to the case of Sobolev spaces of higher order). The following version of the IMS localization formula can be verified by the direct substitution
where is the Laplace on . Rescaling (42) we get
where are real coefficients depending on masses and the second term is relatively compact ims (). We introduce
From (43) it follows that
where we define
The Hamiltonian defined for contains intercluster interactions of the following four clusters , , , , while all interaction cross–terms between these four clusters are contained in . (For some partitions it might happen that one of the four clusters is empty). If we define by the projection operator on the proper symmetry subspace for particles within the cluster then by the HVZ theorem reed (); teschl (); beattie ()
where we define
Note that .
Suppose that satisfies . If then