Mapping Between Charge-Dyon and Position-Dependent Mass Systems∗

Anderson L.de Jesus and Alexandre G.M.Schmidt

Departamento de F´ısica de Volta Redonda,Instituto de Ciˆencias Exatas — Universidade Federal Fluminense,R.Des.Ellis Hermydio Figueira,783,Bloco C,Volta Redonda RJ,CEP 27213-145,Brazil

Abstract The primary purpose of this work is to reproduce the scenario composed of a charge-dyon system utilizing position-dependent effective mass (PDM) background in the non-relativistic and in the relativistic regimes.In the nonrelativistic case we substitute the exact charge-dyon eigenfunction into PDM Schrodinger equation,in the Zhu-Kroemer parametrization,and then solve it for the mass distribution considering M =M(r).Analogously,in the relativistic case we study the Klein-Gordon equation for a position-dependent mass,and in this case,we are able to analytically solve the equation for M =M(r,θ).

Key words: position-dependent mass,Schrodinger equation,Klein-Gordon equation,dyon

1 Introduction

Dyon is a hypothetical particle having electric and magnetic charges proposed by Schwinger in 1969 as a phenomenological alternative to quark model.[1]In his article,Schwinger speculates that hadrons are composed of these dual charged particles,having fractional electric and magnetic charges and following the idea of electroweak interaction with its vectorial bosons exchanging electric charge,he postulated the existence of a vectorial bosonS(strong)of unit magnetic charge intermediating chargeexchange process for the dyon.The study of dyons and isolated magnetic monopoles has received attention over the latest years,in part by the possible explanation of the electric charge quantization by Dirac.[2]In the experimental aspect there is also plenty of research and the physics of magnetic monopoles was reproduced in spin ice system,[3−5]as well as,more recently,in quantum fields,[6]synthetic magnetic fields,[7]creation of a Dirac monopoleantimonopole pair in a spin-1 Bose-Einstein condensate[8]and by the use of metamaterials.[9]

On the other hand,quantum systems with positiondependent effective mass (PDM) have been studied and used in a great number of works and a considerable interest in this subject has grown up over the latest years.The quantum systems with position-dependent effective mass arose initially in the study of transport phenomena in semiconductors of variable,position-dependent chemical composition and the ordering problem was studied by von Roos.[10]PDM systems can also be used to model scattering in abrupt heterostructures,[11]quantum dots[12]and,surprisingly,Dirac magnetic monopoles.The investigation about mapping charge-monopole into PDM systems was initially carried out by the authors in Ref.[13]where an exact mapping was constructed between these systems.In the present work,our goal is to make a similar map between charge-dyon and a PDM system,in other words,we study the possibility of having a PDM charge with the same eigenfunctions and energy eigenvalues of the constant mass charge-dyon system.Thus,we look for a PDM system which could,in principle,be realized in the laboratory and that can simulate another one —the charge dyon system — that in the present day cannot be studied experimentally.The idea of simulating a difficult or even unrealizable[14]quantum system using a controllable condensed matter system is explored in several ways.[15−17]

To determine the mass distribution — which on its turn can also be interpreted as an effective potential —which exactly reproduces the eigenfunctions and energy eigenvalues of the target charge-dyon system we proceed as in Ref.[13].We substitute both the exact eigenfunctions and energy eigenvalues of the charge-dyon system into the Schrodinger equation and solve it not forψbut the massM.In the non-relativistic case,the Hamiltonian bears the so-called ordering problem,and we use the Zhu-Kroemer parametrization[18]which is one of the two parametrizations physically admissible.[19]The differential equation for the massMbecomes non-linear,and we look for simple solutions studying a two-point boundary value problem numerically.Still,in the relativistic problem,using the PDM Klein-Gordon equation,we manage to obtain an exact solution for the mass distribution as a function of spherical coordinates (r,θ).

The outline for our paper is the following: In Sec.2 we derive the mapping equation satisfied by the mass distribution in the non-relativistic case.In Sec.3 we turn to the relativistic case,utilizing the PDM Klein-Gordon equation,and obtain an exact solution form=m(r,θ).Finally,in Sec.4 we conclude the work.

2 Mapping a Non-Relativistic Dyon-Charge System into Position-Dependent Effective Mass Background

In this section we introduce the study of the quantum system composed of a non-relativistic charged point particle in the presence of a dyon,using a position-dependent effective mass background.The physical model can be obtained from an eigenvalue problem,

where the Hamiltonian when the mass is positiondependent reads,

where the effective potentialUeffhas the form,

in some cases the external potentialVextmay be zero.The parametersα,β,andγobey the constrainα+β+γ=−1.

