Next Article in Journal
Forty-Nanometer Plasmonic Lithography Resolution with Two-Stage Bowtie Lens
Next Article in Special Issue
New Insights into Improving the Photovoltaic Performance of Dye-Sensitized Solar Cells by Removing Platinum from the Counter Electrode Using a Graphene-MoS2 Composite or Hybrid
Previous Article in Journal
Anchor Loss Reduction in Micro-Electro Mechanical Systems Flexural Beam Resonators Using Trench Hole Array Reflectors
Previous Article in Special Issue
Structural Control of Nanofibers According to Electrospinning Process Conditions and Their Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electric Field Effects on Curved Graphene Quantum Dots

by
Sergio de-la-Huerta-Sainz
1,
Angel Ballesteros
1 and
Nicolás A. Cordero
1,2,3,*
1
Physics Department, Universidad de Burgos, 09001 Burgos, Spain
2
International Research Center in Critical Raw Materials for Advanced Industrial Technologies (ICCRAM), Unversidad de Burgos, 09001 Burgos, Spain
3
Institute Carlos I for Theoretical and Computational Physics (IC1), 18016 Granada, Spain
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(11), 2035; https://doi.org/10.3390/mi14112035
Submission received: 1 September 2023 / Revised: 28 October 2023 / Accepted: 29 October 2023 / Published: 31 October 2023

Abstract

:
The recent and continuous research on graphene-based systems has opened their usage to a wide range of applications due to their exotic properties. In this paper, we have studied the effects of an electric field on curved graphene nanoflakes, employing the Density Functional Theory. Both mechanical and electronic analyses of the system have been made through its curvature energy, dipolar moment, and quantum regeneration times, with the intensity and direction of a perpendicular electric field and flake curvature as parameters. A stabilisation of non-planar geometries has been observed, as well as opposite behaviours for both classical and revival times with respect to the direction of the external field. Our results show that it is possible to modify regeneration times using curvature and electric fields at the same time. This fine control in regeneration times could allow for the study of new phenomena on graphene.

1. Introduction

There is no doubt that there has been a long-lasting wave of advancements and innovation revolving around graphene, with strange phenomena and potential applications appearing on a yearly basis. Almost every field seems to have an application where the extraordinary properties of graphene (either pristine or oxidised) can shine, from engineering [1,2,3,4] to catalysis [5,6,7,8], medicine [9,10,11,12,13], sensing [14,15,16], optics/electromagnetics [17,18,19,20], hydraulics [21,22], and energy management [23,24,25], to cite a few examples.
In spite of the mechanical properties of graphene being quite remarkable, the electronic ones are by far the main factor responsible for its attractiveness as a research topic and the basis for several groundbreaking applications. The peculiar band structure of graphene, calculated by Wallace many decades ago [26] by an approximated yet still useful method, presents the electrons as massless Dirac quasiparticles [27], allowing them to reach high speeds [28] that result in its incredible electrical conductivity [29,30]. This, alongside the ballistic transport [31,32] observed on non-perturbed systems, makes graphene a great replacement for common conductive materials. Also, due to its near-relativistic behaviour, graphene is often used as an experimental analogue in order to replicate otherwise nearly impossible to measure phenomena, like the Klein paradox [33,34,35], Zitterbewegung [36,37,38,39], or the properties of spacetime with a negative curvature [40,41,42,43].
Perhaps one of its most interesting applications derives from the interaction between graphene and an external electric field. As the band structure is so important for its electronic properties, it is quickly deduced than they can be modified and tuned easily by adjusting the field, allowing feats such as moving the Fermi energy between bands and changing the type and density of charge carriers—the ambipolar electric field effect [44,45]—or widening the gap between valence and conduction bands to obtain an insulating behaviour, if needed. Taking also into account the high mobility of electrons in this material and its nanometric scale, graphene is a very promising alternative material for fast and incredibly small electronic components, such as field-effect transistors [46]; the ballistic transport would ensure a quick response and reduced heating, and the robustness and transparency of the material could extend its use to flexible devices and screens [47,48,49].
A known effect of an electric field on graphene is the appearance of the curvature. When a transversal field is applied, the material bends alongside the direction of the field, as has been both predicted and observed recently [50,51,52,53]. While curved graphene is not a strange view, as it naturally appears in the material itself even when suspended and unperturbed [29,54,55,56,57,58], the effects of the deformation on the electronic behaviour are notable. Intrinsic ripples are known to disturb the otherwise free electron paths [27], and local deformations, known as nanobubbles, can be easily produced in the material, allowing for a localised tuning of electric conductivity and more exotic phenomena, like pseudo-magnetic field generation [59,60]. It has been recently shown that double quantum dots can be created by these pseudo-magnetic fields in nanobubbles and that their quantum states can be manipulated to create a controllable qubit [61].
Thus, our aim in this paper about graphene nanoflakes is to study the effects of an electric field on a spherically deformed system, analysing the electronic spectrum via temporal evolution through the quantum regeneration phenomenon—the partial or total recovery of the initial state of a wavepacket after a certain amount of time, known as the regeneration or revival time. An atomistic method, like the Density Functional Theory (DFT), will serve as the main tool, as it has previously for nanoflakes with spherical geometries [62], as well as those with cylindrical and hyperbolical geometries [63], in the absence of an external field. Some of those results will be used here for comparative purposes.

2. Materials and Methods

As this work is a continuation of our previous research about quantum revivals in graphene, the methodology, as well as the system, will remain the same: a hexagonal graphene nanoflake of 10 benzenic rings per edge, with zig-zag edges and hydrogen passivated; the shape will be achieved through the application of the spherical surface equation, Equation (1), for several values of R between 30 Å and 1000 Å:
z = R 2 x 2 y 2 .
This flake will be subjected to a structural optimisation calculated within the Local Density Approximation (LDA) using the Gaussian 16 package [64] until forces on all atoms are below 0.0045 Hartree/Bohr and displacements below 0.0018 Bohr. While other common functionals, such as B3LYP or Generalised Gradient Approximations (GGAs), can be used for smaller systems, LDA results for graphene systems are of a similar or higher quality [65,66], they require lower computational work, and this functional has been successfully used to study carbon nanostructures [67,68,69,70,71,72,73]. Meta-GGAs provide better accuracy than LDA but the computing time is much higher. Taking into account that some of our calculations needed 1 year of CPU time at the LDA level, using those functionals was out of the question. Therefore, we followed the advice in a recent benchmark paper on ab initio descriptions of carbon nanomaterials using 15 different functionals [74] and used LDA. The basis set will remain 6-31G** (with p and d functions as polarisation aids) [75], as it has provided good quality results on all calculations made so far.
Highly spherical nanobubbles can be achieved with external electric fields [51], so a perfect spherical deformation could be enough for this study; however, for the sake of experimental realism and fair comparisons, two different sets of geometric restrictions will be applied to the nanoflake (Figure 1): (i) maximal restriction, forcing all carbon atoms to lie on the spherical surface and (ii) minimal restriction, forcing only the 12 carbon atoms on the vertices to remain on the spherical surface. The first case may seem purely academic, but it is a first approximation to a possible experimental setup. One method to stabilise graphene for electronic applications is using hexagonal boron nitride (h-BN) as a substrate [76,77,78]. A spherical boron nitride fullerene [79,80,81] could be used to support a spherical graphene nanoflake.
Figure 1. Top (upper row) and side (lower row) views of the two boundary conditions for a spherical nanoflake with initial R = 40 Å . Clear differences between both optimised structures can be seen, especially on the edges. The atoms fixed in the second case are highlighted in blue. Images generated with Gaussview 6 [82].
Figure 1. Top (upper row) and side (lower row) views of the two boundary conditions for a spherical nanoflake with initial R = 40 Å . Clear differences between both optimised structures can be seen, especially on the edges. The atoms fixed in the second case are highlighted in blue. Images generated with Gaussview 6 [82].
Micromachines 14 02035 g001
Regarding the electric field, considering the original flat nanoflake was constructed in the xy plane, the fields applied will be in the direction of the z axis—and thus perpendicular to the unperturbed flake plane—with three intensities (0.0050 a.u., 0.0100 a.u., and 0.0250 a.u.) in two opposite directions (one in the same direction as the flake deformation and the other in the opposite direction) each, labelled as ± 50 , ± 100 , and ± 250 , with the positive sign corresponding to the field in the same direction as the flake deformation. These fields could seem to have an astonishing intensity, considering that a field of 1 atomic unit equals 5 × 10 11   V m 1 , but at the nanometric scale such fields can be achieved by setting the electrodes close enough.

