Next Article in Journal
Recent Advances in Synthesis and Applications of MFe2O4 (M = Co, Cu, Mn, Ni, Zn) Nanoparticles
Next Article in Special Issue
Improvement of Temperature and Optical Power of an LED by Using Microfluidic Circulating System of Graphene Solution
Previous Article in Journal
The Biotin–Avidin Interaction in Biotinylated Gold Nanoparticles and the Modulation of Their Aggregation
Previous Article in Special Issue
Oblique and Asymmetric Klein Tunneling across Smooth NP Junctions or NPN Junctions in 8-Pmmn Borophene
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Screening in Graphene: Response to External Static Electric Field and an Image-Potential Problem

1
Donostia International Physics Center (DIPC), Paseo de Manuel Lardizabal 4, E-20018 San Sebastián, Basque Country, Spain
2
Departamento de Polímeros y Materiales Avanzados: Física, Química y Tecnología, Facultad de Ciencias Químicas, Universidad del País Vasco (UPV-EHU), Apdo. 1072, E-20080 San Sebastián, Basque Country, Spain
3
IKERBASQUE, Basque Foundation for Science, E-48011 Bilbao, Basque Country, Spain
4
Department of Physics, Jack and Pearl Resnick Institute, Bar-Ilan University, Ramat-Gan 52900, Israel
5
Department of Physics and Astronomy, Hunter College of the City University of New York, 695 Park Avenue, New York, NY 10065, USA
*
Author to whom correspondence should be addressed.
Nanomaterials 2021, 11(6), 1561; https://doi.org/10.3390/nano11061561
Submission received: 4 May 2021 / Revised: 4 June 2021 / Accepted: 5 June 2021 / Published: 13 June 2021
(This article belongs to the Special Issue Graphene for Electronics)

Abstract

:
We present a detailed first-principles investigation of the response of a free-standing graphene sheet to an external perpendicular static electric field E. The charge density distribution in the vicinity of the graphene monolayer that is caused by E was determined using the pseudopotential density-functional theory approach. Different geometries were considered. The centroid of this extra density induced by an external electric field was determined as z im = 1.048 Å at vanishing E, and its dependence on E has been obtained. The thus determined z im was employed to construct the hybrid one-electron potential which generates a new set of energies for the image-potential states.

1. Introduction

The numerous properties of graphene have been intensively investigated after its experimental realization. Thousands of papers on this material were published. However, there still remains a simple unanswered question regarding the way in which the induced charge density is distributed in the vicinity of a graphene monolayer when an external electric field is applied to the graphene sheet. This topic was addressed, to some degree, by considering the problem of screening of the electric field induced by point charges in graphite [1,2,3,4,5]. Specifically, the in-plane distribution of the induced charge has been actively discussed [5,6,7,8,9,10]. As for the charge distribution in the direction perpendicular to the plane of carbon atoms, it was considered as being localized on it [5].
The perpendicular charge distribution was studied by considering two- and multi-layer graphene films [11,12,13], though to the best of our knowledge, not for monolayer graphene. Moreover, regarding the question around the location of its center of mass with respect to the carbon atoms position, we are unaware of such work for a graphene film of any thickness. As a matter of fact, this question is important since, for instance, the position of the centroid of the induced density determines the so-called image-plane position z im , (here we define the z axis as pointing in the direction perpendicular to the carbon atoms basal plane) that is a “real position” of a solid surface for many phenomena occurring there. It determines a “physical” position of a metal surface when an external perturbation is applied. This problem was widely studied in the case of metal surfaces. In general, this “real” surface position is different from the spatial localization of the top atomic layer or a geometrical crystal edge, staying towards the vacuum side [14,15,16,17,18].
It is usually assumed for a quasi two-dimensional (2D) system that the excess charge is confined within an infinitesimally thin 2D layer [5,19]. Certainly, this assumption is reasonable if the relevant distance largely exceeds the atomic scale. However, it is critical to take into consideration what occurs on the atomic scale. For instance, if one intends to construct a capacitor by adopting graphene sheets, it would be helpful to determine its “physical size” which defines its electrical properties and may be different from the geometrical distance between two graphene layers. Addionally, determination of the spatial localization of the charge induced by an external electric field can be important in understanding the phenomena occurring in field-effect transistors based on 2D materials [20,21,22].
Knowledge of the position of the center of mass of the induced charge density is important in many fields of surface science. Thus, it determines the reference plane for the image-potential felt by an external charge placed in front of a surface. If this charge is an excited electron with energy below the vacuum level, it can be trapped by this image potential in a state belonging to an infinite Rydberg -like series [23,24]. The members of this series are referred to as image-potential states (IPSs).
In the previous work devoted to the IPSs in graphene, it was assumed [25,26] that z im is located at the carbon atom plane, which seems reasonable owing to the mirror symmetry of the system. Consequently, all the quantum states should be symmetric or anti-symmetric with respect to the z = 0 plane. As a result, a double Rydberg -like series of IPSs was predicted [25] to exist in a free-standing graphene monolayer since two surfaces are separated by a single atomic layer of matter only.
Up to now, IPSs for a free-standing graphene were not studied experimentally. On the other hand, numerous measurements were performed on the graphene supported on various metallic or semiconducting substrates. Usually, the interface distance between the graphene sheet and the surface atomic layer is such that the conventional single Rydberg series of a whole system is observed. Thus, in the graphene/metal systems where the graphene atomic layer is placed closer to the substrate, only a single series of IPSs was observed [27,28,29,30,31,32,33,34,35,36]. Nevertheless, there are cases where the distance separating the graphene and the top surface atomic layer is sufficiently large so as to realize the two lowest members of the graphene double-IPS series. In scanning-tunneling microscopy (STM) measurements, evidence for the Stark-shifted first two members (symmetric and antisymmetric ones) of this series was reported in the Gr/SiC(0001) system [37,38]. These states were also clearly observed in two-photon photoemission spectroscopy experiments [39]. However, in the same system, the splitting of the IPS series was not confirmed in the Ref. [40]. In the very recent experimental paper, the arguments in favor of the splitting were presented [41].
For a description of the IPSs in the graphene/substrate systems, a number of potentials have been developed. Indeed, an accurate description of IPSs is a challenge since the conventional density-functional theory (DFT) calculations do not accurately account for a correct long-range interaction in front of solid surfaces. One of the approaches consists of constructing the nonlocal van der Waals functional [42]. Although it does not yield the correct image potential behavior at long distances away from the 2D sheet, it improves the IPS description. Another input employing a conventional DFT scheme based on the local-density approximation (LDA) consists of the construction of a hybrid potential with the same computational cost. Some others use totally model potentials [26,35,43]. Since the binding energies of IPSs are sensitive to the long-range behavior of an effective potential, a key point is the image-plane position z im with respect to the carbon atom plane. Upon construction of the model potential in the Ref. [43], the fitting procedure gave z im = 0.99 Å. This is significantly different from z im = 0 assumed in other publications [25,26].
Our goal in this work is to determine the z im value for free-standing monolayer graphene from the direct DFT calculations of redistribution of its valence charge density upon application of an external electric field. Subsequently, the thus obtained z im is employed for the construction of a new hybrid “LDA+image−tail” potential. With this potential a new set of binding energies for IPSs is obtained and compared with the previous ones.
The rest of this paper is organized as follows. In Section 2, a brief description of our calculation method and some computational details are given. In Section 3, we present our calculated results. A summary and concluding remarks are presented in Section 4.

