On planar quantum dynamics of a magnetic dipole moment in the presence of electric and magnetic fields

The planar quantum dynamics of a neutral particle with a magnetic dipole moment in the presence of electric and magnetic fields is considered. The criteria to establish the planar dynamics reveal that the resulting nonrelativistic Hamiltonian has a simplified expression without making approximations, and some terms have crucial importance for system dynamics.


Introduction
The Aharonov-Bohm effect [1] has been an usual framework for demonstrating the importance of potentials in quantum mechanics. After its experimental verification [2,3] several other analogs effects were being proposed along the last decades. For example, in Ref. [4], it was shown that a particle with a magnetic moment moving in an electric field acquires a quantum phase. This phase has been observed in a neutron interferometer [5] and in a neutral atomic Ramsey interferometer [6]. In Refs. [7,8] it was verified that a neutral particle with an electric dipole moment which moves in a magnetic field acquires a topological phase. The experimental confirmation of this phase was established in Ref. [9]. In Ref. [10] was proposed a unified and fully relativistic treatment of the interaction of the electric and magnetic dipole moments of a particle with the electromagnetic field. The essence of this study reveals that new force on dipoles are obtained using the non-Abelian nature of this interaction, and new experiments analogous to the Aharonov-Bohm effect to test this interaction are proposed. Since this interaction is a consequence of a nonminimal coupling is also interesting to analyze the consequences of this interaction in other contexts. For instance, it may be of interest to study scattering and bound states of neutral fermions in external electromagnetic fields and access other physical quantities a e-mail: edilbertoo@gmail.com such as energy bound states [11,12] and scattering [13,14]. The dual results of the magnetic dipole interaction for the electric dipole interaction for the system considered in Ref. [10] has been established in Ref. [15], where the phase shift in the interference of a magnetic or electric dipole due to the electromagnetic field is obtained relativistically and nonrelativistically.
In this work, we consider the same system addressed in Ref. [15] but now assuming that the dynamics is purely planar, and derive its equation of motion. This system is a generalization, for example, from that studied in Ref. [12], where only the effects of an electric field has been considered. An interesting feature of this system is that even in the absence of electric field it admits bound and scattering states, which does not occur in the previous work. As an application, we consider the problem of bound state for the case of a magnetic dipole moment interacting with electric and magnetic fields generated by an infinitely long charged solenoid, carrying a magnetic field. In our treatment, we consider the self-adjoint extension method [16], which is appropriate to address any system endowed with a singular Hamiltonian (due to localized fields sources or quantum confinement) [17][18][19][20][21][22][23][24][25][26]. We determine the energy spectrum and wave functions by applying boundary conditions allowed by the system.

