Next Article in Journal
SANS Studies of the Gallium–Indium Alloy Structure within Regular Nanopores
Next Article in Special Issue
Modeling and X-ray Analysis of Defect Nanoclusters Formation in B4C under Ion Irradiation
Previous Article in Journal
Size and Zeta Potential Clicked Germination Attenuation and Anti-Sporangiospores Activity of PEI-Functionalized Silver Nanoparticles against COVID-19 Associated Mucorales (Rhizopus arrhizus)
Previous Article in Special Issue
Magnetic Force Probe Characterizations of Nanoscaled Ferromagnetic Domains: Finite-Element Magnetostatic Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DFT Analysis of Hole Qubits Spin State in Germanium Thin Layer

Computing Center, Far Eastern Branch of the Russian Academy of Sciences, 680000 Khabarovsk, Russia
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(13), 2244; https://doi.org/10.3390/nano12132244
Submission received: 12 May 2022 / Revised: 18 June 2022 / Accepted: 27 June 2022 / Published: 29 June 2022

Abstract

:
Due to the presence of a strong spin–orbit interaction, hole qubits in germanium are increasingly being considered as candidates for quantum computing. These objects make it possible to create electrically controlled logic gates with the basic properties of scalability, a reasonable quantum error correction, and the necessary speed of operation. In this paper, using the methods of quantum-mechanical calculations and considering the non-collinear magnetic interactions, the quantum states of the system 2D structure of Ge in the presence of even and odd numbers of holes were investigated. The spatial localizations of hole states were calculated, favorable quantum states were revealed, and the magnetic structural characteristics of the system were analyzed.

1. Introduction

Currently, there is active study into solid-state semiconductor materials that can be used to develop quantum computers. Loss and DiVincenzo [1] proposed a quantum computation model based on the spins of electrons enclosed in quantum dots. Later, the so-called DiVincenzo criteria were formulated [2]; it imposes certain requirements on the qubit system. According to these criteria, the system must be scalable and well characterized. Before starting any type of computation, the qubit states must be initialized with a reasonable speed, which is essential for quantum error correction. The system must have a sufficiently long coherence time for the qubit states. Reading computational results should be carried out without affecting neighboring qubits and, as a result, the entire quantum computing system. The first successful application of spin qubits in semiconductors was realized in gallium arsenide (GaAs) [3]. However, one of the main drawbacks of Groups III to V elements is the spin decoherence caused by surrounding nuclear spins. At the same time, silicon is characterized by much less hyperfine interactions, since it consists mainly of 28Si atoms with zero nuclear spin that can be isotopically purified. The authors of [4] achieved coherence times on the order of one second for isotopically purified silicon. At the same time, for the gate operation to be fast and fully electrically controllable, a spin–orbit interaction is required, which is absent in the electrons of silicon. Maurand R. et al. [5] showed that holes have the necessary spin–orbit interaction. In theoretical studies [6] it has also been shown that not only filled states of the conduction band but also vacant states of the valence band are promising for the realization of spin qubits. For these reasons, it makes more sense to use structures with holes to produce spin qubit systems. Compared with other semiconductors, germanium has a stronger and more controlled spin–orbit coupling. [7,8,9,10]. Theoretical studies show that, near the Gamma point, the states in the ceiling of the valence band of Ge are well-described by the Luttinger–Kohn Hamiltonian [11,12], whose eigenvalues can be grouped into states with heavy-hole (HH) and light-hole (LH), with values of spin projections on the direction of motion equal to ± 3 / 2 and ± / 2 ,   respectively . There are several candidate materials for the design of quantum calculators. Such materials can be planar Ge/SiGe heterostructures [13,14], germanium hut wires (HW) [15,16], and germanium core–shell nanowires (NW) [17,18]. HW germanium nanowires grown on silicon surfaces are of the greatest interest. It has been shown that, during the anisotropic growth of Ge on a Si(001) substrate, the quantum dot clusters are pulled along the [001] or [010] Si directions [19]. For the HW structure, the spin relaxation and phase mismatch times were measured, and a single qubit spin control operation was performed. It was also shown that Ge {105} facet formation plays a key role in determining the stability and homogeneity of nanowires [10,16,20,21].
However, all studies carried out so far are purely experimental in nature; there are no works that theoretically describe the behavior of qubits in the proposed germanium-based quantum systems. Thus, the literature data on germanium nanowires are mainly devoted to the description of their production technologies and experimental studies of the behavior of quantum dots in them. Our work is devoted to a detailed analysis of quantum states of the two-dimensional structure system of germanium in the presence of even and odd numbers of holes. The spatial localizations of hole states are investigated in detail, advantageous quantum states are identified, and the magnetic structural characteristics of the system are analyzed.

