Welcome to visit Communications in Theoretical Physics,
Quantum Physics and Quantum Information

An analytical study of the magnetic properties of the quantum dot/ring system: Rashba, magnetic field, and electron-electron effects

  • M Shirsefat 1 ,
  • M Servatkhah , 1, 2 ,
  • S Hosseini 1, 2
Expand
  • 1Department of Physics, Marvdasht Branch, Islamic Azad University, Marvdasht, Iran
  • 2Department of Physics, Nanotechnology Research Center, Marvdasht Branch, Islamic Azad University, Marvdasht, Iran

Received date: 2023-12-20

  Revised date: 2024-09-09

  Accepted date: 2024-09-12

  Online published: 2024-12-20

Copyright

© 2025 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing. All rights, including for text and data mining, AI training, and similar technologies, are reserved.

Abstract

The magnetic behavior of a two-electron quantum dot/ring system is analytically studied with electron–electron (e–e) interaction taking into account the Rashba spin–orbit interaction (SOI) and magnetic field. The Jacobi transformation has been employed to separate the Hamiltonian of the system to the center of mass and relative terms. The Schrödinger equation is analytically solved, and energy spectra are obtained. Then, the magnetization and susceptibility are calculated. The magnetization decreases by raising the magnetic field without and with SOI, and also without e-e interaction. Also, the SOI slightly modifies the magnetization of the system without e–e interaction. The susceptibility displays a peak structure as the magnetic field changes from low values to high values. The susceptibility by considering e–e interaction and without the SOI is always negative and its value decreases by rising the magnetic field. The susceptibility displays a transition from diamagnetic to paramagnetic with e–e interaction and SOI.

Cite this article

M Shirsefat , M Servatkhah , S Hosseini . An analytical study of the magnetic properties of the quantum dot/ring system: Rashba, magnetic field, and electron-electron effects[J]. Communications in Theoretical Physics, 2025 , 77(4) : 045106 . DOI: 10.1088/1572-9494/ad7c37

1. Introduction

The physical behaviors of nanostructures are interesting and have potential applications in engineering science. Among the nanostructures, quantum dots (QDs) and quantum rings (QRs) have received much attention. The study of thermodynamic and magnetic properties of QDs and QRs are important challenges in the science of physics [17].
The confinement potentials have a crucial role in investigating the thermodynamic and magnetic properties of quantum dot/ring systems. The available confinement potentials such as the Fu–Wang–Jia potential [8, 9] and its special cases also have a substantial role in the accurate prediction of thermodynamic properties for real substance systems and the successful prediction of equilibrium constants for some important chemical reactions [1014].
In past years, many researchers have tried to study the physical properties of nanostructures analytically [1519]. Also, many authors have tried to investigate the effects of electron–electron (e–e) interaction, magnetic field, and spin–orbit interaction (SOI) on the physical properties of QDs and QRs [2025]. As we know, analytical study is very valuable. Although, due to the complicated form of the e–e interaction and SOI, the analytical study of the physical properties of nanostructures requires considerable theoretical effort. The reader can refer to [2628].
One of the most important physical properties of nanostructures are the magnetic properties. Hitherto, many papers have reported on the magnetic properties of QDs and QRs [2934]. Recently, Shi et al [35] investigated dynamic magnetic properties and the magnetocaloric effect of a decorated square structure. It is to be noted that the magnetic properties of QDs and QRs change under various factors like SOI, impurity, e–e interaction, and magnetic field. Several articles have reported on the magnetic properties of QDs under Rashba SOI and e–e interaction [2932, 3640].
According to the importance of the subject, we have tried to study the magnetic properties of a quantum dot/ring system under a magnetic field, e–e interaction, and Rashba SOI. To reach this goal, we analytically solve the Schrödinger equation and obtain the electron energy levels. Then, we determine magnetization, and magnetic susceptibility under the aforementioned factors.

2. Theory and model

