Search

Superconducting qubits are a key contender for quantum computing elements, but they often face challenges like noise and decoherence from two-level systems (TLS). Tantalum (Ta) qubits are notable for their long T$_1$ coherence times nearing milliseconds, mainly due to fewer TLS, though the cause was unclear. Our research explored this by analyzing the air-oxide interface with density functional theory, particularly comparing Nb oxide (Nb$_2$O$_5$) and Ta oxide (Ta$_2$O$_5$). We discovered that Ta$_2$O$_5$ forms a smoother surface with fewer dangling O atoms and TLS than Nb$_2$O$_5$. The greater atomic mass of Ta also lowers the TLS tunnel splittings below the qubit's operating frequency. Furthermore, using external electric fields or SO$_2$ passivation can significantly reduce TLS on Nb surfaces, potentially improving their coherence times.

Superconducting quantum interference devices (SQUIDs) show great promise as quantum bits (qubits) but continue to be hindered by flux noise. The flux noise power spectra of SQUIDs go as $1/f^\alpha$, where $\alpha$ is the temperature-dependent noise exponent. Experiments find $0.5 \lesssim \alpha \lesssim 1$. Furthermore, experiments find that the noise power spectra versus frequency at different temperatures pivot about or cross at a common point for each SQUID. To try to better understand the results and motivated by experimental evidence that magnetic moments on the surface of SQUIDS produce flux noise, we present the results of our Monte Carlo simulations of various spin systems on 2D lattices. We find that only spin glasses produce $\alpha \sim 1$ at low temperature. We find that aliasing of the noise power spectra at high frequencies can lead to spectral pivoting if it is in proximity to a knee at a slightly lower frequency. We show that the pivot frequency depends on the method of site selection and how often the magnetization is recorded. The spectral pivoting that occurs in our simulations is due to aliasing and does not explain the spectral pivoting of experiments., Comment: 16 pages, 14 figures

Charge noise in semiconducting quantum dots has been observed to have a 1/f spectrum. We propose a model in which a pair of quantum dots are coupled to a 2D bath of fluctuating two level systems (TLS) that have electric dipole moments and that interact with each other, i.e., with the other fluctuators. These interactions are primarily via the elastic strain field. We use a 2D nearest-neighbor Ising spin glass to represent these elastic interactions and to simulate the dynamics of the bath of electric dipole fluctuators in the presence of a ground plane representing metal gates above the oxide layer containing the fluctuators. The interactions between the TLS cause the energy splitting of individual fluctuators to change with time. We calculate the resulting fluctuations in the electric potential at the two quantum dots that lie below the oxide layer. We find that 1/f electric potential noise spectra at the quantum dots and cross correlation in the noise between the two quantum dots are in qualitative agreement with experiment. Our simulations find that the cross correlations decrease exponentially with increasing quantum dot separation., Comment: 14 pages, 11 figures

Silicon quantum dot qubits show great promise but suffer from charge noise with a 1/f^\alpha spectrum, where f is frequency and \alpha \lesssim 1. It has recently been proposed that 1/f^\alpha noise spectra can emerge from a few thermally activated two-level fluctuators in the presence of sub-bath temperature fluctuations associated with a two-dimensional electron gas (2DEG)~\cite{Ahn2021}. We investigate this proposal by doing Monte Carlo simulations of a single Ising spin in a bath with a fluctuating temperature. We find that to obtain noise with a $1/f^\alpha$ spectrum with $alpha \lesssim 1 down to low frequencies, the duration of temperature fluctuations must be comparable to the inverse of the lowest frequency at which the noise is measured. This result is consistent with an analytic calculation in which the fluctuator is a two-state system with dynamics governed by time-dependent switching rates. In this case we find that the noise spectrum follows a Lorentzian at frequencies lower than the inverse of the average duration of the lowest switching rate. We then estimate relaxation times of thermal fluctuations by considering thermal diffusion in an electron gas in a confined geometry. We conclude that temperature fluctuations in a 2DEG sub-bath would require an unphysically long duration to be consistent with experimental measurements of 1/f-like charge noise in quantum dots at frequencies extending well below 1 Hz., Comment: 8 pages, 5 figures