2. Calculational Methods and Details

The band structure of a graphene monolayer in the absence and presence of an external electric field of varying intensity was obtained within the LDA by solving the Kohn-Sham equations employing a home-made band structure computer code [44]. We used norm-conserving Troullier-Martin pseudopotentials to describe the electron-ion interaction for the carbon ions [45]. At the iteration stage, the exchange-correlation potential was taken in the form given in the Refs. [46,47]. For the expansion of the wave functions, a plane-wave basis set with an energy cutoff of 50 Rydberg was employed. In a self-consistent procedure, the summation over wave vectors in the irreducible part of the first Brillouin zone (BZ) was performed over a 48 × 48 × 1 k mesh.
The self-consistent procedure was realized by considering a repeated-slab geometry with the lateral lattice constant of 2.424 Å. The external electric field applied in the direction perpendicular to the graphene plane has no translation symmetry. In order to implement it in the repeated-slab geometry, we added to the Hamiltonian a term corresponding to the extra charge σ ( z ) constant in the x-y plane as shown in Figure 1a. Its z-dependence is defined by a Gaussian with a decay length of 1 a.u. This extra charge was placed at a distance of 10 Å from the graphene plane. In order to ensure the neutrality of the system, the charge + σ was removed from the graphene system. The z variation of the extra potential added to the system is schematically shown in Figure 1a. One can see that in the gap between the graphene and the extra charge position, this potential varies linearly from V g to V e with the E = 2 π σ slope. The problem with such a geometry is that there is a discontinuity in the potential between the left and right sides. In order to employ the repeated-slab geometry, we double the unit cell by mirror reflection of the picture of Figure 1a and establishing the distance between the graphene sheets in 20 Å. The resulting lattice constant in the perpendicular direction is 40 Å. We performed calculations considering the electric fields applied to the graphene sheet ranging from −0.4 to 0.5 V/Å with a step of 0.1 V/Å and keeping the in-plane (1 × 1) geometry for carbon ion positions.
In other sets of calculations, we considered a geometry when the external electric field is applied from both sides of the graphene sheet as shown in Figure 1b. In this case, the unit cell contains only one graphene sheet and a lattice parameter of 20 Å is chosen. This geometry allows us to investigate the scale on which the charge density distribution established in graphene can be considered additively. On the other hand, this geometry is not suitable for the determination of the image-plane position since the resulting system is symmetric by construction. The third geometry considered in this study is schematically presented in Figure 1c. In this case, the two planes charged with + σ and σ are located on each side of the graphene sheet that, in turn, is kept neutral. Since the total induced charge of the graphene is zero, this geometry cannot be used for the determination of the z im position. Nevertheless, the polarization induced in the carbon atom plane by the external field can be represented by two charged planes. In such a way, each surface can be considered as the covers of the different capacitors and charged oppositely.

3. Calculation Results

The electronic structure of graphene around the Fermi level at zero external electric field is presented in Figure 2 by thick black lines. The carbon-derived bonding and antibonding π bands are marked as π and π * , respectively. The two lowest energy bands above the Fermi level characterized by strong expansion into the vacuum are marked as 1 + and 1 . At energies above the vacuum level, one can notice the quantization of the bands representing a free-electron continuum due to the finite size of the vacuum interval. In the same figure, we show how the energy position of all these bands changes when the external electric field of 0.4 V/Å (blue curves) or 0.4 V/Å (red curves) is applied. One can notice that the π and π * bands experience a shift of almost the same magnitude from the bare dispersion upon changing the sign of the electric field. On the contrary, the position of the upper energy bands with a strong expansion of its wave functions into the vacuum side changes differently for opposite signs.

3.1. Electric Field Effects

