pith. sign in

arxiv: 2506.04724 · v4 · submitted 2025-06-05 · ❄️ cond-mat.mes-hall

Buried unstrained germanium channels: a lattice-matched platform for quantum technology

Davide Costa , Patrick Del Vecchio , Karina Hudson , Lucas E. A. Stehouwer , Alberto Tosato , Davide Degli Esposti , Vladimir Calvi , Luca Moreschini
show 3 more authors
Mario Lodari Stefano Bosco Giordano Scappucci
This is my paper

Pith reviewed 2026-05-19 11:51 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall
keywords germanium heterostructurestwo-dimensional hole gasquantum wellsspin-orbit interactionquantum transportlattice-matched barriersquantum technology
0
0 comments X

The pith

Unstrained germanium channels with lattice-matched strained silicon-germanium barriers enable high-mobility two-dimensional hole gases for quantum devices.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper proposes and demonstrates a new platform for quantum technology that uses buried unstrained germanium channels paired with a strained but lattice-matched silicon-germanium barrier. This approach removes the requirement for metamorphic buffers that introduce defects in conventional strained germanium or silicon quantum wells. In a 52 nm thick barrier structure, a two-dimensional hole gas is realized with mobility reaching 1.33×10^5 cm²/Vs and percolation density as low as 1.4×10^10 cm^{-2}. Transport data indicate significant heavy-hole and light-hole mixing through density-dependent mass and g-factor values. The platform offers prospects for strong spin-orbit coupling and integration with superconductors.

Core claim

By forming a heterojunction between unstrained Ge and a lattice-matched ε-SiGe barrier, the authors eliminate metamorphic substrates and achieve a low-disorder 2D hole gas with high mobility of 1.33×10^5 cm²/Vs and percolation density of 1.4(1)×10^10 cm^{-2}. Quantum transport reveals density-dependent in-plane effective mass and out-of-plane g-factor due to heavy-hole-light-hole mixing, with in-plane g-factors twice as large as in strained Ge, as confirmed in quantum point contact measurements.

What carries the argument

The lattice-matched heterojunction between unstrained germanium channel and strained silicon-germanium barrier that confines the two-dimensional hole gas without substrate defects.

If this is right

  • Quantum processors based on this platform could scale more reliably due to reduced defects from avoiding metamorphic buffers.
  • Strong spin-orbit interaction in the unstrained Ge channel may enable faster qubit operations.
  • The platform supports isotopic purification and superconducting pairing for hybrid quantum systems.
  • Measurements show larger in-plane g-factors compared to strained Ge, offering new tuning options for spin qubits.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This structure might simplify growth processes and improve yield in large-scale quantum device fabrication.
  • The heavy-hole-light-hole mixing could provide additional control knobs for qubit manipulation not present in purely strained systems.
  • Further studies could explore combining this with other 2D materials for enhanced functionality in quantum circuits.

Load-bearing premise

The silicon-germanium barrier maintains its full strain and lattice match to the germanium channel over the full 52 nm thickness without relaxation or harmful interface defects.

What would settle it

Direct evidence of strain relaxation in the barrier via structural characterization or a significant drop in hole mobility and rise in percolation density would show the platform does not work as claimed.

Figures

