Next Article in Journal
Analysis of SiC/Si Heterojunction Band Energy and Interface State Characteristics for SiC/Si VDMOS
Next Article in Special Issue
Real-Time Imaging of Plasmonic Concentric Circular Gratings Fabricated by Lens–Axicon Laser Interference Lithography
Previous Article in Journal
Electronic Structure and Transport Properties of Bi2Te3 and Bi2Se3 Single Crystals
Previous Article in Special Issue
Microscale Engineering of n-Type Doping in Nanostructured Gallium Antimonide: AC Impedance Spectroscopy Insights on Grain Boundary Characterization and Strategies for Controlled Dopant Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spin Polarization and Flat Bands in Eu-Doped Nanoporous and Twisted Bilayer Graphenes

1
School of Physics and Engineering, ITMO University, 197101 St. Petersburg, Russia
2
Department of Chemistry, Kyungpook National University, Daegu 41566, Republic of Korea
3
International Research Center of Spectroscopy and Quantum Chemistry (IRC SQC), Siberian Federal University, Svobodniy pr. 79/10, 600041 Krasnoyarsk, Russia
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(10), 1889; https://doi.org/10.3390/mi14101889
Submission received: 10 September 2023 / Revised: 25 September 2023 / Accepted: 28 September 2023 / Published: 30 September 2023
(This article belongs to the Special Issue Feature Papers of Micromachines in 'Materials and Processing' 2023)

Abstract

:
Advanced two-dimensional spin-polarized heterostructures based on twisted (TBG) and nanoporous (NPBG) bilayer graphenes doped with Eu ions were theoretically proposed and studied using Periodic Boundary Conditions Density Functional theory electronic structure calculations. The significant polarization of the electronic states at the Fermi level was discovered for both Eu/NPBG(AA) and Eu/TBG lattices. Eu ions’ chemi- and physisorption to both graphenes may lead to structural deformations, drop of symmetry of low-dimensional lattices, interlayer fusion, and mutual slides of TBG graphene fragments. The frontier bands in the valence region at the vicinity of the Fermi level of both spin-polarized 2D Eu/NPBG(AA) and Eu/TBG lattices clearly demonstrate flat dispersion laws caused by localized electronic states formed by TBG Moiré patterns, which could lead to strong electron correlations and the formation of exotic quantum phases.

1. Introduction