We have examined the way in which the induced charge density profile n ind ( r , E ) varies with the strength and direction of the applied electric field E. Figure 3 reports n ind ( z , E ) obtained by averaging n ind ( r , E ) in the x-y plane for the values of E ranging from −0.4 to 0.5 V/Å. In order to perform a comparison, n ind ( z , E ) is normalized by the amplitude of E. One can see that its shape deviates qualitatively from the total valence density depicted by the green solid line. This can be understood, since the total density is dominated by the σ bands that have a maximum at z = 0 . On the contrary, the induced density is generated mainly by π bands. Additionally, one can observe that the shape of the induced density only slightly depends on the sign and the magnitude of E. In general, we observe that at larger E, the shapes of n ind ( z , E ) are almost the same. However, upon reduction of the E amplitude, the variations in n ind ( z , E ) gradually increase (hardly noticeable in Figure 3). This has consequences in the calculated centroid of the induced charge density versus E, defined as
z im ( E ) = z n ind ( z , E ) d z n ind ( z , E ) d z
and presented in Figure 4. Linear interpolation gives a value of 1.048 Å for z im ( E = 0 ) . Curiously, by constructing a model potential to describe IPSs measured experimentally in graphene/substrate systems, a very close value of 0.99 Å was established for z im in graphene monolayer [43]. A similar value was chosen for the crystal border in graphene in the Ref. [48]. We expect that the value of z im obtained here should not be affected significantly by the presence of the substrate once the valence electronic structure of graphene is not modified strongly by the substrate. In Figure 4, one can notice that upon approaching the E = 0 limit, z im ( E ) starts to deviate from the linear behavior. Moreover, this deviation is different for negative and positive E. In the former case, z im shifts downward, whereas in the latter case it is shifted upward. This can be explained by the fact that with reduction in the magnitude of E, the size of the Fermi surface shrinks and the possible calculation oscillations increase. For comparison, in the insert of Figure 4, we present the way in which z im ( E ) varies with E in a free-standing Al(111) monolayer. Since, for Al, the Fermi surface is large because there are three valence electrons, the deviation from the linear behavior is small.
When we apply an external electric field to the graphene sheet from both sides according to the scheme depicted in Figure 1b, the induced charge density has a symmetric shape owing to the mirror symmetry. Its shape can be reproduced very well by superimposing that of Figure 3 onto a reflection of itself, thereby demonstrating that the response is additive. It means that once one knows how the electronic system of a graphene sheet responds to an external electric field applied from one side, the response to a more complex external perturbing field can be readily evaluated.
Clearly, in the symmetric case (Figure 1b), the calculated centroid of the induced charge density is placed at z = 0 . However, knowing that each side responds independently to external electric fields applied from both the respective sides, the resulting charge distribution can be described in the electrostatic limit by two planes charged with the same signs and located at z = z im and z = z im . We believe that this picture should hold for a case when the valence electronic system is perturbed in a photoemission experiment, for example. In this case, an excited electron is promoted above the Fermi level. If its kinetic energy is lower than the work function, it can be trapped in the discrete IPSs whose number is two times larger than in the conventional Rydberg series of the hydrogen atom [25].
In the case of the geometry described in Figure 1c, the charge redistribution in the neutral graphene caused by placing it inside a capacitor can be represented at the electrostatic level by two charged planes with the opposite signs located at z = z im and z = + z im . Moreover, we found that the shape of the calculated induced density is also reproduced very well by employing the charge density distributions obtained for the positive and negative Es reported in Figure 1a. Since the calculated induced densities and the fitting results are very similar we do not include such a figure.

3.2. Image-Potential States

In our numerical calculations devoted to IPSs, the perpendicular lattice constant was increased up to 80 Å which allowed us to obtain convergent energies for the six lowest-energy members of the series. As it was mentioned previously, IPSs cannot be properly described with the use of conventional DFT calculations since the long-ranged image-potential on the vacuum side is not reproduced correctly. Additionally, the tight-binding methods are not inherently desired for its description [49,50]. Indeed, such states are a result of screening by the valence electron system of an external point charge placed in front of a system. This many-body information is not contained in the one-particle DFT Hamiltonian. In order to overcome this problem, maintaining the computational cost at the DFT level, we constructed a hybrid “LDA+image−tail” potential V ( r ) which replaces the LDA local exchange-correlation potential term V xc ( r ) in the DFT Hamiltonian. At | z | smaller than a certain z o value this potential coincides with V xc ( r ) . For | z | > z o , it has the following form:
V ( r ) = 1 A ( x , y ) · e λ ( x , y ) · | z sgn ( z ) · z im | 4 | z sgn ( z ) · z im | .
The parameters A ( x , y ) and λ ( x , y , ) are defined from the smoothness conditions for V ( r ) and its derivative at the matching planes | z | = z o . In this work, these parameters depend on the x and y coordinates, since V xc ( x , y , z ) still has a small corrugation at the matching plane. The only parameter left is z o which is unknown. In the following, we present results for three values of z o to show the sensitivity of the image-potential state energies to it.
In Figure 5, the thick dashed line shows the hybrid “LDA+image−tail” potential averaged in the x-y plane constructed with z im = 1.048 Å and z o = 1.6 Å. One can see how at distances z larger than z o it evolves from the averaged LDA potential (thin long-dashed line) to the image-potential defined as 1 / 4 ( z z im ) (thin dashed line). Notice that the potentials we construct here and employ for the band structure calculations are symmetric according to the z = 0 plane. Here, we show its behavior for positive z only. For comparison, in Figure 5 by thick dotted line, we show the hybrid potential constructed for z im = 0 and z o = 1.6 Å of the Ref. [25]. One can see that the hybrid potential constructed with z im = 1.048 Å is noticeably lower for z larger than z o . This results in larger binding energies of IPSs. This is confirmed by the values obtained at the center of the BZ as reported in Table 1. One can see that the binding energy of the lowest-energy symmetric 1 + state increases from the 1.47 eV of the Ref. [25] to 1.58 eV here. Almost the same change is experienced by the antisymmetric 1 state. For the states with larger numbers, this shift is notably smaller. Certainly, as n is increased this difference is gradually reduced. With z im = 1.048 Å by employing z o larger than 1.6 Å we encountered a problem with the construction of the hybrid potential. Beyond this value for z o , the two matching conditions for the hybrid potential cannot be fulfilled since the image potential with z im = 1.048 Å is located too far away on the right-hand side of the LDA potential, as seen in Figure 5. Notice that the downward shift of the IPSs is observed over a whole BZ. Nevertheless, this does not significantly affect the interaction of IPSs with the scattering resonances [51] around the K ¯ point.
The approach described above for construction of a hybrid “LDA+image−tail” potential is an adoption of the conventional image-potential picture employed for solid surfaces [52]. It can also be safely applied for sufficiently thick films as well. However, in a film consisting of just a single atomic layer, the situation might be different. In such a system, in a photoemission experiment, an excited electron can occupy a quantum state with the charge density symmetrical with respect to the atom plane, contrary to what occurs for solids where only a single surface is involved. The presence of two independent surfaces was indeed taken into account in our model presented above. Nevertheless, let us consider the situation from another point of view by applying a simple image-potential picture in a different way. In this case, we will replace an excited electron with some spatial charge density distribution by two point charges having 1 2 e located at distances z and z as denoted by A and A , respectively, in Figure 6. In the graphene sheet, these two point charges create the screening charges B and B whose centers of gravity are located at z im and z im , respectively. The spacial arrangement of these screening charges should be such to ensure efficient screening and avoid any charge current. Assuming the distance z is large, let us account for the interaction between the charge A with real charges A , B, and B . The interaction of A with the charge B can be replaced by interaction with its image point charge C with positive sign located at z + 2 z im , like it occurs at a metal surface. The interaction with a point charge A is obviously a Coulomb-like one. However, the interaction with a charge B is not obvious. The space distribution of this part of the total screening charge is such as to screen the point charge A . Therefore, since z is large, for a point charge A, it can be considered as a point charge located at z im . Counting all these three contributions at the first order in 1 / z , the resulting potential takes the form V ( z ) = 1 / 4 | z + z im | , that is, it looks like a charge A interacting with a point charge located at z = z im . A factor of four in the denominator is due to a fractionally charged electron with half its charge representing the A and B charges. Clearly, this latter model might be reasonable for IPSs with high numbers n. However, it may not be good for n = 1 + since the maximum of its wave function is localized [25] around 2 Å, that is, being very close to the z im = 1.048 Å position.
Based on this picture, we constructed a hybrid “LDA+image−tail” potential for which z im is placed at −1.048 Å. This potential with the matching plane at z o = 1.6 Å is shown as the thick solid line in Figure 5. Our calculated energies for the five lowest image-potential states are reported in Table 1. Comparing them with those obtained for z im = 0 and z im = 1.048 Å we observe a significant reduction of the binding energies, especially for n = 1 . Varying z o , we do not encounter problems with construction of the hybrid potential, contrary to the situation with z im = 1.048 Å. For completeness, in Table 1, we report the image-potential energies obtained with z o = 2.1 and 2.6 Å as well. The respective hybrid potentials are reported in the insert of Figure 5. One can see that the effect of the variation in the potential caused by changing z o on the states n = 1 + and n = 1 is substantial. Thus, for the lowest-energy image-potential state the binding energy may vary from 1.58 to 1.27 eV depending on the values of z im and z o . Indeed, one can see how by increasing z o , the value for the state 1 + is approaching the energy of 1.17 eV for the surface state [53,54] obtained in the LDA calculation [25]. However, for the state 1 , it is not the case.
We believe, the measurements of energies of free-standing graphene will provide important information on its screening properties. It may contribute to establishing a detailed picture of what is going on there due to an external perturbation. So far, all the experiments on IPSs were performed on supported graphene. Sensitivity of the image-potential states to the environment where graphene was kept was significant. Our findings point out that they can also provide important information about free-standing graphene screening properties.