In our developments we will use spherical coordinates and in the non-relativistic problem we will use natural units — namely h=me=c=1 and,whereαthe so-called fine-structure constant,namelyα ≈1/137.To determine the effective potentialUefffor a particular system,we need to know the mass functionMand also to select a given ordering,that must obey certain physical considerations proposed by Dutra and Almeida.[19]In this work we will use a parametrization of von Roos Hamiltonian which fulfills the test,that is,the Zhu-Kroemer(ZK)ordering whereα=γ=−1/2,andβ=0.This choice produces an effective potentialUeffgiven by,

The main idea is to construct a map between a system composed of a quantum scalar charged particle of massm0=1 interacting with a dyon located at the origin into a charged particle interacting with a position-dependent effective mass background.In order to obtain such a mapping we substitute the exact wavefunction of the chargedyon system into the PDM Schrodinger equation and solve it for the mass.The target system — a charged particle with position-dependent mass,or in other words,a charged particle interacting with the effective potential(4)— reproduces exactly the behavior of the original chargedyon system.

The eigenfunctions for the charged particle interacting with a dyon are

where the angular part is known as generalized spherical harmonics,[20]

where(cosθ) are the Jacobi polynomials,µ=egis the product of electric and magnetic charges — and it must be such that Dirac quantization condition holds,[2]i.e.,µ=n/2,wherenis an integer number,we consider here,whereQis the electric charge of dyon.The quantum numbers are such thatl=µ,µ+1,...and−l ≤m ≤l.The normalized radial functions reads,

wherek2=−2m0E >0,the function1F1is the confluent hypergeometric function,Nis a radial quantum numberN=0,1,2,...andℓis related to the angular quantum numberlvia

2.1 Non-Linear Differential Equation for the Mass Distribution

Substituting the eigenfunctionψdyoninto the PDM Schrodinger equation with Hamiltonian (2),using the effective potential (4),we obtain

where we used ZK-parametrization.So as to search a differential equation for the mass distributionM(r,θ,ϕ),let us use the following property of theYµlmfunctions,

substituting the above result into Eq.(9) we obtain,

using the confluent differential equation the first four terms yields,

and recalling that the generalized spherical harmonics are eigenfunctions ofLz,we write

We intend to investigate if there exists a solution,so let us simplify the equation takingm+µ=0.This constrain gives an azimuthal symmetry to the system and Eq.(9) becomes,

which is a quite difficult equation.Let us tackle this equation considering a special where the mass depends only on the radial coordinate,consequently∂M/∂θ=0,using the gradient and Laplacian in spherical coordinates,we obtain a non-linear ODE for the mass distribution functionM(r).This differential equation reads,

where the term 4M/rinvolves the producteQ,remembering that we are consideringeQ=−1.In the system composed of a charged particle and a dyon the energy spectrum is discrete,because of the electric chargeQof dyon,so the energy eigenvalues are

namely these are the allowed bound states of the charged particle interacting with a static dyon field.

2.2 Numerical Solutions

Considering that the mass depends only on the radial coordinate,namelyM=M(r),we need to take a particular value ofθ,sayθ0,in order to solve Eq.(15)numerically.Let us choose the initial conditionsM(ri)=M(rf)=1/2 forri=0.22,andrf=1,and letθ0=π/4.We also search for solutions wherem+µ=0 so let us takeµ=l=−m=1.Equation (15) will be numerically solved by theMathematicasoftware — this is an example of the so-called two-point boundary value problem—only for certain values of energy,given by the expression (16).

Fig.1 (Color online) Plot of numerical solutions of Eq.(15) for µ=1 and l=1.Boundary conditions are M(r) using M(ri) = M(rf)=1/2 for ri=0.22 and rf =1,here we fix θ0 =π/4.As close as one approaches the origin,where the dyon is located,the solution diverges.

Fig.2 (Color online) Plot of numerical solutions of Eq.(15) for some values of θ.In these solutions we consider µ=1, l=1,and N=0.Boundary conditions are M(r) using M(ri) = M(rf)=1/2 for ri=0.22 and rf =1,for θ between 10◦,and 140◦approximately,there is a continuum of radial solutions.

Fig.3 (Color online) Effective potentials calculated using solutions of Eq.(15) with M(ri) = M(rf)=1/2 for θ0 = π/4, ri=0.22,and rf=1.The three curves for N=0,1,2 (blue,yellow and green) are overlapped.Parameters are µ=l =1.

The solutions for three values of quantum numberNare presented in Fig.1.The solutions tend to infinity nearriand possess local maximum.The effective potential corresponding to each solution is given by Eq.(3)and the curves for each radial solution are plotted in Fig.2 and for this case it has no local minimum.The radial solutions also depend onθ,this dependence is presented in Fig.3,where we enumerate some values ofθand their respective radial solutions.