Consider a dot/ring system with two electrons under an external magnetic field $B,$ and the Rashba SOI. The total Hamiltonian of the system is expressed as:
$\begin{eqnarray}H={H}_{0}+{H}_{R}+\displaystyle \frac{{e}^{2}}{\left|{{\boldsymbol{r}}}_{2}-{{\boldsymbol{r}}}_{1}\right|}+{H}_{Z}.\end{eqnarray}$
The first term $\left({H}_{0}\right)$ is the Hamiltonian of two electrons without the e–e interaction and the Rashba effect. The second term is the Rashba Hamiltonian, the third term is the e–e interaction, and the fourth term is the Zeeman Hamiltonian. The Hamiltonian of ${H}_{0}$ is given by
$\begin{eqnarray}{H}_{0}=\displaystyle \frac{1}{2{m}^{\ast }}\displaystyle \sum _{j=1}^{2}{\left({{\boldsymbol{P}}}_{j}+e{{\boldsymbol{A}}}_{j}\right)}^{2}+V\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right),\end{eqnarray}$
where ${m}^{* }$ is the electron effective mass, and ${{\boldsymbol{A}}}_{j}\,=\tfrac{B}{2}\left(-{y}_{j},{x}_{j},0\right)=\tfrac{B{r}_{j}}{2}{\boldsymbol{\varphi }}$ is the vector potential. Also, $V\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right)$ is the confining potential with two terms, a parabolic-inverse square (PIS), and a Gaussian term. One can write $V\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right)={V}_{1}\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right)+{V}_{2}\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right).$
The PIS potential is written as [29],
$\begin{eqnarray}{V}_{1}\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right)=\displaystyle \frac{1}{2}{m}^{\ast }{\omega }_{0}^{2}{\left|{{\boldsymbol{r}}}_{2}-{{\boldsymbol{r}}}_{1}\right|}^{2}+\displaystyle \frac{{\hslash }^{2}\xi }{2{m}^{\ast }}\displaystyle \frac{1}{{\left|{{\boldsymbol{r}}}_{2}-{{\boldsymbol{r}}}_{1}\right|}^{2}},\end{eqnarray}$
where ${\omega }_{0}$ is the harmonic frequency, and $\xi $ is a dimensionless parameter. The Gaussian potential is written as:
$\begin{eqnarray}{V}_{2}\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right)=-{V}_{0}\exp \left[-{\left(\displaystyle \frac{\left|{{\boldsymbol{r}}}_{2}-{{\boldsymbol{r}}}_{1}\right|}{{R}_{0}}\right)}^{q}\right],\end{eqnarray}$
where ${V}_{0}$ and ${R}_{0}$ are the depth and range of the potential. We suppose that $q=2$ and $r/{R}_{0}\ll 1,$ hence the Gaussian potential can be expressed by:
$\begin{eqnarray}{V}_{2}\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right)=-\left({V}_{0}-\displaystyle \frac{{V}_{0}}{{R}_{0}^{2}}{\left|{{\boldsymbol{r}}}_{2}-{{\boldsymbol{r}}}_{1}\right|}^{2}\right).\end{eqnarray}$
Combining equations (3) and (5), we have:
$\begin{eqnarray}\begin{array}{c}V\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right)=-{V}_{0}\\ \,+\,\left(\displaystyle \frac{1}{2}{m}^{\ast }{\omega }_{0}^{2}+\displaystyle \frac{{V}_{0}}{{R}_{0}^{2}}\right){\left|{{\boldsymbol{r}}}_{2}-{{\boldsymbol{r}}}_{1}\right|}^{2}+\displaystyle \frac{{\hslash }^{2}\xi }{2{m}^{\ast }}\displaystyle \frac{1}{{\left|{{\boldsymbol{r}}}_{2}-{{\boldsymbol{r}}}_{1}\right|}^{2}}.\end{array}\end{eqnarray}$
The Rashba Hamiltonian is given by
$\begin{eqnarray}{H}_{R}=\displaystyle \frac{{\alpha }_{R}}{\hslash }\displaystyle \sum _{j=1}^{2}\left[{{\boldsymbol{\sigma }}}_{j}\times \left({{\boldsymbol{P}}}_{j}+e{{\boldsymbol{A}}}_{j}\right)\right].\end{eqnarray}$
Here, ${\alpha }_{R}$ is the Rashba coupling parameter, and $\sigma \,=\left({\sigma }_{x},{\sigma }_{y},{\sigma }_{z}\right)$ are the Pauli matrices. The Zeeman term is given by
$\begin{eqnarray}{H}_{Z}=\displaystyle \frac{1}{2}{g}^{\ast }{\mu }_{B}B\displaystyle \sum _{j=1}^{2}{\sigma }_{j,z},\end{eqnarray}$
where ${g}^{* }$ is the Landé factor of the electron, and ${\mu }_{B}$ is the Bohr magneton.
We employ the Jacobi coordinates [33, 35] as follows:
$\begin{eqnarray}\begin{array}{c}{\boldsymbol{R}}=\displaystyle \frac{{{\boldsymbol{r}}}_{1}+{{\boldsymbol{r}}}_{2}}{2},\,{\boldsymbol{r}}=\,{{\boldsymbol{r}}}_{2}-{{\boldsymbol{r}}}_{1},\,\\ {{\boldsymbol{P}}}_{R}={{\boldsymbol{P}}}_{1}+{{\boldsymbol{P}}}_{2},\,{{\boldsymbol{P}}}_{r}=\displaystyle \frac{{{\boldsymbol{P}}}_{2}-{{\boldsymbol{P}}}_{1}}{2}.\end{array}\end{eqnarray}$
Therefore, we obtain the following relations:
$\begin{eqnarray}\begin{array}{c}\displaystyle \frac{1}{2{m}^{\ast }}\displaystyle \sum _{j=1}^{2}{\left({{\boldsymbol{P}}}_{j}+e{{\boldsymbol{A}}}_{j}\right)}^{2}=\,\displaystyle \frac{{{P}_{R}}^{2}}{2M}+\displaystyle \frac{{{P}_{r}}^{2}}{2\mu }+\displaystyle \frac{1}{2}{\omega }_{c}{l}_{zR}\\ \,+\,\displaystyle \frac{1}{2}{\omega }_{c}{l}_{zr}+\displaystyle \frac{1}{8}M{\omega }_{c}^{2}{R}^{2}+\displaystyle \frac{1}{8}\mu {\omega }_{c}^{2}{r}^{2},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{c}V\left({{\boldsymbol{r}}}_{1},{{\boldsymbol{r}}}_{2}\right)=-{V}_{0}\\ \,+\,\left(\displaystyle \frac{1}{2}{m}^{\ast }{\omega }_{0}^{2}+\displaystyle \frac{{V}_{0}}{{R}_{0}^{2}}\right){r}^{2}+\displaystyle \frac{{\hslash }^{2}\xi }{2{m}^{* }}\displaystyle \frac{1}{{r}^{2}},\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{c}{H}_{R}=\displaystyle \frac{{\alpha }_{R}}{\hslash }\displaystyle \sum _{j=1}^{3}\left[{\sigma }_{x}^{\left(j\right)}\left({P}_{y}^{\left(j\right)}+e{A}_{y}^{\left(j\right)}\right)-{\sigma }_{y}^{\left(j\right)}\left({P}_{x}^{\left(j\right)}+e{A}_{x}^{\left(j\right)}\right)\right]\\ \,=\,\displaystyle \frac{{\alpha }_{R}}{\hslash }\left(\begin{array}{cc}0 & {H}_{12}^{R}\\ {H}_{21}^{R} & 0\end{array}\right),\end{array}\end{eqnarray}$
where
$\begin{eqnarray}\begin{array}{c}{H}_{12}^{R}=\displaystyle \sum _{j=1}^{2}\left[\displaystyle \frac{\hslash }{{\rm{i}}}\left(\displaystyle \frac{\partial }{\partial {y}_{j}}+{\rm{i}}\displaystyle \frac{\partial }{\partial {x}_{j}}\right)+\displaystyle \frac{eB}{2}\left({x}_{j}-{\rm{i}}{y}_{j}\right)\right]\\ \,=\,\displaystyle \frac{\hslash }{{\rm{i}}}\left(\displaystyle \frac{\partial }{\partial {R}_{y}}+{\rm{i}}\displaystyle \frac{\partial }{\partial {R}_{x}}\right)+eB\left({R}_{x}-{\rm{i}}{R}_{y}\right),\end{array}\end{eqnarray}$
$\begin{eqnarray}\begin{array}{c}{H}_{21}^{R}=\displaystyle \sum _{j=1}^{2}\left[\displaystyle \frac{\hslash }{{\rm{i}}}\left(\displaystyle \frac{\partial }{\partial {y}_{j}}-{\rm{i}}\displaystyle \frac{\partial }{\partial {x}_{j}}\right)+\displaystyle \frac{eB}{2}\left({x}_{j}+{\rm{i}}{y}_{j}\right)\right]\\ \,=\,\displaystyle \frac{\hslash }{{\rm{i}}}\left(\displaystyle \frac{\partial }{\partial {R}_{y}}-{\rm{i}}\displaystyle \frac{\partial }{\partial {R}_{x}}\right)+eB\left({R}_{x}+{\rm{i}}{R}_{y}\right).\end{array}\end{eqnarray}$
According to equations (8)–(14), the total Hamiltonian is separated into the center of mass and relative coordinate parts, $H={H}_{r}+{H}_{R}.$
$\begin{eqnarray}\begin{array}{c}{H}_{r}{\rm{\Phi }}=-\displaystyle \frac{{\hslash }^{2}}{2\mu }\left[\displaystyle \frac{{{\rm{d}}}^{2}}{{\rm{d}}{r}^{2}}+\displaystyle \frac{1}{r}\displaystyle \frac{{\rm{d}}}{{\rm{d}}r}+\displaystyle \frac{1}{{r}^{2}}\displaystyle \frac{{{\rm{d}}}^{2}}{{\rm{d}}{\theta }^{2}}\right]{\rm{\Phi }}\\ \,+\,\displaystyle \frac{1}{2}{\omega }_{c}{l}_{zr}{\rm{\Phi }}+\displaystyle \frac{{e}^{2}}{r}{\rm{\Phi }}+\displaystyle \frac{{\hslash }^{2}\xi }{2\mu }\displaystyle \frac{1}{{r}^{2}}{\rm{\Phi }}\\ \,-\,{V}_{0}{\rm{\Phi }}+\displaystyle \frac{1}{2}\mu {{\rm{\Omega }}}^{2}{r}^{2}{\rm{\Phi }}={E}_{r}{\rm{\Phi }},\end{array}\end{eqnarray}$
where ${{\rm{\Omega }}}^{2}=\tfrac{{\omega }_{c}^{2}}{4}+{\omega }_{0}^{2}+\tfrac{2{V}_{0}}{\mu {R}_{0}^{2}}.$
$\begin{eqnarray*}\begin{array}{l}{H}_{R}{\rm{\Psi }}=-\displaystyle \frac{{\hslash }^{2}}{2M}\left[\displaystyle \frac{{\partial }^{2}}{\partial {R}_{x}^{2}}+\displaystyle \frac{{\partial }^{2}}{\partial {R}_{y}^{2}}\right]{\rm{\Psi }}+\displaystyle \frac{1}{2}{\omega }_{c}{l}_{zR}{\rm{\Psi }}\\ \,+\,\displaystyle \frac{1}{8}M{\omega }_{c}^{2}{R}^{2}{\rm{\Psi }}+\left(\begin{array}{cc}\displaystyle \frac{1}{2}{g}^{\ast }{\mu }_{B}B & 0\\ 0 & -\displaystyle \frac{1}{2}{g}^{\ast }{\mu }_{B}B\end{array}\right){\rm{\Psi }}\end{array}\end{eqnarray*}$
$\begin{eqnarray}+\,\displaystyle \frac{{\alpha }_{R}}{\hslash }\left(\begin{array}{cc}0 & {H}_{12}^{R}\\ {H}_{21}^{R} & 0\end{array}\right){\rm{\Psi }}={E}_{R}{\rm{\Psi }}.\end{eqnarray}$
To solve equation (15), consider the wave function as ${\rm{\Phi }}\left(r,\theta \right)\,=\tfrac{{{\rm{e}}}^{{\rm{i}}m\theta }}{\sqrt{2\pi }}f\left(r\right),$ and obtain:
$\begin{eqnarray}\begin{array}{c}\displaystyle \frac{{{\rm{d}}}^{2}f}{{\rm{d}}{r}^{2}}+\displaystyle \frac{1}{r}\displaystyle \frac{{\rm{d}}f}{{\rm{d}}r}-\displaystyle \frac{{m}^{2}+\xi }{{r}^{2}}f\\ \,-\,\displaystyle \frac{2\mu }{{\hslash }^{2}}\left[\displaystyle \frac{1}{2}\mu {{\rm{\Omega }}}^{2}{r}^{2}+\displaystyle \frac{{e}^{2}}{r}-{V}_{0}-\displaystyle \frac{m\hslash }{2}{\omega }_{c}-{E}_{r}\right]f=0.\end{array}\end{eqnarray}$
Here, we define the following notations:
$\begin{eqnarray}\begin{array}{c}\displaystyle \frac{2\mu \left({E}_{r}+{V}_{0}\right)}{{\hslash }^{2}}+\displaystyle \frac{m\mu {\omega }_{c}}{\hslash }=\,{\varepsilon }_{r},\,\displaystyle \frac{{\mu }^{2}{{\rm{\Omega }}}^{2}}{{\hslash }^{2}}={\lambda }^{2},\,\displaystyle \frac{2\mu {e}^{2}}{{\hslash }^{2}}=\eta .\end{array}\end{eqnarray}$
Considering the radial wave function as $f\left(r\right)\,={{\rm{e}}}^{-c{r}^{2}-dr}\displaystyle \displaystyle \sum _{n}{a}_{n}{r}^{n+k}$ where $\gamma ,$ and $\delta $ are constant. Inserting $f\left(r\right)$ into equation (17), we have:
$\begin{eqnarray}\begin{array}{c}\displaystyle \displaystyle \sum _{n=0}{a}_{n}\left[\left(n+k\right)\left(n+k-1\right)+\left(n+k\right)-\left({m}^{2}+\xi \right)\right]{r}^{n+k-2}\\ \,+\,\displaystyle \displaystyle \sum _{n=0}{a}_{n}\left[-2d\left(n+k\right)-\eta -d\right]{r}^{n+k-1}\\ \,+\,\displaystyle \displaystyle \sum _{n=0}{a}_{n}\left[-4c\left(n+k\right)+{d}^{2}-4c+{\varepsilon }_{r}\right]{r}^{n+k}\\ \,+\,\displaystyle \displaystyle \sum _{n=0}{a}_{n}\left[4cd\right]{r}^{n+k+1}+\displaystyle \sum _{n=0}{a}_{n}\left[4{c}^{2}-{\lambda }^{2}\right]{r}^{n+k+2}=0.\end{array}\end{eqnarray}$
According to equation (19), the following equations can be obtained:
$\begin{eqnarray}{k}^{2}={m}^{2}+\xi ,\,{\varepsilon }_{r}=4c\left(n+k\right)-{d}^{2}+4c,\end{eqnarray}$
$\begin{eqnarray}4{c}^{2}-{\lambda }^{2}=0,\,2d\left(n+k\right)+\eta +d=0.\end{eqnarray}$
Using equations (18), (20), and (21), we get the following energy eigenvalues:
$\begin{eqnarray}\begin{array}{c}{E}_{r}=\hslash {\rm{\Omega }}\left(n+1+\sqrt{{m}^{2}+\xi }\right)\\ \,-\,\displaystyle \frac{2\mu {e}^{4}}{{\hslash }^{2}{\left(2n+1+2\sqrt{{m}^{2}+\xi }\right)}^{2}}+\displaystyle \frac{1}{2}m\hslash {\omega }_{c}-{V}_{0}.\end{array}\end{eqnarray}$
Now, equation (16) can be written as:
$\begin{eqnarray*}\begin{array}{c}-\displaystyle \frac{{\hslash }^{2}}{2M}\left[\displaystyle \frac{{\partial }^{2}}{\partial {R}_{x}^{2}}+\displaystyle \frac{{\partial }^{2}}{\partial {R}_{y}^{2}}\right]{\rm{\Psi }}+\displaystyle \frac{1}{2}{\omega }_{c}{l}_{zR}{\rm{\Psi }}\\ +\,\displaystyle \frac{1}{8}M{\omega }_{c}^{2}{R}^{2}{\rm{\Psi }}+\left(\begin{array}{cc}\displaystyle \frac{1}{2}{g}^{\ast }{\mu }_{B}B & 0\\ 0 & -\displaystyle \frac{1}{2}{g}^{\ast }{\mu }_{B}B\end{array}\right){\rm{\Psi }}\end{array}\end{eqnarray*}$
$\begin{eqnarray}+\,\displaystyle \frac{{\alpha }_{R}}{\hslash }\left(\begin{array}{cc}0 & {H}_{12}^{R}\\ {H}_{21}^{R} & 0\end{array}\right){\rm{\Psi }}={E}_{R}{\rm{\Psi }}.\end{eqnarray}$
Choosing ${R}_{x}=R\,\cos \,\phi ,$ ${R}_{y}=R\,\sin \,\phi ,$ and $R=\sqrt{{R}_{x}^{2}+{R}_{y}^{2}},$ we have:
$\begin{eqnarray}\begin{array}{c}\displaystyle \frac{\partial }{\partial {R}_{x}}=\,\cos \,\phi \displaystyle \frac{\partial }{\partial R}-\displaystyle \frac{\sin \,\phi }{R}\displaystyle \frac{\partial }{\partial \phi },\\ \displaystyle \frac{\partial }{\partial {R}_{y}}=\,\sin \,\phi \displaystyle \frac{\partial }{\partial R}+\displaystyle \frac{\cos \,\phi }{R}\displaystyle \frac{\partial }{\partial \phi }.\end{array}\end{eqnarray}$
Using the above equation, we can write:
$\begin{eqnarray}\begin{array}{c}{H}_{12}^{R}=\displaystyle \frac{\hslash }{{\rm{i}}}{{e}}^{-{\rm{i}}\phi }\left[{\rm{i}}\displaystyle \frac{\partial }{\partial R}+\displaystyle \frac{1}{R}\displaystyle \frac{\partial }{\partial \phi }\right]+{e}BR{{e}}^{-{\rm{i}}\phi },\\ {H}_{21}^{R}=\displaystyle \frac{\hslash }{{\rm{i}}}{{e}}^{{\rm{i}}\phi }\left[-{\rm{i}}\displaystyle \frac{\partial }{\partial R}+\displaystyle \frac{1}{R}\displaystyle \frac{\partial }{\partial \phi }\right]+{e}BR{{\rm{e}}}^{{\rm{i}}\phi }.\end{array}\end{eqnarray}$
Considering ${\rm{\Psi }}=\left(\begin{array}{c}u\left(R\right){{\rm{e}}}^{{\rm{i}}m\phi }\\ v\left(R\right){{\rm{e}}}^{{\rm{i}}\left(m+1\right)\phi }\end{array}\right)$ and inserting it into equation (25), one has:
$\begin{eqnarray*}\begin{array}{c}\displaystyle \frac{{{\rm{d}}}^{2}u}{{\rm{d}}{R}^{2}}+\displaystyle \frac{1}{R}\displaystyle \frac{{\rm{d}}u}{{\rm{d}}R}+\left({k}^{2}-\displaystyle \frac{{m}^{2}}{{R}^{2}}\right)u\\ \,-\,\displaystyle \frac{{\omega }_{c}^{2}}{4{\hslash }^{2}}{R}^{2}u-\displaystyle \frac{2M{\alpha }_{R}}{\hslash }\left(\displaystyle \frac{{\rm{d}}}{{\rm{d}}R}+\displaystyle \frac{m+1}{R}\right)v\end{array}\end{eqnarray*}$
$\begin{eqnarray}-\,\displaystyle \frac{2M{\alpha }_{R}eB}{\hslash }Rv=0,\,\end{eqnarray}$
$\begin{eqnarray*}\displaystyle \frac{{{\rm{d}}}^{2}v}{{\rm{d}}{R}^{2}}+\displaystyle \frac{1}{R}\displaystyle \frac{{\rm{d}}v}{{\rm{d}}R}+\left({q}^{2}-\displaystyle \frac{{\left(m+1\right)}^{2}}{{R}^{2}}\right)v-\displaystyle \frac{{\omega }_{c}^{2}}{4{\hslash }^{2}}{R}^{2}v\,\end{eqnarray*}$
$\begin{eqnarray}\,-\,\displaystyle \frac{2M{\alpha }_{R}}{\hslash }\left(-\displaystyle \frac{{\rm{d}}}{{\rm{d}}R}+\displaystyle \frac{m}{R}\right)u-\displaystyle \frac{2M{\alpha }_{R}{e}B}{\hslash }Ru=0,\end{eqnarray}$
where,
$\begin{eqnarray}{k}^{2}=-\displaystyle \frac{Mm{\omega }_{c}}{\hslash }-\displaystyle \frac{M{g}^{\ast }{\mu }_{B}B}{{\hslash }^{2}}+\displaystyle \frac{2M{E}_{R}}{{\hslash }^{2}},\end{eqnarray}$
$\begin{eqnarray}{q}^{2}=-\displaystyle \frac{M\left(m+1\right){\omega }_{c}}{\hslash }+\displaystyle \frac{M{g}^{\ast }{\mu }_{B}B}{{\hslash }^{2}}+\displaystyle \frac{2M{E}_{R}}{{\hslash }^{2}}.\end{eqnarray}$
To solve equations (26) and (27), the following solutions are considered,
$\begin{eqnarray}u=A{{\rm{J}}}_{m}\left(kR\right)+B{{\rm{J}}}_{m}\left(qR\right),\end{eqnarray}$
$\begin{eqnarray}v=A{{\rm{J}}}_{m+1}\left(kR\right)+B{{\rm{J}}}_{m+1}\left(qR\right).\end{eqnarray}$
Inserting equations (30) and (31) into (26) and (27). Then, using the recurrence relations, and Bessel differential equation, the energy levels ${E}_{R}$ of the system can be obtained.