The thermal and acoustic properties displayed by a wide variety of glasses at low temperatures are well described by the model of tunneling two level systems (TLS). We review the standard TLS model as well as developments that have occurred since the earlier work. In spite of half a century of effort and impressive successes, mysteries still remain. We discuss some of these longstanding puzzles about glasses at low temperatures and speculate about future directions., Comment: 28 pages, 1 figure, 1 table. To be published as Chapter 5 of the book entitled, "Low-Temperature Thermal and Vibrational Properties of Disordered Solids: A Half-Century of Universal "Anomalies" in Glasses," edited by M. A. Ramos, World Scientific Press

The goal of immunotherapy is to enhance the ability of the immune system to kill cancer cells. Immunotherapy is more effective and, in general, the prognosis is better, when more immune cells infiltrate the tumor. We explore the question of whether the spatial distribution rather than just the density of immune cells in the tumor is important in forecasting whether cancer recurs. After reviewing previous work on this issue, we introduce a novel application of maximum entropy to quantify the spatial distribution of discrete point-like objects. We apply our approach to B and T cells in images of tumor tissue taken from triple negative breast cancer (TBNC) patients. We find that there is a distinct difference in the spatial distribution of immune cells between good clinical outcome (no recurrence of cancer within at least 5 years of diagnosis) and poor clinical outcome (recurrence within 3 years of diagnosis). Our results highlight the importance of spatial distribution of immune cells within tumors with regard to clinical outcome, and raise new questions on their role in cancer recurrence.

We present a novel view of the standard model of tunneling two level systems (TLS) to explain the puzzling universal value of a quantity, $C\sim 3\times 10^{-4}$, that characterizes phonon scattering in glasses below 1 K as reflected in thermal conductivity, ultrasonic attenuation, internal friction, and the change in sound velocity. Physical considerations lead to a broad distribution of phonon-TLS couplings that (1) exponentially renormalize tunneling matrix elements, and (2) reduce the TLS density of states through TLS-TLS interactions. We find good agreement between theory and experiment for a variety of individual glasses., Comment: Resubmission: several typos were corrected and the notation clarified in some places. No major changes

Superconducting qubits are hampered by flux noise produced by surface spins from a variety of microscopic sources. Recent experiments indicated that hydrogen (H) atoms may be one of those sources. Using density functional theory calculations, we report that H atoms either embedded in, or adsorbed on, an a-Al2O3(0001) surface have sizeable spin moments ranging from 0.81 to 0.87 uB with energy barriers for spin reorientation as low as ~10 mK. Furthermore, H adatoms on the surface attract gas molecules such as O2, producing new spin sources. We propose coating the surface with graphene to eliminate H-induced surface spins and to protect the surface from other adsorbates., Comment: 12 pages, 4 figures

A major obstacle to using SQUIDs as qubits is flux noise. We propose that the heretofore mysterious spins producing flux noise could be $O_2$ molecules adsorbed on the surface. Using density functional theory calculations, we find that an $O_2$ molecule adsorbed on an {\alpha}-alumina surface has a magnetic moment of ~1.8 {\mu}B. When the spin is oriented perpendicular to the axis of the O-O bond, the barrier to spin rotations is about 10 mK. Monte Carlo simulations of ferromagnetically coupled, anisotropic XY spins on a square lattice find 1/f magnetization noise, consistent with flux noise in Al SQUIDs., Comment: 22 pages include main text plus supplemental material