Quantum phases [1,2,3] are strongly correlated electronic states that can be formed under specific conditions near 4 K temperature. Even at such low temperatures, any quantum electronic subsystem may be influenced by quantum fluctuations caused by a sudden change in occupation numbers or in the symmetry of electronic states. As a physical parameter is varied, quantum fluctuations can cause a sudden phase transition, like superconductor-Mott insulator transition, which separates two quantum phases that have different electronic symmetries with entirely different transport properties. Quantum phase transitions are sudden and fundamental rearrangements of strongly correlated electronic systems as they evolve from well-defined excitations in Mott insulators to a completely different set of excitations in superconducting phases of high-temperature superconductors [4,5].
The connection between localized electronic states and electronic correlations is a fundamental concept in condensed matter physics and materials science [4,5,6,7,8,9,10,11]. Electron confinement caused by the localization of some electronic states leads to electron–electron correlations with the formation of either superconducting Cooper electron pairs or Mott insulators, in which electron–electron repulsion dominates. The close proximity of the localized non-overlapping electronic states leads to pronounced strong electron correlations, which determine the main physical properties of crystalline lattices, which often exhibit unique and exotic electronic properties, such as high-temperature superconductivity, metal-insulator transitions, and magnetism. For such atomic crystalline lattices, traditional band theory, which assumes weak correlations, breaks down, and more advanced theoretical approaches, like the Hubbard, Anderson, t-J, etc., models, are used to describe the electronic behavior. Localized electronic states, especially when they are strongly localized, can lead to enhanced electronic correlations due to the close proximity and enhanced Coulomb repulsion of the electronic subsystems.
In crystalline lattices which are characterized by flat electronic bands, the electronic states are typically highly localized in real space [12]. This means that electrons within the flat bands tend to be concentrated in specific regions or sites within the crystal lattice. This enhanced electron–electron interaction is a hallmark of correlated electron systems. The strong electron–electron interactions within flat bands could lead to pronounced correlation effects and play a crucial role in determining the material’s properties. For example, the presence of flat bands can enhance the likelihood of electron pairing, which can be a key factor in the emergence of superconductivity or other exotic electronic states.
The unique physical properties of 2D bigraphenes (BGs) [13] make them promising 2D materials for advanced photo- [14], spin [15], and quantum applications [16]. For instance, BG-based superlattices, such as nanoporous bigraphene (NPBG) or twisted bilayer graphenes (TBGs), were introduced and studied theoretically [17,18] and experimentally [19,20] with multiple evidences of remarkable quantum and spin properties.
At particular magic angles with the largest one of around 1.05° [21], the Moiré patterns created by the twisted layers cause the electronic bands of the material to become nearly flat near the Fermi level and cause strong electron–electron interactions and correlations. At magic angles TBGs exhibit challenging strongly correlated [21,22,23,24], Mott insulating [25], superconducting [26], ferromagnetic [27,28,29], and topological [21,30] phases caused by the orbital motion of the valence electrons rather than their spins [31]. Strongly correlated systems [32] may lead to superconductivity and long-range quantum entanglement. Flat valence bands in the vicinity of the Fermi level may determine the main physical properties of low-dimensional lattices since the energy range of the bands is so narrow that the Coulomb interactions dominate over kinetic energy, putting these materials in a strongly correlated regime.
Twisted bilayer graphene at magic angles and at higher temperatures Tc is a strongly correlated Mott insulator [33], which is characterized by a non-conducting ground state produced by strong repulsive interactions between the electrons. During the quantum phase transition at Tc~0 K [1], the insulator turns into a superconductor when a small number of charge carriers are added to the bigraphene lattice. The stacking of graphene sheets allows access to a new family of materials with electronic behavior that is exquisitely sensitive to the atomic alignment between the layers, which affects interlayer electron motion. At very low rotation angles, Moiré patterns are produced by the misaligned lattices with dramatically enlarged unit cells and with completely reconstructed electronic states, which leads to strong correlation coupling between the electrons localized at different graphene layers. New narrow bands emerge at rotation angles below 1.05° when TBG is close to charge neutrality. The electrons in these narrow bands in the vicinity of the Fermi level are found mainly in the regions of the Moiré pattern, in which the atoms are stacked directly above each other, and TBG can be thought of as a triangular lattice of weakly coupled quantum dots with a residual tunnelling of electrons between them. The adding of charge carriers to twisted bilayer graphene leads to a Berezinskii–Kosterlitz–Thouless phase transition with Tc around 1.7 K with the formation of 2D superconducting quantum phase.
Recently, it was shown that superlattices of potential bumps driven by folded holes in nanoporous bigraphene may induce visible spin polarization by breaking the balance between two graphene sublattices, thereby creating unpaired spins [34]. Dai et al., discovered that the graphene sheet’s magnetism only appears in vacancies with unpaired electrons, opening the door to manipulate the magnetism in these two-dimensional materials [35]. As a result, due to the presence of vacancy-related defects in graphene’s lattice, experimental studies were able to discover magnetic moments ranging from 1.0 to 1.5 μ [36,37,38]. It was shown that NPBG may exhibit distinct spin polarization due to intense local electric fields generated by doping graphene lattice by spin non-polarized sp-metals, such as Li, Ca, and Al [17]. It could be a significant step forward since any mechanism that allows for the manipulation of magnetism at the nanoscale is essential for designing novel devices for spin-related applications.
The spin polarization of graphene-based low-dimensional materials can be achieved by the creation of either complex heterostructures or defect lattices based on transition metal dichalcogenides/graphene [39] heterostructures, adatom lattices [40], vacancies [38], or zig-zag edges [41]. Spin polarization is caused by the breakdown of lattice inversion symmetry and the asymmetry of spin-up and spin-down Density of Electronic States (DOS) when one spin state could become more populated than the other [38]. In particular, spin polarization could be induced by the application of a magnetic field [42], which could affect spin transport in graphene as well as the use of rare-earth metals [43].
High magnetoresistance in graphene was discovered in graphene coordinated to EuO and EuS surfaces [44]. It was suggested that the EuO/graphene heterostructures composed of a sufficiently wide and short graphene sheet and gated EuO strips would pave the way to design advanced spintronic devices [44]. Furthermore, it was suggested that coupling graphene to the model magnetic EuS insulator produces a significant (>14 T) magnetic exchange field (MEF) with the potential to reach hundreds of tesla, resulting in an orders-of-magnitude enhancement of the spin signal originating from the Zeeman spin Hall effect [45].
Experimental data demonstrated that both EuS/graphene [45] and EuO/graphene [46] heterostructures generate spin polarization that can function as 2D spin logic and memory devices that would be useful in various advanced quantum applications. A possible synthesis for Eu/bilayer graphene could be achieved by europium deposition on single-layer graphene transferred to a common SiO2/Si substrate. There are two main options for its synthesis: the intercalation of Eu under the graphene sheet or the adsorption of Eu on top of the material [46,47,48,49]. It is well known that Eu intercalation exhibits intricate structural patterns with stripes, compact islands, and channels [50]. It is possible to create ordered superstructures of intercalated Eu under the conditions of Eu distillation, where a high-temperature regime aids in the evaporation of the excess Eu [51,52]. At a high-temperature regime, however, one cannot preclude the formation of Eu nanoclusters between the substrate and the graphene [53]. As a result, the adsorption of Eu on graphene is preferable, which is predicted to proceed under milder conditions. However, this synthetic approach is more difficult since it necessitates that the amount of deposited Eu fit the EuC6 stoichiometry for single-layer graphene.
The main goal of this investigation is theoretical design and the study of complex heterostructures based on nanoporous and twisted bigraphenes with either adsorbed or intercalated Europium ions using Density Functional Theory Periodic Boundary Conditions electronic structure calculations. It was shown that Eu-doping causes spin polarization with a different spatial localization of α and β spin states and fundamental alteration of their band structures with the formation of flat bands in the vicinity of the Fermi level, which may promote strong electronic correlations with the formation of quantum phases, which are promising for advanced spin- and quantum-related applications.

2. Computational Details

The electronic structure calculations of low-dimensional crystalline lattices were performed using the Vienna Ab-initio Simulation Package (VASP) [54,55,56] within Density Functional Theory (DFT) [57,58] and Periodic Boundary Conditions (PBC). Plane-wave basis set coupled with the projector augmented wave (PAW) method [59,60], GGA-PBE [61] functional, and Grimme D3 correction [62] for Van der Waals interaction were used in the study. Taking into account the correlation effects, the simplified form of Hubbard U correction proposed by Dudarev et al. [63,64] was implemented with U = 7.5 eV and J = 0.6 eV [65]. Mönkhorst-Pack Brillouin zone k-point sampling was implemented, and the k-point mesh contained 3 × 3 × 1 k-points along with a, b, and c directions, respectively, being used for carbon-based media optimization. For all elements involved in the electronic structure calculations, PAW potentials were used. For carbon, 4 outer electrons, and for Eu, 8 outer electrons, were treated as valence electrons. For Eu, semi-core f-electrons were treated as core states despite being higher in energy than other valence states. For the density of states calculations, the 6 × 6 × 1 k-point mesh was used. A vacuum interval of 20 Å was set normal to the plane to avoid artificial interactions between adjacent unit cell images. In all calculations, the cut-off energy was equal to 600 eV. During the optimization procedure, the maximum force acting on atoms less than 0.001 eV/Å was used as a stopping criterion for structural minimization.

3. Results and Discussion

3.1. Crystalline Lattices of Proposed Heterostructures

