An effective random matrix theory description is developed for the universal gap fluctuations and the ensemble averaged density of states of chaotic Andreev billiards for finite Ehrenfest time. It yields a very good agreement with the numerical calculation for Sinai-Andreev billiards. A systematic linear decrease of the mean field gap with increasing Ehrenfest time $tau_E$ is observed but its derivative with respect to $tau_E$ is in between two competing theoretical predictions and close to that of the recent numerical calculations for Andreev map. The exponential tail of the density of states is interpreted semi-classically.
We examine the density of states of an Andreev billiard and show that any billiard with a finite upper cut-off in the path length distribution $P(s)$ will possess an energy gap on the scale of the Thouless energy. An exact quantum mechanical calculation for different Andreev billiards gives good agreement with the semi-classical predictions when the energy dependent phase shift for Andreev reflections is properly taken into account. Based on this new semi-classical Bohr-Sommerfeld approximation of the density of states, we derive a simple formula for the energy gap. We show that the energy gap, in units of Thouless energy, may exceed the value predicted earlier from random matrix theory for chaotic billiards.
We present a classical and quantum mechanical study of an Andreev billiard with a chaotic normal dot. We demonstrate that in general the classical dynamics of these normal-superconductor hybrid systems is mixed, thereby indicating the limitations of a widely used retracing approximation. We show that the mixed classical dynamics gives rise to a wealth of wavefunction phenomena, including periodic orbit scarring and localization of the wavefunction onto other classical phase space objects such as intermittent regions and quantized tori.
Comparing the results of exact quantum calculations and those obtained from the EBK-like quantization scheme of Silvestrov et al [Phys. Rev. Lett. 90, 116801 (2003)] we show that the spectrum of Andreev billiards of mixed phase space can basically be decomposed into a regular and an irregular part, similarly to normal billiards. We provide the first numerical confirmation of the validity of this quantization scheme for individual eigenstates and discuss its accuracy.
We demonstrate that the exact quantum mechanical calculations are in good agreement with the semiclassical predictions for rectangular Andreev billiards and therefore for a large number of open channels it is sufficient to investigate the Bohr-Sommerfeld approximation of the density of states. We present exact calculations of the classical path length distribution $P(s)$ which is a non-differentiable function of $s$, but whose integral is a smooth function with logarithmically dependent asymptotic behavior. Consequently, the density of states of rectangular Andreev billiards has two contributions on the scale of the Thouless energy: one which is well-known and it is proportional to the energy, and the other which shows a logarithmic energy dependence. It is shown that the prefactors of both contributions depend on the geometry of the billiards but they have universal limiting values when the width of the superconductor tends to zero.