The quantum mechanics of the Josephson effect is the core ingredient for quantum technologies with superconducting circuits. A new avenue was recently opened in this field by predicting that the Josephson quantum mechanics in the odd parity sector, when a quasiparticle is trapped in an Andreev bound state, is fundamentally different from the conventional one in the even sector. The focus was then on a Josephson junction surrounded by an electromagnetic environment formed of a collection of bosonic modes, including the case of an ohmic environment. Here we consider the distinct case of a superconducting qubit made of a single Josephson junction whose environment reduces to a capacitance. We find a novel structure for the low-lying discrete states in the odd sector, which is altogether different from the one that appears in the even sector. Our study of the bound-state spectrum ranges from the Coulomb-dominated (Cooper pair box) to the Josephson-dominated (transmon) regime. Our prediction could be tested in forthcoming experiments with superconductor/semiconductor/superconductor junctions, which have been studied intensively in recent years, both using nanowires as well as two-dimensional electron gases.

A Josephson junction may be in a stable odd-parity state when a single quasiparticle is trapped in an Andreev bound state. Embedding such junction in an electromagnetic environment gives rise to a special quantum mechanics of the superconducting phase that we investigate theoretically. Our analysis covers several representative cases, from the lifting of the supercurrent quench due to quasiparticle poisoning for a low ohmic impedance of the environment, to a Schmid transition in a current-biased junction that for odd parity occurs at four times bigger critical impedance. For intermediate impedances, the supercurrent in the odd state is higher than in the even one.

In conventional Bardeen-Cooper-Schrieffer superconductors, Cooper pairs of electrons of opposite spin (i.e., singlet structure) form the ground state. Equal-spin triplet pairs (ESTPs), as in superfluid He3, are of great interest for superconducting spintronics and topological superconductivity, yet remain elusive. Recently, odd-parity ESTPs were predicted to arise in (few-)monolayer superconducting NbSe2, from the noncollinearity between the out-of-plane Ising spin-orbit field (due to the lack of inversion symmetry in monolayer NbSe2) and an applied in-plane magnetic field. These ESTPs couple to the singlet order parameter at finite field. Using van der Waals tunnel junctions, we perform spectroscopy of superconducting NbSe2 flakes, of 2-25 monolayer thickness, measuring the quasiparticle density of states (DOS) as a function of applied in-plane magnetic field up to 33 T. In flakes ≲15 monolayers thick the DOS has a single superconducting gap. In these thin samples, the magnetic field acts primarily on the spin (vs orbital) degree of freedom of the electrons, and superconductivity is further protected by the Ising field. The superconducting energy gap, extracted from our tunneling spectra, decreases as a function of the applied magnetic field. However, in bilayer NbSe2, close to the critical field (up to 30 T, much larger than the Pauli limit), superconductivity appears to be more robust than expected from Ising protection alone. Our data can be explained by the above-mentioned ESTPs.

We show that the annihilation dynamics of excess quasiparticles in superconductors may result in the spontaneous formation of large spin-polarized clusters. This presents a novel scenario for spontaneous spin polarization. We estimate the relevant scales for aluminum, finding the feasibility of clusters with total spin S≃10^{4}ℏ that could be spread over microns. The fluctuation dynamics of such large spins may be detected by measuring the flux noise in a loop hosting a cluster.

Recently we predicted that the Andreev bound-state spectrum of four-terminal Josephson junctions may possess topologically protected zero-energy Weyl singularities, which manifest themselves in a quantized transconductance in units of 4e2/h when two of the terminals are voltage biased [R.-P. Riwar, M. Houzet, J. S. Meyer, and Y. V. Nazarov, Nature Commun. 7, 11167 (2016)2041-172310.1038/ncomms11167]. Here, using the Landauer-Büttiker scattering theory, we compute numerically the currents flowing through such a structure in order to assess the conditions for observing this effect. We show that the voltage below which the transconductance becomes quantized is determined by the interplay of nonadiabatic transitions between Andreev bound states and inelastic relaxation processes. We demonstrate that the topological quantization of the transconductance can be observed at voltages of the order of 10-2Δ/e,Δ being the the superconducting gap in the leads.