We have performed extensive molecular dynamics simulations to study noise power spectra of density and potential energy fluctuations of a Lennard-Jones model of a fluid in the supercritical region. Emanating from the liquid-vapor critical point, there is a locus of isobaric specific heat maxima, called the Widom line, which is often regarded as an extension of the liquid-vapor coexistence line. Our simulation results show that the noise power spectrum of the density fluctuations on the Widom line of the liquid-vapor transition exhibits three distinct $1/f^{\gamma}$ behaviors with exponents $\gamma$= 0, 1.2 and 2, depending on the frequency $f$. We find that the intermediate frequency region with an exponent $\gamma \sim$ 1 appears as the temperature approaches the Widom temperature from above or below. On the other hand, we do not find three distinct regions of $1/f^{\gamma}$ in the power spectrum of the potential energy fluctuations on the Widom line. Furthermore, we find that the power spectra of both the density and potential energy fluctuations at low frequency have a maximum on the Widom line, suggesting that the noise power can provide an alternative signature of the Widom line.

Recent experiments implicate spins on the surface of metals as the source of flux noise in SQUIDs, and indicate that these spins are able to relax without conserving total magnetization. We present a model of $1/f$ flux noise in which electron spins on the surface of metals can relax via hyperfine interactions. Our results indicate that flux noise would be significantly reduced in superconducting materials where the most abundant isotopes do not have nuclear moments such as zinc and lead., Comment: 4+ pages, no figures

We propose that stress can decrease the internal friction of amorphous solids, either by increasing the potential barriers of defects, thus reducing their tunneling and thermal activation that produce loss, or by decreasing the coupling between defects and phonons. This stress can be from impurities, atomic bonding constraints, or externally applied stress. Externally applied stress also reduces mechanical loss through dissipation dilution. Our results are consistent with the experiments, and predict that stress could substantially reduce dielectric loss and increase the thermal conductivity., Comment: 9 pages, 7 figures

We study the charge noise $S_Q$ in Josephson qubits produced by fluctuating two level systems (TLS) with electric dipole moments in the substrate. The TLS are driven by an alternating electric field of angular frequency $\Omega$ and electric field intensity $I$. It is not widely appreciated that TLS in small qubits can easily be strongly saturated if $I\gg I_c$, where $I_c$ is the critical electric field intensity. To investigate the effect of saturation on the charge noise, we express the noise spectral density in terms of density matrix elements. To determine the dependence of the density matrix elements on the ratio $I/I_c$, we find the steady state solution for the density matrix using the Bloch-Redfield differential equations. We then obtain an expression for the spectral density of charge fluctuations as a function of frequency $f$ and the ratio $I/I_c$. We find $1/f$ charge noise at low frequencies, and that the charge noise is white (constant) at high frequencies. Using a flat density of states, we find that TLS saturation has no effect on the charge noise at either high or low frequencies., Comment: 13 pages, 9 figures, accepted for publication in Phys. Rev. B

We provide a general scheme to predict and derive the contribution to the noise spectrum of a stochastic sequence of pulses from the distribution of pulse parameters. An example is the magnetization noise spectra of a 2D Ising system near its phase transition. At $T\le T_c$, the low frequency spectra is dominated by magnetization flips of nearly the entire system. We find that both the predicted and the analytically derived spectra fit those produced from simulations. Subtracting this contribution leaves the high frequency spectra which follow a power law set by the critical exponents., Comment: 4 pages, 5 figures. We improved text and included a predicted noise curve in Figure 4. 2 examples from Figure 3 are removed

We show that charge noise $S_Q$ in Josephson qubits can be produced by fluctuating two level systems (TLS) with electric dipole moments in the substrate using a flat density of states. At high frequencies the frequency and temperature dependence of the charge noise depends on the ratio $J/J_c$ of the electromagnetic flux $J$ to the critical flux $J_c$. It is not widely appreciated that TLS in small qubits can easily be strongly saturated with $J/J_c\gg 1$. Our results are consistent with experimental conclusions that $S_Q\sim 1/f$ at low frequencies and $S_Q\sim f$ at high frequencies., Comment: 4 pages

