1. Introduction
Electronic properties of a 2D semiconductor and a semiconductor bulk differs substantially. As the semiconductor film’s thickness reduces to a nanometric range the ratio between number of surface and bulk atoms is increased substantially providing a drastic variation of its electronic properties. Semiconductor quantum well (QW), a layer of nanometric thickness, sandwiched between two layers of a different semiconductor material, today enjoy widespread use in different optoelectronic applications. Progress achieved at last years in manufacturing of nanostructures have allowed for fabricating layers with topologically nontrivial surfaces such as a twisted single crystalline NbSe3 ribbon in a form of a one-sided Mobius strip [
1] or rolled-up quantum wells (RUQWs), in the form of spirals [
2,
3,
4,
5,
6,
7,
8,
9].
RUQW is a structure consisting of semiconductor layer of nanometer thickness curls up after being released from the substrate as a result of the relaxation of elastic stresses. Both sides of the layer are cylindrical surfaces whose cross sections have the form of spirals. The curvature of this structure and its optoelectronic properties depend strongly on the thickness of the layer and the lattice mismatch of the substrate and the layer. Due to existence of a variable curvature along the spiral cross section, electronic properties of RUQWs have characteristic features, which distinguishes them from planar quantum wells and nanotubes. In the last two decades this rolled-up semiconductors have sparked great attentions and been widely used as optical resonators [
10,
11,
12,
13,
14,
15,
16,
17] to realize laser structures [
18,
19,
20], and the alteration of energy band in self-rolled nanomembranes containing
quantum wells (QWs) has been investigated [
21]. Spectral analyses reveal that the small strain change in structures with different curvatures the optical transitions and the light emission intensity of the GaAs QW layer are essentially influenced by the strain evolution. In addition, the first order Stark effect revealed in these structures opens an opportunity for the fabrication of high performance rolled-QW infrared photo-detectors.
Wide scopes for RUQWs of InGaAs/GaAs heterostructures applications as optical materials for quantum information technology are discussed previously in ref. [
22], where it has been found an essential modification of the band structure of InGaAs semiconductor quantum wells by changing their geometrical configuration. Particularly, the photoluminescence demonstrates a strong energy shift of the valence-band states in comparison to flat structures due to an inversion of the heavy-hole with the light-hole states in a rolled-up InGaAs quantum well, which in addition conduces to a strong change in the optical selection rules for the rolled-up quantum wells. Previously, the relation between the photoluminescence characteristics and curvature for rolled-up GaAs quantum wells has been studied in ref. [
23]. It has been revealed that the room-temperature photoluminescence (PL) characteristics of rolled-up RUQWs show dramatic PL intensity enhancement compared to their planar counterparts. Systematic shift of PL peak position as a function of tube curvature, attributed to strain induced band structure change, is established.
In this work, we analyze theoretically an effect of a possible increase of the probability of the carriers tunneling between adjacent spiral wings in the presence of the external electric field, enabled by a special topology of RUQW, that could be a possible source of other novel useful properties. In order to analyze this effect, one should solve a boundary value problem for the Schrodinger equation inside of a domain, which represents a waveguide formed by a semiconductor tape rolled up around z-axis. In the case, when the tape thickness is essentially smaller in comparison to the separations between wings, the motion across the waveguide is considerably faster than along it, that opens a chance for application the well-known adiabatic approximation (AA) to separate the corresponding coordinates of the electron’s motions along and across the spiral waveguide. Such approach has been used earlier in references [
24,
25] for analyzing the energy spectrum of the free electron confined in infinite-barrier RUQW in a form of the Archimedean spiral. It has been found an effective one-dimensional Schrödinger equation that describes in curvilinear coordinates slow motion of the electron inside the waveguide. In this paper, we use a similar formalism to derive an effective wave equation in curvilinear coordinates corresponding to the low-lying energy levels of the electron bound to the shallow donor. Unlike of a similar problem for a free electron, the movements of the electron bound to the donor at Z-direction and along the spiral no longer can be separated and the corresponding effective wave equation is two-dimensional. In our numerical work, to solve it we use the double Fourier series expansion method.
Our calculations for donors with different locations show that energies shiftings in the presence of the external electric field depend heavily on the geometry of the spiral, the strength, and the orientation of the applied field, inducing in some especial cases a permanent dipole moments. We believe that our numerical procedure gives an efficient tool for the study of impurity-related phenomena in RUQWs. The article is organized as follows. In
Section 2, we describe the formalism, in which we solve in natural coordinates for electron and donor wave equation in RUQW and determine low-lying energies as functions of the external electric field. In
Section 3, we present results of calculations and discussion of the Stark effect for shallow donor confined in RUQW. Finally, the conclusion is drawn in
Section 4.
2. Theoretical model
To analyze the features of RUQW’s electronic properties related to their curvature, let us to consider a model of the homogeneous rolled-up QW represented schematically in
Figure 1, flat at z-direction and rolled-up in x-y plane. We assume that the structure along Z-axis is very long and therefore the free electron wave function depends on the z coordinate as in a plane wave. It allows us to reduce the one-particle problem to study of the movement of an electron only inside helical waveguide. In addition, if the QW’s thickness is very small, the electron movement transversely the waveguide is faster than along it, and one can use an advantage of the well-known adiabatic approximation to separate coordinates corresponding to these two degrees of freedom. For this purpose, it is appropriate to describe the RUQW geometry in natural coordinates given by unit vectors (
) shown in
Figure 1. The vector
describes the tangential direction for the curvilinear coordinate
s corresponding to the arc length measured from the beginning of the spiral (
), the vector
describes the normal direction and the corresponding coordinate
h measures the distance from the exterior border of the QW along the normal, while the vector
describes the vertical direction with corresponding to the coordinate
z that gives the distance from the bottom of the structure (
). Here and in what follows, we denote as
the spiral length, as
d and
L the thickness and the height of the QW, respectively. In addition, we assume that in the adiabatic limit
the normal coordinate
and positions of the donor and electron are defined by vectors
and
, respectively. The external electric field is applied in plane XOY along the line that forms the angle
with the X-axis.
We assume that cross-sections of the RUQW’s along XOY plane are spirals, whose equations in polar coordinates correspond approximately to a generalized Archimedean equation
;
. The arc length
s and the curvature radius
of this curve depend on the polar angle as follows:
In our numerical work, we first found the dependence of the polar angle on the arc length,
by using the inverse cubic spline-interpolation procedure. Once such dependence is found, then one can define similar dependencies for the Cartesian and polar coordinates of the electron and donor positions:
Below, we use the effective Bohr radius
as a dimensionless unit of the distance, and the effective Rydberg,
as a dimensionless unit of the energy. In these units, the effective 2D potential energy of the donor confined in a RUQW is given by the following expression [
24,
25,
26]:
The Four terms in Equation (
3) correspond to the contribution to the potential energy from the lateral confinement inside the nanolayer, from the curvature variation along the spiral [
25,
26], from the external electric field and from the electron-donor attraction, respectively. In
Figure 2, we show curves of potential energies along the Archimedean spiral with parameters
,
and with the donor located at the point with polar coordinate
as functions of the arc length, given by the relations (
1)–(
3) for the zero-electric field (solid line) and for the case
(dash line). One can see in
Figure 2 that the curve of the potential energy undergoes an essential change in the presence of the external electric field obtaining a form of oscillation with increasing amplitude. Such transformation of the potential along the spiral can provide a strong polarizability and the formation a giant dipole moment in the ground state of the shallow donor.
To clarify why the external electric field applied the X-axis can induce a permanent giant moment in a donor we present in
Figure 3 some plots, which demonstrate the successive evolution of the potential curves along the spiral at the cross-section
under increasing electric field for the electron (a) and for the donor (b).
The curve of the free electron’s potential in
Figure 3a for the zero electric field case (solid line) is smooth and it has a shallow minimum at the initial point of the spiral (
), where the curvature of the spiral is the largest. In the presence of the external electric field directed toward initial spiral point along X-axis, this potential is transformed to a set of periodically arranged quantum wells with successively increasing depth at the radial direction (dash and dot lines). Bottoms and tops of quantum wells are located at points of corresponding wings with polar coordinates
and
, respectively. As the electric field grows the bottoms of all wells descend, while their heights are increased. It is seen that larger is the electric field, the bigger should be the displacement of the electron from the initial point of the spiral.
Similar potential curves for the donor placed inside RUQW at the position with polar angle are presented in
Figure 3b. As
(solid line), the main deepest minimum is situated at the donor location while other minima with lesser deep are located at points of nearest wings with polar angles
, situated at the same side of spiral. With an increase in the electric field, additional minima appear on both sides of the main minimum at points with polar angles
, at that the minima on the right side are deeper than those on the left. The larger the electric field, the deeper are additional quantum wells (cf. dashed and dotted curves). One can expect that such alteration of the potential curves under increasing electric field could lead to an instability of the donor ground state when the energy a similar state of one of additional minima at right side becomes lower than of the main minimum and the electron jumps from the wing with the donor toward the right side of the spiral. Such transition can be accompanied by an origination of a giant dipole moment and a strong polarizability of the donor. The final 2D wave equation in dimensionless units for the electron inside a two-dimensional surface with the potential energy (
3) has in natural coordinates the following form [
7]:
Here we assume that
, i.e., the donor is located at the middle along the Z-axis of the structure (depreciated surface states). In our numerical work, we solve the boundary value problem (
4) by using the double Fourier series expansion method, in which the wave function is represented as follows:
The functions
present an orthogonal and normalized basis, which satisfy the boundary conditions. Substituting (
5) in (
4) one can reduce the eigenvalue problem (
4) to the following secular equation:
3. Results
In this section, we present results of calculations for the energy spectrum of donors confined in an infinite barrier GaAs RUQW’s cylindrical surface of the height , and the homogeneous thickness , whose cross-section along the plane perpendicular to the axis Z is the finite Archimedean spiral, given in polar coordinates by the relation , with , and , . In addition, we show below results corresponding to the donor position with coordinates , , and for the direction of the applied electric field given by the polar angle .
Energies dependencies on the electric field presented in
Figure 4 demonstrate significant differences of the corresponding curves for the electron and the donor. In the contrast to the electron energy spectrum in
Figure 4a, which is quasi-continuous, the corresponding spectrum of the donor in
Figure 4c is splits in bands with decreasing gaps as the of band number grows typically for of Coulomb-like potentials spectra. Also, one can observe in
Figure 4a,b multiple crossovers of the curves corresponding to the lower energies, which originate a sharp jump of the slope ground state energy dependency for electric field value about 0.1 kV/cm.
Both for the electron and the donor gaps between energies inside the subbands corresponding to states of the electron localized in wings of the spiral, remote from the donor position are very small and therefore, it is more appropriate the description of the energy spectrum in terms of the density of the electronic states (DOS)
. In
Figure 5, we show curves of the DOS for different values of the electric field both for the electron and for the donor. In both cases one can observe a consecutive broadening of the principal peak of the DOS with increasing of the electric field. But for the donor, this broadening is almost symmetrical in contrast to the electron, for which the broadening of the peak is accompanied by its shifting to the red end of the spectrum. In addition, in the case of the electron the increase of the electric field induces the fine structure in the curve of the DOS.
One can also observe in
Figure 4 that ground state energies dependencies on the external electric field
in both cases as for the electron and for the donor have critical points, in which their slopes are changed abruptly. Such a feature in behavior of the curves can be related to a jump of the electron between two adjacent wings of the spiral induced by the external electric field, which should be accompanied by a sharp variation of the dipole moment of the donor
p, defined by the relation
.
In
Figure 6a we display the values of the averaged distance between electron and donor, characterizing the induced dipole moment of donor confined in the rolled-up QW calculated as a function of the external electric field. As the electric field increasing becomes larger than
the electron-donor separation changes its sign and it grows by jump. The explanation of this effect is given in
Figure 6b, where we present the density of the electron distribution before and after the jump. It is seen that before the jump the electron is located mainly at the same wing where is placed the donor, while after the jump it is displaced toward the adjacent wing and it is mainly located at opposite sides of the spiral.
Finally,
Figure 7, we present the energies of the donor of as a function of the angle
between electric field direction and X axis. It is seen that for a weak electric field as the donor dipole moment is small, the lower energies in
Figure 7a almost does not depend on the electric field orientation, because these state are strongly bound to thed donor. The behaviour of the lower levels in
Figure 7b is essentially changed as the electric field is equal to 0.1 kV/cm and the electron-donor separaion is large according to the results presented in
Figure 6a. The energies of the ground and of first excited states become sensitive to the orientation of the electric field in contrasst to the corresponding dependencies in
Figure 7a.
4. Conclusions
Low-dimensional systems based on semiconductors have been the subject of interest for last decades. Quantum properties of these systems offer the possibility of analyzing the physical behavior in the nanoscale; and numerous possible applications. Among other uses, they can be used, for example, as an active surface in heterogeneous catalysis, for the growth of other low-dimensional materials. Donor impurities confined inside a thin semiconductor film of nanometer scale has a spectral properties similar to an artificial two-dimensional hydrogen-like atom studied previously in detail for thin plane films. However, a greater diversity of spectral properties could be expected from donors confined in thin films with curvature. In a particular case of a thin semiconductor film in a form of a spiral, the donor confined in such a structure cannot be no longer considered nor as a two-dimensional hydrogen-like atom neither as a three-dimensional, due to long-range interaction between electron and donor located at different adjacent wings of the spiral. In this work, we propose a simple method, which allows us to analyze effect of impurity doping on the spectral properties of a spirally rolled-up quantum wells in the presence of the external electric field, applied at a direction perpendicular to spiral axis. Considering as an example an Archimedean shape of the spiral, we have investigated theoretically the effect of the curvature and the location of a single shallow donor, the intensity and the orientation of the external electric field on the donor energies. We present novel curves of the electron and donor energies dependencies on the intensity y orientation of the electric field and of the density of the corresponding states. We have found a high sensitivity of all curves of the energies dependencies both for the electron and the donor on the external electric field, which is characterized by abruptly change the dipole moment und the polarizability of the structure, related to the electron translations between adjacent wings of the spiral. Our method should be applicable to a variety of complex quasi-two-dimensional structures with curvature for which more rigorous approaches requires extensive numerical calculations.