The planar Pauli equation
We begin with the Dirac equation in (3 + 1) dimensions [15] which governs the nonrelativistic dynamics of a neutral particle that possesses magnetic dipole moment, in the presence of electric and magnetic fields (h = c = 1) where µ is the magnetic dipole moment, F µν is the electromagnetic tensor whose components are given by F 0i , F i j = −E i , −ε i jk B k , and σ 0 j , σ i j = iα j , −ε i jk Σ k , where Σ k is the spin operator, are the components of the operator σ µν = i[γ µ , γ ν ]/2, which are given in terms of the Dirac matrices. With this notation, it is possible to show that the spin is coupled to the electromagnetic field tensor through the term where E and B are the electric and magnetic field strengths. This result is explicitly calculated in the following representation of the γ-matrices: with σ σ σ = (σ 1 , σ 1 , σ 3 ) being the Pauli matrices. Equation (1) can be written aŝ where the operator is the Dirac Hamiltonian. The non-relativistic limit of Eq.
(3) was established in Ref. [15], and the relevant equation is found to be: where ψ is a two-component spinor, witĥ where µ µ µ = µσ σ σ . Our goal is to analyze the physical implications of the Hamiltonian (6), when we assume that the system dynamics is now planar. This is established as follows. By detaching the third component of Eq. (6), we get If we assume that dynamic is planar, the above Hamiltonian provides an important result, namely, the [p 3 − (µ µ µ × E) 3 ] 2 term leads exactly to the quantity µ 2 E 2 /2M. The planar case is accessed by requiring that p z = z = 0 together with the imposition of the fields should not have the third direction. This question can also be understood when we look at the symmetry under z translations, which allows us to access the solutions of the planar Dirac equation. This type of simplification is in fact manifested only when we assume that the particle moves in the plane. Thus, since A µ = (Φ, A), we write the electric and magnetic fields, respectively, as Fields E and B above are now intrinsically two-dimensional. Note that the square of Eq. (8) gives exactly E 2 1 + E 2 2 for the planar case. Also, the restriction imposed on the potential A reveals that the (µ µ µ · B) i term in Eq. (7) is now identically zero. Now, we can show that the quantity (µ µ µ × E) 3 is given by µ (σ 1 E 2 − σ 2 E 1 ), and the third term of Eq. (7) results Thus, we now can write Eq. (7) aŝ where the magnetic interaction term (µ µ µ · B) 3 gives the only explicit dependence of the spin. In Ref. [15] it was assumed that the charge density ρ = ∇ ∇ ∇·E and also the µ 2 E 2 /2M term (for thermal neutrons) are negligible. In fact, this approximation can only be performed, if we are only interested in the study of phase shift. However, when we want to study the dynamics of the system, such as the scattering and bound states problems, all terms of the equation motion of must be taken into account. This, for example, has been addressed by Hagen [27] to show that there is an exact equivalence between the AB and AC effects for spin-1/2 particles. In these effects, the ∇ ∇ ∇ · E and ∇ ∇ ∇ × A being proportional to a delta function, such terms must now contribute to the dynamics of the system, and can not be neglected. For this reason, since we are dealing with Aharonov-Bohm-like system for spin-1/2 particles, the µ (∇ ∇ ∇ · E) i and (µ µ µ · B) 3 terms in Eq. (11), may not be negligible. Let's clarify this issue. Consider an infinitely long solenoid, carrying a magnetic field B, and with a charge density λ distributed uniformly about it along the z-axis. The electric field and magnetic flux tube (in cylindrical coordinates) generated by this configuration are known to be where φ is the magnetic flux inside the tube, and the vector potential in the Coulomb gauge is We see that the fields E and B are proportional to a δ function. By using Eqs. (12) and (13), Pauli equation (5) is now written aŝ witĥ where φ E = 2µλ and φ B = 2Mµφ . From Eq. (15), we can see that ψ is an eigenfunction of σ z , whose eigenvalues designate by s = ±1, that is, σ z ψ = ±ψ = sψ. Thus, sinceĤ commutes with the operatorsĴ z = −i∂ ϕ + σ z /2, whereĴ z is the total angular momentum operator in the z-direction, we seek solutions of the form with m + 1/2 = ±1/2, ±3/2, . . . , (m ∈ Z). Inserting (17) into Eq. (15), we can extract the radial equation for f m (r) where and Note that, even in the absence of electric field, bound and scattering states are possible. This does not occur, for example, in the system studied in Ref. [12], where the particle interacts only with an electric field. Moreover, if a magnetic field is present, the physical system changes completely. We will see later that this fact directly influences on the expression for the self-adjoint extension parameter and, hence, on the boundary conditions allowed by the operator H 0 . In other words, this has direct implications in the dynamics of the system. This can be seen more easily by studying the signal of φ E −sφ B in Eq. (19), where several possible combinations of φ E , φ B and s gives us the possibilities for the existence of bound and scattering states. As a result of these combinations, we have φ E − sφ B < 0, scattering and bound states, φ E − sφ B > 0, scattering states.
The case φ E = sφ B is not of interest here because it cancels the term that explicitly depends on the spin.