3. Results

The study of the effect of an electric field has been carried out from both mechanical and electronic points of view, the first through the curvature energy and the latter by calculating the temporal evolution of a wavepacket and checking if quantum regeneration occurs. To this end, spherical surfaces of various curvature radii have been used to deform the flake, performing an optimisation on each combination of radius, restriction set, and electric field.

3.1. Curvature Energy

As was the case in our previous works, the mechanical study will use the total energy of each nanoflake, taking the flat configuration as the origin of energies, i.e.,
E = E curved E flat ,
where E curved is the energy of the curved flake and E flat is the energy of the planar configuration (both energies calculated with the same electric field).
Figure 2 shows the curvature energy plotted against 1 / R 2 —which is the Gaussian curvature of the flake for the maximal restriction case—for different field strengths. In order to make figures easier to understand, we have used point markers for results corresponding to maximal restriction and dashed lines for results corresponding to minimal restriction. In this latter case, the lines are merely guides for the eye connecting the results corresponding to consecutive values of 1 / R 2 (we have used the same set of values as for the maximal restriction case). While all curves show that, as expected, flakes with higher curvatures have higher energies, the trends vary depending of the relative direction of the field; for zero fields, the plot has a nearly parabolic shape, with the flat case being the lowest energy geometry; for positive fields, the plots are also parabolic, but the minimal energy is reached for a non-planar structure, which showcases the curvature induced on the material as a transversal field is applied; lastly, negative fields result in an initial rise in energy at lower curvatures, as the field is trying to induce a deformation opposite to the pre-existing one, an effect countered by bending as curvature grows.
Stronger fields lead to analogous behaviours, but the effects are bigger as the field increases. The optimal structure is more curved and less energetic for positive fields, while the initial energy increase appears earlier and is greater for negative fields. However, a strong enough negative field can override the curvature effect and produce a maximum in energy.

3.2. Electric Dipole Moment

Due to its close interaction with external fields, the behaviour of the electric dipole moment has also been studied. From the flake shape and symmetry, a null moment would be expected for the flat case and a non-zero moment on the negative z axis for the spherical ones; this moment would be intensified with positive fields, as the charge distribution that originated it becomes more pronounced and decreases or even flips with negative fields.
In Figure 3, we have plotted the z component of the dipole moment against 1 / R 2 for each of the seven field values. As expected, in general, the dipole along the z direction grows in absolute magnitude with an increasing curvature, since a bigger curvature allows for a greater charge displacement. In the absence of an external electric field (green points), the planar structure ( 1 / R 2 = 0 ) has no global dipole moment due to the symmetry of the electron cloud with respect to the graphene plane.
Figure 4 shows the dependence of one particular case with the field, showing an almost perfectly linear trend, as expected. Without an electric field, the dipole moment is not exactly zero because curvature causes the centre of charge of the electron cloud (negative) to not be at the exact same position as the centre of charge of the atomic nuclei (positive). The first is slightly below the second and this fact creates a very small downwards dipole. With a positive electric field, i.e., the field pointing upwards, atomic nuclei tend to move up (and with them the centre of positive charge) while the electronic cloud tends to move down (and with it the centre of negative charge). This translates into a larger negative dipole moment that increases with the strength of the field. The opposite happens when the applied electric field is negative, i.e., when the field points downwards, atomic nuclei tend to move down, while the electronic cloud tends to move up. This results in a larger positive dipole moment that increases with the field strength.
It is interesting to see the behaviour for negative fields in Figure 3, where in this case the dipole component initially decreases before growing. For positive fields, the charge shift induced by the field and the curvature work together to achieve a greater downwards dipole; for negative fields, they oppose, and as charges can move longer distances in these conditions—the centre and edges are farther away—the curvature achieves a totally opposite effect.
The direction of the dipolar moment suggests a higher electron density on the sphere inner region, something already studied on larger curved carbon nanostructures [30,83], in which the curvature reduces the effective distance between atoms and provokes a rehybridisation on them, with a higher overlap and a consequent energy decrease. While there are works studying the connection between the local charge distribution and graphene rippling [84], these are oriented to a more continuous and extended system; alas, it is not justified to assume a direct connection between our small, static system and a much bigger and more dynamic phenomenon like graphene corrugation.

3.3. Regeneration Times