We present a simple microscopic model to show how fluctuating two--level systems in the Josephson junction tunnel barrier of thickness $L$ can modify the potential energy of the barrier and produce critical current noise spectra. We find low frequency $1/f$ noise that goes as $L^5$. Our values are in good agreement with recent experimental measurements of critical current noise in Al/AlO$_{x}$/Al Josephson junctions. We also investigate the sensitivity of the noise on the nonuniformity of the tunnel barrier., Comment: 4 pages, 3 figures

Dilute dipolar systems in three dimensions are expected to undergo a spin glass transition as the temperature decreases. Contrary to this, we find from Wang-Landau Monte Carlo simulations that at low concentrations $x$, dipoles randomly placed on a cubic lattice with dipolar interactions do not undergo a phase transition. We find that in the thermodynamic limit the ``glass'' transition temperature $T_g$ goes to zero as $1/\sqrt{N}$ where $N$ is the number of dipoles. The entropy per particle at low temperatures is larger for lower concentrations ($x=4.5%$) than for higher concentrations ($x=20%$).

Dielectric loss from two-level states is shown to be a dominant decoherence source in superconducting quantum bits. Depending on the qubit design, dielectric loss from insulating materials or the tunnel junction can lead to short coherence times. We show that a variety of microwave and qubit measurements are well modeled by loss from resonant absorption of two-level defects. Our results demonstrate that this loss can be significantly reduced by using better dielectrics and fabricating junctions of small area $\lesssim 10 \mu \textrm{m}^2$. With a redesigned phase qubit employing low-loss dielectrics, the energy relaxation rate has been improved by a factor of 20, opening up the possibility of multi-qubit gates and algorithms., Comment: shortened version submitted to PRL

Noise and decoherence are major obstacles to the implementation of Josephson junction qubits in quantum computing. Recent experiments suggest that two level systems (TLS) in the oxide tunnel barrier are a source of decoherence. We explore two decoherence mechanisms in which these two level systems lead to the decay of Rabi oscillations that result when Josephson junction qubits are subjected to strong microwave driving. (A) We consider a Josephson qubit coupled resonantly to a two level system, i.e., the qubit and TLS have equal energy splittings. As a result of this resonant interaction, the occupation probability of the excited state of the qubit exhibits beating. Decoherence of the qubit results when the two level system decays from its excited state by emitting a phonon. (B) Fluctuations of the two level systems in the oxide barrier produce fluctuations and 1/f noise in the Josephson junction critical current I_o. This in turn leads to fluctuations in the qubit energy splitting that degrades the qubit coherence. We compare our results with experiments on Josephson junction phase qubits., Comment: 23 pages, Latex, 6 encapsulated postscript figures

Noise is often considered to be a nuisance. Here we argue that it can be a useful probe of fluctuating two level systems in glasses. It can be used to: (1) shed light on whether the fluctuations are correlated or independent events; (2) determine if there is a low temperature glass or phase transition among interacting two level systems, and if the hierarchical or droplet model can be used to describe the glassy phase; and (3) find the lower bound of the two level system relaxation rate without going to ultralow temperatures. Finally we point out that understanding noise due to two level systems is important for technological applications such as quantum qubits that use Josephson junctions., Comment: 15 pages, 4 figures, Latex, to be published in J. Low Temp. Phys. issue in honor of S. Hunklinger

We have performed a Monte Carlo study of a classical three dimensional Coulomb system in which we systematically increase the positional disorder. We start from a completely ordered system and gradually transition to a Coulomb glass. The phase transition as a function of temperature is second order for all values of disorder. We use finite size scaling to determine the transition temperature T_C and the critical exponent nu. We find that T_C decreases and that nu increases with increasing disorder. We also observe changes in the specific heat, the single particle density of states, and the staggered occupation as a function of disorder and temperature., Comment: 14 pages, Latex, 21 encapsulated postscript figures, to be published in Phys. Rev. B

We show that anomalous diffusion can result when the steps of a random walk are not statistically independent. We present an algorithm that counts all the possible paths of particles diffusing on random graphs with arbitrary degree distribution. Using this to calculate the mean square displacement, we show that in sharp contrast to continua, random walks on random graphs can exhibit anomalous behavior and yet have well-defined and predictable properties., Comment: 4 pages, 6 figures

We show that a glass transition, signaled by a peak in the specific heat vs. temperature, can occur because a glassy system that shows no signs of aging progresses so slowly through the energy landscape that the time needed to obtain an accurate estimate of the thermodynamic averages exceeds the observation time. We find that below the glass transition temperature of a three dimensional binary mixture of soft spheres, the specific heat increases with measurement time spans orders of magnitude longer than previously recognized equilibration times., Comment: 4 pages, Latex, 3 encapsulated postscript figures

We introduce a new quantity to probe the glass transition. This quantity is a linear generalized compressibility which depends solely on the positions of the particles. We have performed a molecular dynamics simulation on a glass forming liquid consisting of a two component mixture of soft spheres in three dimensions. As the temperature is lowered (or as the density is increased), the generalized compressibility drops sharply at the glass transition, with the drop becoming more and more abrupt as the measurement time increases. At our longest measurement times, the drop occurs approximately at the mode coupling temperature $T_C$. The drop in the linear generalized compressibility occurs at the same temperature as the peak in the specific heat. By examining the inherent structure energy as a function of temperature, we find that our results are consistent with the kinetic view of the glass transition in which the system falls out of equilibrium. We find no size dependence and no evidence for a second order phase transition though this does not exclude the possibility of a phase transition below the observed glass transition temperature. We discuss the relation between the linear generalized compressibility and the ordinary isothermal compressibility as well as the static structure factor., Comment: 18 pages, Latex, 26 encapsulated postscript figures, revised paper is shorter, to appear in Phys. Rev. E

We show that 1/f noise is produced in a 3D electron glass by charge fluctuations due to electrons hopping between isolated sites and a percolating network at low temperatures. The low frequency noise spectrum goes as \omega^{-\alpha} with \alpha slightly larger than 1. This result together with the temperature dependence of \alpha and the noise amplitude are in good agreement with the recent experiments. These results hold true both with a flat, noninteracting density of states and with a density of states that includes Coulomb interactions. In the latter case, the density of states has a Coulomb gap that fills in with increasing temperature. For a large Coulomb gap width, this density of states gives a dc conductivity with a hopping exponent of approximately 0.75 which has been observed in recent experiments. For a small Coulomb gap width, the hopping exponent approximately 0.5., Comment: 8 pages, Latex, 6 encapsulated postscript figures, to be published in Phys. Rev. B

At low temperatures electron hopping in a three dimensional Coulomb glass produces fluctuations in the single particle density of states and hence in the resistivity. This results in a low frequency resisitivity noise spectrum which goes as (1/f)^(alpha) where alpha is very close to 1. This holds down to extremely low frequencies., Comment: 5 pages, Latex, 1 postscript figures

We introduce a new quantity to probe the glass transition. This quantity is a linear generalized compressibility which depends solely on the positions of the particles. We have performed a molecular dynamics simulation on a glass forming liquid consisting of a two component mixture of soft spheres in three dimensions. As the temperature is lowered, the generalized compressibility drops sharply at the glass transition. Our results are consistent with the kinetic view of the glass transition, but not with an underlying second order phase transition., Comment: 5 pages, Latex, 4 postscript figures, to be published in Phil. Mag. B, some corrections to previous version

We show that thermally activated interstitial and vacancy defects can lead to first order melting of a vortex lattice. We obtain good agreement with experimentally measured melting curve, latent heat, and magnetization jumps for YBCO and BSCCO. The shear modulus of the vortex liquid is frequency dependent and crosses over from zero at low frequencies to a finite value at high frequencies. We also find a small surface tension between the vortex line liquid and the vortex lattice., Comment: 10 pages, Latex, 6 encapsulated postscript figures