2. Computational Details

All calculations of the atomic structures, their total energy, and charge distribution were performed using the Quantum ESPRESSO software package [22]. The ultrasoft, fully relativistic form of the pseudopotential for germanium within the generalized gradient Perdew–Burke–Ernzerhof approximation in the spin–orbit interaction approximation was used. During the relaxation of the atomic structures of the unit cell and the germanium slab, all atoms were given complete freedom. For structural relaxation, we used the BFGS quasi-newton algorithm [22]. Special sets of k-points were used to sample the Brillouin zone. A 6 × 6 × 6 set was used for the Ge unit cell, and a 2 × 6 × 1 set was used for the ultrathin germanium layer. The cut-off energy of the plane waves was 680.28 eV. The values of the interatomic forces, after structural relaxation, did not exceed the value of 0.026 eV/Å. The atomic geometry and distribution of the charge and magnetic characteristics of the structures were analyzed using the Vesta software package [23].

3. Results

After the full atomic relaxation of the unit cell of bulk germanium with Fd-3m symmetry consisting of eight Ge atoms [24], we obtained the following cell parameters: a = b = c = 5.616191 Å (Figure 1a). The experimental data for this structure are a = b = c = 5.657820(5) Å [25]. The quasi-two-dimensional structure of germanium with direction (105) was constructed from a bulk cubic cell with the number of atomic layers equal to three and with 14 Ge atoms. The symmetry of the non-relaxed slab structure corresponded to the monoclinic P21/c group with a unit cell basis equal to a = 14.31850 Å, b = 5.61620 Å and the angle β equal to 128.89° (Figure 1b). After full relaxation, the 2D atomic layer remains in the monoclinic structure, but is transformed into the P2/m symmetry with cell parameters equal to a = 11.4742 Å, b = 4.3403 Å, β = 81.53° (Figure 1c). Thus, we see that, during the relaxation of the atomic structure, the cell parameters shrink, with the parameter a decreasing by almost 20% and the parameter b decreasing by almost 23%. However, since there is a transformation of the angle β, there is a corresponding increase in the interatomic distance from 2.463 Å to 2.660–2.730 Å. A similar atomic structure, in the so-called J-germanium phase, was recently theoretically predicted in another paper [26].
Then, as a model of a hole qubit in a 2D layer of germanium, a hole was created (a lack of one electron). For this purpose, one electron was removed from the structure, resulting in the formation of a hole in its place. Figure 2 shows the hole state distribution for one hole in the 2D germanium structure with an isosurface level equal to 0.004 (marked in light green). These states were determined by the difference in the bulk charge densities. The neutral system charge density was subtracted from the charge density of the charged system (with a hole present). The hole localization states for systems with two and three holes were determined in the same way. The yellow color in the figure represents the states with increasing charge density; these localizations occur as a result of atomic shifts during the hole formation in the structure. This distribution, which characterizes the hole state, is located in the center of the structure. We see from the figure that similar mirror hole states are observed closer to the edge of the structure. These mirror states are located 11.4742 Å away from each other. Thus, such states arise due to the presence of the magnetic space group P2b2/m in the 2D germanium structure under consideration, which imposes conditions for the translation of magnetic states through the two cell parameters 2b along the Y axis {2′010 |0 1 0} [27,28]. Indeed, our calculated structure of germanium with P2/m symmetry and cell parameters equal to a = 11.4742 Å, b = 4.3403 Å, β = 81.53° and shown on Figure 1c can be transformed into a triclinic structure with an elementary basis of seven atoms and parameters equal to a = 4.3403 Å, b = 5.7371 Å, α = 81.53°. Then, the distance of 11.4742 Å to which the magnetic state is translated corresponds exactly to twice the value of the parameter b. The difference between the complete magnetic states for the ±1 spin directions is 0.72 μeV. Thus, we see that, for a single-hole qubit in germanium, the quantum state |1˃ with a spin down direction s = −1 is the most energetically advantageous state compared to the |0˃ state (i.e., with a spin up direction s = +1). The magnetization corresponding to the |1˃ states is more localized on the surface atoms of the 2D germanium layer and the |0˃ states in the center of the structure (Figure 3).
The formation of two holes in the germanium structure leads to the mutual destruction of their magnetic components such that the total magnetization will be zero. Figure 4 shows the localization of hole states for two holes in the 2D germanium structure with an isosurface level equal to 0.004. Thus, an even number of holes does not lead to total magnetization; this result is consistent with the experimental results [20]. With an even number of holes, the most favorable state is the ground singlet state, i.e., with different spin directions.
The formation of three holes in the germanium structure, i.e., an odd number of them, results in the quantum state |0˃, with an upward spin direction s = +1, becoming the most advantageous state at 0.14 μeV compared to the |1˃ state. This is because, with three holes, the first two occupy a favorable ground singlet state in one orbital, and the third must occupy another higher orbital [20]. As a result, there is Coulomb repulsion between the two holes occupying the singlet state and the third [20,29,30], so a system with three holes is easier to transfer between the |0˃ and |1˃ states compared to a system with a single hole. This would require only 0.14 μeV. Figure 5 shows the localization of hole states for three holes in the germanium structure.