Quantum regeneration phenomena will be studied using the electronic spectrum of each nanoflake as a starting point. For that purpose, we define a wavepacket as a linear combination of eigenfunctions u n with varying coefficients c n :
| Ψ ( 0 ) = n = 0 c n | u n .
We can calculate the state of the packet after an arbitrary time t in order to study its time evolution:
| Ψ ( t ) = n = 0 c n | u n e i ħ E n t ,
with E n being the eigenvalue corresponding to the | u n eigenstate. The weights of each eigenfunction will follow a Gaussian distribution defined by level numbers—positive integers for unoccupied states and negative for occupied ones—centred around the LUMO+4 ( n 0 = 5) and a width σ of 0.7, ensuring a wavepacket narrow enough to observe quantum regeneration:
c n = 1 σ π e ( n n 0 ) 2 2 σ 2 .
Using level numbers allows us to use an analytical approximation for the regeneration times, as described in [85]. This approximation begins with a Taylor expansion of the spectrum E ( n ) around the central level n 0 :
E n = E n 0 + E n 0 ( n n 0 ) + 1 2 ! E n 0 ( n n 0 ) 2 + 1 3 ! E n 0 ( n n 0 ) 3 + .
By combining Equations (4) and (6), we reach an expanded evolved state:
| Ψ ( t ) = n = 0 c n | u n e i ħ ( E n 0 + E n 0 ( n n 0 ) + 1 2 ! E n 0 ( n n 0 ) 2 + 1 3 ! E n 0 ( n n 0 ) 3 + ) t ,
in which each term on the exponential function constitutes a temporal scale. This expression allows for the definition of several regeneration times, directly connected to the spectrum derivatives:
T Zb = π ħ | E n 0 | ,
T Cl = 2 π ħ | E n 0 | ,
T Re = 2 π ħ | E n 0 | / 2 ! , and
T Sup = 2 π ħ | E n 0 | / 3 ! .
The main times of interest here will be T Cl , the classical time, and T Re , the revival time, as they will be easier to observe and rationalise. Being dependent on the first and second derivatives, respectively, the classical time gives information about the absolute value of the energy differences and thus the compression of the spectrum, while the revival time informs about the similarity between energy differences and thus about the regularity of the spectrum.
T Zb is related to Zitterbewegung, a relativistic phenomenon only observed on particular wavepackets (cat states), different from ours. Here, it is just a mere phase factor. In contrast, T Sup , the super-revival time, would correspond to a higher scale oscillation on the temporal evolution, which is difficult to detect due to the decreasing value of higher-order derivatives and as such will not be considered either.
The interpolating function selected to calculate the derivatives appearing in these times is a second-degree polynomial, obtained by a least squares fit to the three central levels of the wavepacket. By using three points—the level numbers and their energies—an exact expression for the polynomial can be obtained:
E ( n ) = b a 2 ( n 2 n 0 2 ) + 1 2 ( a + b + 2 a n 0 2 b n 0 ) n + E n 0 ,
with n 0 being the central orbital number, E n 0 its energy, and a and b the absolute energy differences with the lower ( n 0 1 ) and upper ( n 0 + 1 ) levels, respectively. Both first and second derivatives can be analytically computed easily:
E ( n ) = ( b a ) n + 1 2 ( a + b + 2 a n 0 2 b n 0 )
and
E ( n ) = b a .
These expressions will be useful in later sections.
As for the actual temporal evolution, it can be monitored through the autocorrelation function, defined as the overlap of initial ( t = 0 ) and later ( t > 0 ) states:
A ( t ) = Ψ ( 0 ) | Ψ ( t ) = m = 0 n = 0 c m * c n u m | u n e i ħ E n t
Working with orthonormal states in the wavepacket makes it possible to substitute Kronecker’s delta functions for the overlap integrals and cancel all cases where n m :
A ( t ) = m = 0 n = 0 c m * c n e i ħ E n t δ m n = n = 0 | c n | 2 e i ħ E n t .
The physically sound quantity is the squared modulus of the autocorrelation function:
| A ( t ) | 2 = m = 0 n = 0 | c m | 2 | c n | 2 e i ħ ( E n E m ) t .
Taking into account the definition of the cosine function, we can rewrite this expression as follows:
| A ( t ) | 2 = n = 0 | c n | 4 + n = 0 m < n 2 | c m | 2 | c n | 2 cos ( ω m n t ) ,
where ω m n = ( E m E n ) / ħ . This formula relates the temporal evolution to the energy differences between different levels on the packet, rather than the energies themselves. As a sum of cosines, a complete revival is guaranteed (if the spectrum is commensurable) after a long enough time, when the least common multiple of all individual periods is reached.
We can focus on our particular wavepacket to further simplify this expression. While is it true that the wavepacket comprises five levels, the coefficient distribution and their squares on the general formula makes the three central ones (4, 5, and 6) the main contributors, so the general expression can be simplified. We can ignore any terms involving other levels. We can go even further and remove also the ω 46 term, as the distribution is narrow enough to make its contribution barely noticeable, resulting in the following expression:
| A ( t ) | 2 2 | c 4 | 2 | c 5 | 2 cos a ħ t + 2 | c 5 | 2 | c 6 | 2 cos b ħ t ,
where we have removed the first sumatory and proportionality constants, as the periodicity of the function remains the same, and a and b are, like before, the energy differences between level 5 and levels 4 and 6, respectively. Being a sum of cosines, and considering c 4 and c 6 are identical due to the distribution’s symmetry, it can be easily converted to a product of cosines:
| A ( t ) | 2 4 | c 4 | 2 | c 5 | 2 cos a + b 2 ħ t cos a b 2 ħ t .
The two new arguments found here are interesting: the first one is the average of both main energy differences, and the second gives information about how similar they are, which makes them perfect candidates for the origin of both classical and revival times, respectively.
Finally, Figure 5 shows a typical view of | A ( t ) | 2 plotted against time, with a clear pattern of a high-frequency oscillation modulated by a low-frequency one, whose periods correspond to the values of T Cl and T Re , respectively; the former is located at the first maximum of the time evolution and the latter at the first maximum of the enveloping curve. This corresponds to the shape or a general sum (or product) of cosines, corroborating our previous mathematical explanation.
In general, the observation of both times is an easy task, but it is not guaranteed: as both times approximate in value, the interference between them increases, resulting in a larger shift of the first maximum, quickly deviating from the T Cl analytical value, making the visualisation of T Re harder, if not impossible.
Summarising, there are two different ways to calculate regeneration times. The first one is looking for the maxima in the modulus of the autocorrelation function by numerically searching for the first maximum of the function (for T Cl ) or the first maximum of its enveloping curve (for T Re ). The second one is using the analytical expressions given in Equations (9) and (10).

3.3.1. Classical Time

Numerical and analytical values of classical time are plotted against 1 / R 2 in Figure 6 for the fixed surface case and in Figure 7 for the fixed vertices case. Taking the fixed surface as an ideal situation, we can see that, without any field applied, both numerical and analytical values of T Cl grow with 1 / R 2 for low curvatures; this trend does not last, as a heavy deviation between them arises as we progress to more curved flakes, resulting in an overall monotonic growth for analytical T Cl and a seemingly parabolic plot for numerical T Cl . The discrepancy between the two values is attributed to the interference from T Re . Analytical T Cl , being inversely related to the first derivative, shows a compression of the spectrum in the levels studied. The small distortion on the highest positive field is due, as will become clear later, to an important shift on the levels of the packet.
This base trend is intensified with positive fields, with the discrepancy happening at lower curvatures and lower values of T Cl with stronger fields, reaching a limit value at a certain 1 / R 2 . For negative fields, both values start with a decrease—with a good agreement between them—and seem to stabilise with an increasing curvature. These limit values barely change with the field intensity but do with its direction.
The fixed vertices case is similar, though there are sudden changes on the plot due to the shift to a more optimal structure as the curvature grows. However, the main effects are still observable, so previous analysis and conclusions still apply.

3.3.2. Revival Time