4. Conclusions

In this theoretical study we have reported the detailed charge density distribution produced in free-standing graphene by external static electric fields with three geometries. The image-plane position was established. Surprisingly, it is rather large, located at 1.048 Å outside the carbon atom plane. Using this information, we constructed a new potential felt by an electron excited to the image-potential states. We checked several kinds of such a potential, demonstrating sensitivity of the energies of the lowest image-potential states to the details of this potential. It would be of interest to obtain the experimental information, such as from the photoemission spectroscopy, on the image-potential state energies for free-standing graphene. We believe that the experimentally determined image-potential energies will be extremely helpful for development of a more detailed picture for the graphene potential and how it reacts to the external perturbation on the atomic scale.
Our data on z im gives support to the value used for construction of an effective potential in the graphene/substrate systems [43]. Such potentials can be developed for the study of IPSs and interface states in a large class of molecular layers with the π - π interaction similar to graphene [55,56,57,58,59]. Moreover, the information on the image-plane position can be useful for the construction of effective potentials in the systems with more complex geometries like fullerens and nanotubes, where the nearly-free states and the super-atomic orbitals, a subject of intense ongoing research, are inherently linked to IPSs in a flat graphene layer [60,61,62,63,64,65,66,67]. We believe that such a study as ours will not only be restricted to the carbon atoms case, since image-potential states can be realized in many other quasi-2D systems of current interest, like phosphorene, silicene and germanene [68], borophene [69], MXenes [70,71,72], and molecular overlayers on graphene [73].

Author Contributions

Conceptualization, V.M.S., E.K. and G.G.; methodology, V.M.S.; software and calculations, V.M.S.; writing—original draft preparation, V.M.S.; writing—review and editing, V.M.S., E.K., and G.G. All authors have read and agreed to the published version of the manuscript.

Funding

V.M.S. acknowledges support from the Project of the Basque Government for consolidated groups of the Basque University, through the Department of Universities (Q-NANOFOT IT1164-19) and from the Spanish Ministry of Science and Innovation (Grant No. PID2019–105488GB–I00). G.G. would like to acknowledge the support from the Air Force Research Laboratory (AFRL) through Grant No. FA9453-21-1-0046.

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:
2DTwo-dimensional
IPSImage-potential state
STMScanning-tunneling microscopy
DFTDensity-functional theory
LDALocal-density approximation
BZBrillouin zone