Experimentally, the concentration of quasiparticles in gapped superconductors always largely exceeds the equilibrium one at low temperatures. Since these quasiparticles are detrimental for many applications, it is important to understand theoretically the origin of the excess. We demonstrate in detail that the dynamics of quasiparticles localized at spatial fluctuations of the gap edge becomes exponentially slow. This gives rise to the observed excess in the presence of a vanishingly weak nonequilibrium agent.

Topological materials and their unusual transport properties are now at the focus of modern experimental and theoretical research. Their topological properties arise from the bandstructure determined by the atomic composition of a material and as such are difficult to tune and naturally restricted to ≤3 dimensions. Here we demonstrate that n-terminal Josephson junctions with conventional superconductors may provide novel realizations of topology in n-1 dimensions, which have similarities, but also marked differences with existing 2D or 3D topological materials. For n≥4, the Andreev subgap spectrum of the junction can accommodate Weyl singularities in the space of the n-1 independent superconducting phases, which play the role of bandstructure quasimomenta. The presence of these Weyl singularities enables topological transitions that are manifested experimentally as changes of the quantized transconductance between two voltage-biased leads, the quantization unit being 4e²/h, where e is the electric charge and h is the Planck constant.

We study the density of states in disordered s-wave superconductors with a small gap anisotropy. We consider disorder in the form of common nonmagnetic scatterers and pairing-potential impurities, which interact with electrons via an electric potential and a local distortion of the superconducting gap. Using quasiclassical Green functions, we determine the bound-state spectrum at a single impurity and the density of states at a finite concentration of impurities. We show that, if the gap is isotropic, an isolated impurity with suppressed pairing supports an infinite number of Andreev states. With growing impurity concentration, the energy-dependent density of states evolves from a sharp gap edge with an impurity band below it to a smeared BCS singularity in the so-called universal limit. Within one spin sector, pairing-potential impurities and weak spin-polarized magnetic impurities have essentially the same effect on the density of states. We note that, if a gap anisotropy is present, the density of states becomes sensitive to ordinary potential disorder, and the existence of Andreev states localized at pairing-potential impurities requires special conditions. An unusual feature related to the anisotropy is a nonmonotonic dependence of the gap edge smearing on impurity concentration.

We use the Landauer-Büttiker scattering theory for electronic transport to calculate the current cross-correlations in a voltage-biased three-terminal junction with all superconducting leads. At low bias voltage, when charge transport is due to coherent multiple Andreev reflections, we find large cross-correlations compared with their normal-state value. Furthermore, depending on the parameters that characterize the properties of the scattering region between the leads, the cross-correlations can reverse their sign with respect to the case of non-interacting fermionic systems.

We study a short multichannel superconducting junction subject to dc and ac phase biases. The ac modulation changes the occupation of the Andreev bound states formed at the constriction by transitions between bound states and the continuum. In a short junction, the nonequilibrium Andreev bound-state population may relax through processes that conserve parity of the occupation number on the same bound state and processes that do not conserve it. We argue that the parity-conserving processes occur on a much faster time scale. In this case, even a weak driving may lead to a large nonequilibrium quasiparticle population scaling with the number of channels and results in a large deviation of the supercurrent from its equilibrium value. We show that this effect is accompanied by a quasiparticle current which may lead to a measurable charge imbalance in the vicinity of the junction. Furthermore, we study the time evolution of the supercurrent after switching off the ac drive. On a time scale where parity relaxation is negligible, the supercurrent relaxes to a stationary nonequilibrium state. Finally, we briefly outline the regime of ultraweak driving where the ac-induced processes occur on a time scale comparable to that of parity relaxation.