The results obtained for the revival time, which are depicted in Figure 8, showcase two very different behaviours according to the direction of the field. From the initial decrease at zero field conditions, positive fields make it reach lower values at a faster rate; negative fields cause a clear divergence of T Re , reaching a maximum—infinite, ideally—at a certain value of 1 / R 2 , which shifts to a lower curvature as the field intensifies. In contrast to the classical time case, here both numerical and analytical values are very similar (indicating no interference from super-revival time).
As before, the fixed vertices case (Figure 9) merely introduces small jumps to this behaviour, but the general trends remain, which allows for an experimental setup to reach similar results to the ideal case.
By looking at the results presented in Figure 6, Figure 7, Figure 8 and Figure 9, it can be seen that both T Cl and T Re lie in the 0.01–1 ps range (except for the divergences in revival times). Since typical molecular vibrational frequencies vary from less than 10 13 Hz to approximately 10 14 Hz, there could be some interference from atomic nuclei motions. Our calculations make use of the Born–Oppenheimer approximation and do not consider these motions. Nevertheless, taking into account that the mean amplitude of atomic vibrations in graphite at low temperatures is approximately 0.05 Å [86,87] and that the diagonal of our nanoflake is 46.6 Å their effects can be neglected.

3.3.3. Electronic Spectrum

The origin of all observed features for both regeneration times lies within the electronic spectrum. Figure 10 shows how the five levels of the wavepacket change with the external field, taking the central one as the energy origin. A clear effect of the field is the shift of upper and lower levels, reaching an almost perfect degeneracy of n = 5 with n = 4 for big positive fields and a near degeneracy with n = 6 for negative ones. While all eigenvalues change with the field, as expected, the ones of interest experiment a greater effect, which ultimately translates to the regeneration times observed.
Analysing the spectrum at a fixed field value yields similar results, although there are clear differences depending on the field direction, as can be seen in Figure 11. Here, we can also appreciate further effects of both the field and curvature, like the opposite effect of the curvature increase on positive or negative fields, and even some crossings between levels, explaining the distortions on T Cl at high positive fields.
The main consequence is the appearance of the divergence for negative fields. Taking a full view at the spectrum we can deduce that, since it changes from a degeneracy with the upper level to another with the lower one, there is one field value that makes both energy differences equal and thus the spectrum symmetric in this range. From an analytical point of view, and recalling Equations (13) and (14), this would mean a = b , leading to the following results:
T Cl = 2 π ħ a and
T Re = 2 π ħ 0 / 2 ! = ,
giving an infinite T Re , explaining the divergence, and a finite T Cl value. From a numerical perspective, this results in the following:
| A ( t ) | 2 4 | c 4 | 2 | c 5 | 2 cos a ħ t cos 0 2 ħ t = 4 | c 4 | 2 | c 5 | 2 cos a ħ t ,
with a single wave of a 2 π ħ / a period modulated by another one of an infinite period, matching the analytical values.
The same analysis can be performed with the extreme cases, on which the central level becomes degenerate with the upper or lower one, resulting in a or b becoming null. For example, for positive fields ( a = 0 ) ,
E ( n 0 ) = b n 0 + 1 2 ( b 2 b n 0 ) = b 2 ,
T Cl = 2 π ħ b / 2 = 4 π ħ b ,
E ( n 0 ) = b ,
T Re = 2 π ħ b / 2 ! = 4 π ħ b .
Both classical and revival times would reach an equal value of double the classical time corresponding to b, as can be seen on the temporal evolution, again presenting a single frequency. This time, however, it is truly a single oscillation:
| A ( t ) | 2 4 | c 4 | 2 | c 5 | 2 c o s b 2 ħ t c o s b 2 ħ t = = 4 | c 4 | 2 | c 5 | 2 c o s 2 b 2 ħ t ,
resulting in both periods having the same value of 4 π ħ / b and a final temporal evolution comprised of a single wave of half the frequency, which is a squared cosine, as observed. The finite value of T Re is, of course, impossible to observe on the autocorrelation plot due to its lack of modulated superstructure. The main conclusion about this analysis is that the analytical approach is consistent with the numerical value for a narrow enough wavepacket, as long as both half-sum and half-difference—of the two main energy differences, that is—are not close in value.
Finally, all levels contributing to the wavepacket have been studied as a whole in order to check if under certain conditions of field or curvature the spectrum presented a regular distribution, similar to a harmonic oscillator, via the standard deviation of the energy differences between consecutive levels. The divergences observed were a strong hint, although they only account for the central levels. No combination of field and curvature shows any global regularity spanning all five levels, and the most equally spaced case is not regular enough to represent a truly harmonic oscillator.

4. Conclusions

This paper about curved graphene quantum dots expands our previous results by applying an external, perpendicular electric field to a spherical curved nanoflake, in order to complement the results obtained for different Gaussian curvatures without a field.
From a mechanical point of view, positive—towards the deformation—fields yield a curved nanoflake as the most stable configuration, as observed experimentally, while negative—opposing the deformation—fields cause the destabilisation of the system, which is capable of overcoming the bending effect for strong enough fields. Minimally restricted systems achieve further stabilisation, while the general trends remain unchanged. The electric dipole moment shows a clearly linear dependence with the field value, including flipping its direction as the field does the same; its magnitude increases as the deformation grows, with a small balance between opposite phenomena for negative fields.
As for the regeneration times, we have observed a practical equivalency between the same curvature and field. Classical times show initially opposite tendencies but an overall decrease towards a constant value, which is different depending on the field direction. Revival times show a monotonous decrease for positive fields and a divergence for certain values of 1 / R 2 for negative fields, similar to the one previously obtained for the hyperbolic case and also others present in graphene and silicene systems in an external magnetic or electric field, respectively. The connection between our results and these similar looking phenomena, as well as the observation of pseudo-magnetic fields on strained graphene, is yet to be fully justified. Lowering the nanoflake fixation to the minimum—vertices only—mildly distorts these results.
Finally, the analysis of the electronic spectrum provides an explanation for all regeneration times trends: the external field causes an energy shift of the levels in the wavepacket, with the central one reaching degeneracy with the upper or lower level depending on the field intensity and direction. Near-degenerate situations at intense fields give constant, equal classical and revival times; a symmetric level distribution gives rise to the observed divergences. This fine control in regeneration times achieved by adjusting the field and curvature could allow for the study of new phenomena on graphene. For instance, as we pointed out in the Introduction, it has been recently shown that double quantum dots can be created by the pseudo-magnetic fields in nanobubbles and that their quantum states can be manipulated to create a controllable qubit. By changing the curvature and the intensity of the electric field, the properties of these qubits could be optimised.

Author Contributions

Conceptualisation, A.B. and N.A.C.; methodology, N.A.C.; software, S.d.-l.-H.-S. and N.A.C.; validation, S.d.-l.-H.-S. and N.A.C.; formal analysis, A.B. and N.A.C.; investigation, S.d.-l.-H.-S.; writing—original draft preparation, S.d.-l.-H.-S.; writing—review and editing, A.B. and N.A.C.; visualisation, S.d-l-H-S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Regional Government of Castilla y León (Junta de Castilla y León), by the Ministry of Science and Innovation MICIN (grant PGC2018-097831-B-I00), and by the European Union NextGenerationEU (PRTR C17.I1).

Data Availability Statement

The data presented in this study are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
DFTDensity Functional Theory
GGAGeneralised Gradient Approximation
LDALocal Density Approximation
LUMOLower Unoccupied Molecular Orbital