References

  1. Visscher, P.B.; Falicov, L.M. Dielectric screening in a layered electron gas. Phys. Rev. B 1971, 3, 2541–2547. [Google Scholar] [CrossRef]
  2. Pietronero, L.; Strässler, S.; Zeller, H.R.; Rice, M.J. Charge distribution in c direction in lamellar graphite acceptor intercalation compounds. Phys. Rev. Lett. 1978, 41, 763–767. [Google Scholar] [CrossRef]
  3. Pietronero, L.; Strässler, S.; Zeller, H.R. Nonlinear screening in layered semimetals. Solid State Commun. 1979, 30, 399–401. [Google Scholar] [CrossRef]
  4. Safran, S.A.; Hamann, D.R. Electrostatic interactions and staging in graphite intercalation compounds. Phys. Rev. B 1980, 22, 606–612. [Google Scholar] [CrossRef]
  5. DiVincenzo, D.P.; Mele, E.J. Self-consistent effective-mass theory for intralayer screening in graphite intercalation compounds. Phys. Rev. B 1984, 29, 1685–1694. [Google Scholar] [CrossRef]
  6. DiCenzo, S.D.; Wertheim, G.K.; Basu, S.; Fischer, J.E. Charge distribution in potassium graphite. Phys. Rev. B 1981, 24, 2270–2273. [Google Scholar] [CrossRef]
  7. DiCenzo, S.D.; Basu, S.; Wertheim, G.K.; Buchanan, D.N.E.; Fischer, J.E. In-plane charge distribution in potassium-intercalated graphite. Phys. Rev. B 1982, 25, 620–626. [Google Scholar] [CrossRef]
  8. Grunes, L.A.; Ritsko, J.J. Valence and core excitation spectra in K, Rb, and Cs alkali-metal stage-1 intercalated graphite. Phys. Rev. B 1983, 28, 3439–3446. [Google Scholar] [CrossRef]
  9. Peres, N.M.R.; Guinea, F.; Neto, A.H.C. Electronic properties of disordered two-dimensional carbon. Phys. Rev. B 2006, 73, 125411. [Google Scholar] [CrossRef] [Green Version]
  10. Polini, M.; Tomadin, A.; Asgari, R.; MacDonald, A.H. Density functional theory of graphene sheets. Phys. Rev. B 2008, 78, 115426. [Google Scholar] [CrossRef] [Green Version]
  11. Yu, E.K.; Stewart, D.A.; Tiwari, S. Ab initio study of polarizability and induced charge densities in multilayer graphene films. Phys. Rev. B 2008, 77, 195406. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, R.-N.; Dong, G.-Y.; Wang, S.-F.; Fu, G.-S.; Wang, J.-L. Intra- and inter-layer charge redistribution in biased bilayer graphene. AIP Adv. 2016, 6, 035213. [Google Scholar] [CrossRef] [Green Version]
  13. Gao, Y.L.; Okada, S. Carrier distribution control in bilayer graphene under a perpendicular electric field by interlayer stacking arrangements. Appl. Phys. Express 2021, 14, 035001. [Google Scholar] [CrossRef]
  14. Lang, N.D.; Kohn, W. Theory of metal surfaces: Induced surface charge and image potential. Phys. Rev. B 1973, 7, 3541–3550. [Google Scholar] [CrossRef]
  15. Serena, P.A.; Soler, J.M.; Garcia, N. Self-consistent image potential in a metal surface. Phys. Rev. B 1986, 34, 6767–6769. [Google Scholar] [CrossRef] [PubMed]
  16. Inglesfield, J.E. The screening of an electric field at an Al(001) surface. Surf. Sci. 1987, 188, L701–L707. [Google Scholar] [CrossRef]
  17. Eguiluz, A.G.; Hanke, W. Evaluation of the exchange-correlation potential at a metal surface from many-body perturbation theory. Phys. Rev. B 1989, 39, 10433–10436. [Google Scholar] [CrossRef] [PubMed]
  18. Kiejna, A. Image plane position at a charged surface of stabilized jellium. Surf. Sci. 1993, 287–288, 618–621. [Google Scholar] [CrossRef]
  19. McCann, E. Asymmetry gap in the electronic band structure of bilayer graphene. Phys. Rev. B 2006, 74, 161403(R). [Google Scholar] [CrossRef] [Green Version]
  20. 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 atomicall thin carbon films. Science 2004, 306, 666–669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Friori, G.; Bonaccorso, F.; Iannaccone, G.; Palacios, T.; Neumaier, D.; Seabaugh, A.; Banerjee, S.K.; Colombo, L. Electronics based on two-dimensional materials. Nat. Nanotechnol. 2014, 9, 768–779. [Google Scholar] [CrossRef] [PubMed]
  22. Grillo, A.; Di Bartolomeo, A.; Urban, F.; Passacantando, M.; Caridad, J.M.; Sun, J.; Camilli, L. Observation od 2D conduction in ultrathin germanium arsenide field-effect transistors. ACS Appl. Mater. Interfaces 2020, 12, 12998–13004. [Google Scholar] [CrossRef]
  23. Echenique, P.M.; Pendry, J.B. Existence and detection of Rydbeg states at surfaces. J. Phys. C Solid State Phys. 1978, 11, 2065–2075. [Google Scholar] [CrossRef]
  24. Echenique, P.M.; Pendry, J.B. Theory of image states at metal surfaces. Prog. Surf. Sci. 1989, 32, 111–159. [Google Scholar] [CrossRef]
  25. Silkin, V.M.; Zhao, J.; Guinea, F.; Chulkov, E.V.; Echenique, P.M.; Petek, H. Image potential states in graphene. Phys. Rev. B 2009, 80, 121408(R). [Google Scholar] [CrossRef] [Green Version]
  26. de Andres, P.L.; Echenique, P.M.; Niesner, D.; Fauster, T.; Rivacova, A. One-dimensional potential for image-potential states on graphene. New J. Phys. 2014, 16, 023012. [Google Scholar] [CrossRef] [Green Version]
  27. Borca, B.; Barja, S.; Garnica, M.; Sánchez-Portal, D.; Silkin, V.M.; Chulkov, E.V.; Hermanns, C.F.; Hinarejos, J.J.; Vxaxzquez de Parga, A.L.; Arnau, A.; et al. Potential energy landscape for hot electrons in periodically nanostructured graphene. Phys. Rev. Lett. 2010, 105, 036804. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, H.G.; Hu, H.; Pan, Y.; Mao, J.H.; Gao, M.; Guo, H.M.; Du, S.X.; Greber, T.; Gao, H.-J. Graphene based quantum dots. J. Phys. Condens. Matter 2010, 22, 302001. [Google Scholar] [CrossRef] [Green Version]
  29. Niesner, D.; Fauster, T.; Dadap, J.I.; Zaki, N.; Knox, K.R.; Yeh, P.-C.; Bhandari, R.; Osgood, R.M.; Petrović, M.; Kralj, M. Trapping surface electrons on graphene layers and islands. Phys. Rev. B 2012, 85, 081402(R). [Google Scholar] [CrossRef] [Green Version]
  30. Armbrust, N.; Güdde, J.; Jakob, P.; Höfer, U. Time-resolved two-photon photoemission of unoccupied electronic states of periodically rippled graphene on Ru(0001). Phys. Rev. Lett. 2012, 108, 056801. [Google Scholar] [CrossRef]
  31. Nobis, D.; Potenz, M.; Niesner, D.; Fauster, T. Image-potential states of graphene on noble-metal surfaces. Phys. Rev. B 2013, 88, 195435. [Google Scholar] [CrossRef] [Green Version]
  32. Niesner, D.; Fauster, T. Image-potential states and work function of graphene. J. Phys.: Condens. Matter 2014, 26, 393001. [Google Scholar] [CrossRef] [Green Version]
  33. Craes, F.; Runte, S.; Klinkhammer, J.; Kralj, M.; Michely, T.; Busse, C. Mapping image potential states on graphene quantum dots. Phys. Rev. Lett. 2013, 111, 056804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Achilli, S.; Tognolini, S.; Fava, E.; Ponzoni, S.; Drera, G.; Cepek, C.; Patera, L.L.; Africh, C.; del Castillo, E.; Trioni, M.I.; et al. Surface states characterization in the strongly interacting graphene/Ni(111) system. New J. Phys. 2018, 20, 103039. [Google Scholar] [CrossRef]
  35. Lin, Y.; Li, Y.-Z.; Sadowski, J.T.; Jin, W.-C.; Dadap, J.I.; Hybertsen, M.S.; Osgood, R.M., Jr. Excitation and characterization of image potential state electrons on quasi-free-standing graphene. Phys. Rev. B 2018, 97, 165413. [Google Scholar] [CrossRef] [Green Version]
  36. Tognolini, S.; Achilli, S.; Ponzoni, S.; Longetti, L.; Mariani, C.; Trioni, M.I.; Pagliara, S. On- and off-resonance measurement of the Image State lifetime at the graphene/Ir(111) interface. Surf. Sci. 2019, 679, 11–16. [Google Scholar] [CrossRef]
  37. Bose, S.; Silkin, V.M.; Ohmann, R.; Brihuega, I.; Vitali, L.; Michaelis, C.H.; Mallet, P.; Veuillen, J.Y.; Schneider, M.A.; Chulkov, E.V.; et al. Image potential states as a quantum probe of graphene interfaces. New J. Phys. 2010, 12, 023028. [Google Scholar] [CrossRef]
  38. Sandin, A.; Pronschinske, A.; Rowe, J.E.; Dougherty, D.B. Incomplete screening by epitaxial graphene on the Si face of 6H-SiC(0001). Appl. Phys. Lett. 2010, 97, 113104. [Google Scholar] [CrossRef]
  39. Takahashi, K.; Imamura, M.; Yamamoto, I.; Azuma, J.; Kamada, M. Image potential states in monolayer, bilayer, and trilayer epitaxial graphene studied with time- and angle-resolved two-photon photoemission spectroscopy. Phys. Rev. B 2014, 89, 155303. [Google Scholar] [CrossRef]
  40. Gugel, D.; Niesner, D.; Eikhoff, C.; Wagner, S.; Weinelt, M.; Fauster, T. Two-photon photoemission from image-potential states of epitaxial graphene. 2D Mater. 2015, 2, 045001. [Google Scholar] [CrossRef] [Green Version]
  41. Ambrosio, G.; Achilli, S.; Pagliara, S. Resonance intensity of the n = 1 image potential state of graphene on SiC via two-photon photoemission. Surf. Sci. 2021, 703, 121722. [Google Scholar] [CrossRef]
  42. Hamada, I.; Hamamoto, Y.; Morikawa, Y. Image potential states from the van der Waals density functional. J. Phys. Chem. 2017, 147, 044708. [Google Scholar] [CrossRef] [PubMed]
  43. Armbrust, N.; Güdde, J.; Höfer, U. Formation of image-potential states at the graphene/metal interface. New J. Phys. 2015, 17, 103043. [Google Scholar] [CrossRef]
  44. Silkin, V.M.; Chulkov, E.V.; Sklyadneva, I.Y.; Panin, V.E. Self-consistent pseudopotential calcualtion of the aluminum energy spectrum. Soviet Phys. J. 1984, 27, 762–767. [Google Scholar] [CrossRef]
  45. Troullier, N.; Martins, J.L. Efficient pseudopotentials for plane-wave calculations. Phys. Rev. B 1991, 43, 1993–2006. [Google Scholar] [CrossRef]
  46. Ceperley, D.M.; Alder, B.J. Ground state of the electron gas by a stochastic method. Phys. Rev. Lett. 1980, 45, 566–569. [Google Scholar] [CrossRef] [Green Version]
  47. Perdew, J.P.; Zunger, A. Self-interaction correction to density-functional approximations for many-electron systems. Phys. Rev. B 1981, 23, 5048–5079. [Google Scholar] [CrossRef] [Green Version]
  48. Krasovskii, E.E. Ab initio theory of photoemission from graphene. Nanomaterials 2021, 11, 1212. [Google Scholar] [CrossRef]
  49. Kogan, E.; Nazarov, V.U.; Silkin, V.M.; Kaveh, M. Energy bands in graphene: Comparison between the tight-binding model and ab initio calculations. Phys. Rev. B 2014, 89, 165430. [Google Scholar] [CrossRef] [Green Version]
  50. Kogan, E.; Silkin, V.M. Electronic structure of graphene: (Nearly) free electron bands versus tight-binding bands. Phys. Status Solidi B 2017, 254, 1700035. [Google Scholar] [CrossRef] [Green Version]
  51. Nazarov, V.U.; Krasovskii, E.E.; Silkin, V.M. Scattering resonances in two-dimensional crystals with application to graphene. Phys. Rev. B 2013, 87, 041405. [Google Scholar] [CrossRef] [Green Version]
  52. Chulkov, E.V.; Silkin, V.M.; Echenique, P.M. Image potential states on metal surfaces: Binding energies and wave functions. Surf. Sci. 1999, 437, 330–352. [Google Scholar] [CrossRef]
  53. Posternak, M.; Balderischi, A.; Freeman, A.J.; Wimmer, E.; Weinelt, M. Prediction of electronic interlayer states in graphite and reinterpretation of alkali bands in graphite intercalation compounds. Phys. Rev. Lett. 1983, 50, 761–764. [Google Scholar] [CrossRef]
  54. Posternak, M.; Balderischi, A.; Freeman, A.J.; Wimmer, E. Prediction of electronic surface states in layered materials: Graphite. Phys. Rev. Lett. 1984, 52, 863–866. [Google Scholar] [CrossRef]
  55. Tsirkin, S.S.; Zaitsev, N.L.; Nechaev, I.A.; Tonner, R.; Höfer, U.; Chulkov, E.V. Inelastic decay of electrons in Shockley-type metal-organic interface states. Phys. Rev. B 2015, 92, 235434. [Google Scholar] [CrossRef] [Green Version]
  56. Armbrust, N.; Schiller, F.; Güdde, J.; Höfer, U. Model potential for the description of metal/organic interface states. Sci. Rep. 2017, 7, 46561. [Google Scholar] [CrossRef]
  57. Eschmann, L.; Sabitova, A.; Temirov, R.; Tautz, F.S.; Kruger, P.; Rohlfing, M. Electric and thermoelectric transport in graphene and helical metal in finite magnetic fields. Phys. Rev. B 2019, 100, 125155. [Google Scholar] [CrossRef]
  58. Marks, M.; Armbrust, N.; Güdde, J.; Höfer, U. Impact of interface-state formation on the charge-carrier dynamics at organic-metal interfaces. New J. Phys. 2020, 22, 093042. [Google Scholar] [CrossRef]
  59. Stalberg, K.; Shibuta, M.; Höfer, U. Temperature effects on the formation and the relaxation dynamics of metal-organic interface states. Phys. Rev. B 2020, 102, 121401. [Google Scholar] [CrossRef]
  60. Feng, M.; Zhao, J.; Petek, H. Atomlike, hollow-core-bound molecular orbitals of C60. Science 2008, 320, 359–362. [Google Scholar] [CrossRef]
  61. Zhao, J.; Feng, M.; Yang, J.; Petek, H. The superatom states of fullerenes and their hybridization into the nearly free electron bands of fullerites. ACS Nano 2009, 3, 853–864. [Google Scholar] [CrossRef] [PubMed]
  62. Dutton, G.J.; Dougherty, D.B.; Jin, W.; Reutt-Robey, J.E.; Robey, S.W. Superatom orbitals of C60 on Ag(111): Two-photon photoemission and scanning tunneling spectroscopy. Phys. Rev. B 2011, 84, 195435. [Google Scholar] [CrossRef]
  63. Zhao, J.; Zheng, Q.J.; Petek, H.; Yang, J.L. Nonnuclear nearly free electron conduction channels induced by doping charge in nanotube-molecular sheet composites. J. Phys. Chem. A 2014, 118, 7255–7260. [Google Scholar] [CrossRef]
  64. Gumbs, G.; Balassis, A.; Iurov, A.; Fekete, P. Strongly localized image states of spherical graphitic particles. Sci. World J. 2014, 2014, 726303. [Google Scholar] [CrossRef] [Green Version]
  65. Knorzer, J.; Fey, C.; Sadeghpour, H.R.; Schmelcher, P. Control of multiple excited image states around segmented carbon nanotubes. J. Chem. Phys. 2015, 143, 204309. [Google Scholar] [CrossRef] [Green Version]
  66. Johansson, J.O.; Bohl, E.; Campbell, E.E.B. Super-atom molecular orbital excited states of fullerenes. Philos. Trans. R. Soc. A 2016, 374, 20150322. [Google Scholar] [CrossRef]
  67. Shibuta, M.; Yamamoto, K.; Guo, H.L.; Zhao, J.; Nakajima, A. Highly dispersive nearly free electron bands at a 2D-assembled C60 monolayer. J. Phys. Chem. C 2020, 124, 734–741. [Google Scholar] [CrossRef]
  68. Borca, B.; Castenmiller, C.; Tsvetanova, M.; Sotthewes, K.; Rudenko, A.N.; Zandvliet, H.J.W. Image potential states of germanene. 2D Mater. 2020, 7, 035021. [Google Scholar] [CrossRef]
  69. Kong, L.J.; Liu, L.R.; Chen, L.; Zhong, Q.; Cheng, P.; Li, H.; Zhang, Z.H.; Wu, K.H. One-dimensional nearly free electron states in borophene. Nanoscale 2019, 11, 15605–15611. [Google Scholar] [CrossRef]
  70. Khazaei, M.; Ranjbar, A.; Ghorbani-Asi, M.; Arai, M.; Sasaki, T.; Liang, Y.Y.; Yunoki, S. Nearly free electron states in MXenes. Phys. Rev. B 2016, 93, 205125. [Google Scholar] [CrossRef] [Green Version]
  71. Jiang, X.; Kuklin, A.V.; Baev, A.; Ge, Y.Q.; Ågren, H.; Zhang, H.; Prasad, P.N. Two-dimensional MXenes: From morphological to optical, electric, and magnetic properties and applications. Phys. Rep. 2020, 848, 1–58. [Google Scholar] [CrossRef]
  72. Wang, M.Y.; Khazaei, M.; Kawazoe, Y.; Liang, Y.Y. First-principles study of a topological phase transition induced by image potential states in MXenes. Phys. Rev. B 2021, 103, 035433. [Google Scholar] [CrossRef]
  73. Wella, S.A.; Sawada, H.; Kawaguchi, N.; Muttaqien, F.; Inagaki, K.; Hamada, I.; Morikawa, Y.; Hamamoto, Y. Hybrid image potential states in molecular overlayers on graphene. Phys. Rev. Mater. 2017, 1, 061001. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic illustration of three geometries considered in this work of a graphene sheet (solid circles) interacting with an external electric charge uniformly distributed in the x-y plane with density σ at a chosen distance in the z direction.In the geometry (a) this extra charge is placed on the right. In the case (b) the charges of the same signs are located on the left and right sides. Panel (c) illustrates the geometry when the charges of the opposite sings placed on each side.