Initial atomistic models of the heterostructures were developed using (i) AA and AB bigraphenes [5] with either adsorbed or intercalated Eu ions; (ii) AA and AB NPBG [17], and (iii) twisted bilayer graphene had a twisting angle equal to 21.8° (TBG-21; but for the simplicity, it will be referred to simply as TBG, Figure 1). For the sake of comparison, the structure of pristine AA and AB bigraphenes with 0 twisting angle and both Eu-doped/intercalated AA and AB bigraphenes with 0 twisting angle were calculated as well. The Eu atom was placed inside the pore of Eu/NPBG, atop of Eu/TBG, and between carbon-based fragments (Eu/BG) with the consequent structural optimizations. The structural parameters (number of atoms per unit cell, magnetic moments, and unit cell parameters) are displayed in Table 1. The initial symmetry of carbon-based BG, TBG, and NPBG precursors was found to be hexagonal with a = b and γ = 120°, but during structural optimization, a drop of symmetry had been found for some of them. The unit cells of all heterostructures are presented in the Supplementary Materials.
For the sake of comparison, the optimization of the 2D atomic lattices of parent TBG-21 and both AA/NPBG and AB/NPBG was performed using the same PBC GGA-PBE + U PAW D3 approach. Structural optimization revealed small effects of DFT + U corrections on structure and electronic properties of all 2D lattices. The Parent TBG-21 lattice suffered small structural distortions constituting graphene fragments, which drastically increased with the Ei adsorption.
The PBE PBC atomic structure optimization of Eu/BG-AA and Eu/BG-AB experienced drastic cell distortion. In both cases, the lattice symmetry dropped from hexagonal to monoclinic with some slides of graphene fragments (Figure 2) along the xy direction. In both BG (AA) and BG (AB) cases the Eu interlayer led to the bending of the graphene fragments induced by repulsion with Eu ions.
Topological stability of complex low-dimensial crystalline lattices with multiple non-equivalent sublattices is a key issue for the design and study of 2D periodic atomic lattices from both from theoretical [56] and experimental [66] points of view. In contrast to perfectly planar regular hexagonal/triangular symmetry 2D graphene, h-BN, g-C3N4, g-C4N3, MoS2, and MoSe2 lattices [67,68,69] may form considerably complex shapes, like rolls, nanotubes, aperiodic flakes and structure waves, caused by the slight or moderate mutual structural mismatch of structural units, like pentagons, hexagons, heptagons, etc., which constitute complex low-dimensional periodic atomic lattices. Following Topology Conservation Theorem [70] the lattices may release the excess of artificial mechanical stress caused by artificial symmetrical constrains of a perfectly planar lattice implemented by inappropriate application of Periodic Boundary Conditions through bending with out-of-plane expansion.
The lattices of the graphene fragments of Eu/TBG demonstrated significant reconstruction during structural optimization, with the formation of corrugated carbon lattices formed by triangles, pentagons, hexagons, heptagons, octagons, and even nonagons (Figure 3). In particular, Eu/TBG demonstrated a non-planar wave-shaped lattice together with interlayer stitching through covalent carbon–carbon bonds, which prevents the rolling of the heterostructure [70]. The interaction between TBG and Eu-induced structural recombination combined with the out-of-plane distortions is responsible for the reduction in mechanical stress and the relaxation of the fragments of the entire interface.

3.2. Electronic Structure and Spin States of Eu-Doped Bigraphene-Based Heterostructures

The analysis of Total (TDOS) and Partial (PDOS) Density of Electronic States for Eu/BG (Figure 4a) unequivocaly demonstrates its spin non-polarized nature. The valence electronic bands (VB) are formed by BG fragments contributions, while conduction bands (CB) involve a partial contribution of Eu adatom. The band structure (Figure S1) demonstrates the shift of the graphene Dirac cone by 0.4 eV down to the VB region. The DOS analysis of Eu/NPBG (AA) (Figure 4b) demonstrates the presence of distinctive spin polarization at the Fermi level with the absence of spin polarization in the [−8, −0.5] eV and [0.5, 8] eV energy regions. VB are formed by the NPBG (AA) fragment, while CB formation involves the partial contribution of Eu adatom together with the contribution of the NPBG (AA) fragment. The analysis of graphene and Eu partial density of states of Eu/NPBG (AB) (Figure 4c) demonstrates its spin non-polarized nature as well with CB mostly formed by the NPBG (AB) fragment with minor contributions of Eu states, and VB formed just by NPBG (AB). The DOS analysis of Eu/TBG (Figure 4d) demonstrates distinctive spin polarization in the [−0.5, 0.5] eV region, with VB formed by PDOS TBG fragment, and CB formed by PDOSes of TBG and Eu states.
To study the impact of DFT + U corrections, the Uf = 7.5 eV and Jf = 0.6 eV parameters [65] were chosen to treat strong electronic correlations defined by the presence of Eu f-electrons. The results of DOS and band structure calculations (Figures S3 and S4) demonstrated the absence of any impact of Hubbard correction at electron density at the Fermi level due to the small Eu contribution, i.e., the electron density at the Fermi level is mostly formed by carbon-based fragments consisting of only s- and p-shell-occupied electrons.
The band structure of spin-polarized Eu/NPBG (AA) and Eu/TBG heterostructures are presented in Figure 5. Both lattices can be assigned as direct flat band semiconductors because the curvature of the valence bands is smaller than the accuracy of both hybride PBE0 and GGA-PBE functionals, which is equal to 2–3 kcal/mol [71,72]. Band structure analysis proved the spin polarization of the heterostructures only in the vicinity of the Fermi levels, as it was found during DOS analysis (see above). Spin polarization could be seen in the [−0.5, 0.5] eV region for both Eu/NPBG and Eu/TBG. The symmetry of the bands and its comparison between α and β states demonstrate significantly different band dispersion for the spin-up and spin-down electrons. The shape of the second conduction band for the α spin state (−2a) is similar to the shape of the second valence band for the β spin state (2b). This fact can indicate the correspondence of the −2a band to the 2b band, induced by the Eu coordination atop TBG since the DOS at the Fermi level is fully determined by the TBG fragment (Figure 4d). Tuning the electronic structure of graphene to control spectroscopic properties was discussed as well in Refs. [73,74,75].
The band shift could be identified in the case of the Eu/NPBG (AA) interface as well. The first flat valence band of the α spin state (−1a) looks fairly similar to the first flat conduction band of the β spin state (1b), so the −1a band-to-1b band correspondence can occur. The flat nature of the band dispersion law for boundary CB and VB demonstrates their localized nature and the absence of chemical interactions between the states. In case of flat bands, the electronic transitions between VB and CB would not experience phonon relaxation, which could enhance the transition performance.
The spin density distribution for Eu/NPBG (AA) and Eu/TBG are presented in Figure 6. The different spatial localization of α and β spin states is the direct evidence of the spin-polarized nature of the heterostructures. Both α and β spin densities are localized at both the Eu atom and TBG or NPBG AA fragments. In the case of Eu/NPBG (AA), the spin density is spreading to the NPBG pore, inducing the electronic redistribution of the NPBG fragment (yellow for the α channel and cyan for the β one).
In particular, for Eu/TBG, the spin density is fully localized at the TBG fragment (Figure 6c,d). Neither α (yellow) nor β (cyan) spin densities could be found at the Eu atom, while in the case of the TBG fragment, spin density was distributed for the entire unit cell. The length of C-Eu bonds (2.65 Å for Eu/NPBG (AA) and 2.64 Å for the Eu/TBG) directly indicates the absence of strong C-Eu covalent bonding, while the spin polarization of the fragments could be promoted simply by the exchange interactions between the states.