2.1. Susceptibility and magnetization

The partition function is expressed by,
$\begin{eqnarray}Z=\displaystyle \displaystyle \sum _{s}{{\rm{e}}}^{-\beta {E}_{s}},\end{eqnarray}$
where ${E}_{s}$ is the energy eigenvalue, $\beta =1/{k}_{B}T,$ and ${k}_{{\rm{B}}}$ is Boltzmann’s constant. The magnetization, and susceptibility are given by,
$\begin{eqnarray}M=\displaystyle \frac{1}{Z}\displaystyle \displaystyle \sum _{s}\left(-\displaystyle \frac{\partial {E}_{s}}{\partial B}\right){{\rm{e}}}^{-\beta {E}_{s}},\end{eqnarray}$
$\begin{eqnarray}\chi =\displaystyle \frac{1}{Z}\displaystyle \displaystyle \sum _{s}\left(-\displaystyle \frac{{\partial }^{2}{E}_{s}}{\partial {B}^{2}}\right){{\rm{e}}}^{-\beta {E}_{s}}.\end{eqnarray}$

3. Numerical results

We have chosen GaAs to show the numerical results. The following parameters have been employed: ${m}^{\ast }=0.067\,{m}_{0},$ ${g}^{\ast }=-0.44,$ ${\alpha }_{R}$ = 5.4 meV-nm, $\hslash {\omega }_{0}$ = 0.019 meV, and ${V}_{0}$ = 300 meV.
In figure 1(a), a schematic of the quantum dot/quantum ring (QD/QR) structure is shown. This system is composed of a QD separated from a surrounding concentric QR by a tunneling barrier ${V}_{0}.$ Figure 1(b) displays the curve of the actual single-particle potential energy as a function of distance with arbitrary values [32].
Figure 1. Schematic diagram of QD/QR structure; (b) the shape of the actual single-particle potential energy, with arbitrary values.
Figure 2 displays the reduced energy $\left(E/\hslash {\omega }_{0}\right)$ of the two-electron system versus the magnetic field, and without Rashba SOI. Two curves correspond to the ground and first excited states. We know that in the absence of a magnetic field, the energy levels of various states within an atomic orbital are identical. If a magnetic field is applied, the states split according to the Zeeman effect and take a great variety of values. The amount of splitting is entirely dependent on the strength of the field applied. As can be seen, the energy increases as the magnetic field enhances. In the absence of Rashba effect, as the magnetic field increases, the system receives more energy and the energy values ​​increase. According to the obtained relations, the energy of the system depends on the magnetic field value. We can see from figure 2 that a singlet → triplet transition at a critical value of ${B}_{c}$ in the presence of e–e interaction.
Figure 2. Reduced energy $\left(E/\hslash {\omega }_{0}\right)$ versus the magnetic field with e–e interaction and without SOI.
Figure 3 demonstrates reduced magnetization $\left(M/{\mu }_{B}^{* }\right)$ of the two-electron system versus reduced magnetic field $\left({\omega }_{0}/{\omega }_{c}\right)$ without the e–e interaction at $T=1\,{\rm{K}}.$ The plotted curves correspond to the presence and absence of Rashba SOI. It can be seen that the magnetization increases by enhancing $B$ without and with SOI. We know that at low temperatures ($T=1\,{\rm{K}}$ ), only the lower energy states are more likely to be occupied. From the figure, it is observed that the SOI slightly modifies the magnetization of the system. When $B$ rises, the degeneracy is lifted due to the influence of the field on the Landau level and Zeeman splitting.
Figure 3. Reduced magnetization $\left(M/{\mu }_{B}^{* }\right)$ versus the reduced magnetic field $\left({\omega }_{0}/{\omega }_{c}\right)$ without e–e interaction with T = 1 K.
In figure 4, the susceptibility is plotted versus $B$ for different temperatures with and without e–e interaction. As can be seen, the susceptibility rises by increasing the magnetic field, and it has a peak structure. The susceptibility reduces with rising temperature because the system disorder increases. In the presence of e–e interaction, the maximum value of the curve moves to lower magnetic fields. By considering e–e interaction, the system receives more addition energy (third term in equation (1)). When the susceptibility reaches a maximum, the system disorder decreases. We see also the width and height of the plotted curves depend on the temperature. In the absence of SOI, the system order depends on temperature, magnetic field, and e–e interaction. There is a competition between three factors. Since the used temperatures are low, the main competition is between e–e interaction and magnetic field.
Figure 4. The susceptibility versus the magnetic field without SOI.
Figure 5 presents the susceptibility versus temperature without and with e–e interaction for $B$ = 0.5 T. As can be seen, the susceptibility rises by increasing the temperature without and with e–e interaction. With rising the temperature, it seems that a negative saturation occurs. This means that the system obtains a specific order. By considering e–e interaction, more addition energy enters to the system and susceptibility ($\chi =-\tfrac{{\partial }^{2}F}{\partial {B}^{2}}$ ) decreases in the negative direction. It is to be noted that the curves in figures 4 and 5 have been plotted without Rashba SOI.
Figure 5. The susceptibility versus the temperature without SOI for $B$ = 0.5 T.
In figure 6, the susceptibility has been presented versus temperature without and with e–e interaction taking into account Rashba SOI at zero magnetic field. As can be seen, the susceptibility displays a transition from diamagnetic to paramagnetic at a special temperature. The maximum value of the curve moves to lower temperatures in the presence of e–e interaction. Also, the width and height of the curves change with e–e interaction. By considering the SOI, there is an energy splitting as ${E}_{+}$ and ${E}_{-}.$ The occupation probability of each energy branch depends on the temperature, and e–e interaction. By considering e–e interaction, the system obtains more order at low temperatures and the susceptibility shows a maximum at low temperature.
Figure 6. The susceptibility versus the temperature with SOI for $B$ = 0.
Figure 7 demonstrates the susceptibility versus the temperature in the presence and absence of e–e interaction for $B=1\,{\rm{T}}.$ The curves were plotted in the presence and absence of the Rashba SOI. There is a paramagnetic peak in the absence of e–e interaction. The maximum value of the curve with SOI is lower than that without SOI. By increasing the temperature, the thermodynamic fluctuations enhance and the system shows a diamagnetic behavior with and without e–e interaction and Rashba SOI. We also observe that the susceptibility has negative and positive values. This means that the system shows diamagnetic and paramagnetic behaviors. The regions of diamagnetic and paramagnetic depend on the values of the temperature.
Figure 7. The susceptibility versus the temperature with SOI for $B$ = 1 T.
In figure 8, we plotted energies of the two low-lying states of a two-electron QD as a function of field $B$ for ${V}_{0}$ = 20 meV and ${R}_{0}$ = 10 nm. Our results have been compared with the Gaussian potential model [28]. As the magnetic field increases, the energy of the state $l$ = 0 increases but the energy of the states with non-vanishing $l$ decreases. This leads to a sequence of different ground states with $l$ = 0, −1, −2, −3, −4, … with an increase in magnetic field. An interesting phenomenon of energy level crossing happens because the coulomb energy gets smaller and smaller with increasing $\left|l\right|$ along with the average distance between the electrons. In the presence of the e–e interaction, the energy is higher because the e–e interaction is supposed to increase the total energy.
Figure 8. Energies of the two low-lying states of a two-electron QD as a function of magnetic field for ${V}_{0}$ = 20 meV and ${R}_{0}$ = 10 nm. Our results have been compared with the Gaussian potential model [28].