3 Mapping a Relativistic Dyon-Charge System into Position-Dependent Effective Mass Background - The Klein-Gordon Equation

In the present section we accomplish the mapping between a relativistic charge-dyon system and positiondependent mass system.This mapping is carried out using the Klein-Gordon equation introducing a positiondependent mass.The Klein-Gordon equation for constant mass reads,

in the relativistic treatment we are going to use natural units where h=me=c=1 andthe d’Alembertian operator reads,

and let us assume the mass to be dependent on spherical coordinates (r,θ,ϕ),that is,m=m(r,θ,ϕ).Taking the Laplacian in spherical coordinates we obtain the following expression to the PDM Klein-Gordon equation,

following the same methodology of the previous sections we need to substitute the exact solution for the relativistic charge-dyon system into PDM Klein-Gordon equation.This solution is well-known,see for instance,[20]and it reads,

whereωis the eigenvalue of energy (in present section we use the discrete part of the spectrum),FNℓ(r) is given in Eq.(7) andYµlm(θ,ϕ) are the generalized spherical harmonics.Replacing the above exact solution into Eq.(19) and solving for the mass yields,

utilizing Eq.(10) and considering the azimuthal symmetry,m+µ=0,we can set up the following expression,

The second step is to deal with the radial part.The radial equation for the relativistic charge-dyon system via Klein-Gordon equation is the confluent differential equation,

using this equation to eliminate the terms involving the derivatives ofFNℓ(r) and using also the relativistic relationk2=ω2we can obtain the exact mass distribution for the mapping between a charge-dyon system and PDM system in the relativistic regime.So our exact mass distribution reads,

wherem0is scalar particle mass interacting relativistically with the dyon located at the origin.So what we intend to do is to reproduce this scenario in a PDM background,keeping in mind there is no real dyon,but instead a relativistic scalar charged particle having a positiondependent mass,which can also be interpreted as a constant mass particle interacting with an effective potentialUeff.The discrete energy spectrum for this bound system has the following form,

whereN=0,1,2,...is the quantum number for energy discrete levels.The relation between the angular quantum numbersℓandlis,

It is important to remember that the dyon shares the same problem the Dirac magnetic monopole has: the vector potential has a singularity which gives rise to the wellknown Dirac string.[2]It is either located atθ=0 or atθ=πsince the vector potential has a pole at one of these points.Wu and Yang[21]correctly interpreted the problem and presented a solution introducing two vector potentials related by a gauge transformation and thus eliminated the string.In Fig.4(a)we plot the mass distribution which would produce the charge-dyon system with the Dirac string located atθ=π.Note that the mass increases asθapproaches the string.In fact it diverges forθ →π.In Fig.4(b) we plot the mass distribution following the Wu-Yang construction,namely,we introduce a vector potential valid for 0≤θ ≤π/2+αthat has a singularity atθ=π,and another vector potential valid forπ/2−α ≤θ ≤πand singular atθ=0.Both of them are related by a gauge transformation.The first one produces a mass distribution (24) and the second an analogous function where instead of (1−cosθ)2we have(1+cosθ)2.So,a position-dependent electric charge such thatm=m(r,θ) behaves as it would be interacting with a dyon located at the origin.

Fig.4 (Color online) Plot of analytical solution of Eq.(19) given by Eq.(24).In Fig.4(a) we plot the mass distribution which would produce the pure Dirac magnetic monopole.Observe that the mass increases as one approaches the Dirac string which is located at θ = π.In Fig.4(b) we plot the mass distribution which would produce the Dirac magnetic monopole without the string using the well-known Wu-Yang construction.[21] The above position-dependent mass reproduces exactly the eigenfunction of a relativistic charge-dyon system withµ=1, l =0, N =1,and eQ=−1.

4 Conclusion

We introduced a mapping between a charge-dyon and a PDM system in non-relativistic and relativistic regimes,in the former we replaced the exact solution (7) into PDM Schrodinger equation to obtain the non-linear partial differential equation (14) for the mass distribution,so considering a radial dependence only,we solved this equation numerically and plotted the results.The treatment for the relativistic case is quite similar: instead of using PDM Schrodinger equation we tackled the PDM Klein-Gordon equation (19).The ordering problem does not occur because of the absence of product between the mass and momentum operators,the last case is then more straightforward and gives us an analytical expression for the mass distributionm(r,θ),the plot of this solution is shown in Fig.4.In Table 1,we present a summary of our results.This approach could lead to an experimental model,using a condensed matter system,which could simulate a charge-dyon system.

Table 1 Summary of the mapping between charge-dyon and PDM systems in the non-relativistic and relativistic regimes,respectively.