4. Conclusions

In this study, the key Eu-doped low-dimensional 21.8° twisted bilayer graphene and both AA and AB nanoporous bigraphenes were theoretically proposed and studied using the PW PBC DFT approach. Based on electronic structure calculations, it was shown that Eu coordination on top of both the bigraphene and twisted bilayer graphene induces the glide of graphene fragments and drastic structural deformation for Eu/TBG. In accordance with the Topology Conservation Theorem, corrugated TBG graphene fragments form structural waves coupled with interlayer graphene fusion, which would stabilize the local structure and the long-order topology of 2D graphene lattice. The electronic structure calculations revealed that both Eu/NPBG (AA) and Eu/TBG are spin-polarized direct flat band semiconductors with a spin density mostly localized at graphene fragments. The localization of uncompensated spins highlights the leading role of carbon fragments in the spin polarization of the heterostructures. Heterostructure DOSes clearly demonstrate the similar distribution of electronic density of both valence and conductivity bands in the vicinity of the Fermi levels. The total density of states of the valence bands is mostly formed by carbon contributions, while conduction bands have some contributions to the Eu partial electronic states. The distinctive spin polarization and flat bands of the proposed graphene-based low-dimensional quantum materials may promote advanced spin- and quantum-related applications.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/mi14101889/s1. Figure S1. The band structure of Eu/BG heterostructure. Spin-up (α) bands are depicted in black, and spin-down (β) bands are depicted in purple. Figure S2. The band structure of Eu/NPBG heterostructures. Spin-up (α) bands are depicted in black, and spin-down (β) bands are depicted in purple. Figure S3. The band structure of Eu/TBG heterostructures calculated at GGA PBE (left panel) and GGA PBE+U (right panel) levels of theory. Spin-up (α) bands are depicted in black, and spin-down (β) bands are depicted in purple. Figure S4. Total density of states (TDOS) of Eu/TBG calculated using DFT (black curve), DFT+U with f-electons treated by Dudarev method (green curve), and DFT+U with d-electons treated by Dudarev method (purple curve). Crystal data are provided for: Eu/BG (AA), Eu/BG (AB), Eu/NPBG (AA), Eu/NPBG (AB), Eu/TBG.

Author Contributions

Conceptualization, P.V.A.; Methodology, P.V.A.; Software, I.A.M.; Validation, S.P.P.; Formal analysis, I.A.M., S.P.P. and P.V.A.; Investigation, I.A.M., G.T.O. and P.V.A.; Data curation, I.A.M. and G.T.O.; Writing—original draft, I.A.M. and P.V.A.; Writing—review & editing, S.P.P. and P.V.A.; Project administration, P.V.A.; Funding acquisition, S.P.P. and P.V.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Research Foundation of the Republic of Korea, grant NRF 2021R1A2C1010455; ITMO Fellowship program (Iu.A.M.); and by Project No. FSRZ-2023-0006 of the Ministry of Science and High Education of the Russian Federation (IRC SQC team).

Data Availability Statement

All data available in the text of the paper, SI Section and on request to the authors.

Acknowledgments

