Thermalization slowing down in multidimensional Josephson junction networks

We characterize thermalization slowing down of Josephson junction networks in one, two, and three spatial dimensions for systems with hundreds of sites by computing their entire Lyapunov spectra. The ratio of Josephson coupling EJ to energy density h controls two different universality classes of thermalization slowing down, namely, the weak-coupling regime, EJ/h -> 0, and the strong-coupling regime, EJ/h -> infinity. We analyze the Lyapunov spectrum by measuring the largest Lyapunov exponent and by fitting the rescaled spectrum with a general ansatz. We then extract two scales: the Lyapunov time (inverse of the largest exponent) and the exponent for the decay of the rescaled spectrum. The two universality classes, which exist irrespective of network dimension, are characterized by different ways the extracted scales diverge. The universality class corresponding to the weak-coupling regime allows for the coexistence of chaos with a large number of near-conserved quantities and is shown to be characterized by universal critical exponents, in contrast with the strong-coupling regime. We expect our findings, which we explain using perturbation theory arguments, to be a general feature of diverse Hamiltonian systems.