References

  1. Cao, C.F.; Yu, B.; Guo, B.F.; Hu, W.J.; Sun, F.N.; Zhang, Z.H.; Li, S.N.; Wu, W.; Tang, L.C.; Song, P.; et al. Bio-inspired, sustainable and mechanically robust graphene oxide-based hybrid networks for efficient fire protection and warning. Chem. Eng. J. 2022, 439, 134516. [Google Scholar] [CrossRef]
  2. Chen, B.; Ding, J.; Bai, X.; Zhang, H.; Liang, M.; Zhu, S.; Shi, C.; Ma, L.; Liu, E.; Zhao, N.; et al. Engineering Pocket-Like Graphene–Shell Encapsulated FeS2: Inhibiting Polysulfides Shuttle Effect in Potassium-Ion Batteries. Adv. Funct. Mater. 2022, 32, 2109899. [Google Scholar] [CrossRef]
  3. Su, M.; Han, G.; Gao, J.; Feng, Y.; He, C.; Ma, J.; Liu, C.; Shen, C. Carbon welding on graphene skeleton for phase change composites with high thermal conductivity for solar-to-heat conversion. Chem. Eng. J. 2022, 427, 131665. [Google Scholar] [CrossRef]
  4. Baskakov, S.A.; Baskakova, Y.V.; Kabachkov, E.N.; Dvoretskaya, E.V.; Krasnikova, S.S.; Korepanov, V.I.; Michtchenko, A.; Shulga, Y.M. On the State of Graphene Oxide Nanosheet in a Polyurethane Matrix. Nanomaterials 2023, 13, 553. [Google Scholar] [CrossRef]
  5. Lee, K.; Lim, J.; Lee, M.J.; Ryu, K.; Lee, H.; Kim, J.Y.; Ju, H.; Cho, H.S.; Kim, B.H.; Hatzell, M.C.; et al. Structure-controlled graphene electrocatalysts for high-performance H2O2 production. Energy Environ. Sci. 2022, 15, 2858–2866. [Google Scholar] [CrossRef]
  6. Kim, J.; Kim, S.; Jung, E.; Mok, D.H.; Paidi, V.K.; Lee, J.; Lee, H.S.; Jeoun, Y.; Ko, W.; Shin, H.; et al. Atomic Structure Modification of Fe–N–C Catalysts via Morphology Engineering of Graphene for Enhanced Conversion Kinetics of Lithium–Sulfur Batteries. Adv. Funct. Mater. 2022, 32, 2110857. [Google Scholar] [CrossRef]
  7. Vasseghian, Y.; Le, V.T.; Joo, S.W.; Dragoi, E.N.; Kamyab, H.; Chelliapan, S.; Klemeš, J.J. Spotlighting graphene-based catalysts for the mitigation of environmentally hazardous pollutants to cleaner production: A review. J. Clean. Prod. 2022, 365, 132702. [Google Scholar] [CrossRef]
  8. Wang, Y.; Wu, D.; Lv, P.; He, B.; Li, X.; Ma, D.; Jia, Y. Theoretical insights into the electroreduction of nitrate to ammonia on graphene-based single-atom catalysts. Nanoscale 2022, 14, 10862–10872. [Google Scholar] [CrossRef]
  9. Afifi, M.; El-Naggar, M.E.; Muhammad, S.; Alghamdi, N.; Wageh, S.; Abou Taleb, M.F.; Mostafa, M.S.; Salem, S.; El-Tantawy, I. Compositional Adjusting and Antibacterial Improvement of Hydroxyapatite/Nb2O5/Graphene Oxide for Medical Applications. J. Inorg. Organomet. Polym. Mater. 2022, 32, 2160–2172. [Google Scholar] [CrossRef]
  10. Ashrafizadeh, M.; Saebfar, H.; Gholami, M.H.; Hushmandi, K.; Zabolian, A.; Bikarannejad, P.; Hashemi, M.; Daneshi, S.; Mirzaei, S.; Sharifi, E.; et al. Doxorubicin-loaded graphene oxide nanocomposites in cancer medicine: Stimuli-responsive carriers, co-delivery and suppressing resistance. Expert Opin. Drug Deliv. 2022, 19, 355–382. [Google Scholar] [CrossRef]
  11. Biru, E.I.; Necolau, M.I.; Zainea, A.; Iovu, H. Graphene Oxide–Protein-Based Scaffolds for Tissue Engineering: Recent Advances and Applications. Polymers 2022, 14, 1032. [Google Scholar] [CrossRef] [PubMed]
  12. Mousavi, S.M.; Hashemi, S.A.; Kalashgrani, M.Y.; Omidifar, N.; Bahrani, S.; Vijayakameswara Rao, N.; Babapoor, A.; Gholami, A.; Chiang, W.H. Bioactive Graphene Quantum Dots Based Polymer Composite for Biomedical Applications. Polymers 2022, 14, 617. [Google Scholar] [CrossRef] [PubMed]
  13. Srimaneepong, V.; Skallevold, H.E.; Khurshid, Z.; Zafar, M.S.; Rokaya, D.; Sapkota, J. Graphene for Antimicrobial and Coating Application. Int. J. Mol. Sci. 2022, 23, 499. [Google Scholar] [CrossRef] [PubMed]
  14. Purwidyantri, A.; Ipatov, A.; Domingues, T.; Borme, J.; Martins, M.; Alpuim, P.; Prado, M. Programmable graphene-based microfluidic sensor for DNA detection. Sens. Actuators B Chem. 2022, 367, 132044. [Google Scholar] [CrossRef]
  15. Shangguan, Q.; Chen, Z.; Yang, H.; Cheng, S.; Yang, W.; Yi, Z.; Wu, X.; Wang, S.; Yi, Y.; Wu, P. Design of Ultra-Narrow Band Graphene Refractive Index Sensor. Sensors 2022, 22, 6483. [Google Scholar] [CrossRef]
  16. He, S.; Zhang, Y.; Gao, J.; Nag, A.; Rahaman, A. Integration of Different Graphene Nanostructures with PDMS to Form Wearable Sensors. Nanomaterials 2022, 12, 950. [Google Scholar] [CrossRef]
  17. Tang, B.; Guo, Z.; Jin, G. Polarization-controlled and symmetry-dependent multiple plasmon-induced transparency in graphene-based metasurfaces. Opt. Express 2022, 30, 35554–35566. [Google Scholar] [CrossRef]
  18. Ye, Z.; Wu, P.; Wang, H.; Jiang, S.; Huang, M.; Lei, D.; Wu, F. Multimode tunable terahertz absorber based on a quarter graphene disk structure. Results Phys. 2023, 48, 106420. [Google Scholar] [CrossRef]
  19. Lai, R.; Shi, P.; Yi, Z.; Li, H.; Yi, Y. Triple-Band Surface Plasmon Resonance Metamaterial Absorber Based on Open-Ended Prohibited Sign Type Monolayer Graphene. Micromachines 2023, 14, 953. [Google Scholar] [CrossRef]
  20. Wei, J.; Zang, Z.; Zhang, Y.; Wang, M.; Du, J.; Tang, X. Enhanced performance of light-controlled conductive switching in hybrid cuprous oxide/reduced graphene oxide (Cu2O/rGO) nanocomposites. Opt. Lett. 2017, 42, 911–914. [Google Scholar] [CrossRef]
  21. Hu, Y.; Yao, H.; Liao, Q.; Lin, T.; Cheng, H.; Qu, L. The promising solar-powered water purification based on graphene functional architectures. EcoMat 2022, 4, e12205. [Google Scholar] [CrossRef]
  22. Ren, Z.; Zhang, H.; Liu, N.; Lei, D.; Zhang, Q.; Su, T.; Wang, L.; Su, J.; Gao, Y. Self-powered 2D nanofluidic graphene pressure sensor with Serosa-Mimetic structure. EcoMat 2023, 5, e12299. [Google Scholar] [CrossRef]
  23. Liu, J.B.; Gong, H.S.; Ye, G.L.; Fei, H.L. Graphene oxide-derived single-atom catalysts for electrochemical energy conversion. Rare Met. 2022, 41, 1703–1726. [Google Scholar] [CrossRef]
  24. Yang, L.; Su, Q.; Si, B.; Zhang, Y.; Zhang, Y.; Yang, H.; Zhou, X. Enhancing bioenergy production with carbon capture of microalgae by ultraviolet spectrum conversion via graphene oxide quantum dots. Chem. Eng. J. 2022, 429, 132230. [Google Scholar] [CrossRef]
  25. Yang, P.; Yang, X.; Liu, W.; Guo, R.; Yao, Z. Graphene-based electrocatalysts for advanced energy conversion. Green Energy Environ. 2022, 8, 1265–1278. [Google Scholar] [CrossRef]
  26. Wallace, P.R. The Band Theory of Graphite. Phys. Rev. 1947, 71, 622–634. [Google Scholar] [CrossRef]
  27. Katsnelson, M.I.; Geim, A.K. Electron scattering on microscopic corrugations in graphene. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2008, 366, 195–204. [Google Scholar] [CrossRef]
  28. Bolotin, K.; Sikes, K.; Jiang, Z.; Klima, M.; Fudenberg, G.; Hone, J.; Kim, P.; Stormer, H. Ultrahigh electron mobility in suspended graphene. Solid State Commun. 2008, 146, 351–355. [Google Scholar] [CrossRef]
  29. Ishigami, M.; Chen, J.H.; Cullen, W.G.; Fuhrer, M.S.; Williams, E.D. Atomic structure of graphene on SiO2. Nano Lett. 2007, 7, 1643–1648. [Google Scholar] [CrossRef]
  30. Castro Neto, A.H.; Guinea, F.; Peres, N.M.; Novoselov, K.S.; Geim, A.K. The electronic properties of graphene. Rev. Mod. Phys. 2009, 81, 109–162. [Google Scholar] [CrossRef]
  31. Borunda, M.F.; Hennig, H.; Heller, E.J. Ballistic versus diffusive transport in graphene. Phys. Rev. B 2013, 88, 125415. [Google Scholar] [CrossRef]
  32. Chu, Z.; He, L. Origin of room-temperature single-channel ballistic transport in zigzag graphene nanoribbons. Sci. China Mater. 2015, 58, 677–682. [Google Scholar] [CrossRef]
  33. Katsnelson, M.I.; Novoselov, K.S.; Geim, A.K. Chiral tunnelling and the Klein paradox in graphene. Nat. Phys. 2006, 2, 620–625. [Google Scholar] [CrossRef]
  34. Bai, C.; Zhang, X. Klein paradox and resonant tunneling in a graphene superlattice. Phys. Rev. B 2007, 76, 075430. [Google Scholar] [CrossRef]
  35. Van Duppen, B.; Peeters, F.M. Klein paradox for a pn junction in multilayer graphene. Europhys. Lett. 2013, 102, 27001. [Google Scholar] [CrossRef]
  36. Breit, G. An Interpretation of Dirac’s Theory of the Electron. Proc. Natl. Acad. Sci. USA 1928, 14, 553–559. [Google Scholar] [CrossRef]
  37. Schrödinger, E. Über die kräftefreie Bewegung in der Relativistischen Quantenmechanik; Sitzungsberichte der Preussischen Akademie der Wissenschaften, Physikalisch-Mathematische Klasse: Berlin, Germany, 1930; pp. 418–428. [Google Scholar]
  38. Rusin, T.M.; Zawadzki, W. Theory of electron Zitterbewegung in graphene probed by femtosecond laser pulses. Phys. Rev. B-Condens. Matter Mater. Phys. 2009, 80, 045416. [Google Scholar] [CrossRef]
  39. Zawadzki, W.; Rusin, T.M. Zitterbewegung (trembling motion) of electrons in semiconductors: A review. J. Phys. Condens. Matter 2011, 23, 143201. [Google Scholar] [CrossRef]
  40. Iorio, A. Graphene: QFT in curved spacetimes close to experiments. J. Phys. Conf. Ser. 2013, 442, 012056. [Google Scholar] [CrossRef]
  41. Iorio, A.; Lambiase, G. Quantum field theory in curved graphene spacetimes, Lobachevsky geometry, Weyl symmetry, Hawking effect, and all that. Phys. Rev. D 2014, 90, 025006. [Google Scholar] [CrossRef]
  42. Iorio, A. Curved spacetimes and curved graphene: A status report of the Weyl symmetry approach. Int. J. Mod. Phys. D 2015, 24, 1530013. [Google Scholar] [CrossRef]
  43. Gallerati, A. Negative-curvature spacetime solutions for graphene. J. Phys. Condens. Matter 2021, 33, 135501. [Google Scholar] [CrossRef]
  44. Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Zhang, Y.; Dubonos, S.V.; Grigorieva, I.V.; Firsov, A.A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666–669. [Google Scholar] [CrossRef] [PubMed]
  45. Obeng, Y.; Srinivasan, P. Graphene: Is it the future for semiconductors? An overview of the material, devices, and applications. Electrochem. Soc. Interface 2011, 20, 47–52. [Google Scholar] [CrossRef]
  46. Lone, S.; Bhardwaj, A.; Pandit, A.K.; Gupta, S.; Mahajan, S. A Review of Graphene Nanoribbon Field-Effect Transistor Structures. J. Electron. Mater. 2021, 50, 3169–3186. [Google Scholar] [CrossRef]
  47. Briggs, B.D.; Nagabhirava, B.; Rao, G.; Geer, R.; Gao, H.; Xu, Y.; Yu, B. Electromechanical robustness of monolayer graphene with extreme bending. Appl. Phys. Lett. 2010, 97, 223102. [Google Scholar] [CrossRef]
  48. Li, X.; Lau, S.P.; Tang, L.; Ji, R.; Yang, P. Multicolour light emission from chlorine-doped graphene quantum dots. J. Mater. Chem. C 2013, 1, 7308. [Google Scholar] [CrossRef]
  49. Jo, J.W.; Lee, J.U.; Jo, W.H. Graphene-based electrodes for flexible electronics. Polym. Int. 2015, 64, 1676–1684. [Google Scholar] [CrossRef]
  50. Wang, Z.; Philippe, L.; Elias, J. Deflection of suspended graphene by a transverse electric field. Phys. Rev. B-Condens. Matter Mater. Phys. 2010, 81, 155405. [Google Scholar] [CrossRef]
  51. Georgiou, T.; Britnell, L.; Blake, P.; Gorbachev, R.V.; Gholinia, A.; Geim, A.K.; Casiraghi, C.; Novoselov, K.S. Graphene bubbles with controllable curvature. Appl. Phys. Lett. 2011, 99, 2011–2014. [Google Scholar] [CrossRef]
  52. Osváth, Z.; Lefloch, F.; Bouchiat, V.; Chapelier, C. Electric field-controlled rippling of graphene. Nanoscale 2013, 5, 10996. [Google Scholar] [CrossRef] [PubMed]
  53. Mangum, J.M.; Harerimana, F.; Gikunda, M.N.; Thibado, P.M. Mechanisms of spontaneous curvature inversion in compressed graphene ripples for energy harvesting applications via molecular dynamics simulations. Membranes 2021, 11, 516. [Google Scholar] [CrossRef] [PubMed]
  54. Meyer, J.C.; Geim, A.K.; Katsnelson, M.I.; Novoselov, K.S.; Booth, T.J.; Roth, S. The structure of suspended graphene sheets. Nature 2007, 446, 60–63. [Google Scholar] [CrossRef] [PubMed]
  55. Meyer, J.C.; Geim, A.K.; Katsnelson, M.I.; Novoselov, K.S.; Obergfell, D.; Roth, S.; Girit, C.; Zettl, A. On the roughness of single- and bi-layer graphene membranes. Solid State Commun. 2007, 143, 101–109. [Google Scholar] [CrossRef]
  56. Fasolino, A.; Los, J.H.; Katsnelson, M.I. Intrinsic ripples in graphene. Nat. Mater. 2007, 6, 858–861. [Google Scholar] [CrossRef] [PubMed]
  57. Carlsson, J.M. Graphene: Buckle or break. Nat. Mater. 2007, 6, 801–802. [Google Scholar] [CrossRef]
  58. Deng, S.; Berry, V. Wrinkled, rippled and crumpled graphene: An overview of formation mechanism, electronic properties, and applications. Mater. Today 2016, 19, 197–212. [Google Scholar] [CrossRef]
  59. Ramezani Masir, M.; Moldovan, D.; Peeters, F. Pseudo magnetic field in strained graphene: Revisited. Solid State Commun. 2013, 175–176, 76–82. [Google Scholar] [CrossRef]
  60. Kang, D.H.; Sun, H.; Luo, M.; Lu, K.; Chen, M.; Kim, Y.; Jung, Y.; Gao, X.; Parluhutan, S.J.; Ge, J.; et al. Pseudo-magnetic field-induced slow carrier dynamics in periodically strained graphene. Nat. Commun. 2021, 12, 5087. [Google Scholar] [CrossRef]
  61. Park, H.C.; Han, J.; Myoung, N. A strain-engineered graphene qubit in a nanobubble. Quantum Sci. Technol. 2023, 8, 025012. [Google Scholar] [CrossRef]
  62. De-la Huerta-Sainz, S.; Ballesteros, A.; Cordero, N.A. Quantum Revivals in Curved Graphene Nanoflakes. Nanomaterials 2022, 12, 1953. [Google Scholar] [CrossRef] [PubMed]
  63. De-la Huerta-Sainz, S.; Ballesteros, A.; Cordero, N.A. Gaussian Curvature Effects on Graphene Quantum Dots. Nanomaterials 2023, 13, 95. [Google Scholar] [CrossRef] [PubMed]
  64. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision B.01; Gaussian Inc.: Wallingford, CT, USA, 2016.
  65. Girifalco, L.A.; Hodak, M. Van der Waals binding energies in graphitic structures. Phys. Rev. B 2002, 65, 125404. [Google Scholar] [CrossRef]
  66. Hod, O. Graphite and hexagonal boron-nitride have the same interlayer distance. Why? J. Chem. Theory Comput. 2012, 8, 1360–1369. [Google Scholar] [CrossRef]
  67. Khantha, M.; Cordero, N.A.; Molina, L.M.; Alonso, J.A.; Girifalco, L.A. Interaction of lithium with graphene: An ab initio study. Phys. Rev. B 2004, 70, 125422. [Google Scholar] [CrossRef]
  68. Cordero, N.A.; Alonso, J.A. The interaction of sulfuric acid with graphene and formation of adsorbed crystals. Nanotechnology 2007, 18, 485705. [Google Scholar] [CrossRef]
  69. Khantha, M.; Cordero, N.A.; Alonso, J.A.; Cawkwell, M.; Girifalco, L.A. Interaction and concerted diffusion of lithium in a (5,5) carbon nanotube. Phys. Rev. B 2008, 78, 115430. [Google Scholar] [CrossRef]
  70. Cordero, N.A.; Alonso, J.A. Interaction of Surfactants Containing a Sulfuric Group with a (5,5) Carbon Nanotube. J. Phys. Chem. C 2010, 114, 17249–17256. [Google Scholar] [CrossRef]
  71. Ayala, I.G.; Cordero, N.A.; Alonso, J.A. Surfactant effect of sulfuric acid on the exfoliation of bilayer graphene. Phys. Rev. B 2011, 84, 165424. [Google Scholar] [CrossRef]
  72. Ayala, I.G.; Cordero, N.A. Interaction of sodium bisulfate with mono- and bi-layer graphene. J. Nanoparticle Res. 2012, 14, 1071. [Google Scholar] [CrossRef]
  73. Uchida, K.; Furuya, S.; Iwata, J.I.; Oshiyama, A. Atomic corrugation and electron localization due to Moiré patterns in twisted bilayer graphenes. Phys. Rev. B 2014, 90, 155451. [Google Scholar] [CrossRef]
  74. López-Carballeira, D.; Polcar, T. Ab initio description of nanodiamonds: A DFT and TDDFT benchmark. Diam. Relat. Mater. 2020, 108, 107959. [Google Scholar] [CrossRef]
  75. Francl, M.M.; Pietro, W.J.; Hehre, W.J.; Binkley, J.S.; Gordon, M.S.; DeFrees, D.J.; Pople, J.A. Self-consistent molecular orbital methods. XXIII. A polarization-type basis set for second-row elements. J. Chem. Phys. 1982, 77, 3654–3665. [Google Scholar] [CrossRef]
  76. Dean, C.R.; Young, A.F.; Meric, I.; Lee, C.; Wang, L.; Sorgenfrei, S.; Watanabe, K.; Taniguchi, T.; Kim, P.; Shepard, K.L.; et al. Boron nitride substrates for high-quality graphene electronics. Nat. Nanotechnol. 2010, 5, 722–726. [Google Scholar] [CrossRef] [PubMed]
  77. Banszerus, L.; Janssen, H.; Otto, M.; Epping, A.; Taniguchi, T.; Watanabe, K.; Beschoten, B.; Neumaier, D.; Stampfer, C. Identifying suitable substrates for high-quality graphene-based heterostructures. 2D Mater. 2017, 4, 025030. [Google Scholar] [CrossRef]
  78. Martini, L.; Mišeikis, V.; Esteban, D.; Azpeitia, J.; Pezzini, S.; Paletti, P.; Ochapski, M.W.; Convertino, D.; Hernandez, M.G.; Jimenez, I.; et al. Scalable High-Mobility Graphene/hBN Heterostructures. ACS Appl. Mater. Interfaces 2023, 15, 37794–37801. [Google Scholar] [CrossRef]
  79. Pokropivny, V.V.; Skorokhod, V.V.; Oleinik, G.S.; Kurdyumov, A.V.; Bartnitskaya, T.S.; Pokropivny, A.V.; Sisonyuk, A.G.; Sheichenko, D.M. Boron Nitride Analogs of Fullerenes (the Fulborenes), Nanotubes, and Fullerites (the Fulborenites). J. Solid State Chem. 2000, 154, 214–222. [Google Scholar] [CrossRef]
  80. Batista, R.J.; Mazzoni, M.S.; Chacham, H. A theoretical study of the stability trends of boron nitride fullerenes. Chem. Phys. Lett. 2006, 421, 246–250. [Google Scholar] [CrossRef]
  81. Chkhartishvili, L. Boron nitride nanosystems of regular geometry. J. Phys. Conf. Ser. 2009, 176, 012014. [Google Scholar] [CrossRef]
  82. Dennington, R.; Keith, T.A.; Millam, J.M. GaussView, Version 6; Semichem Inc.: Shawnee Mission, KS, USA, 2016.
  83. Lenosky, T.; Gonze, X.; Teter, M.; Elser, V. Energetics of negatively curved graphitic carbon. Nature 1992, 355, 333–335. [Google Scholar] [CrossRef]
  84. Kim, E.A.; Castro Neto, A.H. Graphene as an electronic membrane. Europhys. Lett. 2008, 84, 57007. [Google Scholar] [CrossRef]
  85. Robinett, R.W. Quantum wave packet revivals. Phys. Rep. 2004, 392, 1–119. [Google Scholar] [CrossRef]
  86. Kellett, E.; Jackets, B.; Richards, B. A study of the amplitude of vibration of carbon atoms in the graphite structure. Carbon 1964, 2, 175–183. [Google Scholar] [CrossRef]
  87. Kelly, B. Thermal vibration amplitudes of carbon atoms in the graphite lattice parallel to the basal planes. J. Nucl. Mater. 1970, 34, 189–192. [Google Scholar] [CrossRef]