This study was supported by the National Research Foundation of the Republic of Korea, grant NRF 2021R1A2C1010455. Iu. M. expresses gratitude to the Information Technology Centre, Novosibirsk State University, for providing access to a supercomputer. Iu. M. expresses gratitude to the ITMO Fellowship program for their financial support. The IRC SQC team was supported by Project No. FSRZ-2023-0006 of the Ministry of Science and High Education of the Russian Federation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sachdev, S. Quantum Phase Transitions; Cambridge University Press: Cambridge, UK, 2011; ISBN 9780511973765. [Google Scholar]
  2. Gregg, J. The Ehrenfest Classification of Phase Transitions: Introduction and Evolution. Arch. Hist. Exact Sci. 1998, 53, 51–81. [Google Scholar]
  3. Dobrosavljevic, V.; Trivedi, N.; Valles, J.M. (Eds.) Conductor-Insulator Quantum Phase Transitions; Oxford University Press: Oxford, UK, 2012. [Google Scholar]
  4. Avramov, P.V.; Ovchinnikov, S.G. Effects of strong electron correlations in X-ray and electron spectra of High-Tc superconductors. Phys. Solid State 2000, 42, 788–809. [Google Scholar] [CrossRef]
  5. Avramov, P.V.; Ovchinnikov, S.G. Strong Electron Correlation Effects in X-ray and Photoelectron Spectra of High-Temperature Superconductors. J. Struct. Chem. 1999, 40, 108–151. [Google Scholar] [CrossRef]
  6. Avramov, P.V.; Ovchinnikov, S.G. The strong electron correlation effects in XAFS spectra of HTSC cuprates. J. Phys. IV 1997, 7, C2-183–C2-185. [Google Scholar] [CrossRef]
  7. Avramov, P.V.; Ovchinnikov, S.G. Underestimation of forbidden gap in the electronic spectra of HTSC. Phys. Solid State 1997, 39, 389–391. [Google Scholar] [CrossRef]
  8. Avramov, P.V.; Ovchinnikov, S.G.; Gavrichkov, V.A.; Ruzankin, S.P. The Theory of X-ray Absorption Spectra of Strongly Correlated Copper Oxides. Phys. C 1997, 278, 94–106. [Google Scholar] [CrossRef]
  9. Avramov, P.V.; Ovchinnikov, S.G. Non-diagram transitions in polarized X-ray absorption CuL3-spectra of HTSC. Phys. Solid State 1996, 38, 1761–1768. [Google Scholar]
  10. Avramov, P.B.; Ovchinnikov, S.G. Influence of strong electron correlations on the form of the X-ray CuK Absorption spectra of La2-xSrxCuO4. J. Exp. Theor. Phys. 1995, 81, 811–816. [Google Scholar]
  11. Ovchinnikov, S.G.; Avramov, P.V. Manifestation of strong electron correlation effects in X-ray and X-ray electron spectra of copper in High-temperature superconductors. Phys. Solid State 1995, 37, 1405–1409. [Google Scholar]
  12. Kittel, C. Introduction to Solid State Physics, 8th ed.; Willey: Hoboken, NJ, USA, 2005; ISBN 0-471-41526-X. [Google Scholar]
  13. Ghiasi, T.S.; Kaverzin, A.A.; Dismukes, A.H.; de Wal, D.K.; Roy, X.; van Wees, B.J. Electrical and thermal generation of spin currents by magnetic bilayer graphene. Nat. Nanotechnol. 2021, 16, 788–794. [Google Scholar] [CrossRef]
  14. Das, S.; Pandey, D.; Thomas, J.; Roy, T. The role of graphene and other 2D materials in solar photovoltaics. Adv. Mater. 2019, 31, 1802722. [Google Scholar] [CrossRef] [PubMed]
  15. Afzal, A.M.; Khan, M.F.; Nazir, G.; Dastgeer, G.; Aftab, S.; Akhtar, I.; Seo, Y.; Eom, J. Gate modulation of the spin-orbit interaction in bilayer graphene encapsulated by WS2 films. Sci. Rep. 2018, 8, 3412. [Google Scholar] [CrossRef] [PubMed]
  16. Oyeniyi, G.T.; Melchakova, I.A.; Chernozatonskii, L.A.; Avramov, P.V. Nanodiamond islands confined between two graphene sheets as perspective 2D quantum materials. Carbon 2022, 196, 1047–1053. [Google Scholar] [CrossRef]
  17. Melchakova, I.A.; Tenev, T.G.; Vitanov, N.V.; Tchaikovskaya, O.N.; Chernozatonskii, L.A.; Yakobson, B.I.; Avramov, P.V. Extreme structure and spontaneous lift of spin degeneracy in doped perforated bilayer graphenes. Carbon 2022, 192, 61–70. [Google Scholar] [CrossRef]
  18. Lee, B.; Kang, J. Electric-Field-Tunable Bandgaps in the Inverse-Designed Nanoporous Graphene/Graphene Heterobilayers. Adv. Electron. Mater. 2022, 8, 2200252. [Google Scholar] [CrossRef]
  19. Wang, S.; Dai, S.; Jiang, D. Continuously tunable pore size for gas separation via a bilayer nanoporous graphene membrane. ACS Appl. Nano Mater. 2018, 2, 379–384. [Google Scholar] [CrossRef]
  20. Shao, H.; Wu, Y.-C.; Lin, Z.; Taberna, P.-L.; Simon, P. Nanoporous carbon for electrochemical capacitive energy storage. Chem. Soc. Rev. 2020, 49, 3005–3039. [Google Scholar] [CrossRef]
  21. Choi, Y.; Kim, H.; Peng, Y.; Thomson, A.; Lewandowski, C.; Polski, R.; Zhang, Y.; Arora, H.S.; Watanabe, K.; Taniguchi, T. Correlation-driven topological phases in magic-angle twisted bilayer graphene. Nature 2021, 589, 536–541. [Google Scholar] [CrossRef]
  22. Xie, Y.; Lian, B.; Jäck, B.; Liu, X.; Chiu, C.-L.; Watanabe, K.; Taniguchi, T.; Bernevig, B.A.; Yazdani, A. Spectroscopic signatures of many-body correlations in magic-angle twisted bilayer graphene. Nature 2019, 572, 101–105. [Google Scholar] [CrossRef]
  23. Jiang, Y.; Lai, X.; Watanabe, K.; Taniguchi, T.; Haule, K.; Mao, J.; Andrei, E.Y. Charge order and broken rotational symmetry in magic-angle twisted bilayer graphene. Nature 2019, 573, 91–95. [Google Scholar] [CrossRef]
  24. Kerelsky, A.; McGilly, L.J.; Kennes, D.M.; Xian, L.; Yankowitz, M.; Chen, S.; Watanabe, K.; Taniguchi, T.; Hone, J.; Dean, C. Maximized electron interactions at the magic angle in twisted bilayer graphene. Nature 2019, 572, 95–100. [Google Scholar] [CrossRef] [PubMed]
  25. Liao, Y.-D.; Xu, X.-Y.; Meng, Z.-Y.; Kang, J. Correlated insulating phases in the twisted bilayer graphene. Chin. Phys. B 2021, 30, 017305. [Google Scholar] [CrossRef]
  26. Cao, Y.; Fatemi, V.; Fang, S.; Watanabe, K.; Taniguchi, T.; Kaxiras, E.; Jarillo-Herrero, P. Unconventional superconductivity in magic-angle graphene superlattices. Nature 2018, 556, 43–50. [Google Scholar] [CrossRef]
  27. Alavirad, Y.; Sau, J. Ferromagnetism and its stability from the one-magnon spectrum in twisted bilayer graphene. Phys. Rev. B 2020, 102, 235123. [Google Scholar] [CrossRef]
  28. Saito, Y.; Ge, J.; Rademaker, L.; Watanabe, K.; Taniguchi, T.; Abanin, D.A.; Young, A.F. Hofstadter subband ferromagnetism and symmetry-broken Chern insulators in twisted bilayer graphene. Nat. Phys. 2021, 17, 478–481. [Google Scholar] [CrossRef]
  29. Lin, J.-X.; Zhang, Y.-H.; Morissette, E.; Wang, Z.; Liu, S.; Rhodes, D.; Watanabe, K.; Taniguchi, T.; Hone, J.; Li, J.I.A. Spin-orbit–driven ferromagnetism at half moiré filling in magic-angle twisted bilayer graphene. Science 2022, 375, 437–441. [Google Scholar] [CrossRef] [PubMed]
  30. Park, M.J.; Jeon, S.; Lee, S.; Park, H.C.; Kim, Y. Higher-order topological corner state tunneling in twisted bilayer graphene. Carbon 2021, 174, 260–265. [Google Scholar] [CrossRef]
  31. Sharpe, A.L.; Fox, E.J.; Barnard, A.W.; Finney, J.; Watanabe, K.; Taniguchi, T.; Kastner, M.A.; Goldhaber-Gordon, D. Evidence of orbital ferromagnetism in twisted bilayer graphene aligned to hexagonal boron nitride. Nano Lett. 2021, 21, 4299–4304. [Google Scholar] [CrossRef]
  32. Balents, L.; Dean, C.R.; Efetov, D.K.; Young, A.F. Superconductivity and strong correlations in moiré flat bands. Nat. Phys. 2020, 16, 725–733. [Google Scholar] [CrossRef]
  33. Cao, Y.; Fatemi, V.; Demir, A.; Fang, S.; Tomarken, S.L.; Luo, J.Y.; Sanchez-Yamagishi, J.D.; Watanabe, K.; Taniguchi, T.; Kaxiras, E.; et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 2018, 556, 80–84. [Google Scholar] [CrossRef]
  34. Nkou, F.B.S.; Douma, D.H.; Boungou, B.R.M.; Raji, A.T.; Kenmoe, S.; Mabiala, B. Theoretical insights into magnetization in graphene containing single and interacting nanoporous defects. Phys. E Low-Dimens. Syst. Nanostructures 2021, 128, 114564. [Google Scholar] [CrossRef]
  35. Dai, X.Q.; Zhao, J.H.; Xie, M.H.; Tang, Y.N.; Li, Y.H.; Zhao, B. First-principle study of magnetism induced by vacancies in graphene. Eur. Phys. J. B 2011, 80, 343–349. [Google Scholar] [CrossRef]
  36. Zhang, Y.; Gao, F.; Gao, S.; He, L. Tunable magnetism of a single-carbon vacancy in graphene. Sci. Bull. 2020, 65, 194–200. [Google Scholar] [CrossRef] [PubMed]
  37. Zhang, Y.; Li, S.-Y.; Huang, H.; Li, W.-T.; Qiao, J.-B.; Wang, W.-X.; Yin, L.-J.; Bai, K.-K.; Duan, W.; He, L. Scanning tunneling microscopy of the π magnetism of a single carbon vacancy in graphene. Phys. Rev. Lett. 2016, 117, 166801. [Google Scholar] [CrossRef]
  38. Douma, D.H.; Poaty, L.T.; Nkou, F.B.S.; Nianga-Obambi, P.S.; Raji, A.T.; Malonda-Boungou, B.R.; M’Passi-Mabiala, B. Influence of Nitrogen and Boron on the magnetization of nanoporous graphene: A first-principle investigation. Comput. Condens. Matter. 2022, 30, e00603. [Google Scholar] [CrossRef]
  39. Sukhanova, E.V.; Kvashnin, D.G.; Popov, Z.I. Induced spin polarization in graphene via interactions with halogen doped MoS2 and MoSe2 monolayers by DFT calculations. Nanoscale 2020, 12, 23248–23258. [Google Scholar] [CrossRef]
  40. González-Herrero, H.; Gómez-Rodríguez, J.M.; Mallet, P.; Moaied, M.; Palacios, J.J.; Salgado, C.; Ugeda, M.M.; Veuillen, J.-Y.; Yndurain, F.; Brihuega, I. Atomic-scale control of graphene magnetism by using hydrogen atoms. Science 2016, 352, 437–441. [Google Scholar] [CrossRef]
  41. Blackwell, R.E.; Zhao, F.; Brooks, E.; Zhu, J.; Piskun, I.; Wang, S.; Delgado, A.; Lee, Y.-L.; Louie, S.G.; Fischer, F.R. Spin splitting of dopant edge state in magnetic zigzag graphene nanoribbons. Nature 2021, 600, 647–652. [Google Scholar] [CrossRef]
  42. Khan, M.E.; Wali, Q.; Aamir, M.; Kim, Y.-H. Spin transport properties of carbon nanotubes by ferromagnetic zigzag triangular defects: A first-principles study. Mater. Today Commun. 2022, 32, 104074. [Google Scholar] [CrossRef]
  43. Wolf, M.J.; Sürgers, C.; Fischer, G.; Beckmann, D. Spin-polarized quasiparticle transport in exchange-split superconducting aluminum on europium sulfide. Phys. Rev. B 2014, 90, 144509. [Google Scholar] [CrossRef]
  44. Zhai, X.; Blanter, Y.M. Proximity-induced diversified magnetic states and electrically controllable spin polarization in bilayer graphene: Towards layered spintronics. Phys. Rev. B 2022, 106, 75425. [Google Scholar] [CrossRef]
  45. Wei, P.; Lee, S.; Lemaitre, F.; Pinel, L.; Cutaia, D.; Cha, W.; Katmis, F.; Zhu, Y.; Heiman, D.; Hone, J. Strong interfacial exchange field in the graphene/EuS heterostructure. Nat. Mater. 2016, 15, 711–716. [Google Scholar] [CrossRef] [PubMed]
  46. Averyanov, D.V.; Sokolov, I.S.; Tokmachev, A.M.; Parfenov, O.E.; Karateev, I.A.; Taldenkov, A.N.; Storchak, V.G. High-temperature magnetism in graphene induced by proximity to EuO. ACS Appl. Mater. Interfaces 2018, 10, 20767–20774. [Google Scholar] [CrossRef] [PubMed]
  47. Swartz, A.G.; Odenthal, P.M.; Hao, Y.; Ruoff, R.S.; Kawakami, R.K. Integration of the ferromagnetic insulator EuO onto graphene. ACS Nano 2012, 6, 10063–10069. [Google Scholar] [CrossRef]
  48. Klinkhammer, J.; Förster, D.F.; Schumacher, S.; Oepen, H.P.; Michely, T.; Busse, C. Structure and magnetic properties of ultra thin textured EuO films on graphene. Appl. Phys. Lett. 2013, 103, 131601. [Google Scholar] [CrossRef]
  49. Klinkhammer, J.; Schlipf, M.; Craes, F.; Runte, S.; Michely, T.; Busse, C. Spin-polarized surface state in euo (100). Phys. Rev. Lett. 2014, 112, 16803. [Google Scholar] [CrossRef]
  50. Schumacher, S.; Förster, D.F.; Rösner, M.; Wehling, T.O.; Michely, T. Strain in epitaxial graphene visualized by intercalation. Phys. Rev. Lett. 2013, 110, 86111. [Google Scholar] [CrossRef]
  51. Huttmann, F.; Klar, D.; Atodiresei, N.; Schmitz-Antoniak, C.; Smekhova, A.; Martinez-Galera, A.J.; Caciuc, V.; Bihlmayer, G.; Blügel, S.; Michely, T. Magnetism in a graphene-4 f− 3 d hybrid system. Phys. Rev. B 2017, 95, 75427. [Google Scholar] [CrossRef]
  52. Schumacher, S.; Huttmann, F.; Petrović, M.; Witt, C.; Förster, D.F.; Vo-Van, C.; Coraux, J.; Martinez-Galera, A.J.; Sessi, V.; Vergara, I. Europium underneath graphene on Ir (111): Intercalation mechanism, magnetism, and band structure. Phys. Rev. B 2014, 90, 235437. [Google Scholar] [CrossRef]
  53. Anderson, N.A.; Hupalo, M.; Keavney, D.; Tringides, M.C.; Vaknin, D. Intercalated europium metal in epitaxial graphene on SiC. Phys. Rev. Mater. 2017, 1, 54005. [Google Scholar] [CrossRef]
  54. Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558–561. [Google Scholar] [CrossRef] [PubMed]
  55. Kresse, G.; Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metalamorphous- semiconductor transition in germanium. Phys. Rev. B 1994, 49, 14251–14269. [Google Scholar] [CrossRef] [PubMed]
  56. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B-Condens. Matter Mater. Phys. 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  57. Hohenberg, P.; Kohn, W. Inhomogeneous electron gas. Phys. Rev. 1964, 136, B864. [Google Scholar] [CrossRef]
  58. Kohn, W.; Sham, L.J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, A1133. [Google Scholar] [CrossRef]
  59. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [PubMed]
  60. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  61. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef]
  62. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef]
  63. Dudarev, S.L.; Botton, G.A.; Savrasov, S.Y.; Humphreys, C.J.; Sutton, A.P. Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study. Phys. Rev. B 1998, 57, 1505. [Google Scholar] [CrossRef]
  64. Anisimov, V.I.; Zaanen, J.; Andersen, O.K. Band theory and Mott insulators: Hubbard U instead of Stoner I. Phys. Rev. B 1991, 44, 943. [Google Scholar] [CrossRef]
  65. An, J.M.; Barabash, S.V.; Ozolins, V.; Van Schilfgaarde, M.; Belashchenko, K.D. First-principles study of phase stability of Gd-doped EuO and EuS. Phys. Rev. B-Condens. Matter Mater. Phys. 2011, 83, 064105. [Google Scholar] [CrossRef]
  66. Sayyad, M.; Qin, Y.; Kopaczek, J.; Gupta, A.; Patoary, N.; Sinha, S.; Benard, E.; Davis, A.; Yumigeta, K.; Wu, C.-L.; et al. Strain Anisotropy Driven Spontaneous Formation of Nanoscrolls from 2D Janus Layers. arXiv 2023, arXiv:2306.00162. [Google Scholar] [CrossRef]
  67. Kuklin, A.; Ågren, H.; Avramov, P.V. Structural stability of single-layer PdSe2 with pentagonal puckered morphology and its nanotubes. Phys. Chem. Chem. Phys. 2020, 22, 8289–8295. [Google Scholar] [CrossRef]
  68. Avramov, P.; Demin, V.; Luo, M.; Choi, C.H.; Sorokin, P.B.; Yakobson, B.; Chernozatonskii, L. Translation Symmetry Breakdown in Low-Dimensional Lattices of Pentagonal Rings. J. Phys. Chem. Lett. 2015, 6, 4525–4531. [Google Scholar] [CrossRef] [PubMed]
  69. Avramov, P.V.; Fedorov, D.G.; Sorokin, P.B.; Sakai, S.; Entani, S.; Ohtomo, M.; Matsumoto, Y.; Naramoto, H. Intrinsic Edge Asymmetry in Narrow Zigzag Hexagonal Heteroatomic Nanoribbons Causes their Subtle Uniform Curvature. J. Phys. Chem. Lett. 2012, 3, 2003–2008. [Google Scholar] [CrossRef]
  70. Avramov, P.V.; Kuklin, A.V. Topological and Quantum Stability of Low-Dimensional Crystalline Lattices with Multiple Nonequivalent Sublattices. New J. Phys. 2022, 24, 103015. [Google Scholar] [CrossRef]
  71. Torbjörnsson, M.; Ryde, U. Comparison of the accuracy of DFT methods for reactions with relevance to nitrogenase. Electron. Struct. 2021, 3, 034005. [Google Scholar] [CrossRef]
  72. Bogojeski, M.; Vogt-Maranto, L.; Tuckerman, M.E.; Müller, K.R.; Burke, K. Quantum chemical accuracy from density functional approximations via machine learning. Nat. Commun. 2020, 11, 5223. [Google Scholar] [CrossRef]
  73. 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]
  74. 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] [PubMed]
  75. Chen, Z.; Cai, P.; Wen, Q.; Chen, H.; Tang, Y.; Yi, Z.; Wei, K.; Li, G.; Tang, B.; Yi, Y. Graphene Multi-Frequency Broadband and Ultra-Broadband Terahertz Absorber Based on Surface Plasmon Resonance. Electronics 2023, 12, 2655. [Google Scholar] [CrossRef]