4. Conclusions

The atomic and electronic structure of a 2D germanium layer with a crystallographic direction (105) was studied in detail. Spatial localizations of the hole qubit states were investigated, advantageous quantum states were identified, and an analysis of the magnetic structural characteristics of the system was given. It is shown that, for a nanoscale germanium layer with a thickness of 0.27 nm, its atomic structure is transformed into a structure with a monoclinic spatial group with P2/m symmetry. This transformation is accompanied by an increase in the interatomic distance. We analyzed the quantum states of the hole qubits in the system in the presence of an even and odd number of holes. The results show that, for a single hole, the advantageous quantum state at 0.72 μeV is the |1˃ state, with a spin down direction s = −1, compared to the |0˃ state, with a spin up direction s = +1. An even number of holes in the system does not result in full magnetization. The formation of three holes causes the quantum state |0˃, with a spin up direction s = +1, to become the most advantageous state at 0.14 μeV compared to the state |1˃, so a system with three holes is easier to transfer between the quantum states |0˃ and |1˃ compared to a system with one hole. The paper shows that hole qubits are characterized by the condition of the translation of their magnetic states through two cell parameters 2b along the Y axis. We are confident that our theoretical results will be relevant and promising for use by technologists and experimentalists in the design and study of quantum computing systems based on hole qubits.

Author Contributions