Figure 2. Curvature energy vs. 1 / R 2 plots for different values of field intensity and direction—with positive fields in red tones and negative ones in blue tones—with markers and dotted lines for total and minimal restriction, respectively.
Figure 2. Curvature energy vs. 1 / R 2 plots for different values of field intensity and direction—with positive fields in red tones and negative ones in blue tones—with markers and dotted lines for total and minimal restriction, respectively.
Micromachines 14 02035 g002
Figure 3. z component value of the electric dipole moment vector against 1 / R 2 for all seven field cases.
Figure 3. z component value of the electric dipole moment vector against 1 / R 2 for all seven field cases.
Micromachines 14 02035 g003
Figure 4. z component of the electric dipole against external field, for a nanoflake with R = 150 Å. Because of the relatively small variation within the same field, only a single case is displayed here, for clarity.
Figure 4. z component of the electric dipole against external field, for a nanoflake with R = 150 Å. Because of the relatively small variation within the same field, only a single case is displayed here, for clarity.
Micromachines 14 02035 g004
Figure 5. Temporal evolution plot for a fixed surface nanoflake, with R = 150 Å in a −50 field, as an exemplary case. Dotted lines represent the analytical values calculated for both T Cl and T Re .
Figure 5. Temporal evolution plot for a fixed surface nanoflake, with R = 150 Å in a −50 field, as an exemplary case. Dotted lines represent the analytical values calculated for both T Cl and T Re .
Micromachines 14 02035 g005
Figure 6. Classical time against 1 / R 2 for all fixed surface calculations, with markers for numerical values and dotted lines for analytical ones.
Figure 6. Classical time against 1 / R 2 for all fixed surface calculations, with markers for numerical values and dotted lines for analytical ones.
Micromachines 14 02035 g006
Figure 7. Classical time against 1 / R 2 for all fixed vertices calculations, with markers for numerical values and dotted lines for analytical ones.
Figure 7. Classical time against 1 / R 2 for all fixed vertices calculations, with markers for numerical values and dotted lines for analytical ones.
Micromachines 14 02035 g007
Figure 8. Revival time against 1 / R 2 for all fixed surface calculations, with markers for numerical values and dotted lines for analytical ones.
Figure 8. Revival time against 1 / R 2 for all fixed surface calculations, with markers for numerical values and dotted lines for analytical ones.
Micromachines 14 02035 g008
Figure 9. Revival time against 1 / R 2 for all fixed vertices calculations, with markers for numerical values and dotted lines for analytical ones.
Figure 9. Revival time against 1 / R 2 for all fixed vertices calculations, with markers for numerical values and dotted lines for analytical ones.
Micromachines 14 02035 g009
Figure 10. Eigenvalues for all five levels contributing to the wavepacket against external field, for a nanoflake with R = 100 Å.
Figure 10. Eigenvalues for all five levels contributing to the wavepacket against external field, for a nanoflake with R = 100 Å.
Micromachines 14 02035 g010
Figure 11. Eigenvalues for all five levels at fields +250 (left) and −250 (right) plotted against Gaussian curvature. All data correspond to fixed surface nanoflakes.
Figure 11. Eigenvalues for all five levels at fields +250 (left) and −250 (right) plotted against Gaussian curvature. All data correspond to fixed surface nanoflakes.
Micromachines 14 02035 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

de-la-Huerta-Sainz, S.; Ballesteros, A.; Cordero, N.A. Electric Field Effects on Curved Graphene Quantum Dots. Micromachines 2023, 14, 2035. https://doi.org/10.3390/mi14112035

AMA Style

de-la-Huerta-Sainz S, Ballesteros A, Cordero NA. Electric Field Effects on Curved Graphene Quantum Dots. Micromachines. 2023; 14(11):2035. https://doi.org/10.3390/mi14112035

Chicago/Turabian Style

de-la-Huerta-Sainz, Sergio, Angel Ballesteros, and Nicolás A. Cordero. 2023. "Electric Field Effects on Curved Graphene Quantum Dots" Micromachines 14, no. 11: 2035. https://doi.org/10.3390/mi14112035

APA Style

de-la-Huerta-Sainz, S., Ballesteros, A., & Cordero, N. A. (2023). Electric Field Effects on Curved Graphene Quantum Dots. Micromachines, 14(11), 2035. https://doi.org/10.3390/mi14112035

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop