Abstract
Parasitic twolevel tunnelling systems originating from structural material defects affect the functionality of various microfabricated devices by acting as a source of noise. In particular, superconducting quantum bits may be sensitive to even single defects when these reside in the tunnel barrier of the qubit’s Josephson junctions, and this can be exploited to observe and manipulate the quantum states of individual tunnelling systems. Here, we detect and fully characterize a system of two strongly interacting defects using a novel technique for highresolution spectroscopy. Mutual defect coupling has been conjectured to explain various anomalies of glasses, and was recently suggested as the origin of lowfrequency noise in superconducting devices. Our study provides conclusive evidence of defect interactions with full access to the individual constituents, demonstrating the potential of superconducting qubits for studying material defects. All our observations are consistent with the assumption that defects are generated by atomic tunnelling.
Introduction
Superconducting quantum bits^{1} have recently achieved a breakthrough by demonstrating excellent gate fidelities and long coherence times in a fully scalable architecture^{2}, placing the realization of an integrated quantum computing chip within reach. The solidstate approach, however, bears the burden that the material of the quantum device itself may host parasitic defects that give rise to twolevel systems (TLS) acting as a sparse decohering bath.
First signatures of coherent TLS in phase qubits were found in spectroscopy data, where observed avoided level crossings manifest the defects’ twolevel quantum character^{3,4}. Often, these defects show longer coherence times than the qubit itself^{5}, and thus might be useful as quantum memories^{6} and resources for quantum algorithms^{7,8}. Phase qubits were used in several attempts to identify the physical origin of those TLS, for example, by obtaining statistics on frequencies and coupling strengths^{9}, estimating their density^{4}, measuring the temperature dependence of their coherence times^{5} or verifying theoretical models describing their origin^{10}.
The possibility of a direct interaction between TLS has been invoked in the past to explain the line width broadening and spectral diffusion of ultrasonically excited ensembles of TLS in glasses^{11,12} as well as various other lowtemperature properties of disordered solids^{13,14}. TLS are furthermore a widely accepted model to explain noise in superconducting circuits, and mutual TLS coupling was recently suggested as the origin of the lowfrequency noise observed in microwave resonators^{15}.
Here, we report the first clear experimental evidence of two coherently interacting TLS residing in the tunnel barrier of a Josephson junction (JJ). The data are obtained with a new technique for highresolution defect spectroscopy that exploits the tunability of TLS by mechanical strain and their strong coupling to a superconducting qubit. To characterize the coupled defect system in more detail, we build on this technique and perform coherent twophoton spectroscopy that directly reveals the TLS’ coupling strengths and independent parameters. Interpretation of the measurement based on atomic tunnelling systems fully accounts for all observations.
Results
Atomic tunnelling systems
To explain the microscopic origin of the TLS in superconducting electronics, several models^{16,17,18,19,20} have been proposed. However, all experimental results obtained so far, including the recent demonstration that the energy of the TLS is tunable by static mechanical strain^{21}, are readily explained assuming that they originate from atomic tunnelling systems. As in the wellstudied model describing the lowtemperature thermal, dielectric and acoustic properties of disordered solids^{22,23,24}, it is assumed that some atoms or small groups of atoms are able to tunnel between two energetically almost equivalent sites within the disordered oxide material of the device. These systems give rise to twolevel excitations in a wide energy range of up to the order E≈k_{B}·1 K or ≈h·20 GHz. In bulk disordered solids, TLS are found in large numbers but, in contrast to their counterparts present in superconducting qubits, cannot be addressed individually.
According to the tunnelling model, an atomic tunnelling system is described as a particle in a doublewell potential as shown in Fig. 1a. The energies of the two wells differ by the asymmetry ε and the tunnelling amplitude between them is denoted as Δ, resulting in a level splitting of the two eigenstates given by . Tunnelling systems couple to the environment predominantly by variation δε of their asymmetry energy with δε depending linearly on strain fields and, if the tunnelling entity moves a nonzero charge, as well on electric field—the latter serving as an apparent explanation for the observed coupling of the TLS to the qubit circuit. A variation δΔ of the tunnelling amplitude induced by strain or external fields is generally believed to be negligible^{12,24}. We have recently verified the linear strain dependence of ε and the corresponding hyperbolic variation of the energy splitting E by tracking individual TLS with a phase qubit while bending the chip circuit with a piezo actuator^{21}. In our setup, sketched in Fig. 1b, an applied piezo voltage V_{p} results in variable strain fields in the order of 10^{−6}/V.
Defect spectroscopy
In this study, we detect and analyse TLS using a superconducting qubit. These devices rely on JJs as nonlinear circuit elements, which are realized as two superconducting films separated by a thin, insulating tunnelling barrier, consisting of a 2–3nmthick structurally disordered layer of aluminium oxide. A sketch of the employed phase qubit^{25} including measurement and manipulation circuitry is shown in Fig. 1c. The qubit’s level splitting and their population are controlled by externally applied flux bias and resonant microwave pulses, respectively, and a DCSQUID is used for qubit readout.
In order to trace the energies of individual TLS while applying strain to the qubit chip, we use a spectroscopy scheme based on the pulse protocol depicted in Fig. 2a. The qubit is first biased at a frequency far away from the intended spectroscopy region and excited by a resonant microwave πpulse. Applying appropriate flux bias, it is then tuned to the probing frequency f_{h} where it resides for the holding time τ. If at this frequency the qubit is in resonance with a certain TLS, the excitation is shared between the systems^{6,9}. This results in coherent oscillations that effectively swap the quantum states of the two systems at a frequency determined by their coupling strength as shown in Fig. 2b. A change δP in the qubit excitation probability, measured after the interaction time τ, thus reveals the presence of a TLS. We chose τ to be about half the qubit’s life and coherence times T_{1} and T_{2} (here, both ≈100 ns) to reach a compromise between the loss of signal due to qubit relaxation and the sensitivity needed for detecting also weakly coupled TLS. We repeat this procedure for a range of probing frequencies f_{h} and vary the mechanical strain applied to the sample.
By plotting the change in qubit population as a function of both applied strain and probing frequency, we obtain defect spectra like the one shown in Fig. 2d. Dark traces indicate the resonance frequencies of individual TLS, which are tuned by strain as expected for atomic tunnelling systems. Some TLS have a tunnelling energy Δ that falls within the frequency range accessible by the qubit (~6.5–10 GHz for our sample), while also their asymmetry energy ε is tuned through zero in the investigated strain range. Accordingly, for those TLS, we can clearly observe the hyperbolic strain dependence of their resonance frequencies around minima given by Δ/h. We note that the distribution of TLS resonances changes completely once the sample is warmed to room temperature; see Supplementary Fig. 1 for various examples. This can be explained by a modification of the atomic configuration changing the TLS environment locally, and also by an offset in the applied strain due to thermal dilatation of the sample fixture. As long as the temperature is kept below ~10–20 K, the properties of the majority of TLS remain constant over several months of measurements.
For TLS that are strongly coupled to the qubit, the chosen interaction time τ may exceed the duration of one swap operation. Since the latter also depends on the detuning, a fringelike interference pattern occurs around the traces; see Supplementary Fig. 2. More details on this effect and other artefacts in such defect spectra are discussed in Supplementary Note 1.
Mutually coupled TLS
In the defect spectroscopy example shown in Fig. 2d, two different effects of interactions between TLS can be identified: (i) the resonance frequency of certain TLS is observed to switch between two values in a random, telegraphsignallike fashion, with a switching frequency that depends on the applied strain. This effect is explained by assuming that the observed TLS couples nonresonantly to an incoherent defect that fluctuates between its positions. (ii) Much less frequently observed are nonhyperbolic strain responses in the TLS resonance frequency as well as level splittings that are the typical signature of two resonantly coupled coherent quantum systems.
In the remainder of this article, we discuss a particularly clear manifestation of a system of two coupled TLS, whose defect spectroscopy signature revealed the Sshaped trace with avoided level crossings shown in Fig. 3. To explain this spectroscopic feature, we construct a model of two interacting defects denoted as ‘TLS1’ and ‘TLS2’ (see also inset to Fig. 3). Qualitatively speaking, we observe that TLS1 couples only weakly to the applied external strain so that the steplike increase of its resonance frequency is almost exclusively owing to its nonresonant coupling to TLS2. In the region shown, the doublewell potential of TLS2 is straintuned through its symmetry point ε_{2}(V_{p})=0 at V_{p}≈ −14 V such that the probability density of its ground state is shifted gradually from one potential well to the other. The accompanied shift in atomic positions is mediated via internal strain or electric field to TLS1, which responds by modifying its energy splitting, here from 6 to 6.8 GHz. Eventually, the applied strain tunes TLS2 into resonance with TLS1 and their coherent interaction gives rise to the associated and observed level repulsions.
In the following, we outline our theory of the coupled defect system and show how it can be fully characterized by analysing strainspectroscopic data. The Hamiltonian of a single TLS is written as:
where Δ_{i} is the tunnelling energy and the asymmetry energy ε_{i}(V_{p}) depends linearly on external strain, that is, voltage V_{p} of the piezo drive. Here and in the following, we use the tilde to distinguish operators such as the Paulimatrices in the eigenbasis from those in the localized basis σ_{j}. The energy splitting in the diagonal basis is . The transformation to the diagonal basis corresponds to a rotation about the angle ξ_{i}(V_{p}), defined by tan ξ_{i}=Δ_{i}/ε_{i} with sin ξ_{i}=Δ_{i}/E_{i} and cos ξ_{i}=ε_{i}/E_{i}. For example, the operators σ_{z,i}, whose eigenvalues identify the particle positions, transform as
We write the Hamiltionian of two coupled TLS as H_{T}=H_{1}+H_{2}+H_{12}, with the interaction term,
Within the tunnelling model, the defect’s mutual coupling parameter g comprises their electric dipole interaction as well as a strainmediated elastic contribution. Our spectroscopic data does not allow us to distinguish between these coupling mechanisms. However, since TLS interact with the qubit only electrically, we can determine the projections of the TLS’s electric dipole moments onto the field in the qubit junction independently. Details of this analysis are included in the Supplementary Note 2.
In the diagonal basis, the interaction mentioned in equation (3) consists of four terms that are combinations of and for each TLS (Supplementary Note 3). However, for explaining the observed Sshaped signal (Fig. 3), only two terms are relevant. The energy shift of TLS1, that is, the amplitude of the ‘S’, is given by the longitudinal coupling component which stems from the term . The transversal component results in exchange coupling and defines the size of the level repulsion . The remaining two parts of and , yield only minor energy shifts and can be neglected to first order. However, in our numerical fits to the spectrum, we take the full interaction equation (3) into account.
Twophoton swap spectroscopy
To fully explore the coherently coupled defect system with even higher precision, we focus on the region near the right anticrossing where both TLS are straintuned into resonance. This results in the fourlevel energy spectrum illustrated in Fig. 4a, which we map out by performing microwave swap spectroscopy on the system prepared in different entangled states. For this, we follow the sequence sketched in Fig. 4b. The qubit is first prepared in its excited state and then tuned to either the lower or upper branch of the avoided level crossing, realizing a swap operation with the entangled state or , respectively (Fig. 4a). Directly afterwards, the JJ qubit is again excited and tuned through a lower frequency range in order to find the transition that brings the system of coupled TLS to the fully excited state .
Data in the upper panel of Fig. 4c were obtained by measuring the qubit after the first swap operation, thus indicating the avoided level crossing as in our usual defect spectroscopy protocol in order to calibrate the preparation of the chosen entangled state. The complete sequence results in the data shown in the lower panel of Fig. 4c, where we plot the difference in qubit population between two experiments in which the TLS system was prepared in either one of the two entangled state by adjusting the frequency of the first swap pulse. This experiment clearly reveals the transitions to the fully excited state in excellent agreement with theory (solid lines in Fig. 4c). Moreover, we directly obtain the energies E_{1} and E_{2} of the two unperturbed TLS as well as the longitudinal and transversal coupling strengths as indicated in Fig. 4a. A detailed analysis of this experiment is contained in Supplementary Note 4.
The spectrum calculated from our theoretical model, shown in the inset of Fig. 3, allows one to fit all system parameters and reproduces the data with high accuracy (Supplementary Fig. 3). The interTLS coupling strength g= −872 MHz is calculated from its components and , obtained in swap spectroscopy for the TLS tuned into resonance at V_{p}=7 V as explained above.
Together with the energy splitting E_{i} of the individual TLS at their resonance, we can determine their straindependent mixing angles ξ_{i} to fully characterize the system. The obtained TLS parameters including their dependence on applied piezo voltage are summarized in Table 1, more details on this evaluation are given in Supplementary Note 3.
A final remark concerns the probability of finding two TLS spaced closely enough to expect an interaction of similar strength as in our experiment. Considering only electric dipole interaction and assuming for both TLS an electric dipole moment of d~1 eÅ, which is consistent with the results of our work and agrees with recent observations and theory^{4,10,19,20}, we can make a rough estimate of the maximal distance between the two interacting TLS to be on the order of 5 nm (Supplementary Note 4). When distributing the ~50 TLS visible in the accessible frequency range of ~6–9 GHz evenly onto the area of the 1 μm^{2} large qubit junction, on average each TLS occupies an area of ~140 × 140 nm^{2}. Thus, although one may expect an increased TLS density at interfaces, observing two coherently interacting TLS in JJs is indeed a rare case.
Discussion
The experimental techniques presented here provide a novel spectroscopic view onto the bath of sparse material defects by accessing the quantum states of individual TLS and small coupled systems with a superconducting qubit while straintuning their internal degree of freedom. This lays the ground for further experiments to clarify the microscopic origin of TLS, which is vital for the advancement of various kinds of nanofabricated devices whose functionality is hampered by defects. So far, the atomic tunnelling model readily explains all effects observed with TLS in tunnel junctions and, in particular, the here demonstrated strong coupling of TLS to mechanical strain and their interaction with both coherent as well as randomly fluctuating defects. Our results open way to detailed testing of the 50year old tunnelling model on the basis of individual TLS, for example, by performing defect spectroscopy for statistical analyses of the TLS distribution and by studying the strain dependence of TLS coherence.
Methods
Superconducting qubit sample
The phase qubit sample^{25} used in this work was fabricated in the group of J. M. Martinis at University of California, Santa Barbara. The qubit junction had an area of about 1 μm^{2}, fabricated using aluminium as electrode material and its thermally grown oxide as a tunnel barrier. All data have been obtained at a sample temperature of about 35 mK. The mechanical strain was controlled by bending the sample chip with a piezo transducer^{21}.
Additional information
How to cite this article: Lisenfeld, J. et al. Observation of directly interacting coherent twolevel systems in an amorphous material. Nat. Commun. 6:6182 doi: 10.1038/ncomms7182 (2015).
References
 1
Clarke, J. & Wilhelm, F. K. Superconducting quantum bits. Nature 453, 1031–1042 (2004).
 2
Barends, R. et al. Superconducting quantum circuits at the surface code threshold for fault tolerance. Nature 508, 500–503 (2014).
 3
Simmonds, R. W., Lang, K. M., Hite, D. A., Pappas, D. P. & Martinis, J. M. Decoherence in Josephson qubits from junction resonances. Phys. Rev. Lett. 93, 077003–077006 (2004).
 4
Martinis, J. M. et al. Decoherence in Josephson qubits from dielectric loss. Phys. Rev. Lett. 95, 210503–210506 (2005).
 5
Lisenfeld, J. et al. Measuring the temperature dependence of individual twolevel systems by direct coherent control. Phys. Rev. Lett. 105, 230504–230507 (2010).
 6
Neeley, M. et al. Process tomography of quantum memory in a Josephsonphase qubit coupled to a twolevel state. Nat. Phys. 4, 523–526 (2008).
 7
Müller, C., Shnirman, A. & Weides, M. P. T1echo sequence—protecting the state of a qubit in the presence of coherent interaction. Phys. Rev. A 86, 032335–032342 (2012).
 8
Grabovskij, G. J. et al. Entangling microscopic defects via a macroscopic quantum shuttle. New J. Phys. 13, 063015–063024 (2011).
 9
Shalibo, Y. et al. Lifetime and coherence of twolevel defects in a Josephson junction. Phys. Rev. Lett. 105, 177001–177004 (2010).
 10
Cole, J. H. et al. Quantitative evaluation of defectmodels in superconducting phase qubits. Appl. Phys. Lett. 97, 252501–252503 (2010).
 11
Arnold, W. & Hunklinger, S. Experimental evidence for the direct interaction between twolevel systems in glasses at very low temperatures. Solid State Commun. 17, 883 (1972).
 12
Black, J. L. & Halperin, B. I. Spectral diffusion, phonon echoes, and saturation recovery in glasses at low temperatures. Phys. Rev. B 16, 2879–2895 (1977).
 13
Hunklinger, S. & Enns, C. inSeries of Directions in Condensed Matter Physics vol. 17, (ed. Boolchand P. 499–551World Scientific (2000).
 14
Burin, A. L., Natelson, D., Osheroff, D. D. & Kagan, Y. inTunneling Systems in Amorphous and Crystalline Solids (ed. Esquinazi P. Springer (1998).
 15
Burnett, J. et al. Experimental evidence for the direct interaction between twolevel systems in glasses at very low temperatures. Nat. Commun. 5, 4119 (2014).
 16
Martin, I., Bulaevskii, L. & Shnirman, A. Tunneling spectroscopy of twolevel systems inside a Josephson junction. Phys. Rev. Lett. 95, 127002 (2005).
 17
Constantin, M. & Yu, C. C. Microscopic model of critical current noise in Josephson junctions. Phys. Rev. Lett. 99, 207001–207004 (2007).
 18
De Sousa, R., Whaley, K. B., Hecht, T., Von Delft, J. & Wilhelm, F. K. Microscopic model of critical current noise in Josephsonjunction qubits: Subgap resonances and Andreev bound states. Phys. Rev. B 80, 094515–094523 (2009).
 19
DuBois, T. C., Per, M. C., Russo, S. P. & Cole, J. H. Delocalized oxygen as the origin of twolevel defects in Josephson junctions. Phys. Rev. Lett. 110, 077002–077006 (2013).
 20
Agarwal, K., Martin, I., Lukin, M. D. & Demler, E. Polaronic model of twolevel systems in amorphous solids. Phys. Rev. B 87, 144201–1442014 (2013).
 21
Grabovskij, G. J., Peichl, T., Lisenfeld, J., Weiss, G. & Ustinov, A. V. Strain tuning of individual atomic tunneling systems detected by a superconducting qubit. Science 338, 232–234 (2012).
 22
Phillips, W. A. Tunneling states in amorphous solids. J. Low Temp. Phys. 7, 351–360 (1972).
 23
Anderson, P. W., Halperin, B. I. & Varma, C. Anomalous lowtemperature thermal properties of glasses and spin glasses. Philos. Mag. 25, 1–9 (1972).
 24
Phillips, W. A. Twolevel states in glasses. Rep. Prog. Phys. 50, 1657–1708 (1987).
 25
Steffen, M. et al. State tomography of capacitively shunted phase qubits with high fidelity. Phys. Rev. Lett. 97, 050502–050505 (2005).
Acknowledgements
We would like to thank M. Ansmann and J. M. Martinis (University of California, Saint Barbara) for providing the sample that we measured in this work. We thank E. Demler, T. DuBois, C. Enss, A. Heimes, M.D. Lukin, J.E. Mooij, A. Shnirman and A. Würger for their fruitful discussions. Jared H. Cole and Clemens Müller acknowledge the support of the Australian Research Council and the RMIT Foundation through an International Research Exchange Fellowship. This work was supported by the Deutsche Forschungsgemeinschaft DFG and the State of BadenWürttemberg through the DFG Center for Functional Nanostructures (CFN) as well as by the EU project SOLID. We acknowledge support by Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of Karlsruhe Institute of Technology. Partial support by the Ministry for Education and Science of Russian Federation under contract no. 11.G34.31.0062 and in the framework of Increase Competitiveness Program of the National University of Science and Technology MISIS under contract no. K22014025 is gratefully acknowledged.
Author information
Affiliations
Contributions
The experiments were conceived by J.L., G.J.G., G.W. and A.V.U., and performed by J.L. and G.J.G. Theory and simulations were done by C.M., J.H.C. and G.J.G. The article was written by J.L., G.J.G., G.W. and C.M. with contributions from all authors.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures 14, Supplementary Notes 15, and Supplementary References. (PDF 14052 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Lisenfeld, J., Grabovskij, G., Müller, C. et al. Observation of directly interacting coherent twolevel systems in an amorphous material. Nat Commun 6, 6182 (2015). https://doi.org/10.1038/ncomms7182
Received:
Accepted:
Published:
Further reading

Acoustic spectral holeburning in a twolevel system ensemble
npj Quantum Information (2021)

Quantum sensors for microscopic tunneling systems
npj Quantum Information (2021)

Engineering highcoherence superconducting qubits
Nature Reviews Materials (2021)

Resolving the positions of defects in superconducting quantum bits
Scientific Reports (2020)

Amplitude and frequency sensing of microwave fields with a superconducting transmon qudit
npj Quantum Information (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.