4. Conclusions

The effects of e–e interaction, Rashba SOI, and magnetic field on the susceptibility of GaAs quantum dot/ring system have been theoretically studied. According to the results, it is deduced that the magnetic field, Rashba SOI, and e–e interaction have important roles in the magnetic behavior of the system. The system can show both diamagnetic and paramagnetic behaviors with a proper selection of system parameters. The transition from paramagnetic to diamagnetic depends on temperature. It is also shown that e–e interaction causes singlet–triplet transitions in our system. The position of peak in the susceptibility has been shown to depend on the temperature and also the strength of the e–e interaction. This system can be of interest from the point of view of magnetic device applications. In future, the thermodynamic and also the magnetocaloric properties of the system in the presence of SOI can be studied.

Competing interests

The authors declare that there is no competing interests that might beperceived to influence the results and/or discussion reported in this paper.

Funding

The authors declare that no funds, grants, or other support were received during the preparation of this manuscript

Author contributions

In this work M Servatkhah proposed and conceived the main idea as the supervisor (Conceptualization). M Shirsefat developed the theory and conducted the formal analysis. Programming has been done by M Shirsefat and S. Hosseini (Advisor). M Servatkhah and S Hosseini verified the analytical methods and investigated the findings of this work (Methodology and validation). The original draft of the paper was written by M Servatkhah, M Shirsefat and S Hosseini, all authors discussed the results and contributed to the review and editing of the final corrected manuscript. The whole process has been done under the supervision of M Servatkhah. The descriptions are accurate and agreed by all authors.

Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

