We present a numerical approach which allows the solving of Bethe equations whose solutions define the eigenstates of Gaudin models. By focusing on a new set of variables, the canceling divergences which occur for certain values of the coupling strength no longer appear explicitly. The problem is thus reduced to a set of quadratic algebraic equations. The required inverse transformation can then be realized using only linear operations and a standard polynomial root finding algorithm. The method is applied to Richardsons fermionic pairing model, the central spin model and generalized Dicke model.
In this work, we generalize the numerical approach to Gaudin models developed earlier by us to degenerate systems showing that their treatment is surprisingly convenient from a numerical point of view. In fact, high degeneracies not only reduce the number of relevant states in the Hilbert space by a non negligible fraction, they also allow to write the relevant equations in the form of sparse matrix equations. Moreover, we introduce a new inversion method based on a basis of barycentric polynomials which leads to a more stable and efficient root extraction which most importantly avoids the necessity of working with arbitrary precision. As an example we show the results of our procedure applied to the Richardson model on a square lattice.
We introduce a new family of deep neural network models. Instead of specifying a discrete sequence of hidden layers, we parameterize the derivative of the hidden state using a neural network. The output of the network is computed using a black-box differential equation solver. These continuous-depth models have constant memory cost, adapt their evaluation strategy to each input, and can explicitly trade numerical precision for speed. We demonstrate these properties in continuous-depth residual networks and continuous-time latent variable models. We also construct continuous normalizing flows, a generative model that can train by maximum likelihood, without partitioning or ordering the data dimensions. For training, we show how to scalably backpropagate through any ODE solver, without access to its internal operations. This allows end-to-end training of ODEs within larger models.
We report development and microwave characterization of rf SQUID (Superconducting QUantum Interference Device) qubits, consisting of an aluminium-based Josephson junction embedded in a superconducting loop patterned from a thin film of TiN with high kinetic inductance. Here we demonstrate that the systems can offer small physical size, high anharmonicity, and small scatter of device parameters. The hybrid devices can be utilized as tools to shed further light onto the origin of film dissipation and decoherence in phase-slip nanowire qubits, patterned entirely from disordered superconducting films.
In this note, we discuss implications of the results obtained in [MTV4]. It was shown there that eigenvectors of the Bethe algebra of the quantum gl_N Gaudin model are in a one-to-one correspondence with Fuchsian differential operators with polynomial kernel. Here, we interpret this fact as a separation of variables in the gl_N Gaudin model. Having a Fuchsian differential operator with polynomial kernel, we construct the corresponding eigenvector of the Bethe algebra. It was shown in [MTV4] that the Bethe algebra has simple spectrum if the evaluation parameters of the Gaudin model are generic. In that case, our Bethe ansatz construction produces an eigenbasis of the Bethe algebra.