Conceptualization, methodology, investigation, and writing: A.C.; calculation, investigation, analysis, writing, review, and editing: M.A.; investigation, visualization, review, and editing: M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Russian Science Foundation of the Russian Federation (project No. 22-23-01186).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Loss, D.; DiVincenzo, D.P. Quantum Computation with Quantum Dots. Phys. Rev. A 1998, 57, 120–126. [Google Scholar] [CrossRef] [Green Version]
  2. DiVincenzo, D.P. The Physical Implementation of Quantum Computation. Fortschr. Phys. 2000, 48, 771–783. [Google Scholar] [CrossRef] [Green Version]
  3. Koppens, F.H.L.; Buizert, C.; Tielrooij, K.J.; Vink, I.T.; Nowack, K.C.; Meunier, T.; Kouwenhoven, L.P.; Vandersypen, L.M.K. Driven Coherent Oscillations of a Single Electron Spin in a Quantum Dot. Nature 2006, 442, 766–771. [Google Scholar] [CrossRef] [PubMed]
  4. Muhonen, J.T.; Dehollain, J.P.; Laucht, A.; Hudson, F.E.; Kalra, R.; Sekiguchi, T.; Itoh, K.M.; Jamieson, D.N.; McCallum, J.C.; Dzurak, A.S.; et al. Storing Quantum Information for 30 Seconds in a Nanoelectronic Device. Nat. Nanotech. 2014, 9, 986–991. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Maurand, R.; Jehl, X.; Kotekar-Patil, D.; Corna, A.; Bohuslavskyi, H.; Laviéville, R.; Hutin, L.; Barraud, S.; Vinet, M.; Sanquer, M.; et al. A CMOS Silicon Spin Qubit. Nat. Commun. 2016, 7, 13575. [Google Scholar] [CrossRef]
  6. Kloeffel, C.; Loss, D. Prospects for Spin-Based Quantum Computing in Quantum Dots. Annu. Rev. Condens. Matter Phys. 2013, 4, 51–81. [Google Scholar] [CrossRef] [Green Version]
  7. Higginbotham, A.P.; Larsen, T.W.; Yao, J.; Yan, H.; Lieber, C.M.; Marcus, C.M.; Kuemmeth, F. Hole Spin Coherence in a Ge/Si Heterostructure Nanowire. Nano Lett. 2014, 14, 3582–3586. [Google Scholar] [CrossRef] [Green Version]
  8. Marcellina, E.; Hamilton, A.R.; Winkler, R.; Culcer, D. Spin-Orbit Interactions in Inversion-Asymmetric Two-Dimensional Hole Systems: A Variational Analysis. Phys. Rev. B 2017, 95, 075305. [Google Scholar] [CrossRef] [Green Version]
  9. Watzinger, H. Ge Hut Wires—From Growth to Hole Spin Resonance. Ph.D. Thesis, Faculty of the Graduate School of the Institute of Science and Technology Austria, Klosterneuburg, Austria, 2018. [Google Scholar] [CrossRef]
  10. Scappucci, G.; Kloeffel, C.; Zwanenburg, F.A.; Loss, D.; Myronov, M.; Zhang, J.-J.; De Franceschi, S.; Katsaros, G.; Veldhorst, M. The Germanium Quantum Information Route. Nat. Rev. Mater. 2021, 6, 926–943. [Google Scholar] [CrossRef]
  11. Luttinger, J.M.; Kohn, W. Motion of Electrons and Holes in Perturbed Periodic Fields. Phys. Rev. 1955, 97, 869–883. [Google Scholar] [CrossRef]
  12. Winkler, R.; Culcer, D.; Papadakis, S.J.; Habib, B.; Shayegan, M. Spin Orientation of Holes in Quantum Wells. Semicond. Sci. Technol. 2008, 23, 114017. [Google Scholar] [CrossRef] [Green Version]
  13. Aggarwal, K.; Hofmann, A.; Jirovec, D.; Prieto, I.; Sammak, A.; Botifoll, M.; Martí-Sánchez, S.; Veldhorst, M.; Arbiol, J.; Scappucci, G.; et al. Enhancement of Proximity-Induced Superconductivity in a Planar Ge Hole Gas. Phys. Rev. Res. 2021, 3, L022005. [Google Scholar] [CrossRef]
  14. Jirovec, D.; Hofmann, A.; Ballabio, A.; Mutter, P.M.; Tavani, G.; Botifoll, M.; Crippa, A.; Kukucka, J.; Sagi, O.; Martins, F.; et al. A Singlet-Triplet Hole Spin Qubit in Planar Ge. Nat. Mater. 2021, 20, 1106–1112. [Google Scholar] [CrossRef] [PubMed]
  15. Vukusic, L. Charge Sensing and Spin Relaxation Times of Holes in Ge Hut Wires. Ph.D. Thesis, Faculty of the Graduate School of the Institute of Science and Technology Austria, Klosterneuburg, Austria, 2018. [Google Scholar] [CrossRef]
  16. Gao, F.; Wang, J.; Watzinger, H.; Hu, H.; Rančić, M.J.; Zhang, J.; Wang, T.; Yao, Y.; Wang, G.; Kukučka, J.; et al. Site-Controlled Uniform Ge/Si Hut Wires with Electrically Tunable Spin–Orbit Coupling. Adv. Mater. 2020, 32, 1906523. [Google Scholar] [CrossRef] [Green Version]
  17. Lauhon, L.J.; Gudiksen, M.S.; Wang, D.; Lieber, C.M. Epitaxial Core–Shell and Core–Multishell Nanowire Heterostructures. Nature 2002, 420, 57–61. [Google Scholar] [CrossRef]
  18. Kodambaka, S.; Tersoff, J.; Reuter, M.C.; Ross, F.M. Germanium Nanowire Growth Below the Eutectic Temperature. Science 2007, 316, 729–732. [Google Scholar] [CrossRef] [PubMed]
  19. Zhang, J.J.; Katsaros, G.; Montalenti, F.; Scopece, D.; Rezaev, R.O.; Mickel, C.; Rellinghaus, B.; Miglio, L.; De Franceschi, S.; Rastelli, A.; et al. Monolithic Growth of Ultrathin Ge Nanowires on Si(001). Phys. Rev. Lett. 2012, 109, 085502. [Google Scholar] [CrossRef] [Green Version]
  20. Katsaros, G.; Kukučka, J.; Vukušić, L.; Watzinger, H.; Gao, F.; Wang, T.; Zhang, J.-J.; Held, K. Zero Field Splitting of Heavy-Hole States in Quantum Dots. Nano Lett. 2020, 20, 5201–5206. [Google Scholar] [CrossRef]
  21. Kukucka, J. Implementation of a Hole Spin Qubit in Ge Hut Wires and Dispersive Spin Sensing. Ph.D. Thesis, Faculty of the Graduate School of the Institute of Science and Technology Austria, Klosterneuburg, Austria, 2020. [Google Scholar] [CrossRef]
  22. Giannozzi, P.; Andreussi, O.; Brumme, T.; Bunau, O.; Buongiorno Nardelli, M.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Cococcioni, M.; et al. Advanced Capabilities for Materials Modelling with Quantum ESPRESSO. J. Phys. Condens. Matter. 2017, 29, 465901. [Google Scholar] [CrossRef] [Green Version]
  23. Momma, K.; Izumi, F. VESTA 3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 1272–1276. [Google Scholar] [CrossRef]
  24. Aroyo, M.I. (Ed.) International Tables for Crystallography: Space-Group Symmetry, 2nd ed.; International Union of Crystallography: Chester, UK, 2016; Volume A, ISBN 978-0-470-97423-0. [Google Scholar]
  25. Hom, T.; Kiszenik, W.; Post, B. Accurate Lattice Constants from Multiple Reflection Measurements. II. Lattice Constants of Germanium Silicon, and Diamond. J. Appl. Crystallogr. 1975, 8, 457–458. [Google Scholar] [CrossRef]
  26. Fan, Q.; Peng, H.; Zhang, W.; Yu, X.; Yun, S. Physical Properties of Group 14 Elements in P2/m Phase. J. Solid State Chem. 2022, 305, 122641. [Google Scholar] [CrossRef]
  27. Gallego, S.V.; Tasci, E.S.; de la Flor, G.; Perez-Mato, J.M.; Aroyo, M.I. Magnetic Symmetry in the Bilbao Crystallographic Server: A Computer Program to Provide Systematic Absences of Magnetic Neutron Diffraction. J. Appl. Crystallogr. 2012, 45, 1236–1247. [Google Scholar] [CrossRef]
  28. Perez-Mato, J.M.; Gallego, S.V.; Tasci, E.S.; Elcoro, L.; de la Flor, G.; Aroyo, M.I. Symmetry-Based Computational Tools for Magnetic Crystallography. Annu. Rev. Mater. Res. 2015, 45, 217–248. [Google Scholar] [CrossRef]
  29. Hanson, R.; Kouwenhoven, L.P.; Petta, J.R.; Tarucha, S.; Vandersypen, L.M.K. Spins in Few-Electron Quantum Dots. Rev. Mod. Phys. 2007, 79, 1217–1265. [Google Scholar] [CrossRef] [Green Version]
  30. Kouwenhoven, L.P.; Austing, D.G.; Tarucha, S. Few-Electron Quantum Dots. Rep. Prog. Phys. 2001, 64, 701–736. [Google Scholar] [CrossRef]
