Two-dimensional electron gases (2DEGs) with high mobility, engineered in semiconductor heterostructures host a variety of ordered phases arising from strong correlations, which emerge at sufficiently low temperatures. The 2DEG can be further controlled by surface gates to create quasi-one dimensional systems, with potential spintronic applications. Here we address the long-standing challenge of cooling such electrons to below 1 mK, potentially important for identification of topological phases and spin correlated states. The 2DEG device was immersed in liquid 3He, cooled by the nuclear adiabatic demagnetization of copper. The temperature of the 2D electrons was inferred from the electronic noise in a gold wire, connected to the 2DEG by a metallic ohmic contact. With effective screening and filtering, we demonstrate a temperature of 0.9 ± 0.1 mK, with scope for significant further improvement. This platform is a key technological step, paving the way to observing new quantum phenomena, and developing new generations of nanoelectronic devices exploiting correlated electron states.
Cooling electrons into the microkelvin temperature range is of interest both for practical purposes and fundamental studies, but current demonstrations are limited to small, specific devices. Here, the authors achieve sub-millikelvin temperatures in a large-area, two-dimensional electron gas.
RCUK | Engineering and Physical Sciences Research Council (EPSRC)https://doi.org/10.13039/501100000266EP/K004077/1LevitinLev V.EC | EU Framework Programme for Research and Innovation H2020 | H2020 Priority Excellent Science | H2020 European Research Council (H2020 Excellent Science - European Research Council)https://doi.org/10.13039/100010663824109LevitinLev V.publisher-imprint-nameNature Portfoliovolume-issue-count1issue-article-count666issue-toc-levels0issue-pricelist-year2022issue-copyright-holderThe Author(s)issue-copyright-year2022article-contains-esmYesarticle-numbering-styleUnnumberedarticle-registration-date-year2022article-registration-date-month1article-registration-date-day17article-toc-levels0toc-levels0volume-typeRegularjournal-productNonStandardArchiveJournalnumbering-styleUnnumberedarticle-grants-typeOpenChoicemetadata-grantOpenAccessabstract-grantOpenAccessbodypdf-grantOpenAccessbodyhtml-grantOpenAccessbibliography-grantOpenAccessesm-grantOpenAccessonline-firstfalsepdf-file-referenceBodyRef/PDF/41467_2022_Article_28222.pdfpdf-typeTypesettarget-typeOnlinePDFissue-typeRegulararticle-typeOriginalPaperjournal-subject-primaryScience, Humanities and Social Sciences, multidisciplinaryjournal-subject-secondaryScience, Humanities and Social Sciences, multidisciplinaryjournal-subject-secondaryScience, multidisciplinaryjournal-subject-collectionScience (multidisciplinary)open-accesstrueIntroduction
Two-dimensional electron gases (2DEGs), created at a GaAs–AlGaAs heterojunction and grown by molecular-beam epitaxy (MBE), have been the building block for the study of low-dimensional physics over the past few decades. This requires the high crystal perfection of MBE growth and the use of modulation doping (see ref. 1 for a review of the state of the art).
The lowest possible electron carrier temperature coupled to the highest sample quality is key to discovering and elucidating new ground states in such systems. The observation of the integer quantum Hall effect2 and shortly thereafter the fractional quantum Hall effect (FQHE)3 has led to extensive studies of two-dimensional electron systems in perpendicular magnetic field B. The lowest Landau level (LL), with filling factors ν < 1 (here ν = nh/eB is determined by the carrier density n2D, Planck’s constant h and elementary charge e) for the lowest spin branch is described by composite fermion (CF) theory4. It exhibits a rich variety of quantum states: fractional quantum Hall states for odd fractions, while ν = 1/2 is a composite fermion Fermi sea5. A key question is the role of residual CF interactions in giving rise to exotic states6. A Wigner solid is found for small ν, which also exhibits re-entrance with a maximum melting temperature of 50 mK above the ν = 1/3 FQHE state7. Interactions are also believed to play a role at other fractional fillings such as ν = 4/11, which has a measured activation energy of 5 mK8.
The nature of quantum states in the second LL is in many cases not resolved, and this 2D electron system is of particular interest for study to lower temperatures than hitherto achieved. We discuss a few examples here. Even denominator states forming in the second Landau level are thought to arise from pairing of composite fermions. Confirmation of exact quantization of the 5/2 state first discovered by Willett et al.9 required electron temperatures as low as 8 mK10. Nevertheless the topology of the ν = 5/2 ground state (Pfaffian, anti-Pfaffian or other) remains an open question. Observations are highly sensitive to disorder11,12. It has been shown that a quantum phase transition to a quantum Hall nematic phase can be induced by applied pressure in both the 5/2 and 7/2 state13. The importance of such topological ground states is reinforced by the proposal that they offer a route to topological quantum computing using non-Abelian anyons14. A non-Abelian phase is also predicted for the 12/5 state.
Many of these states have small energy gaps. The 12/5 state has a gap of order 30 mK, which exhibits non-trivial tilt dependence of the magnetic field. Studies of the energy gaps at other fractional fillings15,16 found that the gaps of the ν = 2 + 6/13 (10 mK) and 2+2/5 (80 mK) states do not follow the CF hierarchy, leading to the proposal of a parton state17. The 2+3/8 state also has a reported gap of only 10 mK. Furthermore, there are indications of novel topological order in the upper spin branch of the second LL, at fillings 3+1/3 (37 mK) and 3+1/5 (104 mK)18.
New states arising from strong correlations in the low-dimensional electron system are also expected in low magnetic fields. Electrostatic confinement of the 2DEG to form 1D wires and zero-dimensional (0D) quantum dots has led to studies of the many-body physics of Luttinger liquids19,20 and the Kondo effect21,22. Furthermore the conduction electrons in GaAs-based 2DEGs are expected to couple to the nuclear spins, undergoing a collective ferromagnetic transition mediated by RKKY interactions which is expected at electron temperatures in the mK regime23. The hyperfine interaction is also predicted24 to create new topological phases in a semiconductor 1D wire coupled to a superconductor and magnetic moments. A transition from metal- to insulator-like behaviour with decreasing temperature has been reported in a 2D hole gas in GaAs25, potentially related to many-body localization. Recently, Wigner crystal-like zigzag chains have been discovered in 1D wires at low electron densities26. At lower temperatures there is potential for quantum communication through 1D spin chains, coupling high fidelity quantum dot qubits27.
Advances in cryogen-free technologies now make low millikelvin temperature platforms widely accessible28, with prospects of commercial solutions in the microkelvin range29,30. However, cooling electrons in semiconductor devices is challenging. In situ cooling via the electron–phonon coupling is severely limited by its strong temperature dependence, while the electrons are susceptible to heat leaks arising from electromagnetic noise. Electromagnetic shielding and ultra-sensitive measurement schemes with low or no electrical excitation31–33 are common strategies for minimising the heat leak. Cooling via the electrical leads has been widely implemented10,33–36, commonly employing immersion in liquid helium. For 2D electrons, this has so far been limited, for many years, to 4 mK, see ref. 34 and references therein. Primary thermometers within samples confirmed 2D electron temperatures of 6 mK33,36.
In response to the challenge to cool nano-electronic devices to ultra-low temperatures (ULT), new approaches have been introduced which focus on refrigeration by nuclear demagnetisation of each lead connected to the sample37. More recently this has been extended by incorporating nuclear refrigeration elements into the device itself. For this on-chip cooling the Coulomb blockade device used is itself a thermometer, providing an accurate measurement of the electron temperature in mesoscopic metallic islands as low as 0.5 mK38,39.
In contrast, and complementary to this approach, the focus of our work is the cooling of a relatively large area (8 mm2) 2DEG. Our approach to cooling is motivated by the requirement for flexibility to cool a wide range of devices in different sample environments. We report an ultimate electron temperature in the 2DEG of 0.9 ± 0.1 mK, cooled in a 3He immersion cell, achieved after improvements in electromagnetic shielding. The 2DEG device is immersed in liquid 3He and cooled through metallic ohmic contacts, which are coupled to the 3He via compact heat exchangers made of sintered Ag powder. Cooling of devices by the immersion cell technique, as compared to on-chip cooling, is remarkably versatile due to modularity of the design. The cooling is provided by a physically remote copper nuclear adiabatic demagnetization stage module, so the magnetic field applied to the sample can be independently adjusted. A broad class of devices can potentially be cooled in this way, and in some cases the 3He liquid pressure can be used as an experimental control parameter40.
The effective cooling of electrons in the 2DEG relies on both reducing the resistance (and hence thermal resistance via the Wiedemann-Franz law) of the 2DEG and ohmic contacts, and reducing the heat leak to the device. We fabricated a high-mobility GaAs-based 2DEG with sub-1 Ω AuNiGe ohmic contacts, see Fig. 1a, b, and reached a 1 fW level of heat leak to the 2DEG, see Fig. 1c, through designing a metal cell that is tight to photons over a wide frequency range, combined with extensive filtering of all electrical lines, see Fig. 2.
Cooling electrons in a 2DEG in 3He immersion cells.
a Schematic of the sample: the 2DEG is created at a GaAs/AlGaAs heterojunction, and electrical connections are provided by AuNiGe ohmic contacts. b Top view of the 4 mm × 4 mm device cooled to 1 mK by the I± and V± ohmic contacts, also used to measure the electrical resistances in the device. A noise thermometer (NT) is connected to the NT ohmic contact to probe the electron temperature Te. c Progressive cooling of the 2DEG device in plastic (d) and metal (e) immersion cells equipped with different low-pass filters. Te is estimated from Eq. (4), assuming that the 3He bath is at temperature Tbath = 0.3 mK.
Experimental platform for cooling low-dimensional electrons.
a Schematic diagram of the metal immersion cell. Silver sinter on the silver base thermalises the liquid 3He to the nuclear demagnetisation refrigerator. The copper lid holds the 2DEG device, the sinter heat exchangers for cooling the I± and V± lines, the 2DEG noise thermometer (NT), and a red light-emitting diode (LED) for illuminating the 2DEG. The NT and its heater were connected using wedge-bonded resistive Au wires and superconducting Al wires acting as thermal breaks62. Another NT is mounted on the nuclear stage (NS) to measure the 3He bath temperature Tbath. The base and the lid sealed with In comprise a He-tight and photon-tight enclosure, complete with a Ag sinter plug in the filling capillary. All incoming electrical connections are fitted with low-pass filters of three types: b discrete RC filters with R = 1 kΩ and C = 10 nF, c insulating Cu powder filters63, and, d Ag epoxy filters with 1 m-long NbTi wires embedded in conducting Ag epoxy. RF-tight coaxial connections were used between these filters and the walls of the immersion cell. Ag epoxy filters with multiple wires, compatible with SQUID NT readout, were used on the SQUID lines. e The inside of the copper lid of the cell, with the Au and Al wires highlighted and the LED shown schematically. The NT is housed in a Nb shield, shown here with the top removed.
A key feature of the experiment was the choice of thermometer. We use a current sensing noise thermometer (NT) as an external module attached to the 2DEG via an ohmic contact. A SQUID current sensor is used to read out the voltage fluctuations across a gold wire41. Due to the fundamental Nyquist relation the NT operates over 5 orders of magnitude in temperature with a single point calibration. It can be configured as a primary thermometer42, and as such underpins the recent redefinition of the Kelvin.
The noise thermometer dissipates no power in the device, but we must take account of both the inevitable residual heat leak into the thermometer, and the potential parallel cooling channel by direct coupling to the 3He in which it is immersed. We have constructed a thermal model which demonstrates that at the lowest temperatures the NT accurately reflects the electron temperature in the 2DEG. The ability to alter the resistance of the 2DEG in situ by optical illumination, due to persistent photoconductivity, was used to experimentally validate the model.
ResultsExperimental setup
Figure 2a is a schematic of the metal immersion cell, which combines a 3He bath with a photon-tight environment. Filtered connections to room-temperature electronics are provided to probe electrical resistances in the device and to operate heaters in thermal transport experiments. Two NTs are connected to the 2DEG and the demagnetisation refrigerator. Illumination by a red light-emitting diode (LED) facing the 2DEG allows us to increase the 2D carrier density and mobility43. Preliminary studies were conducted in a plastic cell shown in Fig. 1d, which was based on ref. 34. See Supplementary Information (SI) for further details.
Electrical transport measurements
Figure 3a shows that the resistance of both the 2DEG and ohmic contacts exhibit a weak temperature dependence. These transport measurements were performed down to 1 mK with no discernible heating using conventional room temperature electronics without a screened room; this demonstrates the effectiveness of the filtering. A change in the resistance was observed between 0.4 and 0.8 K, that has been identified in similar samples44 to be a superconducting transition in the ohmic contacts.
Transport measurements of 2DEG and ohmic contacts.
a Sheet resistance R□ of the 2DEG (red) and the resistance including 2DEG and two ohmic contacts (cyan) measured down to 1 mK. Weak temperature dependence of the resistance is observed in both the 2DEG and the ohmic contacts below the superconducting transition in the contacts which occurs at Tc = 0.6 ± 0.2 K. The electron temperature Te is determined from Eq. (4). b Thermal conductance GΣ between 2DEG NT and 3He bath before (cyan) and after (red) illuminating the 2DEG, modelled as a sum of two parallel channels (c): Ge via electrons in the 2DEG and ohmic contacts, sensitive to illumination, and GHe due to immersion of the gold wires in liquid helium. Assuming the Wiedemann-Franz (WF) law applies in the 2DEG and ohmic contacts, after illumination the electron channel GeWF=L0T/RNTtoIV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{\,{{{{\rm{e}}}}}}^{{{{{\rm{WF}}}}}\,}={L}_{0}T/{R}_{{{{{{{{\rm{NT}}}}}}}}\,{{{{{{{\rm{to}}}}}}}}\,{{{{{{{\rm{IV}}}}}}}}}$$\end{document} alone exceeds the measured GΣ. We attribute the violation of the WF law to superconductivity within the ohmic contacts. The effective electrical resistance Re is obtained by fitting GΣ(T) in (b) to Eq. (1). d In the partially-superconducting ohmic contact model the thermal conductance at T ≪ Tc is determined by the non-superconducting channel Rc*\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${R}_{c}^{* }$$\end{document} according to the WF law. e Response of the 2DEG NT to heating the 2DEG, the non-local thermal conductance GN, is obtained analogously to GΣ and fitted to Eq. (2). Dashed lines show the X4T4 term in Eq. (2), unimportant below TNT = 3 mK. f A lumped-element thermal model described by Eqs. (3) and (4), in which Q̇e\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}$$\end{document} is assumed to be applied in the centre of the 2DEG, splitting Ge into GNT and GIV, each including the 2DEG and the ohmic contacts. g The ratio GΣ/GN describes the fraction of Q̇e\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}$$\end{document} that flows into the NT. The T → 0 limit α of GΣ/GN governs GNT and GIV in (f). Error bars in b, e, g represent s.d.
Thermal transport measurements and thermal model
The use of an external thermometer in our setup requires an investigation of the heat flow in the system in order to relate the electron temperature in the 2DEG Te to the 2DEG NT temperature TNT. We consider two sources of heat: Q̇NT=Q̇NT0+Q̇NTJ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{NT}}}}}}={\dot{Q}}_{\,{{{{\rm{NT}}}}}\,}^{0}+{\dot{Q}}_{\,{{{{\rm{NT}}}}}}^{{{{{\rm{J}}}}}\,}$$\end{document} applied to the NT and Q̇e=Q̇e0+Q̇eJ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}={\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}+{\dot{Q}}_{{{{{\rm{e}}}}}}^{{{{{\rm{J}}}}}\,}$$\end{document} uniformly generated within the 2DEG device. Each includes a residual heat leak (Q̇NT0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{\,{{{{\rm{NT}}}}}\,}^{0}$$\end{document} and Q̇e0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}$$\end{document}); additionally, to study the thermal response Joule heating Q̇NTJ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{\,{{{{\rm{NT}}}}}}^{{{{{\rm{J}}}}}\,}$$\end{document} is applied using the NT heater and Q̇eJ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}^{{{{{\rm{J}}}}}\,}$$\end{document} by driving a current between I+ and I− ohmic contacts.
Figure 3b shows the thermal conductance GΣ(TNT)=dQ̇NT/dTNT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{{{\Sigma}}}({T}_{{{{{\rm{NT}}}}}})=d{\dot{Q}}_{{{{{\rm{NT}}}}}}/d{T}_{{{{{\rm{NT}}}}}}$$\end{document} between the 2DEG NT and the 3He bath, as inferred from measurements of TNT as a function of Q̇NT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{NT}}}}}}$$\end{document} at a constant 3He bath temperature Tbath. Figure 3c illustrates the two parallel channels of GΣ = Ge + GHe, Ge via electrons in the 2DEG and ohmic contacts, and the Kapitza boundary conductance GHe between the NT assembly and the 3He bath45. The illumination of the 2DEG increases Ge, allowing us to separate it from GHe. To reduce GHe the heat exchangers in the immersion cell were plated with approx. 30 μmol/m24He coverage before loading 3He46.
In the initial stage of this research we assumed that both the 2DEG and the ohmic contacts obey the Wiedemann-Franz (WF) law, ubiquitous in electronic transport in the T → 0 limit47–52. The WF law predicts that GeWF=L0T/RNTtoIV\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{\,{{{{{{\rm{e}}}}}}}}^{{{{{{{\rm{WF}}}}}}}\,}={L}_{0}T/{R}_{{{{{{{{\rm{NT}}}}}}}}\,{{{{{{{\rm{to}}}}}}}}\,{{{{{{{\rm{IV}}}}}}}}}$$\end{document}, where RNT to IV = 27 Ω (6 Ω) is the electrical resistance from the NT ohmic contact to I± and V± connected in parallel, measured before (after) illumination, where L0=π2kB2/3e2=2.44×10−8WΩ/K2\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${L}_{0}={\pi }^{2}{k}_{{{{{{{{\rm{B}}}}}}}}}^{2}/3{e}^{2}=2.44\times 1{0}^{-8}\,{{{{{{{\rm{W}}}}}}}}{{\Omega }}/{{{\rm{K}}}}^{2}$$\end{document} is the Lorenz number47. Figure 3b shows that GΣ<GeWF\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{{{\Sigma }}} \; < \; {G}_{\,{{{{{{\rm{e}}}}}}}}^{{{{{{{\rm{WF}}}}}}}}$$\end{document} after illumination, a clear violation of the WF law that we attribute to superconductivity53 in the ohmic contacts44. To the leading order the change in GΣ due to illumination is ΔGΣ = ΔGe ∝ T; to explain this observation we propose a model of partially-superconducting ohmic contacts, see Fig. 3d, such that well below the superconducting transition temperature Tc, Ge = L0T/Re is described by a resistance Re ≥ RNT to IV. The data are well described byGΣ(T)=L0T/Re+K2T2+K3T3,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{{{\Sigma }}}(T)={L}_{0}T/{R}_{{{{{\rm{e}}}}}}+{K}_{2}{T}^{2}+{K}_{3}{T}^{3},$$\end{document}with Re = 59 Ω (27 Ω) before (after) illumination and illumination-independent GHe = K2T2 + K3T3. Similarly the response to Q̇e\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}$$\end{document} was characterised in terms of the non-local thermal conductance GN(TNT)=dQ̇e/dTNT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{{{{{\rm{N}}}}}}({T}_{{{{{\rm{NT}}}}}})=d{\dot{Q}}_{{{{{\rm{e}}}}}}/d{T}_{{{{{\rm{NT}}}}}}$$\end{document}, see Fig. 3e, which was found to followGN(T)=GΣ(T)+X4T4/α,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{{{{{\rm{N}}}}}}(T)=\left({G}_{{{\Sigma }}}(T)+{X}_{4}{T}^{4}\right)/\alpha ,$$\end{document}with an illumination-independent term X4, but with different values of α before and after illumination. The X4T4 term describes an additional cooling mechanism in the 2DEG or ohmic contacts, insignificant below TNT = 3 mK, where X4T4 ≪ GΣ(T). In this low temperature regime, Eq. (2) reduces to GN = GΣ/α where α is the fraction of Q̇e\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}$$\end{document} that flows towards the NT ohmic contact, as shown in Fig. 3f, g. ThenQ̇NT+αQ̇e=∫TbathTNTGΣ(T)dT=L02ReTNT2−Tbath2+K23TNT3−Tbath3+K34TNT4−Tbath4.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{NT}}}}}}+\alpha {\dot{Q}}_{{{{{\rm{e}}}}}}= \,\int\nolimits_{{T}_{{{{{\rm{bath}}}}}}}^{{T}_{{{{{\rm{NT}}}}}}}\,\,{G}_{{{\Sigma }}}(T)\,dT=\frac{{L}_{0}}{2{R}_{{{{{\rm{e}}}}}}}\left({T}_{\,{{{{\rm{NT}}}}}}^{2}-{T}_{{{{{\rm{bath}}}}}\,}^{2}\right)\\ \, +\frac{{K}_{2}}{3}\left({T}_{\,{{{{\rm{NT}}}}}}^{3}-{T}_{{{{{\rm{bath}}}}}\,}^{3}\right)+\frac{{K}_{3}}{4}\left({T}_{\,{{{{\rm{NT}}}}}}^{4}-{T}_{{{{{\rm{bath}}}}}\,}^{4}\right).$$\end{document}The combined heat leak Q̇NT0+αQ̇e0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{\,{{{{\rm{NT}}}}}\,}^{0}+\alpha {\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}$$\end{document} can be inferred from the difference between TNT and Tbath in the absence of Joule heating. We assume Tbath to be equal to the fridge temperature, as justified in SI. By replacing the gold thermal link between the NT ohmic contact and the NT with an aluminium one, we measured Q̇NT0=0.08±0.02\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{\,{{{{\rm{NT}}}}}\,}^{0}=0.08\pm 0.02$$\end{document} fW in a separate experiment (see SI), allowing us to obtain Q̇e\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}$$\end{document} directly, see Fig. 4a. Before illumination by the LED the heat leak was found to be Q̇e0=0.7±0.1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}=0.7\pm 0.1$$\end{document} fW. After illumination there was a dramatic increase of Q̇e0,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0},$$\end{document} followed by a slow relaxation, consistent with slow recombination processes in the heterostructure54. On a time scale of a month, a new stable level Q̇e0=1.9±0.4\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}=1.9\pm 0.4$$\end{document} fW was reached, higher than before the illumination.
Heat leak to the electrons in 2DEG and inferred electron temperature before (blue) and after (red) illumination.
a The heat leak Q̇e0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}$$\end{document} to the 2DEG, determined from Eq. (3). Error bars and thick lines through the data represent s.d. Low-excitation AC transport measurements (filled triangles) do not increase the heat leak significantly over the value it takes when no room-temperature electronics are connected to the experiment (open triangles). After illumination with the LED the heat leak increases dramatically and then gradually decays, stabilising a month later; the orange band that highlights the exponential decay is a guide to the eye. b, c The triangles represent the measurements of TNT vs Tbath corresponding to the Q̇e0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}$$\end{document} shown in (a); the solid lines are fits to Eq. (3). The dotted lines show the electron temperature Te in the 2DEG estimated from Eq. (4). The first 40 days after illumination are excluded from (c) due to evolving Q̇e0\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}$$\end{document}.
To estimate the minimum electron temperature in the 2DEG we use the simple model shown in Fig. 3f, in which Q̇e\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}$$\end{document} is a point source of heat in the middle of the device and additional cooling associated with the X4T4 term in Eq. (2) is ignored:Te2=(1−α)Tbath2+αTNT2+2α(1−α)ReQ̇e/L0.\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${T}_{\,{{{{\rm{e}}}}}\,}^{2}=(1-\alpha )\,{T}_{\,{{{{\rm{bath}}}}}\,}^{2}+\alpha \,{T}_{\,{{{{\rm{NT}}}}}\,}^{2}+2\alpha (1-\alpha )\,{R}_{{{{{\rm{e}}}}}}{\dot{Q}}_{{{{{\rm{e}}}}}}/{L}_{0}.$$\end{document}Figures 4b and c show the inferred electron temperature Te, demonstrating that we have cooled the electrons to 1 mK, despite the unexpected hindrance posed by the superconductivity in the ohmic contacts.
Discussion
We demonstrate that at 1 mK the dominant cooling of the 2D electrons is by the ohmic contacts, which have a thermal conductance of Ge ~ 10−9 W/K2 × T. The electron–phonon coupling in high-mobility 2DEGs has been observed55 to follow Q̇e-ph=61\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e-ph}}}}}}=61\,$$\end{document} eV/s K5×N(Te5−Tph5)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${}^{5}\times N({T}_{\,{{{{\rm{e}}}}}}^{5}-{T}_{{{{{\rm{ph}}}}}\,}^{5})$$\end{document} at phonon temperatures of Tph = 0.3−0.5 K, where N is the number of electrons. For our device this leads to a thermal conductance Ge-ph=dQ̇e-ph/dTe~10−6\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{{{{{\rm{e-ph}}}}}}=d{\dot{Q}}_{{{{{\rm{e-ph}}}}}}/d{T}_{{{{{\rm{e}}}}}} \sim 1{0}^{-6}$$\end{document} W/K5 × T4. Assuming this high-temperature power law extrapolates to ULT, the cross-over Ge = Ge-ph occurs near 100 mK and at 1 mK Ge-ph/Ge = 10−6, so cooling via electron–phonon coupling alone would only achieve Te = 20 mK (70 mK) for Tph = 0.3 mK and Q̇e0~1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0} \sim 1$$\end{document} fW (500 fW).
The X4T4 term in Eq. (2) points towards an additional cooling channel competing with Ge above 3 mK, such as enhanced Ge-ph ≫ 10−6 W/K5 × T4 at ULT; alternatively Ge(T), GNT(T) and GIV(T), see Fig. 3f, may deviate from the G ∝ T behaviour, if the effective contact resistance Rc*\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${R}_{c}^{* }$$\end{document} (see Fig. 3d) is not constant at ULT. Ignoring these effects in Eq. (4) potentially overestimates Te for the preliminary experiments shown in Fig. 1c, but the ultimate electron temperature we report and the qualitative observation of the reduction of the heat leak Q̇e\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}$$\end{document} due to the electromagnetic shielding are robust. In addition to the X4T4 term, a more detailed thermal model, beyond the scope of this work, should take into account the distribution of Q̇e\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{\rm{e}}}}}}$$\end{document} across the 2DEG and 2D and 3D heat flow in the 2DEG and ohmic contacts.
We note that after illumination the thermal resistances are dominated by the ohmic contacts, since Re ≫ RNT to IV. Therefore the thermal conductance through the I± and V± contacts in parallel, GIV = Ge/α, (see Fig. 3f), is only (1 − α)/α = 2.1 times higher than that through the NT ohmic, GNT = Ge/(1 − α), despite being 15 times larger in circumference. This suggests that the large I± ohmics were rendered thermally inactive due to low thermal conductance along them, as the gold wires were bonded 0.5 mm away from their front edges, Fig. 2e, to prevent damage to the 2DEG adjacent to these ohmics. The larger value of (1 − α)/α = 3.2 before illumination is consistent with the I± ohmics being more active when the 2DEG resistance is higher.
A natural question to ask is whether even lower electron temperatures can be achieved. First, the heat leak may be reduced by optimising the device geometry (for some sources Q̇e0∝area\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}\propto \,{{\mbox{area}}}\,$$\end{document}) and further filtering. We recognise that our NT readout scheme limited the filters between the NT and the SQUID to a two-way Ag epoxy filter with a high cut-off frequency. Improved electron cooling is expected in experiments compatible with heavy filtering of every measurement line.
Another approach is to improve the ohmic contacts. Their thermal conductance can be increased if the superconductivity is suppressed with a magnetic field or by a change of recipe. To optimise the performance of the large partially-superconducting ohmic contacts used in this work, a thick gold film could be evaporated on the annealed top surface or a dense network of gold wires could be bonded along the front edge of the contacts. We estimate electron temperatures of 0.4–0.6 mK, if the above steps are combined with the optimised fridge performance (see SI).
Important future steps include characterising and mitigating the heating associated with surface gates used to define mesoscopic samples, measuring the electron temperature in the device directly33,56, and extending the techniques presented here to high magnetic fields. Based on the thermal measurements at 6–100 mK33, the quantum Hall states will cool into the microkelvin regime in our environment. The immersion cooling can efficiently thermalise a variety of degrees of freedom in condensed matter systems, offering a promising path to control the decoherence problem in superconducting electronics57.
In conclusion, the demonstrated cooling of a large area two-dimensional electron system to below 1 mK constitutes a significant technical breakthrough in quantum nanoelectronics. Fundamental studies of lower dimensional electron devices such as quantum wires and quantum dots in the microkelvin regime can be achieved if the heat leak to the 2DEG, generated by the surface gates, is kept at the fW level. Likewise, the addition of a strong magnetic field, which is straightforward to implement in our modular design, coupled with samples of the highest quality, is poised to contribute to the understanding of novel ground states of two-dimensional electron systems at fractional LL fillings, arising from CF interactions. Thus the microkelvin regime is opened up for the study of a rich and diverse array of strongly correlated quantum systems, with high potential for future discovery.
Methods2DEG sample
The 4 mm × 4 mm sample shown in Fig. 1b was fabricated using wafer W476; details of its MBE growth and the subsequent fabrication of AuNiGe ohmic contacts are given elsewhere44, together with the critical field measurements of the superconducting state below 1 K. The as-grown 2D carrier density and mobility are n2D ≈ 2 × 1011 cm−2 and μ ≈ 1 × 106 cm2/Vs giving a calculated sheet resistance of R□ ≈ 31 Ω. After illumination these quantities are 3.3 × 1011 cm−2, 3 × 106 cm2/Vs and 6 Ω. Illumination was performed at 1.5 K with the immersion cell evacuated. The LED was driven with currents up to 2 mA for 1 min, by which time R□ saturated.
Filters
The Ag epoxy filters, shown in Fig. 2d, combine features of several designs32,58,59. Metre-long superconducting NbTi wires were coiled around threaded sterling silver tubes and encapsulated in conducting Ag epoxy. In the resulting lossy coaxial the dissipation occurs in the outer conductor which is in direct metallic contact with the refrigerator, ensuring good thermalisation. These filters with a 100 MHz cut-off were combined with the Cu powder filters with 2 MHz cut-off, and the 16 kHz low-pass RC filters shown in Fig. 2b. To prevent high-frequency leaks the Ag epoxy filters were connected to the immersion cell via semi-rigid coaxial cables with threaded connectors.
A version of the Ag epoxy filter using a NbTi twisted pair (instead of a single wire) embedded in Ag epoxy was inserted between 2DEG NT and its SQUID sensor. This design was chosen for having zero DC resistance and for inducing sufficiently low amounts of noise in the SQUID so as to not mask the NT signal. This filter was mounted directly on top of the immersion cell and sealed to the cell wall with In to maintain a photon-tight enclosure. Another multi-line Ag epoxy filter at 1.5 K was inserted between the SQUID sensor and its room temperature electronics. All remaining filters were mounted at the mixing chamber plate of the dilution refrigerator.
Noise thermometry
Both NTs were read out with integrated 2-stage SQUID current sensors60 and were calibrated against a primary magnetic field fluctuation thermometer at 10–100 mK61.
Electrical transport measurements
were performed using a mains-powered Stanford Research Systems SR-124 analogue lock-in amplifier using its internal oscillator. Below 20 mK measurement currents between 1 and 10 nA were used to ensure sub-fW Joule heating.
Thermal measurements
The thermal conductance GΣ(TNT)=dQ̇NT/dTNT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{{{\Sigma }}}({T}_{{{{{\rm{NT}}}}}})=d{\dot{Q}}_{{{{{\rm{NT}}}}}}/d{T}_{{{{{\rm{NT}}}}}}$$\end{document} was inferred from pairs of measurements of TNT = Ta,b at different levels of Q̇NTJ=Q̇a,b\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{\,{{{{\rm{NT}}}}}}^{{{{{\rm{J}}}}}\,}={\dot{Q}}_{a,b}$$\end{document} at constant Tbath as GΣ(Ta+Tb)/2=Q̇b−Q̇a/Tb−Ta\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{{{\Sigma }}}\left(({T}_{a}+{T}_{b})/2\right)=\left({\dot{Q}}_{b}-{\dot{Q}}_{a}\right)/\left({T}_{b}-{T}_{a}\right)$$\end{document}, Fig. 5a. Identical techniques were used for measuring GN(T). Below 3 mK the 2DEG NT responded slowly, Fig. 5b, and Tbath could not be kept constant sufficiently long due to the single-shot operation of the demagnetisation cooling, so Ta was inferred from measurements before and after Tb, see Fig. 5c. Since similar slow relaxation was observed when NT was isolated from the 2DEG, see Fig. 5b, we conclude that its origin is the large heat capacity of the NT itself, and all other elements of the system thermalise faster. This implies that when dTNT/dt = 0, see Fig. 5d, the entire system is in a steady state, allowing the stationary Eq. (3) to be used to extract the heat leak. The slow NT response practically limits the thermal measurements to TNT > 1.2 mK.
Measurements of the thermal conductance and heat leak.
a The thermal conductance GΣ was inferred from the step-wise measurement of TNT(Q̇NTJ)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${T}_{{{{{\rm{NT}}}}}}({\dot{Q}}_{\,{{{{\rm{NT}}}}}}^{{{{{\rm{J}}}}}\,})$$\end{document} when Tbath was kept constant. b The slow response of TNT to heaters and Tbath below 3 mK required a modified GΣ measurement shown in (c). d Heat leak was inferred via stationary Eq. (3) which applies when dTNT/dt = 0.
In the early experiments shown in Fig. 1c, TNT was too high to ignore the X4T4 term in GN. Here the thermal model given by Eqs. (3) and (4) is inaccurate. The heat leak was inferred from Q̇e0=∫TbathTNTGN(T)dT\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\dot{Q}}_{{{{{{{{\rm{e}}}}}}}}}^{0}=\int\nolimits_{{T}_{{{{{\rm{bath}}}}}}}^{{T}_{{{{{\rm{NT}}}}}}}{G}_{{{{{\rm{N}}}}}}(T)\,dT$$\end{document}. To account for the use of different 2DEG devices of the same design and the lack of 4He plating GN(T) was measured separately in each experiment, see SI for further information.
Acknowledgements
We thank P. Bamford, R. Elsom, I. Higgs and J. Taylor for mechanical support, and V. Antonov for wire bonding. Cu powder filters designed and fabricated by A. Iagallo and M. Venti were instrumental in the early stages of this work. We acknowledge fruitful discussions with C. Ford, R. Haley, P. Meeson, P. See and J. Waldie. This research was supported by EPSRC Programme grant EP/K004077/1 and the EU H2020 European Microkelvin Platform EMP, Grant No. 824109.
Author contributions
J.S., A.D.C., J.N., A.J.C. and L.V.L. conceived the experiments; the wafer was grown by D.A.R. and I.F.; A.D.C., G.C. and J.T.N. designed and fabricated the 2DEG device; A.D.C. and L.V.L. designed and constructed the plastic immersion cell, with input from J.N., A.J.C., J.T.N. and J.S. L.V.L. and H.vd.V. designed and constructed the metal immersion cell and filters and performed ultra-low temperature measurements in both cells, with input from J.N., A.J.C., J.T.N. and J.S.; S.D. and L.V.L. measured the heat leak to the 2DEG; T.T., M.L., A.J.C. and J.T.N. characterised the superconductivity in the ohmic contacts; L.V.L. analysed the data; L.V.L., J.S. and J.T.N. wrote the manuscript with critical input from all authors.
Peer reviewPeer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Data availability
The thermal conductance, electrical resistance, heat leak and thermal time constant data obtained in this work are available in Figshare at https://doi.org/10.6084/m9.figshare.17057063.
Competing interests
The authors declare no competing interests.
ReferencesUmanskyVMBE growth of ultra-low disorder 2DEG with mobility exceeding 35 × 106 cm2/Vs2009311165816612009JCrGr.311.1658U1:CAS:528:DC%2BD1MXktFShtbY%3DKlitzingKVDordaGPepperMNew method for high-accuracy determination of the fine-structure constant based on quantized Hall resistance1980454944971980PhRvL..45..494KTsuiDCStormerHLGossardACTwo-dimensional magnetotransport in the extreme quantum limit198248155915621982PhRvL..48.1559T1:CAS:528:DyaL38XktFGhu7s%3DJainJKComposite-fermion approach for the fractional quantum Hall effect1989631991989PhRvL..63..199J1:STN:280:DC%2BC2sfos1ahug%3D%3D10040805KamburovDWhat determines the Fermi wave vector of composite fermions?20141131968012014PhRvL.113s6801K1:CAS:528:DC%2BC2cXitFGjurrM25415915JainJKComposite fermion theory of exotic fractional quantum Hall effect2015639622015ARCMP...6...39JMaMKThermal and quantum melting phase diagrams for a magnetic-field-induced Wigner solid20201250366012020PhRvL.125c6601M1:CAS:528:DC%2BB3cXhsFCntbnK32745416PanWBaldwinKWWestKWPfeifferLNTsuiDCFractional quantum Hall effect at Landau level filling ν = 4/112015910413012015PhRvB..91d1301PWillettRObservation of an even-denominator quantum number in the fractional quantum Hall effect198759177617791987PhRvL..59.1776W1:CAS:528:DyaL1cXhsVChsw%3D%3D10035326PanWExact quantization of the even-denominator fractional quantum Hall state at ν = 5/2 Landau level filling factor199983353035331999PhRvL..83.3530P1:CAS:528:DyaK1MXmslOmsro%3DMrossDFOregYSternAMargalitGHeiblumMTheory of disorder-induced half-integer thermal Hall conductance20181210268012018PhRvL.121b6801M1:CAS:528:DC%2BC1MXltFSgsbc%3D30085751WangCVishwanathAHalperinBITopological order from disorder and the quantized Hall thermal metal: possible applications to the ν = 5/2 state2018980451122018PhRvB..98d5112W1:CAS:528:DC%2BC1MXltVGqsLY%3DSchreiberKElectron–electron interactions and the paired-to-nematic quantum phase transition in the second Landau level20189171:CAS:528:DC%2BC1cXhtF2jsL7KNayakCSimonSHSternAFreedmanMDas SarmaSNon-Abelian anyons and topological quantum computation200880108311592008RvMP...80.1083N24437221:CAS:528:DC%2BD1cXhtFCku73K1205.81062KumarACsáthyGAManfraMJPfeifferLNWestKWNonconventional odd-denominator fractional quantum Hall states in the second Landau level20101052468082010PhRvL.105x6808K1:STN:280:DC%2BC3M7htlyitQ%3D%3D21231551ShinglaVFinite-temperature behavior in the second Landau level of the two-dimensional electron gas2018972411052018PhRvB..97x1105S1:CAS:528:DC%2BC1MXltVCrsr0%3DBalramACFractional quantum Hall effect at ν = 2 + 6/13: The parton paradigm for the second Landau level20181211866012018PhRvL.121r6601B1:CAS:528:DC%2BC1MXltV2js7s%3D30444400KleinbaumEKumarAPfeifferLNWestKWCsáthyGAGap reversal at filling factors 3 + 1/3 and 3 + 1/5: Towards novel topological order in the fractional quantum Hall regime20151140768012015PhRvL.114g6801K25763967AuslaenderOMTunneling spectroscopy of the elementary excitations in a one-dimensional wire20022958258282002Sci...295..825A1:CAS:528:DC%2BD38XhtFamtrc%3D11823634JompolYProbing spin-charge separation in a Tomonaga-Luttinger liquid20093255976012009Sci...325..597J1:CAS:528:DC%2BD1MXptFOjtrw%3D19644117PotokRMRauIGShtrikmanHOregYGoldhaber-GordonDObservation of the two-channel Kondo effect20074461671712007Natur.446..167P1:CAS:528:DC%2BD2sXisFygtb0%3D1734484917344849IftikharZTwo-channel Kondo effect and renormalization flow with macroscopic quantum charge states20155262332362015Natur.526..233I1:CAS:528:DC%2BC2MXhs1Chu7jO2645005626450056SimonPBrauneckerBLossDMagnetic ordering of nuclear spins in an interacting two-dimensional electron gas2008770451082008PhRvB..77d5108S1175.82062BrauneckerBSimonPInterplay between classical magnetic moments and superconductivity in quantum one-dimensional conductors: Toward a self-sustained topological Majorana phase20131111472022013PhRvL.111n7202B24138267HuangJXiaJSTsuiDCPfeifferLNWestKWDisappearance of metal-like behavior in GaAs two-dimensional holes below 30 mK2007982268012007PhRvL..98v6801H17677869HoS-CImaging the zigzag Wigner crystal in confinement-tunable quantum wires20181211068012018PhRvL.121j6801H1:CAS:528:DC%2BC1MXltV2ju7o%3D30240231BoseSQuantum communication through spin chain dynamics: an introductory overview20074813302007ConPh..48...13B1:CAS:528:DC%2BD2sXnsVCrt7s%3DBateyGIntegration of superconducting magnets with cryogen-free dilution refrigerator systems2009497277342009Cryo...49..727B1:CAS:528:DC%2BD1MXhsVaisLzOBateyGA microkelvin cryogen-free experimental platform with integrated noise thermometry2013151130342013NJPh...15k3034BTodoshchenkoIKaikkonenJ-PBlaauwgeersRHakonenPJSavinADry demagnetization cryostat for sub-millikelvin helium experiments: Refrigeration and thermometry2014850851062014RScI...85h5106T1:STN:280:DC%2BC2M%2FmvVOgtA%3D%3DHashisakaMNoise measurement system at electron temperature down to 20 mK with combinations of the low pass filters2009800961052009RScI...80i6105HSchellerCPSilver-epoxy microwave filters and thermalizers for millikelvin experiments20141042111062014ApPhL.104u1106SIftikhar, Z. et al. Primary thermometry triad at 6 mK in mesoscopic circuits. Nat. Commun.7, 12908 (2016).SamkharadzeNIntegrated electronic transport and thermometry at millikelvin temperatures and in strong magnetic fields2011820539022011RScI...82e3902S1:STN:280:DC%2BC3MjgvVykug%3D%3DXiaJUltra-low-temperature cooling of two-dimensional electron gas20002804914922000PhyB..280..491X1:CAS:528:DC%2BD3cXit1ajtb4%3DNicolíGQuantum dot thermometry at ultra-low temperature in a dilution refrigerator with a 4He immersion cell2019901139012019RScI...90k3901NClarkASchwarzwälderKBandiTMaradanDZumbühlDMethod for cooling nanostructures to microkelvin temperatures2010811039042010RScI...81j3904C1:STN:280:DC%2BC3cbkvV2mug%3D%3DJonesATProgress in cooling nanoelectronic devices to ultra-low temperatures20202017728022020JLTP..201..772J1:CAS:528:DC%2BB3cXhtFSgsrnI332398287679351SarsbyMYurttagülNGeresdiA500 microkelvin nanoelectronics20201114922020NatCo..11.1492S1:CAS:528:DC%2BB3cXlvFCqurs%3D321983827083907LaneJRIntegrating superfluids with superconducting qubit systems20201010123362020PhRvA.101a2336L1:CAS:528:DC%2BB3cXosVGnt7k%3DCaseyACurrent sensing noise thermometry: A fast practical solution to low temperature measurement20141757647752014JLTP..175..764C1:CAS:528:DC%2BC2cXksVSitrc%3DShibaharaAPrimary current-sensing noise thermometry in the millikelvin regime2016374201500542016RSPTA.37450054S1:STN:280:DC%2BC28jjtlOrtg%3D%3DKlemJPersistent photoconductivity in (Al,Ga)As/GaAs modulation doped structures: Dependence on structure and growth temperature198354521452171983JAP....54.5214K1:CAS:528:DyaL3sXltlaisrk%3DBeauchampCBSuperconductivity in AuNiGe ohmic contacts to a GaAs-based high mobility two-dimensional electron gas20201171621042020ApPhL.117p2104B1:CAS:528:DC%2BB3cXitFehu7zJNakayama, T. in Progress in Low Temperature Physics. Brewer, D. (ed.) Vol. 12, chap. 3, 115–194 (Elsevier, 1989).HuYStecherGJGramilaTJRichardsonRCMagnetic coupling in thermal-boundary resistance between thin silver films and liquid 3He in the millikelvin regime199654R9639R96421996PhRvB..54.9639H1:CAS:528:DyaK28Xmt1yktrc%3DSommerfeldAZur elektronentheorie der metalle auf grund der Fermischen statistik1928471321928ZPhy...47....1S1:CAS:528:DyaB1cXhtlSqtQ%3D%3D54.0987.06SymeRTKellyMJPepperMDirect measurement of the thermal conductivity of a two-dimensional electron gas19891337533801989JPCM....1.3375SKumarGSPrasadGPohlROExperimental determinations of the Lorenz number199328426142721993JMatS..28.4261K1:CAS:528:DyaK3sXmtlymtr4%3DChiattiOQuantum thermal conductance of electrons in a one-dimensional wire2006970566012006PhRvL..97e6601C1:STN:280:DC%2BD28rpsVWkuw%3D%3D17026125Jezouin, S. et al. Quantum limit of heat flow across a single electronic channel. Sciencehttps://science.sciencemag.org/content/early/2013/10/02/science.1241912 (2013).MaciáEThermal conductivity in complex metallic alloys: beyond Wiedemann-Franz law2009792451122009PhRvB..79x5112MZavaritskii, N. V. Investigation of the thermal properties of superconductors ii. JETP34,https://www.osti.gov/biblio/4298568 (1958).LinJYDissanayakeABrownGJiangHXRelaxation of persistent photoconductivity in Al0.3Ga0.7As199042585558581990PhRvB..42.5855L1:CAS:528:DyaK3cXmt1alsLs%3DAppleyardNJNichollsJTSimmonsMYTribeWRPepperMA thermometer for the 2D electron gas using 1D thermopower199881349134941998PhRvL..81.3491A1:CAS:528:DyaK1cXmvVCisLs%3DKleinbaumEShinglaVCsáthyGSQUID-based current sensing noise thermometry for quantum resistors at dilution refrigerator temperatures2017880349022017RScI...88c4902Kde GraafSEDirect identification of dilute surface spins on Al2O3: Origin of flux noise in quantum circuits20171180577032017PhRvL.118e7703D28211716BluhmHMolerKADissipative cryogenic filters with zero dc resistance2008790147032008RScI...79a4703BTancrediGSchmidlinSMeesonPJCryogenic coaxial microwave filters2014850261042014RScI...85b6104T1:STN:280:DC%2BC2crgtlWgtA%3D%3DDrungDHighly sensitive and easy-to-use SQUID sensors2007176997042007ITAS...17..699D1:CAS:528:DC%2BD2sXpslelsrY%3DKirsteAEngertJA SQUID-based primary noise thermometer for low-temperature metrology2016374201500502016RSPTA.37450050K1:STN:280:DC%2BC28jjtlOqtw%3D%3DVenturaGThermal conductivity of the normal and superconducting Al/Si 1% alloy1998615765811998NuPhS..61..576VLukashenkoAUstinovAVImproved powder filters for qubit measurements2008790147012008RScI...79a4701L1:STN:280:DC%2BD1c%2FptV2iuw%3D%3DSupplementary information
Supplementary Information
Peer Review File
Supplementary information
The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-28222-x.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.