Figure 1. The atomic lattices of bigraphene-based heterostructures, namely Eu/BG (a), Eu/NPBG(AA) (b), Eu/NPBG(AB) (c), and Eu/TBG, with a (d) side view. Eu atoms are depicted in purple; carbon atoms are depicted in brown.
Figure 1. The atomic lattices of bigraphene-based heterostructures, namely Eu/BG (a), Eu/NPBG(AA) (b), Eu/NPBG(AB) (c), and Eu/TBG, with a (d) side view. Eu atoms are depicted in purple; carbon atoms are depicted in brown.
Micromachines 14 01889 g001
Figure 2. The top view of atomic lattices of Eu/BG (AA) (top) and Eu/BG (AB) (bottom) heterostructures. Eu atoms are depicted in purple, carbon atoms are depicted in brown.
Figure 2. The top view of atomic lattices of Eu/BG (AA) (top) and Eu/BG (AB) (bottom) heterostructures. Eu atoms are depicted in purple, carbon atoms are depicted in brown.
Micromachines 14 01889 g002
Figure 3. The structure and symmetry of Eu/TBG unit cell. Eu atom is depicted in purple, carbon atoms are depicted in brown, the triangular fragment is depicted in yellow, pentagonal fragments are depicted in red, the heptagonal fragment is depicted in purple, the octagonal fragment is depicted in blue, and the nanogonal fragment is depicted in green.
Figure 3. The structure and symmetry of Eu/TBG unit cell. Eu atom is depicted in purple, carbon atoms are depicted in brown, the triangular fragment is depicted in yellow, pentagonal fragments are depicted in red, the heptagonal fragment is depicted in purple, the octagonal fragment is depicted in blue, and the nanogonal fragment is depicted in green.
Micromachines 14 01889 g003
Figure 4. (a) Total (TDOS) and partial (PDOS) density of states of Eu coordinated to bigraphene (Eu/BG), (b) AA nanoporous bigraphene (Eu/NPBG (AA), (c) AB nanoporous bigraphene (Eu/NPBG (AB), and (d) twisted bilayer graphene 21 (AA nanoporous bigraphene, Eu/TBG). TDOSes are depicted in black, and Eu PDOSes are depicted in red.
Figure 4. (a) Total (TDOS) and partial (PDOS) density of states of Eu coordinated to bigraphene (Eu/BG), (b) AA nanoporous bigraphene (Eu/NPBG (AA), (c) AB nanoporous bigraphene (Eu/NPBG (AB), and (d) twisted bilayer graphene 21 (AA nanoporous bigraphene, Eu/TBG). TDOSes are depicted in black, and Eu PDOSes are depicted in red.
Micromachines 14 01889 g004
Figure 5. The band structure of Eu/TBG (left panel) and Eu/NPBG (AA) (right panel) heterostructures. Spin-up (α) bands are depicted in black, and spin-down (β) bands are depicted in purple.
Figure 5. The band structure of Eu/TBG (left panel) and Eu/NPBG (AA) (right panel) heterostructures. Spin-up (α) bands are depicted in black, and spin-down (β) bands are depicted in purple.
Micromachines 14 01889 g005
Figure 6. (a,b) Top and side views of spin density distribution for Eu/NPBG (AA), respectively. (c,d) Top and side views of spin density distribution for Eu/TBG. Spin-up (α) density is depicted in yellow, and spin-down (β) density is depicted in cyan. Carbon atoms are depicted in brown and Eu atom is depicted in purple.
Figure 6. (a,b) Top and side views of spin density distribution for Eu/NPBG (AA), respectively. (c,d) Top and side views of spin density distribution for Eu/TBG. Spin-up (α) density is depicted in yellow, and spin-down (β) density is depicted in cyan. Carbon atoms are depicted in brown and Eu atom is depicted in purple.
Micromachines 14 01889 g006
Table 1. The structural parameters (Å) and magnetic moments (µB) of bigraphene-based heterostructures doped with Eu ions.
Table 1. The structural parameters (Å) and magnetic moments (µB) of bigraphene-based heterostructures doped with Eu ions.
Number of Atoms Per Cell, NatMagnetic Moment Per Cell, µBCell Parameter a, ÅCell Parameter b, Å
Eu/BG (AA)1970.0017.2617.26
Eu/BG (AB)1970.0017.2517.26
Eu/TBG1954.1416.9316.90
Eu/NPBG (AA)1330.7414.7414.74
Eu/NPBG (AB)1370.0014.7214.72
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

Melchakova, I.A.; Oyeniyi, G.T.; Polyutov, S.P.; Avramov, P.V. Spin Polarization and Flat Bands in Eu-Doped Nanoporous and Twisted Bilayer Graphenes. Micromachines 2023, 14, 1889. https://doi.org/10.3390/mi14101889

AMA Style

Melchakova IA, Oyeniyi GT, Polyutov SP, Avramov PV. Spin Polarization and Flat Bands in Eu-Doped Nanoporous and Twisted Bilayer Graphenes. Micromachines. 2023; 14(10):1889. https://doi.org/10.3390/mi14101889

Chicago/Turabian Style

Melchakova, Iu. A., G. T. Oyeniyi, S. P. Polyutov, and P. V. Avramov. 2023. "Spin Polarization and Flat Bands in Eu-Doped Nanoporous and Twisted Bilayer Graphenes" Micromachines 14, no. 10: 1889. https://doi.org/10.3390/mi14101889

APA Style

Melchakova, I. A., Oyeniyi, G. T., Polyutov, S. P., & Avramov, P. V. (2023). Spin Polarization and Flat Bands in Eu-Doped Nanoporous and Twisted Bilayer Graphenes. Micromachines, 14(10), 1889. https://doi.org/10.3390/mi14101889

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