1
Sakr M R 2014 In-plane electron g-factor anisotropy in nanowires due to the spin–orbit interaction Physica E 64 68 71

DOI

2
Gong M, Duan K, Li C F, Magri R, Narvaez G A, He L 2008 Electronic structure of self-assembled InAs/InP quantum dots: comparison with self-assembled InAs/GaAs quantum dots Phys. Rev. B 77 045326

DOI

3
Pershin Y V, Nesteroff J A, Privman V 2004 Effect of spin–orbit interaction and in-plane magnetic field on the conductance of a quasi-one-dimensional system Phys. Rev. B 69 121306

DOI

4
Srouji R, Sakr M R 2015 Thermopower of a ballistic nanowire subjected to an in-plane magnetic field in the presence of Rashba and Dresselhaus interactions Physica E 68 210 214

DOI

5
Pal S, Ghosh M, Duque C A 2019 Impurity related optical properties in tuned quantum dot/ring systems Philos. Magaz. 99 2457

DOI

6
Hashemi P, Servatkhah M, Pourmand R 2021 The effect of rashba spin–orbit interaction on optical far-infrared transition of tuned quantum dot/ring systems Opt. Quant. Electron. 53 567

DOI

7
Hashemi P, Servatkhah M, Pourmand R 2022 Influence of polaron on second and third harmonic generation of tuned quantum dot/ring systems Opt. Commun. 506 127551

