Abstract
A solution to the problem of a hydrogenic atom in a homogeneous dielectric medium with a concentric spherical cavity using the oscillator representation method (ORM) is presented. The results obtained by the ORM are compared with a known exact analytic solution. The energy levels of the hydrogenic atom in a spherical cavity exhibit a shallowdeep instability as a function of the cavity radius. The sharpness of the transition depends on the value of the dielectric constant of the medium. The results of the ORM agree well with the results obtained by the analytic solution when the shallowdeep transition is not too sharp (i.e., when the dielectric constant is not too large) for all values of the cavity radius. The ORM results in the zeroth order approximation diverge significantly in the region of the shallowdeep transition (i.e., for the values of the radius where the shallowdeep transition occurs) when the dielectric constant is high and as a result the transition is sharp. Even for the sharp transition, the ORM results again agree very well with the analytic results at least for the ground state when a commonly used approximation in the ORM is removed. The ORM methodology for the cavity model presented in this article can potentially be used for twoelectron systems in a quantum dot.
1 Introduction
We consider a hydrogenlike atom of nuclear charge
e
in a medium of dielectric constant
ϵ
with a concentric spherical cavity of radius
r
0
around the nucleus and free space dielectric constant
ϵ
0
. The following potential represents such a system:
(1)
V
(
r
)
=
−
e
2
ϵ
0
r
+
(
1
−
ϵ
0
ϵ
)
e
2
ϵ
0
r
0
θ
(
r
0
−
r
)
−
e
2
ϵ
r
θ
(
r
−
r
0
)
,
where
θ
(
x
)
is the Heaviside function which is 0 for
x
< 0 and 1 for
x
≥ 0.
Chaudhuri and Coon [1] treated a more general version of this problem in which the effective mass of the electron is also different in the medium outside the cavity and provided an exact analytic solution. Von Roos [2] pointed out that the Hamiltonian with positiondependent mass is not Hermitian and hence should be abandoned. Still it has been widely and successfully used in semiconductor physics particularly in quantum well problems. Considerable interest and work on the positiondependent effective mass have continued (see ref. [3] and references therein). It is not clear how to treat the positiondependent effective mass in the oscillator representation method (ORM) and it is not so important for our purposes either. So we have not included the positiondependent effective mass in the problem treated in this article. It should also be noted that the second term in
V
(
r
)
was added [1] to ensure the continuity of the potential at
r
0
although it is not an absolute requirement that the potential must be continuous. While this term makes the potential continuous, the electric field is still discontinuous.
It is clear from a simple inspection of the potential that in the limits of
r
0
→
∞
and
r
0
→
0
,
the energy levels are identical to those of a hydrogen atom in free space and in the medium of dielectric constant
ϵ
, respectively. Thus, the energy undergoes a transition from a shallow to a deep level in some range of
r
0
as it increases from 0 to
∞
. The magnitude of the transition is proportional to
ϵ
2
. As discussed in ref. [1], the reason for this transition is as follows. If the free space wavefunction, or more accurately the probability function, fits almost entirely inside the cavity, then the effect of the medium outside the cavity is small. If, on the other hand, the cavity size is such that significant portions of the free space wavefunction fall both inside and outside the cavity, then there is a competition to squeeze in and to stretch out the wavefunction by the two regions resulting in an intermediate energy value. The sharpness of the transition increases as the value of
ϵ
increases. This transition is akin to a secondorder phase transition.
In Section 2, we reproduce the solution to the Schrödinger equation (SE) with the potential in equation (1) provided in ref. [1]. In Sections 3 and 4, we outline the general principles of the ORM developed by Dineykhan and Efimov [4] and present the results of the application of the ORM to the current problem, respectively. In Section 5, we compare and discuss the results obtained by the two methods for different parameters.
2 Cavity model: analytic solution
The SE with the potential given in equation (1) is spherically symmetric. Consequently, the wavefunction can be separated in spherical coordinates and written as
(2)
Ψ
ν
l
m
=
N
ν
l
m
R
ν
l
(
r
)
Y
l
m
(
θ
,
ϕ
)
,
where
N
ν
l
m
is the normalization constant,
R
ν
l
is the radial part of the wavefunction, and
Y
l
m
s are the spherical harmonics. The radial part of the SE in dimensionless form is written as:
(3a)
ρ
⁎
d
2
d
ρ
⁎
2
+
(
2
l
+
2
−
ρ
⁎
)
d
d
ρ
⁎
+
(
ν
⁎
−
l
−
1
)
ρ
⁎
×
u
ν
l
(
ρ
⁎
)
=
0
for
r
<
r
0
,
(3b)
ρ
d
2
d
ρ
2
+
(
2
l
+
2
−
ρ
)
d
d
ρ
+
(
ν
−
l
−
1
)
ρ
×
v
ν
l
(
ρ
)
=
0
for
r
>
r
0
,
where
ρ
⁎
=
2
r
/
(
ν
⁎
a
B
)
,
ρ
=
2
r
ϵ
0
/
(
ν
ϵ
a
B
)
, and
a
B
=
ϵ
0
ℏ
2
m
e
2
is the Bohr radius. The radial wavefunction
R
(
r
)
has been written in terms of
u
and
v
in the two regions
(4)
R
(
r
)
=
R
<
(
r
)
=
A
u
ν
l
ρ
⁎
l
exp
(
−
ρ
⁎
/
2
)
for
r
<
r
0
R
>
r
=
B
v
ν
l
ρ
l
exp
(
−
ρ
/
2
)
for
r
>
r
0
.
The constants A and B are obtained by matching the wavefunctions at
r
=
r
0
and from the normalization condition. The energy eigenvalue
E
is related to the effective principal quantum number
ν
by
(5)
E
=
−
m
e
4
ϵ
2
ℏ
2
1
2
ν
2
,
and
ν
⁎
is related to
ν
by
(6)
1
2
ν
⁎
2
=
ϵ
0
ϵ
2
1
2
ν
2
+
1
−
ϵ
0
ϵ
a
B
r
0
.
Equations 3(a) and 3(b) are confluent hypergeometric equations, each of which has two independent solutions. To obtain a physical solution, the solutions of equations 3(a) and 3(b) are selected so that they are regular at
r
=
0
and
r
=
∞
,
respectively. Chaudhuri et al. [1] provided the following energy eigenvalue condition by matching the logarithmic derivative of the radial wavefunctions corresponding to the two selected solutions appropriate for the two regions:
(7)
1
−
−
ν
+
l
+
1
l
+
1
Φ
(
l
+
2
−
ν
⁎
,
2
l
+
3
;
ρ
0
⁎
)
Φ
(
l
+
1
−
ν
⁎
,
2
l
+
2
;
ρ
0
⁎
+
ν
⁎
ν
ϵ
0
ϵ
1
−
2
ρ
0
Ψ
(
−
ν
−
l
−
1
,
−
2
l
−
1
;
ρ
0
Ψ
(
−
ν
−
l
,
−
2
l
;
ρ
0
)
=
0
.
The functions
Φ
and
Ψ
are confluent hypergeometric functions related to the confluent hypergeometric series _{1}F_{1} and _{2}F_{0}, respectively. While the series representation of the function
Φ
(see equation (15) in ref. [1]) is convergent for all finite values of
ρ
⁎
, the series representation of
Ψ
has zero radius of convergence. However, an integral representation of
Ψ
(see equation (16) in ref. [1]) that satisfies the convergence requirements is used to calculate the function. The eigenvalues
ν
for a given value of the angular momentum quantum number
l
can be calculated by locating the zeroes of the expression in equation (7).
3 ORM
Dineykhan and Efimov [4] developed the ORM arising from ideas and methods of the quantum field theory. Using the ORM they calculated the binding energies of a number of systems with various types of potentials including the Coulomb and powerlaw potentials, exponentially screened Coulomb potential, logarithmic potential [4,5], and a twoelectron quantum dot in a magnetic field [6]. The ORM results for the Coulomb and power law, the exponentially screened Coulomb, and logarithmic potentials agree very well with the results obtained by variational numerical methods. Amin and ElAsser have applied the ORM to calculate the energy spectrum of hydrogenlike atoms in a van der Waals potential [7].
The first key step in the ORM is a transformation of the variables in the SE such that the wavefunction takes a Gaussian asymptotic form. Schrödinger in a paper [8] on solving eigenvalue problems by factorization pointed out the existence of such a transformation in which the Kepler problem is transformed into an oscillator problem in four dimensions. The modified SE in the new expanded space having the Gaussian asymptotic solution exhibit oscillator behavior at large distances.
In the next steps, the canonical variables (coordinate and momentum) in the transformed space are represented in terms of creation and annihilation operators
a
^{†} and
a
, and the Hamiltonian is written in terms of normal ordered products over
a
^{†} and
a
. The normal order operation (also known as Wick’s transform) is an operation in which all the creation operators
a
^{†} are moved to the left and all the annihilation operators
a
are moved to the right. The normal order operation is frequently used in quantum field theory. Wick’s transform :
q
n
: yields the nth order Hermite polynomial in
q
which, of course, is the harmonic oscillator wavefunction apart from the exponential term. Wick’s transform is in fact used in quantum field theory to eliminate infinity arising from the zeropoint energy. See ref. [9] for an exposition on Wick’s calculus. The pure oscillator part with some yet unknown frequency, ω, is extracted from the Hamiltonian written in the form
H
=
H
0
+
H
I
+
ε
0
, where
H
0
=
ω
a
†
a
. In addition, a requirement is imposed that the interaction part,
H
I
,
does not contain terms quadratic in the canonical variables so that they are completely absorbed in the oscillator part. This requirement leads to the following condition:
(8)
∂
ε
0
∂
ω
=
0
,
which determines
ω
, the oscillator frequency. This condition is known as the oscillator requirement condition (ORC).
Another interesting aspect of the ORM is that the dimension of the hyperspace can be a variational parameter and as a result can be noninteger. Even though an integer dimension is used to derive the energy equations, the dimension appears in the end results just as a parameter and thus can be varied to obtain the energy minimum.
Dineykhan and Efimov [4] used the following transformations that provide the Gaussian asymptotic wavefunction behavior:
(9)
r
=
q
2
ρ
,
(10)
r
ψ
(
r
)
=
q
a
ϕ
(
q
)
.
The radial part of the SE in the new variable
q
in the expanded space is then obtained as follows:
(11)
−
1
2
∂
2
∂
q
2
+
d
−
1
q
∂
∂
q
+
W
l
(
q
;
E
)
ϕ
q
=
ε
(
E
)
ϕ
(
q
)
,
where
(12)
W
l
(
q
;
E
)
=
−
K
(
l
,
ρ
,
d
)
2
q
2
+
4
ρ
2
(
q
2
)
(
2
ρ
−
1
)
(
V
(
q
2
)
−
E
)
,
(13)
d
=
2
a
−
2
ρ
+
2
,
and
(14)
K
(
l
,
ρ
,
d
)
=
1
4
[
(
d
−
2
)
2
−
4
ρ
2
(
2
l
+
1
)
2
]
.
It should be noted here that
d
is the dimension of the hyperspace, the energy
E
is now incorporated in the new potential
W
l
,
and
q
2
=
q
i
q
i
summed over the repeated index
i
from 1 to the dimension
d
. As mentioned earlier,
d
is assumed to be an integer for the calculations of
ε
0
,
ε
,
W
l
, until at the end when
ρ
and
d
can be treated as noninteger variational parameters. We also note that the SE is written in dimensionless form by using the length unit as
ℏ
2
/
m
e
2
and the energy unit as
m
e
4
/
ℏ
2
.
The energy spectrum
E
n
l
of the original system is obtained from the radial excitation spectrum
ε
[
n
r
]
of the Hamiltonian of equation (11)
(15)
H
(
E
n
l
)
ϕ
[
n
r
]
(
q
)
=
ε
[
n
r
]
ϕ
[
n
r
]
(
q
)
,
(
n
r
=
0
,
1
,
2
,
…
)
,
and it is determined by
(16)
ε
[
n
r
]
(
E
n
l
)
=
ε
[
n
r
]
(
l
,
ρ
,
d
;
E
)
=
0
,
where
n
r
is the radial quantum number.
If the potential
V
(
r
)
does not have a repulsive character as
r
0
→
0
,
K
(
l
,
ρ
,
d
)
is typically chosen to be zero and that leads to the equation for the dimension,
(17)
d
=
4
ρ
l
+
2
ρ
+
2
.
The oscillator representation is then obtained by writing the Hamiltonian
H
in the form
(18)
H
=
1
2
(
p
2
+
ω
2
q
2
)
+
W
l
q
−
1
2
ω
2
q
2
,
and by introducing the usual creation and annihilation operators,
a
†
and
a
, respectively, in terms of the canonical variables
q
_{ j } and
p
_{ j } (
j = 1, 2,...,
d), as:
(19)
q
j
=
1
2
ω
a
j
+
a
j
†
,
p
j
=
1
i
ω
2
a
j
−
a
j
†
.
The creation and annihilation operators satisfy the standard commutation relation
(20)
a
j
,
a
k
†
=
δ
j
k
,
j
,
k
=
1
,
2
,
…
d
.
After normal ordering the products over
a
^{†} and
a
, in which all the creation operators are moved to the left and the annihilation operators are moved to the right as if they commute and some manipulation the ORM Hamiltonian is obtained as
(21)
H
=
H
0
+
H
I
+
ε
0
,
where
(22)
H
0
=
ω
a
j
†
a
j
,
(23)
H
I
=
∫
−
∞
∞
d
k
2
π
d
W
˜
(
k
)
exp
(
−
k
2
4
ω
)
:
e
2
i
(
k
q
)
:
,
and
(24)
ε
0
=
d
ω
4
+
∫
−
∞
∞
d
k
2
π
d
W
˜
(
k
)
exp
−
k
2
4
ω
.
As mentioned earlier, the symbol :*: represents the normal ordering of the products over
a
^{†} and
a
. In
k
q
=
k
j
q
j
as well as in
k
2
and
q
2
summation over repeated indices are assumed and the notation
e
2
x
=
e
x
−
1
−
x
−
1
2
x
2
is introduced. The function
W
˜
is the Fourier transform of
W
(
q
)
in
d
dimension and given as
(25)
W
˜
(
k
)
=
∫
−
∞
∞
(
d
ξ
)
d
W
l
(
ξ
)exp
(
i
k
ξ
)
,
k
ξ
=
k
j
ξ
j
.
The energy spectrum is then obtained by calculating the contribution of the interaction part
H
I
of the Hamiltonian in the perturbation approach. In the zeroth order approximation, the energy spectrum is determined by
(26)
2
n
r
ω
+
n
r

H
I

n
r
+
ε
0
=
0
.
See ref. [4] for the details of the calculation of the radial eigenstates

n
r
and the matrix elements
n
r

H
I

n
r
. In this article, we will consider only the zeroth order approximation in
H
I
.
4 Cavity model: ORM solution
In this section, we will apply the ORM to the cavity model defined by the potential in equation (1). We note that the calculation of the energy spectrum using the wavefunction matching method is a natural one when the potential is discontinuous. However, to apply the ORM to the cavity model it may be useful for the calculations of various terms to write the potential in a continuous form as follows:
(27)
V
(
r
)
=
−
1
r
−
(
B
−
1
)
r
0
θ
(
r
0
−
r
)
−
(
B
−
1
)
r
θ
(
r
−
r
0
)
,
(28)
θ
(
x
)
=
P
(
τ
)
e
−
i
x
τ
=
lim
σ
→
0
+
∫
−
∞
+
∞
e
−
i
x
τ
τ
+
i
σ
d
τ
,
where
B
=
(
ϵ
0
/
ϵ
)
. The second and the third terms can be treated as exponential and screened Coulomb potentials, respectively, with complex exponentials by using equation (
28) for the Heaviside function. The final results are obtained by applying the operation
P
(
τ
)
over the functions of
(
i
τ
)
.
Now we make the typical approximation in which there is no repulsive character in the potential by setting
K
(
l
,
ρ
,
d
)
=
0
and then
d
is given by equation (17). With this approximation, we obtain after some manipulations the pertinent ORM equations for the cavity model potential given in equation (27):
(29)
ε
0
=
a
0
ω
−
4
E
ρ
2
ω
2
ρ
−
1
Γ
d
2
+
2
ρ
−
1
Γ
d
2
−
a
1
ω
ρ
−
1
−
a
2
ω
2
ρ
−
1
γ
l
d
2
+
2
ρ
−
1
,
ω
r
0
1
ρ
−
a
3
ω
ρ
−
1
γ
u
d
2
+
ρ
−
1
,
ω
r
0
1
/
ρ
,
(30)
n
r

H
I

n
r
=
−
a
1
B
T
(
ρ
−
1
)
ω
ρ
−
1
−
4
E
ρ
2
ω
2
ρ
−
1
Γ(
d
/
2
+
2
ρ
−
1
)
×
T
(
2
ρ
−
1
)
−
a
4
T
12
(
ρ
),
where
(31a)
a
0
=
d
4
,
a
1
=
4
ρ
2
Γ(
d
/
2
+
ρ
−
1
)
Γ(
d
/
2
)
,
a
2
=
4
(
B
−
1
)
ρ
2
Γ(
d
/
2
+
2
ρ
−
1
)
r
0
Γ(
d
/
2
)
,
a
3
=
4
(
B
−
1
)
ρ
2
Γ(
d
/
2
+
ρ
−
1
)
Γ(
d
/
2
)
,
and
a
4
=
4
(
B
−
1
)
ρ
2
Γ(
d
/
2
)
(31b)
T
(
x
)
=
Γ
(
d
/
2
)
Γ
(
n
r
+
1
)
Γ
(
d
/
2
+
n
r
)
x
(
x
−
1
)
Γ(
d
/
2
+
n
r
+
1
)
θ
(
n
r
−
1
)
Γ
(
d
/
2
+
2
)
Γ
3
(
n
r
)
+
∑
t
=
2
n
r
∑
p
=
t
2
t
2
2
t
−
p
.
(
x
)
p
.
(
p
)
t
.
(
t
)
2
t
−
p
×
Γ
(
d
/
2
+
p
+
n
r
−
t
)
Γ
(
p
+
1
)
Γ
(
t
+
1
)
Γ
(
d
/
2
+
p
)
Γ
3
(
n
r
−
t
+
1
)
Γ
(
p
−
t
+
1
)
Γ
(
2
t
−
p
+
1
)
;
(31c)
(
x
)
p
=
x
(
x
−
1
)
(
x
−
2
)
…
(
x
−
p
+
1
)
;
(31d)
T
12
(
ρ
)
=
−
ρ
ω
r
0
1
/
ρ
d
/
2
+
ρ
−
1
ω
ρ
−
1
×
∑
k
=
0
∞
(
−
1
)
k
ω
r
0
1
/
ρ
k
T
(
−
d
/
2
−
k
)
(
d
/
2
+
2
ρ
−
1
+
k
)
(
d
/
2
+
ρ
−
1
+
k
)
Γ
(
k
+
1
)
;
and
γ
l
and
γ
u
are incomplete gamma functions defined as
(31e)
γ
l
(
a
,
x
)
=
1
Γ
(
a
)
∫
0
x
t
a
−
1
e
−
t
d
t
,
γ
u
(
a
,
x
)
=
1
Γ
(
a
)
∫
x
∞
t
a
−
1
e
−
t
d
t
.
In order to obtain the energy spectrum we use the ORC in equation (8) with
ε
0
and
n
r

H
I

n
r
given in equations (29) and (30). We solve the ORC equation for
E
and plug in the expression for
E
in equations (26) and (30) to obtain the equation for
ω
as follows:
(32)
1
+
T
(
2
ρ
−
1
)
2
ρ
a
0
+
n
r
(
2
ρ
−
1
)
ρ
ω
−
a
1
2
1
+
B
(
2
ρ
−
1
)
ρ
T
(
ρ
−
1
)
−
(
ρ
−
1
)
ρ
T
(
2
ρ
−
1
)
+
a
3
2
γ
u
d
/
2
+
ρ
−
1
,
ω
r
0
1
/
ρ
1
−
(
ρ
−
1
)
ρ
T
(
2
ρ
−
1
)
ω
1
−
ρ
+
a
2
(
2
ρ
−
1
)
2
ρ
γ
l
d
/
2
+
2
ρ
−
1
,
ω
r
0
1
/
ρ
T
(
2
ρ
−
1
)
ω
1
−
2
ρ
−
a
4
(
2
ρ
−
1
)
2
ρ
T
12
(
ρ
)
=
0
.
With a solution for
ω
obtained from equation (32) the energy spectrum
E
is then obtained from equations (8) and (30) as
(33)
E
(
n
r
,
l
;
ρ
)
=
−
Γ(
d
/
2
)
ω
2
ρ
4
ρ
2
(
2
ρ
−
1
)
Γ(
d
/
2
+
2
ρ
−
1
)
×
a
0
+
a
1
(
ρ
−
1
)
ω
−
ρ
+
a
2
(
2
ρ
−
1
)
γ
l
d
/
2
+
2
ρ
−
1
,
ω
r
0
1
/
ρ
ω
−
2
ρ
+
a
3
(
ρ
−
1
)
γ
u
d
/
2
+
ρ
−
1
,
ω
r
0
1
/
ρ
ω
−
ρ
.
The energy spectrum given by equation (33) is still a function of the parameter
ρ
. The energy spectrum
E
(
n
r
,
l
)
is obtained by minimizing
E
(
n
r
,
l
;
ρ
)
with respect to
ρ
as
(34)
E
(
n
r
,
l
)
=
E
(
n
r
,
l
;
ρ
min
)
=
min
ρ
E
(
n
r
,
l
;
ρ
)
,
in which
ω
for a given value of
ρ
is determined by finding the solution to equation (
32).
An inspection of equations (31b) and (31d) indicates that
T
(
0
)
=
T
(
1
)
=
T
12
=
0
for
n
r
=
0
states (i.e., 1s, 2p, 3d…). Therefore, there is no contribution from the interaction part of the Hamiltonian,
H
I
,
for these states.
In the limiting case of
B
=
1
, the solution of equation (32) yields
ω
=
2
/
n
, where
n
=
n
r
+
l
+
1
is the principal quantum number, and the energy spectrum appropriately reduces to that of a free space hydrogen atom,
E
(
n
r
,
l
)
=
−
1
2
n
2
. In the limiting case of
r
0
→
0
,
ω
=
2
B
/
n
, and the energy spectrum correctly reduces to that of a hydrogen atom in a medium of dielectric constant
ϵ
, i.e.,
E
(
n
r
,
l
)
=
−
B
2
2
n
2
. In the limiting case of
r
0
→
∞
, it can be shown that
T
(
0
)
=
T
(
1
)
=
T
12
=
0
as the parameter
ρ
tends to 1 and the energy spectrum reduces to that of the free space hydrogen atom.
In some situations, fixing the parameter
ρ
to a value
1
may be sufficient. Noting that
T
(
0
)
=
T
(
1
)
=
0
when
ρ
=
1
, equations (32) and (33) reduce to the simpler forms:
(35)
ω
−
2
n
−
2
n
[
γ
l
(
2
l
+
2
,
r
0
ω
)
+
B
γ
u
(
2
l
+
2
,
r
0
ω
)
]
−
2
n
(
B
−
1
)
Γ
(
2
l
+
2
)
T
12
(
1
)
=
0
,
(36)
E
n
l
=
−
ω
2
8
1
+
8
(
B
−
1
)
r
0
ω
2
γ
l
(
2
l
+
3
,
r
0
ω
)
,
and
(37)
T
12
(
1
)
=
−
(
ω
r
0
)
2
l
+
2
∑
k
=
0
∞
(
−
1
)
k
×
(
ω
r
0
)
k
T
(
−
2
l
−
2
−
k
)
(
2
l
+
3
+
k
)
(
2
l
+
2
+
k
)
Γ
(
k
+
1
)
.
As mentioned earlier,
T
12
=
0
for
n
r
=
0
states (i.e., 1s, 2p, 3d…).
5 Results and discussion
First, we consider the ground state, i.e., the 1s state for which
n
=
1
,
n
r
=
0
,
and
l
=
0
.
In Figure 1(a–d), we have plotted
ν
as a function of
r
0
for a set of parameters
B
(0.088, 0.25, 0.5, and 0.75) obtained by the two methods presented earlier for the 1s state. For the ORM, we have computed the values of
E
1
s
for the
ρ
=
1
case by using equations (35–37). As noted earlier,
T
(
ρ
)
=
T
(
2
ρ
−
1
)
=
T
12
(
ρ
)
=
0
even when
ρ
≠
1
for the
1
s
state. The values of
ν
are then obtained by using equation (5). We have computed the values of
ν
for the analytic solution directly by solving equation (7) for
ν
. The value of
B
=
0.088
corresponds to silicon for which the dielectric constant
ϵ
=
11.4
ϵ
0
.
As seen from the analytic results the shallow to deep transition of the energy level occurs near
r
0
=
1
, where the wavefunction is maximum. Evidently, the ORM results agree very well with the analytic results for all values of
r
0
when the dielectric constant of the medium outside the cavity does not differ significantly from that of the free space cavity. However, the ORM results diverge significantly in the neighborhood of
r
0
where the energy level undergoes shallow to deep transition when the dielectric constant of the medium outside the cavity differs significantly from that of the free space cavity. For the realistic case of silicon (
B
=
0.088
)
the disagreement is stark.
Now in order to eliminate the
ρ
=
1
restriction we computed the
1
s
energy level
E
1s
(
ρ
min
)
by finding the minima of the energy with respect to
ρ
following equations (32), (33), and (34). The dashed line in Figure 2 represents
E
1s
(
ρ
min
)
as a function of
r
0
. The value of
ρ
min
was found to vary between 1 and 2.14. As expected, the values of
ρ
min
differ from 1 only in the transition region. With the minimization of the energy with respect to the parameter
ρ
, or equivalently, the dimension
d
, the ORM results agree very well with the exact analytic results including in the transition region even for the case of sharp transition (
B
=
0.088
).
For all states for which
n
r
=
0
that includes the ground state (1s), the contribution from the interaction part
n
r

H
I

n
r
is zero since
T
(
x
)
=
T
12
=
0
. The first state for which the contribution from the interaction part is nonzero is the 2s state. We present the results in Figure 3(a and b) for the values of
B
=
0.088
and 0.5.
The solid lines in Figure 3(a and b) represent the analytic results. The 2s radial probability function has two lobes with maxima at
r
=
0.76
and at
r
=
5.24
.
Consequently, there are two sharp transitions close to these two values of
r
0
with a plateau in between.
The ORM results with
ρ
=
1
diverge widely in the transition region for the 2s state. So we have shown the ORM results only with
ρ
=
ρ
min
that minimizes the energy. The value of
ρ
min
was found to vary between 1 and 2.45 (i.e.,
d
=
4
to
6.9
)
for
B
=
0.088
and from 1 to 1.28 (i.e.,
d
=
4
to
4.56
)
for
B
=
0.5
.
Even with the minimization with respect to
ρ
,
the ORM result for the sharp transition case (
B
=
0.088
) differs significantly from the analytic result in the transition region. The origin of the discrepancy for this case can be attributed to the extremely shallow zerocrossing of the
ω
function given in equation (32), where the value of the
ω
function remains very close to zero for a wide range of
ω
values. A small change in the
ω
function would move the zerocrossing from the left to the right of the wide shallow region of
ω
. A small error in the
ω
function thus would change the calculated energy significantly. This indicates that higher order perturbation terms in the interaction part may be necessary to obtain better ORM results for sharp transition cases for states with
n
r
≥
1
.
From this analysis, it appears that at least for the ground state in which there is no contribution from the interaction part of the Hamiltonian,
H
I
, the ORM provides excellent results when the energy is minimized with respect to the parameter,
ρ
. On the other hand, for states where the interaction part is nonzero, such as the 2s state, the ORM results are not so accurate in the transition region when the transition is sharp. It appears that for these cases higher order perturbation terms in the interaction contribution are necessary for more accurate results.
Dineykhan and Nazmitdinov [6] provided analytical results for a twoelectron quantum dot in a magnetic field with a harmonic oscillator potential for the quantum dot using the ORM. The problem of a quantum dot in a magnetic field was treated using numerical method [10,11] and analytical method with some approximations [12]. A singlet to triplet state transition was predicted for the ground state in a quantum dot as a function of the applied magnetic field. These transitions were experimentally observed by Ashoori et al. [13]. The harmonic oscillator potential has been used [10,11,13] to represent the confining potential in a quantum dot. However, a quantum dot having a finite size, a finite potential should provide a better model. Ashoori et al. conjectured that the discrepancy of their experimental values of the magnetic fields at which the singlet to triplet transition occurs with the theoretical values to the strictly harmonic confining potentials used in the theoretical calculations. Recently, Chaudhuri [14] obtained analytic expressions for the energy levels of a 2electron system in a 2dimensional quantum dot modeled with a finite Gaussian potential and subjected to a magnetic field using an ORM approach like the one presented in this article. In the appropriate limits, the results are shown to match very well with previous analytical and numerical results. Using the expressions, the magnetic field at which the spinsinglet to spintriplet ground state crossing occurs is calculated. The calculated value is closer to the experimental value compared with the infinite harmonic potential results in previous theoretical models.
References
[1] Chaudhuri S, Coon DD. Analytic solution of the cavity model for interstitial impurities in semiconductors. Phys Rev B. 1984;30:3338–43. Search in Google Scholar
[2] von Roos O. Positiondependent effective masses in semiconductor theory. Phys Rev B. 1983;27:7547–52. Search in Google Scholar
[3] Abdalla MS, Eleuch H. AIP Adv. 2016;6(055011):1. Search in Google Scholar
[4] Dineykhan M, Efimov GV. Rep Math Phys. 1995;36:287. Search in Google Scholar
[5] Dineykhan M, Efimov GV. The oscillator representation and the stability of threebody Coulomb systems. FewBody Syst. 1994;16:59–90. Search in Google Scholar
[6] Dineykhan M, Nazmitdinov RG. Twoelectron quantum dot in a magnetic field: analytical results. Phy. Rev B. 1997;55(20):13707–14. Search in Google Scholar
[7] Amin ME, ElAsser MA. An analytical treatment of the energy spectrum of hydrogenlike atoms perturbed by a generalized van der Waals potential. Braz J Phys. 2009;39(2):301–5. Search in Google Scholar
[8] Schrödinger E. Proc R Irish Acad. 1941;46:183. Search in Google Scholar
[9] Wurm A, Berg M. Wick calculus. Am J Phys. 2008;76(1):65–72. Search in Google Scholar
[10] Maksym PA, Chakrabarty T. Quantum dots in a magnetic field: role of electronelectron interactions. Phys Rev Lett. 1990;65:108–11. Search in Google Scholar
[11] Wagner M, Merkt U, Chaplik AV. Spinsinglet–spintriplet oscillations in quantum dots. Phys Rev B. 1992;45:1951–4. Search in Google Scholar
[12] Taut M. Two electrons in a homogeneous magnetic field: particular analytical solutions. J Phys A. 1994;27:1045–55. Search in Google Scholar
[13] Ashoori RC, Stormer HL, Weiner JSPLN, Baldwin KW, West KW. Nelectron ground state energies of a quantum dot in magnetic field. Phys Rev Lett. 1993;71:613–6. Search in Google Scholar
[14] Chaudhuri S. Twoelectron quantum dot in a magnetic field: Analytic solution for finite potential model. Phys E Low Dimens Syst Nanostruct. 2021;128:114571. Search in Google Scholar
Received: 20200128
Revised: 20201002
Accepted: 20201022
Published Online: 20210305
© 2021 Sid Chaudhuri, published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.