Figure 1. Schematic illustration of three geometries considered in this work of a graphene sheet (solid circles) interacting with an external electric charge uniformly distributed in the x-y plane with density σ at a chosen distance in the z direction.In the geometry (a) this extra charge is placed on the right. In the case (b) the charges of the same signs are located on the left and right sides. Panel (c) illustrates the geometry when the charges of the opposite sings placed on each side.
Nanomaterials 11 01561 g001
Figure 2. Electronic band structure of graphene when E = 0 (thick black lines), 0.4 V/Å (thin red lines), and −0.4 V/Å (thin blue lines) obtained with application of the geometry of Figure 1a. The Fermi level, E F , is placed at zero energy. The position of the vacuum level, E vac , is shown for the zero field. The π and π * bands are marked by corresponding symbols. The two unoccupied lowest energy states around the Γ ¯ point with strong localization in the vacuum are marked as 1 + and 1 according to the Ref. [25].
Figure 2. Electronic band structure of graphene when E = 0 (thick black lines), 0.4 V/Å (thin red lines), and −0.4 V/Å (thin blue lines) obtained with application of the geometry of Figure 1a. The Fermi level, E F , is placed at zero energy. The position of the vacuum level, E vac , is shown for the zero field. The π and π * bands are marked by corresponding symbols. The two unoccupied lowest energy states around the Γ ¯ point with strong localization in the vacuum are marked as 1 + and 1 according to the Ref. [25].
Nanomaterials 11 01561 g002
Figure 3. The valence charge density of graphene averaged in the x-y plane (green thick line) and induced charge densities generated by an applied electric field E for the color-coded values shown in the insets. The induced density for E = θ × 0.1 V/Å is normalized by the value of | θ | . The origin of the z direction is taken as the carbon atom position. The image plane position z im at 1.048 Å is marked by vertical arrow (the positive value is due to the application of the electric field from the right side).
Figure 3. The valence charge density of graphene averaged in the x-y plane (green thick line) and induced charge densities generated by an applied electric field E for the color-coded values shown in the insets. The induced density for E = θ × 0.1 V/Å is normalized by the value of | θ | . The origin of the z direction is taken as the carbon atom position. The image plane position z im at 1.048 Å is marked by vertical arrow (the positive value is due to the application of the electric field from the right side).
Nanomaterials 11 01561 g003
Figure 4. Dependence of the image plane position z im ( E ) in graphene versus the electric field E amplitude (red circles and solid line). The linear interpolation is shown by a blue dashed line. The insert shows the way in which z im ( E ) depends on E in the case of an Al(111) monolayer.
Figure 4. Dependence of the image plane position z im ( E ) in graphene versus the electric field E amplitude (red circles and solid line). The linear interpolation is shown by a blue dashed line. The insert shows the way in which z im ( E ) depends on E in the case of an Al(111) monolayer.
Nanomaterials 11 01561 g004
Figure 5. The LDA potential averaged in the x-y plane as a function of the z distance is shown by the thin dashed line. Hybrid “LDA+image−tail” potentials for z im = 0, 1.048, −1.048 Å with the matching plane at z 0 =1.6 a.u. are presented as thick dotted, dashed, and solid lines, respectively. The corresponding bare image potentials are shown by thin dotted, dashed, and solid lines, respectively. Insert: Hybrid “LDA+image−tail” potentials constructed for z im = −1.048 Å with the matching planes z 0 = 1.6, 2.1, and 2.6 Å are represented by thick solid, dashed, and dashed-dotted lines, respectively.
Figure 5. The LDA potential averaged in the x-y plane as a function of the z distance is shown by the thin dashed line. Hybrid “LDA+image−tail” potentials for z im = 0, 1.048, −1.048 Å with the matching plane at z 0 =1.6 a.u. are presented as thick dotted, dashed, and solid lines, respectively. The corresponding bare image potentials are shown by thin dotted, dashed, and solid lines, respectively. Insert: Hybrid “LDA+image−tail” potentials constructed for z im = −1.048 Å with the matching planes z 0 = 1.6, 2.1, and 2.6 Å are represented by thick solid, dashed, and dashed-dotted lines, respectively.
Nanomaterials 11 01561 g005
Figure 6. Schematic illustration of the charges created in the vicinity of a graphene monolayer. The carbon atom plane is located at z = 0 . The charge distribution in an image-potential state with the centers of gravity on the positive and negative sides according to the graphene plane located at z and z are represented by two point charges as shown by blue circles A and A . The positive charge density generated in graphene in response to this external perturbation is represented by red areas B and B centered at z im and z im , respectively. As a result, the negative point charge A interacts with its own image charge C, negative charge A , and positive charge B .
Figure 6. Schematic illustration of the charges created in the vicinity of a graphene monolayer. The carbon atom plane is located at z = 0 . The charge distribution in an image-potential state with the centers of gravity on the positive and negative sides according to the graphene plane located at z and z are represented by two point charges as shown by blue circles A and A . The positive charge density generated in graphene in response to this external perturbation is represented by red areas B and B centered at z im and z im , respectively. As a result, the negative point charge A interacts with its own image charge C, negative charge A , and positive charge B .
Nanomaterials 11 01561 g006
Table 1. Binding energies (in eV) of the image-potential states in graphene obtained with the hybrid potentials constructed with z o  = 1.6 Å and  z im  = 1.048 and z im  = 0 Å. In the case of z im placed at −1.048 Å the energies are obtained for three values of the matching plane position z o . Last line presents the values of the states obtained in the LDA calculation [25].
Table 1. Binding energies (in eV) of the image-potential states in graphene obtained with the hybrid potentials constructed with z o  = 1.6 Å and  z im  = 1.048 and z im  = 0 Å. In the case of z im placed at −1.048 Å the energies are obtained for three values of the matching plane position z o . Last line presents the values of the states obtained in the LDA calculation [25].
z im ( Å ) z o ( Å ) 1 + 1 2 + 2 3 +
1.0481.61.580.840.290.210.12
01.61.470.720.250.190.11
−1.0481.61.430.640.210.160.10
2.11.300.520.190.150.11
2.61.270.490.200.150.10
LDA 1.170.25---
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Silkin, V.M.; Kogan, E.; Gumbs, G. Screening in Graphene: Response to External Static Electric Field and an Image-Potential Problem. Nanomaterials 2021, 11, 1561. https://doi.org/10.3390/nano11061561

AMA Style

Silkin VM, Kogan E, Gumbs G. Screening in Graphene: Response to External Static Electric Field and an Image-Potential Problem. Nanomaterials. 2021; 11(6):1561. https://doi.org/10.3390/nano11061561

Chicago/Turabian Style

Silkin, Vyacheslav M., Eugene Kogan, and Godfrey Gumbs. 2021. "Screening in Graphene: Response to External Static Electric Field and an Image-Potential Problem" Nanomaterials 11, no. 6: 1561. https://doi.org/10.3390/nano11061561

APA Style

Silkin, V. M., Kogan, E., & Gumbs, G. (2021). Screening in Graphene: Response to External Static Electric Field and an Image-Potential Problem. Nanomaterials, 11(6), 1561. https://doi.org/10.3390/nano11061561

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