DOI

8
Fu K X, Wang M, Jia C S 2019 Improved five-parameter exponential-type potential energy model for diatomic molecules Commun. Theor. Phys. 103 71

DOI

9
Araújo J P, Ballester M Y 2021 A comparative review of 50 analytical representation of potential energy interaction for diatomic systems: 100 years of history Int. J. Quant. Chem. 121 e26808

DOI

10
Ding Q C, Jia C S, Wang C W, Peng X L, Liu J Y, Zhang L H, Jiang R, Zhu S Y, Yuan H, Tang H X 2023 Unified non-fitting formulation representation of thermodynamic properties for diatomic substances J. Mol. Liq. 371 121088

DOI

11
Liang D C, Zeng R, Wang C W, Ding Q C, Wei L S, Peng X L, Liu J Y, Yu J, Jia C S 2022 Prediction of thermodynamic properties for sulfur dioxide J. Mol. Liq. 352 118722

DOI

12
Liu G H, Ding Q C, Wang C W, Jia C S 2022 Unified explicit formulations of thermodynamic properties for the gas NO2, and gaseous BF2 and AlCl2 radicals Chem. Phys. Lett. 830 140788

DOI

13
Liu G H, Ding Q C, Wang C W, Jia C S 2023 Unified non-fitting explicit formulation of thermodynamic properties for five compounds J. Mol. Struct. 1294 136543

DOI

14
Ding Q C, Wang C W, Peng X L, Lu Y Y, Liu G H, Li G, Wei L S, Chen J Q, Yuan H, Jia C S 2024 Unified analytical formulation regarding thermodynamic properties subject to gaseous metal hydroxides J. Mol. Struct. 1298 137075

DOI

15
Chakraborty T, Manaselyan A, Barseghyan M, Laroze D 2018 Controllable continuous evolution of electronic states in a single quantum ring Phys. Rev. B 97 041304 (R)

DOI

16
Ghazaryan A, Manaselyan A, Chakraborty T 2016 Signatures of Majorana fermions in an elliptical quantum ring Phys. Rev. B 93 245108

DOI

17
Khordad R 2018 Simultaneous effects of electron-electron interactions, Rashba spin–orbit interaction and magnetic field on susceptibility of quantum dots J. Magn. Magn. Mater. 449 510

DOI

18
Kirak M 2021 Magnetic and thermodynamic properties of GaAs quantum dot confined by parabolic-inverse square plus gaussian potential J. Magn. Magn. Mater. 536 167481

DOI

19
Khordad R, Vaseghi B 2019 Magnetic properties in three electrons under Rashba spin–orbit interaction and magnetic field Int. J. Quant. Chem. 119 e25994