We show that the time development of the Coulomb gap in a Coulomb glass can involve very long relaxation times due to electron rearrangement and hopping. We find that an applied magnetic field reduces the rate of electron hopping and, hence, Coulomb gap formation. These results are consistent with recent conductance experiments on thin semiconducting and metallic films., Comment: 4 pages, Latex, 3 encapsulated postscript figures

This paper has been withdrawn by the authors. It has been superceded by an improved investigation: ``Thermodynamic Generalized Compressibility in a Glass Forming Liquid,'' H. M. Carruzzo and C. C. Yu, cond-mat/0004320., Comment: withdrawn; superceded by improved investigation:cond-mat/0004320

Vortex lattices in the high temperature superconductors undergo a first order phase transition which has thus far been regarded as melting from a solid to a liquid. We point out an alternative possibility of a two step process in which there is a first order transition from an ordinary vortex lattice to a soft vortex solid followed by another first order melting transition from the soft vortex solid to a vortex liquid. We focus on the first step. This premelting transition is induced by vacancy and interstitial vortex lines. We obtain good agreement with the experimental transition temperature versus field, latent heat, and magnetization jumps for YBCO and BSCCO., Comment: revised version replaces 9705092, 5 pages, Latex, 2 postscript figures, defect line wandering is included, 2 step melting is proposed

We have studied a Kondo hole in a one-dimensional Kondo insulator at half-filling using a density matrix formulation of the numerical renormalization group. The Kondo hole introduces midgap states. The spin density introduced by the hole is localized in the vicinity of the hole. It resides primarily in the f-spins for small exchange coupling $J$ and in the conduction spins for large $J$. We present results on the spin gap, charge gap, and neutral gap. For small $J$, the spin gap is smaller than the charge gap, while for large $J$, the spin gap is larger than the charge gap. The presence of the Kondo hole reduces RKKY interactions as can be seen in the staggered susceptibility., Comment: Revtex + postscript figures uuencoded

We study an Anderson impurity in a semiconducting host using the density matrix renormalization group technique. We use the $U=0$ one--dimensional Anderson Hamiltonian at half filling as the semiconducting host since it has a hybridization gap. By varying the hybridization of the host, we can control the size of the semiconducting gap. We consider chains with 25 sites and we place the Anderson impurity (with $U>0$) in the middle of the chain. We dope the half--filled system with one hole and we find two regimes: when the hybridization of the impurity is small, the hole density and the spin are localized near the impurity. When the hybridization of the impurity is large, the hole and spin density are spread over the lattice. Additional holes avoid the impurity and are extended throughout the lattice. Away from half--filling, the semiconductor with an impurity is analogous to a double well potential with a very high barrier. We also examine the chemical potential as a function of electron filling, and we find that the impurity introduces midgap states when the impurity hybridization is small., Comment: Revtex + figures uuencoded

In order to better understand Kondo insulators, we have studied both the symmetric and asymmetric Anderson lattices at half-filling in one dimension using the density matrix formulation of the numerical renormalization group. We have calculated the charge gap, spin gap and quasiparticle gap as a function of the repulsive interaction U using open boundary conditions for lattices as large as 24 sites. We find that the charge gap is larger than the spin gap for all U for both the symmetric and asymmetric cases. RKKY interactions are evident in the f-spin-f-spin correlation functions at large U in the symmetric case, but are suppressed in the asymmetric case as the f-level approaches the Fermi energy. This suppression can also be seen in the staggered susceptibility and it is consistent with neutron scattering measurements in CeNiSn., Comment: 32 pages, Latex file with Postcript figures.

We investigate whether or not Kondo insulators undergo a magnetic field induced metal-insulator transition in one dimension at half filling using both a density matrix formulation of the numerical renormalization group and a new analytical approach. Contrary to expectations, the quasiparticle gap never vanishes at any field and no metal-insulator transition is found. We discuss generalizing our result to three dimensions and to the asymmetric Anderson lattice. Our result is consistent with recent experiments., Comment: 4 pages (using REVTEX), 2 figures in postscript. Submitted to Physical Review Letters