Figures reproduced from arXiv: 2506.04724 by Alberto Tosato, Davide Costa, Davide Degli Esposti, Giordano Scappucci, Karina Hudson, Luca Moreschini, Lucas E. A. Stehouwer, Mario Lodari, Patrick Del Vecchio, Stefano Bosco, Vladimir Calvi.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Layer schematic of the Ge/SiGe strained-barrier heterostructure and gate stack (left) and simulated band-edges [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Hall density [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Simulation of the first four HH energy levels at [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Landau fan diagram with longitudinal resistivity [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
read the original abstract

Strained Ge ($\epsilon$-Ge) and strained Si ($\epsilon$-Si) buried quantum wells have enabled advanced spin-qubit quantum processors. However, in the absence of suitable lattice-matched substrates, $\epsilon$-Ge and $\epsilon$-Si are deposited on defective, metamorphic SiGe substrates, which may impact device performance and scaling. Here an alternative platform is introduced, based on the heterojunction between unstrained Ge and a lattice-matched strained SiGe ($\epsilon$-SiGe) barrier, eliminating the need for metamorphic buffers altogether. In a structure with a 52-nm-thick $\epsilon$-SiGe barrier, a low-disorder two-dimensional hole gas is demonstrated with a high-mobility of 1.33$\times$10$^5$ cm$^2$/Vs and a low percolation density of 1.4(1)$\times$10$^1$$^0$ cm$^-$$^2$. Quantum transport shows that holes confined in the buried unstrained Ge channel have a strong density-dependent in-plane effective mass and out-of-plane $g$-factor, pointing to a significant heavy-hole$-$light-hole mixing in agreement with theory. Measurements of Zeeman spin-split levels in quantum point contacts further highlight this character, showing a two-fold larger in-plane $g$-factor in Ge than in $\epsilon$-Ge. The prospect of strong spin-orbit interaction, isotopic purification, and of hosting superconducting pairing correlations make this platform appealing for fast quantum hardware and hybrid quantum systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The manuscript introduces an alternative platform for quantum technologies based on a heterojunction between an unstrained Ge channel and a lattice-matched strained SiGe (ε-SiGe) barrier, eliminating metamorphic buffers. In a 52-nm-thick ε-SiGe structure, a low-disorder 2D hole gas is reported with mobility 1.33×10^5 cm²/Vs and percolation density 1.4(1)×10^{10} cm^{-2}. Quantum transport data show density-dependent in-plane effective mass and out-of-plane g-factor due to heavy-hole–light-hole mixing, with quantum point contact measurements indicating a two-fold larger in-plane g-factor than in strained Ge.

Significance. If the strain state and low-disorder claims hold, the platform offers a meaningful advance for scalable spin-qubit and hybrid quantum systems by removing defective metamorphic layers while retaining competitive mobility and percolation values. The g-factor and mass observations provide additional insight into hole states in Ge and support prospects for strong spin-orbit interaction.

major comments (1)
  1. [Growth and structural characterization (or equivalent methods section)] The central claim that the ε-SiGe barrier remains fully pseudomorphic to the unstrained Ge channel over 52 nm (eliminating relaxation and dislocations) is load-bearing for the lattice-matched advantage over metamorphic buffers. No reciprocal-space maps, XRD, or TEM data are presented to confirm the in-plane lattice constant or absence of misfit/threading dislocations; transport metrics alone (mobility and percolation density) are consistent with low disorder but do not exclude partial relaxation whose defects lie outside the hole wavefunction.
minor comments (1)
  1. [Abstract] The percolation density is written as 1.4(1)×10^1$$^0 in the abstract; this appears to be a typesetting error and should be corrected to 1.4(1)×10^{10} cm^{-2}.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive overall assessment. We address the single major comment below.

read point-by-point responses

  1. Referee: [Growth and structural characterization (or equivalent methods section)] The central claim that the ε-SiGe barrier remains fully pseudomorphic to the unstrained Ge channel over 52 nm (eliminating relaxation and dislocations) is load-bearing for the lattice-matched advantage over metamorphic buffers. No reciprocal-space maps, XRD, or TEM data are presented to confirm the in-plane lattice constant or absence of misfit/threading dislocations; transport metrics alone (mobility and percolation density) are consistent with low disorder but do not exclude partial relaxation whose defects lie outside the hole wavefunction.

    Authors: We agree that direct structural confirmation of the pseudomorphic state would strengthen the central claim. The ε-SiGe barrier composition was chosen to lattice-match the Ge channel by design, and the measured mobility of 1.33×10^5 cm²/Vs together with the percolation density of 1.4(1)×10^{10} cm^{-2} are among the best reported for Ge-based 2D hole gases, which is difficult to reconcile with a high density of misfit or threading dislocations intersecting the channel. Nevertheless, transport metrics alone cannot rigorously exclude relaxation whose defects lie outside the hole wavefunction. In the revised manuscript we will add high-resolution X-ray diffraction reciprocal-space maps around the (224) reflection to directly confirm that the in-plane lattice constant of the ε-SiGe barrier matches that of the Ge channel and that no detectable relaxation occurs over the 52 nm thickness. If available, cross-sectional TEM images will also be included to corroborate the absence of dislocations in the relevant region. revision: yes

Circularity Check

0 steps flagged

No circularity: purely experimental platform demonstration

full rationale

The paper reports growth of a heterostructure and direct transport measurements (mobility 1.33×10^5 cm²/Vs, percolation density 1.4(1)×10^10 cm^-2, g-factors) on a buried unstrained Ge channel with ε-SiGe barrier. No derivations, equations, or first-principles predictions are presented that could reduce to inputs by construction. Claims rest on observed data rather than any fitted parameter renamed as prediction or self-citation chain. The strain-state assumption is an unverified experimental precondition but does not create a circular derivation loop; the work is self-contained as an experimental report.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

No free parameters, invented entities or ad-hoc axioms are identifiable from the abstract. The work rests on standard assumptions of epitaxial growth, 2D carrier confinement and effective-mass theory in group-IV semiconductors.

axioms (1)
  • domain assumption Standard models of hole confinement, effective mass and g-factor in germanium heterostructures apply without modification.
    Invoked implicitly when interpreting density-dependent mass and g-factor as evidence of heavy-light hole mixing.

pith-pipeline@v0.9.0 · 5841 in / 1280 out tokens · 66493 ms · 2026-05-19T11:51:14.796243+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Strain engineering of Andreev spin qubits in Germanium

    cond-mat.mes-hall 2026-04 unverdicted novelty 6.0

    Compressive strain suppresses spin splitting in germanium Josephson junctions while tensile or unstrained heterostructures enable GHz-scale splittings for Andreev spin qubits via enhanced spin-orbit effects.

  2. Tailoring Germanium Heterostructures for Quantum Devices with Machine Learning

    cond-mat.mes-hall 2026-04 unverdicted novelty 6.0

    Localized strained silicon spikes in unstrained Ge channels, optimized via multi-objective Bayesian optimization, enhance spin-orbit interaction by up to three orders of magnitude and improve quantum-dot spin qubit qu...

Reference graph

Works this paper leans on

44 extracted references · 44 canonical work pages · cited by 2 Pith papers

  1. [1]

    N. P. De Leon, K. M. Itoh, D. Kim, K. K. Mehta, T. E. Northup, H. Paik, B. S. Palmer, N. Samarth, S. Sangtawesin, and D. W. Steuerman, Science 372, eabb2823 (2021)

    work page 2021

  2. [2]

    Burkard, T

    G. Burkard, T. D. Ladd, A. Pan, J. M. Nichol, and J. R. Petta, Reviews of Modern Physics 95, 025003 (2023)

    work page 2023

  3. [3]

    Prada, P

    E. Prada, P. San-Jose, M. W. A. De Moor, A. Geresdi, E. J. H. Lee, J. Klinovaja, D. Loss, J. Nygård, R. Aguado, and L. P. Kouwenhoven, Nature Reviews Physics 2, 575 (2020)

    work page 2020

  4. [4]

    Scappucci, C

    G. Scappucci, C. Kloeffel, F. A. Zwanenburg, D. Loss, M. Myronov, J.-J. Zhang, S. De Franceschi, G. Kat- saros, and M. Veldhorst, Nature Reviews Materials 6, 926 (2021) . 6

    work page 2021

  5. [5]

    Scappucci, P

    G. Scappucci, P. J. Taylor, J. R. Williams, T. Ginley, and S. Law, MRS Bulletin 46, 596 (2021)

    work page 2021

  6. [6]

    Paquelet Wuetz, D

    B. Paquelet Wuetz, D. Degli Esposti, A.-M. J. Zwerver, S. V. Amitonov, M. Botifoll, J. Arbiol, A. Sammak, L. M. K. Vandersypen, M. Russ, and G. Scappucci, Na- ture Communications 14, 1385 (2023)

    work page 2023

  7. [7]

    Moutanabbir, S

    O. Moutanabbir, S. Assali, A. Attiaoui, G. Daligou, P. Daoust, P. D. Vecchio, S. Koelling, L. Luo, and N. Ro- taru, Advanced Materials 36, 2305703 (2024)

    work page 2024

  8. [8]

    P. W. Deelman, L. F. Edge, and C. A. Jackson, MRS Bulletin 41, 224 (2016)

    work page 2016

  9. [9]

    Corley-Wiciak, C

    C. Corley-Wiciak, C. Richter, M. H. Zoellner, I. Zait- sev, C. L. Manganelli, E. Zatterin, T. U. Schülli, A. A. Corley-Wiciak, J. Katzer, F. Reichmann, W. M. Klesse, N. W. Hendrickx, A. Sammak, M. Veldhorst, G. Scap- pucci, M. Virgilio, and G. Capellini, ACS Applied Ma- terials & Interfaces 15, 3119 (2023)

    work page 2023

  10. [10]

    Corley-Wiciak, M

    C. Corley-Wiciak, M. Zoellner, I. Zaitsev, K. Anand, E. Zatterin, Y. Yamamoto, A. Corley-Wiciak, F. Re- ichmann, W. Langheinrich, L. Schreiber, C. Manganelli, M. Virgilio, C. Richter, and G. Capellini, Physical Re- view Applied 20, 024056 (2023)

    work page 2023

  11. [11]

    L. E. A. Stehouwer, A. Tosato, D. Degli Esposti, D. Costa, M. Veldhorst, A. Sammak, and G. Scappucci, Applied Physics Letters 123, 092101 (2023)

    work page 2023

  12. [12]

    Exploiting epitaxial strained germanium for scaling low noise spin qubits at the micron-scale,

    L. E. A. Stehouwer, C. X. Yu, B. v. Straaten, A. Tosato, V. John, D. D. Esposti, A. Elsayed, D. Costa, S. D. Oost- erhout, N. W. Hendrickx, M. Veldhorst, F. Borsoi, and G. Scappucci, “Exploiting epitaxial strained germanium for scaling low noise spin qubits at the micron-scale,” (2025), arXiv:2411.11526

    work page arXiv 2025

  13. [13]

    People, Physical Review B 34, 2508 (1986)

    R. People, Physical Review B 34, 2508 (1986)

    work page 1986

  14. [14]

    Sammak, D

    A. Sammak, D. Sabbagh, N. W. Hendrickx, M. Lodari, B. Paquelet Wuetz, A. Tosato, L. Yeoh, M. Bollani, M. Virgilio, M. A. Schubert, P. Zaumseil, G. Capellini, M. Veldhorst, and G. Scappucci, Advanced Functional Materials 29, 1807613 (2019)

    work page 2019

  15. [15]

    Lodari, N

    M. Lodari, N. W. Hendrickx, W. I. L. Lawrie, T.-K. Hsiao, L. M. K. Vandersypen, A. Sammak, M. Veldhorst, and G. Scappucci, Materials for Quantum Technology 1, 011002 (2021)

    work page 2021

  16. [16]

    Matthews and A

    J. Matthews and A. Blakeslee, Journal of Crystal Growth 27, 118 (1974)

    work page 1974

  17. [17]

    People and J

    R. People and J. C. Bean, Applied Physics Letters 47, 322 (1985)

    work page 1985

  18. [18]

    J. C. Bean, Science 230, 127 (1985)

    work page 1985

  19. [19]

    M. M. Alam, Y. Wagatsuma, K. Okada, Y. Hoshi, M. Ya- mada, K. Hamaya, and K. Sawano, Applied Physics Ex- press 12, 081005 (2019)

    work page 2019

  20. [20]

    Method for manufacturing a single het- erojunction semiconductor device and such a single het- erojunction semiconductor device,

    G. Scappucci, A. Tosato, M. F. Russ, L. E. A. Stehouwer, and A. Sammak, “Method for manufacturing a single het- erojunction semiconductor device and such a single het- erojunction semiconductor device,” (2024), international publication under the Patent Cooperation Treaty (PCT)

    work page 2024

  21. [21]

    Friesen, M

    M. Friesen, M. A. Eriksson, and S. N. Coppersmith, Applied Physics Letters 89, 202106 (2006)

    work page 2006

  22. [22]

    Degli Esposti, L

    D. Degli Esposti, L. E. A. Stehouwer, O. Gül, N. Samkharadze, C. Déprez, M. Meyer, I. N. Meijer, L. Tryputen, S. Karwal, M. Botifoll, J. Arbiol, S. V. Amitonov, L. M. K. Vandersypen, A. Sammak, M. Veld- horst, and G. Scappucci, npj Quantum Information 10 (2024), 10.1038/s41534-024-00826-9

    work page doi:10.1038/s41534-024-00826-9 2024

  23. [23]

    M. H. Zoellner, M.-I. Richard, G. A. Chahine, P. Zaum- seil, C. Reich, G. Capellini, F. Montalenti, A. Marze- galli, Y.-H. Xie, T. U. Schülli, M. Häberlen, P. Storck, and T. Schroeder, ACS Applied Materials & Interfaces 7, 9031 (2015)

    work page 2015

  24. [24]

    L. A. Tracy, E. H. Hwang, K. Eng, G. A. Ten Eyck, E. P. Nordberg, K. Childs, M. S. Carroll, M. P. Lilly, and S. Das Sarma, Physical Review B 79, 235307 (2009)

    work page 2009

  25. [25]

    Fogelholm, Journal of Physics C: Solid State Physics 13, L571 (1980)

    R. Fogelholm, Journal of Physics C: Solid State Physics 13, L571 (1980)

    work page 1980

  26. [26]

    X. Mi, T. M. Hazard, C. Payette, K. Wang, D. M. Zajac, J. V. Cady, and J. R. Petta, Physical Review B 92, 035304 (2015)

    work page 2015

  27. [27]

    Costa, L

    D. Costa, L. E. A. Stehouwer, Y. Huang, S. Martí- Sánchez, D. Degli Esposti, J. Arbiol, and G. Scappucci, Applied Physics Letters 125, 222104 (2024)

    work page 2024

  28. [28]

    Lodari, A

    M. Lodari, A. Tosato, D. Sabbagh, M. A. Schubert, G. Capellini, A. Sammak, M. Veldhorst, and G. Scap- pucci, Physical Review B 100, 041304 (2019)

    work page 2019

  29. [29]

    L. A. Terrazos, E. Marcellina, Z. Wang, S. N. Copper- smith, M. Friesen, A. R. Hamilton, X. Hu, B. Koiller, A. L. Saraiva, D. Culcer, and R. B. Capaz, Physical Review B 103, 125201 (2021)

    work page 2021

  30. [30]

    Lodari, O

    M. Lodari, O. Kong, M. Rendell, A. Tosato, A. Sam- mak, M. Veldhorst, A. R. Hamilton, and G. Scappucci, Applied Physics Letters 120, 122104 (2022)

    work page 2022

  31. [31]

    Data Repository for

    D. Costa, “Data Repository for "unstrained ge channels: a lattice-matched platform for semiconductor quantum technology",” (2025)

    work page 2025

  32. [32]

    Elsayed, M

    A. Elsayed, M. M. K. Shehata, C. Godfrin, S. Kubicek, S. Massar, Y. Canvel, J. Jussot, G. Simion, M. Mongillo, D. Wan, B. Govoreanu, I. P. Radu, R. Li, P. Van Dorpe, and K. De Greve, npj Quantum Information 10, 1 (2024)

    work page 2024

  33. [33]

    Bosco, M

    S. Bosco, M. Benito, C. Adelsberger, and D. Loss, Phys- ical Review B 104, 115425 (2021)

    work page 2021

  34. [34]

    Hole-spin qubits in germanium beyond the single-particle regime,

    A. Secchi, G. Forghieri, P. Bordone, D. Loss, S. Bosco, and F. Troiani, “Hole-spin qubits in germanium beyond the single-particle regime,” (2025), arXiv:2505.02449

    work page arXiv 2025

  35. [35]

    Hole spin qubits in unstrained Germa- nium layers,

    L. Mauro, M. J. Rodriguez, E. A. Rodriguez-Mena, and Y.-M. Niquet, “Hole spin qubits in unstrained Germa- nium layers,” (2025), arXiv:2506.04977. Supplemental Material for "Buried unstrained germanium channels: a lattice-matched platform for quantum technology" Davide Costa, Patrick Del Vecchio, Karina Hudson, Lucas E. A. Stehouwer, Alberto Tosato, Davide...

    work page arXiv 2025

  36. [36]

    The extracted root mean square surface roughness values are ∼ 0

    The lattice-matched nature of the heterojunction is confirmed by the absence of cross-hatch pattern, which is in turn very visible on the surface of the heterostructure featuring the strained Ge QW. The extracted root mean square surface roughness values are ∼ 0. 4 nm and ∼ 1. 8 nm , respectively. (a) (b) Supplementary Fig. 1. 3D AFM of (a) the Ge/SiGe str...

    work page

  37. [37]

    45 × 10− 3 (compressive strain), respectively

    99 × 10− 3 (tensile strain) and − 7. 45 × 10− 3 (compressive strain), respectively. The larger strain in the SiGe top barrier of the heterojunction arises from the lower Ge content ( 0. 8) with respect to the SiGe barrier ( 0. 83) defining the strained Ge quantum wells. Moreover, the strain map of the Ge quantum well shows signatures of the cross-hatch pat...

    work page

  38. [38]

    a a |mð = |mðm

    suggests that all its eigenvectors |φ ð can be written as |φ ð =     cH+ |m − 1ð cη+ |mð cη− |m + 1ð cH− |m + 2ð     , (11) where the cı are scalars (one scalar for each H + subband, and so on) and the |mð are the Fock states of the Landau levels operators, i.e. a a |mð = |mðm. By projecting the Hamiltonian (

    work page

  39. [39]

    (12) This Hamiltonian can be diagonalized to extract the set of scalars cÄ Ã for any m g − 2

    onto the eigenvectors ( 11), we find a block diagonal Hamiltonian that involves only the integer m: ˜Hmg 1 =     EH +(B) Eη +(B) Eη − (B) EH − (B)     + ³ 0 l2 B      2ΓH ∥ (m − 1) √ 2lBTx√ m − 2µ √ m(m + 1) 0 2Γη ∥m √ 2lBTη√ m + 1 − 2µ √ (m + 1)(m + 2) 2Γη ∥(m + 1) √ 2lBTx √ m + 2 2ΓH ∥ (m + 2)      . (12) This Hamiltonian can be diagona...

    work page

  40. [40]

    We remark that this procedure is independent of the gauge

    is diagonalized directly without any conditions on cÄ Ã. We remark that this procedure is independent of the gauge. We truncate the system to include only the levels from m = − 2 to m = 11 and vary B. This gives a total of 2 × 11 + 3 = 25 magnetic field-dependent spin-polarized Landau levels for each HH subband in the z-direction, which we now label by n g...

    work page

  41. [41]

    Del Vecchio and O

    P. Del Vecchio and O. Moutanabbir, Phys. Rev. B 110, 045409 (2024)

    work page 2024

  42. [42]

    C. G. Van de Walle and R. M. Martin, Phys. Rev. B 34, 5621 (1986)

    work page 1986

  43. [43]

    Winkler, Spin-orbit Coupling Effects in Two-Dimensional Electron and Hole Systems , Vol

    R. Winkler, Spin-orbit Coupling Effects in Two-Dimensional Electron and Hole Systems , Vol. 191 (Springer, 2003)

    work page 2003

  44. [44]

    Winkler, M

    R. Winkler, M. Merkler, T. Darnhofer, and U. Rössler, Phys. Rev. B 53, 10858 (1996) . 11

    work page 1996