DOI

20
Ghanbari A, Khordad R, Taghizadeh F, Nasirizadeh I, Edet C O, Ali N 2022 Impurity effect on thermal properties of tuned quantum dot/ring systems Chem. Phys. Lett. 806 140000

DOI

21
Khordad R, Ghanbari A, Abbasi K, Ghaffaripour A 2023 Harmonic generation of tuned quantum dots including impurity effects J. Comput. Electron. 22 260

22
Khordad R, Mohammadi S A 2022 Analytical study of absorption threshold frequency of tuned quantum dot/ring system Physica B 647 434181

DOI

23
Khordad R, Mohammadi S A 2023 Simultaneous effects of pressure, temperature, and external magnetic field on absorption threshold frequency of tuned quantum dot/ring systems: an analytical study J. Comput. Electron. 22 641

DOI

24
El-Bakkari K, Iqraoun E, Jaouane M, Arraoui R, Sali A, Fakkahi1 A, Ed-Dahmouny A, Azmi H 2023 Magnetic field effect on the linear and nonlinear refractive index and optical absorption of a donor in a GaAs quantum ring Phys. Scr. 98 085102

DOI

25
Hashemi P, Servatkhah M, Pourmand R 2023 Analytical study of the Stark shift effect on the tuned quantum dot/ring systems Pramana-J. Phys. 97 124

DOI

26
Kul E B, Polat M, Guclu A D 2020 Electronic and magnetic properties of graphene quantum dots with two charged vacancies Solid State Commun. 322 114096

DOI

27
Al Shorman M M, Nammas F S, Haddad H, Shukri A A 2018 Heat capacity and entropy of two electrons quantum dot in a magnetic field with parabolic interaction Chin. J. Phys. 56 1057

DOI

28
Sharma H K, Boda A, Boyacioglu B, Chatterjee A 2019 Electronic and magnetic properties of a two-electron Gaussian GaAs quantum dot with spin-Zeeman term: a study by numerical diagonalization J. Magn. Magn. Matter. 469 171

DOI

29
Lafaurie L G, Suaza Y A, Laroze D, Gutiérrez W, Marín J H 2024 Thermodynamic properties of an artificial molecule quantum rings: Geometric and external field effects Physica B 679 415786

DOI

30
Lakaal K, Kria M, El Hamdaoui J, Prasad V V, Nautiyal V V, El-Yadri M, Pérez L M, Laroze D, Feddi E 2022 Polaronic corrections on magnetization and thermodynamic properties of electron-electron in 2D systems with Rashba spin–orbit coupling J. Magn. Magn. Mater. 551 169042

DOI

31
Bejan D, Stan C, Niculescu E C 2018 Optical properties of an elliptic quantum ring: eccentricity and electric field effects Opt. Mater. 78 207

DOI

32
Biborski A, Kądzielawa A P, Gorczyca-Goraj A, Zipper E, Maśka M M, Spałek J 2016 Dot-ring nanostructure: rigorous analysis of many-electron effects Sci Rep. 6 29887

DOI

33
Hassanabadi H, Hamzavi M, Zarrinkamar S, Rajabi A A 2010 Quadratic and coulomb terms for the spectrum of a three-electron quantum dot Few-Body Syst. 48 53

DOI

34
Hassanabadi H, Rahimov H, Zarrinkamar S 2012 Analytical treatment of a three-electron-quantum dot under rashba spin–orbit interaction Few-Body Syst. 52 87

DOI

35
Shi K L, Si N, Jiang W 2024 Dynamic magnetic properties and magnetocaloric effect of a decorated square structure Chin. J. Phys. In press

DOI

36
Goudarzi M, Karamdel J, Hassanabadi Sh Zorriasatein H 2022 Investigation of two-electron quantum dot in a magnetic field by using of quasi-exact-solvable method Solid State Commun. 344 114669

DOI

37
Bhattacharyya K, Debnath D, Chatterjee A 2023 Rashba effect on finite temperature magnetotransport in a dissipative quantum dot transistor with electronic and polaronic interactions Sci Rep. 13 5500

DOI

38
Ibrahim N, Saa M A, Shaer A, Elsaid M 2021 Effects of Rashba spin–orbit interaction and topological defect on the magnetic properties of an electron confined in a 2D quantum dot effects of Rashba spin–orbit interaction and topological defect on the magnetic properties of an electron confined in a 2D quantum dot J. Taibah Univ. Sci. 15 1210

39
Elsaid M K, Alia A A, Shaer A 2020 Rashba spin–orbit interaction effects on thermal and magnetic properties of parabolic GaAs quantum dot in the presence of donor impurity under external electric and magnetic fields Chin. J. Phys. 66 335

DOI

40
Chaudhuri S 2021 Two-electron quantum dot in a magnetic field: Analytic solution for finite potential model Physica E 128 114571

DOI

Outlines

/