1. Introduction
The modern technological state of the art permits the fabrication of materials endowed with new and astounding properties, able to disclose physical effects not seen before. Among these new materials, reduced-dimensionality nano and mesoscopic objects are gaining increasing interest in fundamental as well as applied research. They are characterised by the feature that at least one of their spatial sizes is much smaller that the others; the energy required to populate the levels pertaining to the relatively small dimension is much larger than that required for the other dimensions so that the small dimension can be disregarded in the mathematical description of the physical object; therefore, the system can be safely studied by means of a lower number of degrees of freedom. The conceptual interest in the possibility of experimentally exploring systems with a smaller number of degrees of freedom is evident but the technological repercussions have been shown to be of great importance too.
Among these new physical objects, one can include quantum dots, nano and meso rings and nano tubes, all deserving special consideration. Quantum dots are small clumps of semiconductor, confined within a convenient substrate and possessing a spectrum of discrete energy levels that can be designed almost at will to fulfil particular requirements and for this characteristics are pictorially referred to as artificial atoms; they can be seen as zero dimensionality systems. Particular mention must be devoted to the rich variety of allotropes of the carbon atom as some of them can be used to produce reduced dimensionality systems. For example, C60 fullerene is a spherical arrangement of carbon atoms, situated at the corners of hexagons and pentagons, that can be described as a zero- or a two-dimensional object.
On the other hand, graphene is a flat monolayer of carbon atoms placed at the corners of a pattern of hexagons of side
, with
the Bohr radius, and resembling a honeycomb; it is one of the most interesting materials for fundamental and applied research. At energy lower than 2 eV, owing to the particular symmetry of the lattice, the valence and conduction bands slightly overlap [
1,
2], producing a linear dependence of the electron energy,
E, upon the modulus of the kinetic momentum,
:
where
is the Fermi speed of the electrons in the layer,
c the speed of the light and
a dimensionless numerical constant. The linear dispersion is typical of the photon, which allows the description of the charge carriers as massless, free fermions moving with speed
[
1,
3,
4] and gives the opportunity of studying zero-mass charges which, nevertheless, can be acted upon by electromagnetic forces; in this way, a new phenomenon with no classical counterpart is open to investigation. Even if
with
, the average electron speed in the hydrogen ground state, one can see that the description of the charge carriers by using the Schrödinger equation is impossible and the Dirac equation must be used after the substitution
.
In the Dirac electron theory, the dual presence of a positive–negative energy solution, or electron–positron pair, is inherent and the Klein paradox makes the confinement of relativistic fermions by means of a repulsive potential barrier larger than
quite complicated, it would act as a chasm for the negative energy solution of the Dirac equation. Actually, in the presence of a large enough barrier, the positive energy solution of the Dirac equation disappears from the classical permitted zone of the space and makes a transition into the negative energy state solution by appearing in the classical forbidden zone of the space [
5].
In graphene, the condition of massless fermions is valid at low energy and for an infinite sheet, although the latter condition can be safely overcome by using quite large fragments of film. In spite of the finding that the formalism is quite similar, in graphene, actually, there is no room for the electron–positron interpretation of the positive and negative energy solution and the two solutions are to be interpreted as states in the valence and conduction bands. Nevertheless, the similarity of the theory makes a phenomenon similar to the Klein paradox present also in graphene, where the electron, moving in the conduction band of the film, in impinging upon a potential barrier, appears in the region of the step in the valence band [
6]. Thus, several confinement schemes, based on magnetised pieces of graphene [
2,
7,
8,
9,
10,
11], or on the presence of electric and magnetic fields [
12,
13,
14,
15,
16], have been devised.
Graphene is the strongest known material [
17] and endowed with a quite a large charge mobility [
18]: because of these and other properties, it is one of the most promising materials for applicative research in Medicine, Chemistry, Informatics and Physics. Among the many applications, we mention the use in efficient solar cells [
19], which could help the increasing demand for electrical energy with a quite low environmental impact, and as a frequency multiplier of laser light in the well-known high harmonic generation [
20,
21,
22]. Considering graphene, most interesting is the finding that, by suitable manipulation, it is possible to obtain, from the basic honeycomb structure, different systems possessing different dimensionality. Actually, graphene sheets can be rolled to form quite long and slender cylinders, dubbed nano tubes, and cones with lengths several orders of magnitude larger than the radius and, by cutting from tubes, nano rings can also be fabricated with unidimensional properties. Also by cutting a graphene sheet, two dimensional annuli and disks can be prepared. Another fascinating development is the realisation of elaborate three-dimensional graphene structures inspired by origami [
23,
24,
25].
Today, computer technology prompts the design and fabrication of meso- and nano-sized electronic components for information storage and as the fastest logic circuital elements [
26,
27,
28,
29]. However, saturation in the performance of traditional silicon devices is presently observed so that the quest for new material is active; in this context, graphene and other carbon allotropes, such as C
60, fullerene and carbon tubes, are natural candidates for their quick response to external fields and high-frequency operation. Therefore, the comprehension of the behaviour of carbon-based material as C
60 molecules, quantum rings and annuli acted upon by static or dynamic electromagnetic fields is of paramount importance. Far from being of hindrance to the theoretical approach, the quite large size of these molecules may create symmetries which, judiciously exploited, produce deep simplification in the theoretical approach to the problem and often leads to analytical solutions or crucial insight into the physics of the problem [
30].
This paper deals with a disk of graphene of radius
crossed by a static magnetic field
. The quite large radius
of the disk makes the dispersion law in Equation (
1) and the massless condition of the charge carriers suitable approximations. Moreover, the large size makes the inevitable irregularity and roughness at the disk border insignificant as it is well known that armchair or zigzag edges give different spectra in graphene ribbons [
1,
31]. We analytically solve the Dirac equation for the charge carriers by finding the eigenstates and relative energy as a function of the magnetic field and numerically evaluate them. As a result, we find that the magnetic field has a strong confinement action on the electron and that relatively large orbital angular momentum states are supported by the disk. These states are confined in different regions and might be used to store several bits of information; since large angular momenta states of Bose–Einstein condensate, vortices, confined in spherical surfaces and rings, are currently actively investigated [
32,
33,
34,
35], graphene disks can be suitable candidates for studying vortices states with simple enough equipment.
2. Theory
Several equivalent representations of the Dirac equation can be given that are related by a unitary transformation. The commonest forms are the standard representation that is mostly useful in weak relativistic approximations and the chiral representation which is mostly useful in ultrarelativistic regimes, when the energy of the particle is much larger than the rest energy. For reasons that will soon be clear, in what follows, the chiral representation is used. Thus, the Dirac equation for a fermion of charge
q and mass
in the presence of a magnetic field in the chiral representation is
where
is the triplet of the Pauli matrices,
the canonical momentum operator (denoted by a tilde below the symbol),
ℏ is the reduced Planck constant,
the vector potential,
,
t denotes the time, and
and
are the time-dependent left and right two-component Weyl spinors (denoted by a dot below the symbol), respectivly, forming the four component Dirac spinor in the chiral representation.
and
depict the positive and negative energy states of the fermion. Care must be taken in handling
because, in spite of the look, under a spatial rotation it does not transform like a vector; it must be regarded as a convenient notational device that accounts for the mixing of the spinor components. From Equation (
2), one sees that the mass couples the left and right spinors. In our model, however, the mass and the related rest energy are zero so that the equations for the two components uncouple:
and differ only for a sign; thus, the positive and negative energy solutions behave independently. For the sake of brevity, only the first equation, with the superscript
dropped out, is studied here. For the Dirac charge (after the default substitution
),
assumes the form of a time-dependent Schrödinger equation,
with the Hamiltonian,
One recognises in the operator
the kinetic momentum operator. The Hamiltonian is related to the helicity operator. As mentioned in
Section 1, in general, for studying the properties of ultrarelativistic electrons, Equation (
3) is used approximately. In the case of graphene, for the vanishing of the mass, this use is necessary and offers the opportunity to investigate relativistic effects without making use of large accelerators.
2.1. Total Angular Momentum
In this paper, we consider that the static magnetic field
generated by the vector potential,
pierces a disk of radius
whose centre
lays at the origin of the
-plane,
. The particular geometry of the system suggests that the third component of the total angular momentum is a constant and that to seek common eigenstates of the Hamiltonian and this component of the angular momentum operators is convenient. The calculations to find the expressions for the operators and to obtain their equation of motion require quite an amount of algebra and may be rather long. In what follows, we confine ourselves to outlining the significant steps and the interested reader may find details in [
5,
36].
Thus, after long but straightforward use of commutator algebra, the equation
for the total angular momentum
(in units of
ℏ) is arrived at. One can see that the time derivative of
is orthogonal to the magnetic field, so that the component of the total angular momentum along the magnetic field,
is constant and admits common eigenstates with the Hamiltonian. Here,
is the polar angle. Moreover, as discussed in
Section 1, in materials with reduced dimensionality, the active electron can be described as confined on the disk to be considered as a two-dimensional object.
We begin by considering the eigenstates of the angular momentum. In circular variables, the eigenvalue equation is
, with
the radial variable and
the angular momentum eigenvalue. The two components of the time-independent spinor
correspond to the two projections of the angular momentum along the magnetic field. After splitting
into its components by setting
and sustituting into Equation (
9), one obtains
equivalent to the two coupled differential equations for the angular part:
with the solution
where
m is the quantum number of the orbital angular momentum,
Let us stress that and not m is the conserved quantity. Actually, it is a straightforward task to check, by direct calculations, that is not eigenstate of .
In order to obtain the form of the equation for the radial parts, many tedious steps are required [
36]; the result is
with the prime denoting differentiation with respect to
. Surely, the subscript
m had to be added to the wave function. Actually, notational rigorousness would require that the same index should flag the radial wave functions
and
but symbolism would become heavy and cumbersome and requires refraining from full coherence.
The eigenenergy equation is
where the energy
E, common to both the spin components, and the radial functions
and
must be found simultaneously. The two coupled first-order equations are not of immediate solution and substitution of
from the second into the first equation is convenient. The procedure leads to a second-order differential equation and is similar to the one that leads from the first-order Maxwell equation to the second-order differential equation for the potentials. The result is
To work with a dimensionless variable, it is convenient to define
where
has the dimension of length square. Thus, the parameter
introduces in the problem a length scale, inversely proportional to
, and measures the confinement range induced by the field. In terms of the new variable (
18), Equation (
17) for
thus reads
where
with
the sign of the charge.
The substitution,
leads to the known confluent hypergeometric equation (CHE):
The standard form for the CHE is
admitting two independent solutions,
and
, which can be written as infinite series. The expression of
gives few problems. The expansion of
which is well behaved (as a limit) for integer
b [
37], is not straightforwardly numerically approached because of the presence of the
and of the
function of negative integer that can lead to divergences, the latter are quite difficult to handle. Indeed, the one of interest here is located in the hard subset of the problem because, in our model,
b is always an integer.
Looking at the expression (
23), one determines the parameters of the confluent hypargeometric functions (CHFs) in terms of the relevant physical quantities:
Again, we stress that the feature that b is a positive integer be carefully addressed during the numerical evaluation of .
To summarize, the radial function
has two independent solutions:
and
which correspond to the spin-up component of the spinor.
The analytic evaluation of the radial functions of the spin down is possible starting from these expressions but requires a lot of algebra of the CHF and is detailed in
Appendix A. The final result is
and
In
and
, the presence of the term
is rather interesting as it acts in different ways for positive and negative
m. Actually, the study of a positive charge (
) and positive
m would flow into a tremendous simplification of the formulas as one piece of the wave functions would vanish. The general radial wave functions are given by a linear combination of these functions:
where
and
are constants to be found from the boundary conditions and normalisation condition.
For a fixed value of m, one observes that a can be maintained constant by keeping and this would give the dependence of the eigenenergies from the magnetic field; this is, nevertheless, of little use as is dependent upon .
2.2. Boundary Conditions
In this paper, we consider a disk with a sharp edge. Actually, no real fragment of graphene can have a sharp border and studies of graphene ribbons have shown that the form of the edges is relevant for the determination of the states of the charges [
31]. Again, this approximation can be made weak by considering a disk with a large enough radius
so that roughness is not relevant.
The fulfilment of the boundary conditions permits the determination of the eigenstates and eigenenergies (the eigenvalues of the total angular momentum were already found in Equation (
14))
To avoid long final expressions, let us introduce the constants
and
so that the boundary conditions assume the form of a two-variable, homogeneous linear algebraic equation:
where the unknown coefficients
and
may be determined only if the determinant of the coefficients vanishes:
Remembering that
E is present in the expression of
a, one can see that the values of
E that fulfil the condition in Equation (
40) are the allowed values of the energies. The problem is therefore fully analytically solved:
The value of
is found by the normalisation of the wave functions:
2.3. Helicity
As it quite straightforward to see by direct calculation, the Hamiltonian eigenfunctions are not eigenstates of
and
and the spin projections on the
-are not individually constant but the helicity is constant. The helicity operator
can be defined as
so that the helicity evaluated on the energy eigenstates is constant. Generally, the helicity gives the projection of the spin on the kinetic momentum. In our model, both
and
lie on the
-plane (see Equation (
7)), so only the projection of the spin in this plane is relevant. The finding that in the present model of disk, the Hamiltonian
shares the eigenstates with the helicity is helpful feature as helicity is the only degree of freedom for massless particles. One can expect that, the orbital angular momentum is associated with a rotation around an axis; the fact that here the eigenstates of the energy are not eigenstates of
and that
m is not a proper quantum number, hinders the possibility of giving a transparent meaning to vortex states.
3. Results and Conclusions
In what follows, we adopt atomic units (au) (
); thus, the Bohr radius
is the unit of length, the unit of energy is 27.21 eV and
au
T. The radius of the disk is
= 10,000
, the large value ensuring that the electron can be treated as massless; the magnetic field is taken in the interval
au. To find the values of
E that solve Equation (
40) is a complicated enough task because, as comes from their power expansion, the CHFs are especially sensitive to the value of their parameters; nevertheless, we have checked our results against numerical precision and found no mistakes.
From the expression of the parameters and of the wave functions, we observe that the energy always appears through the ratio
and therefore we use this ratio as the relevant parameter and, generally, call it energy. Through the parameter
a in Equation (
27), the energy enters the equations with the ratio
, showing an intrinsic weak dependence of the eigenvalues upon
. The parameters of the CHF are quite sensitive to the value of the magnetic field
and this dependence is rather crucial at low enough
. Thus, actually, the results are stabler and more reliable at large enough
. For any value of
m, there are several eigenstates fulfilling the boundary conditions and these eiegenstates are numbered with the integer
n; thus, to label the eigenenergies and, consequently the relative eigenstates, the second subscript
n must be introduced as
. For
,
n counts the number of nodes of the eigenstates between the centre and the disk edge. For
, such an appealing correspondence is not realised and the number of nodes can be
. We have performed numerical calculations to determine the eigenvalues and eigenstates of the Dirac electron in the disk for different values of
,
m and
n; in all the calculations, we have considered a negative charge:
.
Figure 1A shows the eigenengies of the Dirac electron versus the magnetic field for
and for different values of the orbital angular momentum,
m. For non-negative
m,
increases with
and
m and observe the mentioned weak dependence of
on
. Relegated in the lower left corner of the figure, the curve for
pops up with its considerably small energy and decreasing behaviour.
Figure 1B shows the eigenenergies of the electron versus
n at
au for different value of the orbital angular momentum. The eigenvalues increase almost linearly with
n and can be extrapolated with confidence.
The meaning of the integer
n can be understood from
Figure 2 where we show few eigenstates at
and different values of
n.
Figure 2A displays the case of spin up. The states peak always at the centre of the disk and present
n nodes (the node at the disk edge is there by construction). The case for the spin down is shown in panel
Figure 2B. Now, the states have a node at the origin and present again
n nodes between border and edge. The presence of the nodes gives an expected explanation of the growth in the energy as shown in
Figure 1B since the nodes have the effect of decreasing the available area of the disk, thus confining the charge in a smaller area.
Actually, the description is not so clear-cut for negative
m. In
Figure 3, the case for
is shown. For spin up,
n again gives the number of nodes between the centre and edge of the disk. In the spin-down case, on the other hand, the eigenstates show
nodes. To stress is that states with the same
n (in
Figure 3A,B) correspond to the same eigenvalue.
Figure 2 suggests a possibility of obtaining a spatial separation of the two-spin polarization. For example, for
, the spin-up eigenstates, in
Figure 2A, peak at the centre of the disk while the spin-down states, in
Figure 2B, have a node at the same point. For the case
, shown in
Figure 3, the situation is reversed and the spin-down states peak at the centre. However, to stress is that the relevant quantity is
, (see Equation (
43)), which emphasizes the mid to the detriment of the innermost (where
) and peripheral (where
) parts of the disk. For example, in
Figure 4, we show
(
) for
au,
and
whose relative eigenstates are shown in
Figure 2. However, still from
Figure 4, one can see that a certain amount of external control for the location of the peak of the two-spin polarization is possible by using the quantum numbers
n and
m as knobs.
In
Figure 5, we show the eigenstates of the Dirac electron versus
for
au and for several values of the orbital angular momentum
for both the cases with spin up and down. In all the curves,
and thus the wave function shows no nodes. A few interesting features emerge. As to be expected, states with a quite large value of the orbital angular momentum are pushed toward the outer region of the disk as if a centrifugal force is acting on the particle even if it is massless. Recently the study of trapped Bose–Einstein condensates (BEC) led to the discussion of quantised vortices [
32] which generally can be defined as states possessing a well-defined orbital quantum number,
. Such states and their control to design devices for reading and writing in computers are today an extensive research activity [
38]. Actually, the possibility of stacking or arranging rings and disks in a coplanar geometry to let them communicate via their magnetic field is appealing [
39]; from this point of view, the use of graphene annuli and disks is particularly attractive because modern technology permits the preparation of multi-layered graphene sheets.
One can associate these concept of energetically rotating states possessing a relatively large orbital angular momentum to a vortex. In relativistic quantum theory for a free particle, the orbital and the spin angular momenta are not individually constant, so we are led to interpret the state as a vortex having vorticity number. Moreover, states with large and small enough values of m overlap quite little and a way to associate different bits of information to them should be envisaged.
Figure 6 clarifies the role of the nodes of the eigenstates for negative
m with the curves taken for
and the states have different number of nodes in spite of their same energy. Summarizing, for
, the degenerate eigenstates can be labelled by the two quantum numbers,
m and
n, with
n counting the number of nodes of the eigenstates; for
, on the other hand,
n can only set the order of the states.
In
Section 1, we mentioned the difficulty, encountered in relativistic quantum theory, of confining electrons by means of a potential barrier.
This difficulty arises from the point that the negative energy states see the barrier as a well and are, then, pulled to climb the cliff. In graphene, the Klein paradox is present but appears as a transition between conduction and valence band of the lattice; however, confinement can be obtained by means of a magnetic field. The curves in
Figure 7 show few states for
and spin up, for sake of simplicity. It is clear enough that the magnetic field has a relevant importance in confining the charges; In general, for a relatively large magnetic field, the wave function does not approach the edge of the disk. Here, one can see how the Klein paradox may be controlled. The paradox lies on the basis that a positive energy solution of the Dirac equation (the electron in the conduction band) cannot be confined by a repulsive potential since in reaching the potential it makes a transition toward a negative energy solution (the electron in the valence band). However, the magnetic field forbids the approach of the particle to the border of the disk and the consequent interaction with the potential; we consider this one of the relevant results of this paper. It is essentially to stress that the necessary magnetic field is rather large but the physical situation is not ordinary in describing massless charges; the use of a smaller disk is even detrimental for the confinement as a larger field is required. Actually, the energy of the electrons increases by decreasing the radius
and this outcome can be exploited if wanted. As a side effect, it can be noticed that a stretching force is exerted by the charge on the disk [
36] wich increases with the confinement and, consequently, with the magnetic field. Such a stretching may become detectable in the near future.
In conclusion, in this paper, we study the wave function of the charge carriers in a considerably large disk of graphene in the presence of a static magnetic field orthogonal to the plane of the disk. Thanks to the particular symmetry of the lattice, the valence electrons of the carbon atoms can be treated as having null mass, like the photon, but moving with a speed much slower that the average speed of the electron in the hydrogen ground state. In spite of this observation, because of the massless condition, the charges can be described only by means of the relativistic Dirac theory with its intrinsic positive–negative energy dualism. Thus, the graphene provides a handy way to observe massless particles that can be acted upon by electromagnetic fields.
The analysis shows that the magnetic field acts as a confining tool of the charges and, perhaps, it may provide the only confinement parameter to defeat the Klein paradox. The eigenstates and related energy can be labelled by the integer quantum numbers
m and
n. The key result is that the position of the peak of the eigenstates can be manipulated using the quantum numbers
m and
n as control knobs and can be made to range from the inner to the outer part of the disk. In this way, by linear superposition it is possible, in a graphene disk, to construct states with different peaks. By associating a bit of information to any eigenstate of the superposition, it is possible to store many bits of information in the disk. Moreover, the disk can support states with high values of the orbital angular momentum
m, thus providing a way to observe states with large enough vortex charge. In this perspective, annuli and disks of graphene have been extensively studied both from the analytical and the fully numerical point of view to investigate their use as quantum dots by varying the edge confinement potential [
12,
13,
15,
40].
Recapitulating, in this paper, we obtained the wave function of the states of the electrons in a magnetised disk of graphene in strict analytical form both for the up and down spin components and labelled them according to two quantum numbers for the orbital angular momentum and for the nodes of the wave functions. These quantum numbers can be used as a knobs for controlling the position of the peaks of the states in the disk in the context of storing bits of information. We showed that the system can support states with quite large angular momentum and suggested that the relativistic treatment does not permit the definition of vortex states as endowed with a precise orbital angular momentum. However, the total angular momentum is a proper quantum number and allows a redefinition of the vortex state. A key new contribution is the evaluation of the states for several values of the magnetic field and proved that confinement may be achieved by means of a magetic field. All of these characteristics make the graphene fragments a quite particular object for basic and applicative research.