Physical regularization and the bound states problem
In this section, we study the dynamics of the system in all space, including the r = 0 region. We consider the problem of bound states. To this end, we use the self-adjoint extension method in the treatment. As is well known, if the Hamiltonian has a singularity point, as is the case of the Hamiltonian in Eq. (19), we must verify that it is self-adjoint in the region of interest. Even though H † 0 = H 0 , their domains could be different. This is the crucial point in our study. The operator H 0 , with domain D(H 0 ), is self-adjoint if D(H † 0 ) = D(H 0 ) and H † 0 = H 0 . However, for this to be established, we must find the deficiency subspaces, with dimensions n + and n − , respectively, called deficiency indices of H 0 [16]. We also know of this theory that a necessary and sufficient condition for H 0 being essentially selfadjoint is that its deficiency indices n + = n − = 0. On the other hand, if n + = n − ≥ 1 the operator H 0 has an infinite number of self-adjoint extensions parametrized by the unitary operators U : N + → N − . With these ideas in mind, we now decompose the Hilbert space H = L 2 (R 2 ) with respect to the angular momentum H = H r ⊗ H ϕ , where H r = L 2 (R + , rdr) and H ϕ = L 2 (S 1 , dϕ), with S 1 denoting the unit sphere in R 2 . The operator −∂ 2 ϕ is known to be essentially self-adjoint in L 2 (S 1 , dϕ). By using the unitary operator [18] V : L 2 (R + , rdr) → L 2 (R + , dr), given by the operator H 0 reads which is essentially self-adjoint for (m − sφ E ) ≥ 1, while for (m − sφ E ) < 1, it admits a one-parameter family of selfadjoint extensions [16]. To characterize this family, we follow the recipe based in boundary conditions given in Ref. [28]. Basically, the boundary condition is a match of the logarithmic derivatives of the zero-energy solutions for Eq. (18) and the solutions for the problem H 0 plus self-adjoint extension. Then, following [28], we temporarily forget the δfunction potential and find the boundary conditions allowed for H 0 . Next, we substitute the problem in Eq. (18) by plus self-adjoint extensions. Here, f ζ is labeled by the parameter ζ of the self-adjoint extension, which is related to the behavior of the wave function at the origin. In order for the H 0 to be a self-adjoint operator in H r , its domain of definition has to be extended by the deficiency subspace, which is spanned by the solutions of the eigenvalue equation where k 2 0 ∈ R is introduced for dimensional reasons. Since H † 0 = H 0 , the only square integrable functions which are solutions of Eq. (29) are the modified Bessel functions of second kind with Im √ ±i > 0. By studying Eq. (30), we verify that it is square integrable only in the range m − sφ E ∈ (−1, 1). In this interval, nevertheless, the Hamiltonian (20) is not self-adjoint. The dimension of such deficiency subspace is (n + , n − ) = (1, 1). So, we have two situations for m − sφ E , i.e., To address both cases of Eq. (31), we write Eq. (30) as Equation (32) allows us to identify the domain of H † 0 as So, to extend the domain D(H 0 ) and make it equal to D(H † 0 ) and therefore to turn H 0 self-adjoint, we get Equation (34) establishes the following result. For each value of ζ , we have a possible domain for D(H 0,ζ ), but are the physical parameters of the problem that will select a particular value of it. The Hilbert space is now specified by [16] f ζ (r) = f m (r) where f m (r), with f m (0) =ḟ m (0) = 0 (ḟ ≡ d f /dr), is the regular wave function and the parameter ζ ∈ [0, 2π) represents a choice for the boundary condition. For each ζ , we have a possible domain for H 0 and the physical situation is the factor that will determine the value of ζ [20,[29][30][31][32][33][34]. Thus, to find a fitting for ζ compatible with the physical situation, a physically motivated form for the magnetic field is preferable for the regularization of the δ -function. This is accomplished by replacing (14) with [35] eA = With this modification, the delta function in Eq. (19) is now regularized as δ (r − a)/a. A remarkable feature of this regularization is that, although the functional structure of δ (r)/r and δ (r − a)/a are quite different, we are free to use any form of potential once that the specific details of the regularization model can be shown to be irrelevant provided that only the contribution is independent of angle and has no δ -function contribution at the origin [35]. It should also be mentioned that the δ (r − a)/a potential is one dimensional and well defined contrary to the two dimensional δ (r)/r. Now, we are in the position to determine a fitting value for ζ . To do so, we consider the zero-energy solutions for f 0 (r) with the regularization, and for f ζ ,0 (r) without the δ function, respectively, i.e., The value of ζ is determined by the boundary condition By integrating Eq. (37) from 0 to a and noting that the behavior of f 0 as a → 0 is f 0 ∼ r |m−sφ E | , the left-hand side of Eq. (39) is found to be To calculate the right-hand side of Eq. (39), we need to use the asymptotic behavior for K ν (z) in the limit z → 0, given by The substitution of Eq. (41) into Eq. (35) leads to with and Γ (±) = Γ (1 ± |m − sφ E |) was defined for purposes of simplification. Inserting (40) and (42) in (39), we obtain which gives us the parameter ζ in terms of the physics of the problem, i.e., the correct behavior of the wave functions when r → 0. As promised above, let us now derive determine the bound states for H 0 . In order systems have a bound state, its energy must be negative, so that in Eq. (28), k is a pure imaginary, i.e., k = iκ, with κ = √ −2ME, where E < 0 is the bound state energy. Then, with the substitution k → iκ, we have The above equation is the modified Bessel equation whose general solution is given by Since these solutions belong to D(H ζ ,0 ), it is of the form (35) for some ζ selected from the physics of the problem.
So, we substitute (46) into (35) and use (41) to calculate the left-hand side of Eq. (39). After these manipulations, we find the relation Solving the above equation for E, we find the energy spectrum Notice that there is no arbitrary parameter in the above equation. Also, to ensure that the energy is a real quantity, we must establish that This inequality is satisfied if Because of the condition that |m − sφ E | < 1, it is sufficient to consider |φ E − sφ B | ≥ 1. A necessary condition for a δ function to generate an attractive potential, which is able to support bound states, is that the coupling constant (φ E − sφ B ) must be negative. Thus, the existence of bound states requires So, it seems that we must have in such way that the flux and the spin must be parallel, and consequently, a minimum value for |φ B | and |φ E | is established.

Conclusions
We have analyzed the planar quantum dynamics of a magnetic dipole moment in the presence of electric and magnetic fields. We have shown that the initial Hamiltonian system (Eq. (6)) reduces to a planar form (Eq. (11)) without making any approximations. As an application, we have considered the bound state problem for the case of a magnetic dipole moment interacting with electric and magnetic fields generated by an infinitely long solenoid, carrying a magnetic field, and with a charge density distributed uniformly about it along the z-axis. The self-adjoint extension approach was used to determine the bound states of the particle in terms of the physics of the problem, in a very consistent way and without any arbitrary parameter.