Figure 1. Structure of bulk germanium Fd-3m (a), original, unrelaxed structure of the 2D P21/c layer with direction (105) cut from the bulk structure (b), relaxed structure of the 2D P2/m layer (c).
Figure 1. Structure of bulk germanium Fd-3m (a), original, unrelaxed structure of the 2D P21/c layer with direction (105) cut from the bulk structure (b), relaxed structure of the 2D P2/m layer (c).
Nanomaterials 12 02244 g001
Figure 2. The localization of hole states (one hole) in the germanium with quantum state |1˃, with a spin down direction s = −1, is the most energetically advantageous state by 0.72 μeV compared with the state |0˃, with a spin up s = +1. The directions of the downward spins s = −1 are indicated by red arrows.
Figure 2. The localization of hole states (one hole) in the germanium with quantum state |1˃, with a spin down direction s = −1, is the most energetically advantageous state by 0.72 μeV compared with the state |0˃, with a spin up s = +1. The directions of the downward spins s = −1 are indicated by red arrows.
Nanomaterials 12 02244 g002
Figure 3. Distribution of the difference in magnetization in the 2D layer of germanium in the presence of a single hole. Yellow indicates states corresponding to |0˃ and blue indicates states corresponding to |1˃.
Figure 3. Distribution of the difference in magnetization in the 2D layer of germanium in the presence of a single hole. Yellow indicates states corresponding to |0˃ and blue indicates states corresponding to |1˃.
Nanomaterials 12 02244 g003
Figure 4. Localization of hole states for two holes in the 2D germanium.
Figure 4. Localization of hole states for two holes in the 2D germanium.
Nanomaterials 12 02244 g004
Figure 5. Localization of hole states for three holes in the 2D germanium.
Figure 5. Localization of hole states for three holes in the 2D germanium.
Nanomaterials 12 02244 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chibisov, A.; Aleshin, M.; Chibisova, M. DFT Analysis of Hole Qubits Spin State in Germanium Thin Layer. Nanomaterials 2022, 12, 2244. https://doi.org/10.3390/nano12132244

AMA Style

Chibisov A, Aleshin M, Chibisova M. DFT Analysis of Hole Qubits Spin State in Germanium Thin Layer. Nanomaterials. 2022; 12(13):2244. https://doi.org/10.3390/nano12132244

Chicago/Turabian Style

Chibisov, Andrey, Maxim Aleshin, and Mary Chibisova. 2022. "DFT Analysis of Hole Qubits Spin State in Germanium Thin Layer" Nanomaterials 12, no. 13: 2244. https://doi.org/10.3390/nano12132244

APA Style

Chibisov, A., Aleshin, M., & Chibisova, M. (2022). DFT Analysis of Hole Qubits Spin State in Germanium Thin Layer. Nanomaterials, 12(13), 2244. https://doi.org/10.3390/nano12132244

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