Reduced basis emulation of pairing in finite systems

Abstract

In recent years, reduced basis methods (RBMs) have been adapted to the many-body eigenvalue problem and they have been used, largely in nuclear physics, as fast emulators able to bypass expensive direct computations while still providing highly accurate results. This work is meant to show that the RBM is an efficient and accurate emulator for the strong correlations induced by the pairing interaction in a variety of finite systems like ultrasmall superconducting grains, interacting topological superfluids and mesoscopic hybrid superconductor-semiconductor devices, all of which require an expensive, beyond-mean-field, particle-number conserving description. These systems are modelled by the number-conserving Richardson pairing Hamiltonian and its appropriate generalizations. Their ground state is solved for exactly using the Density Matrix Renormalization Group. The reduced basis is assembled iteratively from a small number of exact ground state vectors, well-chosen from across the relevant parameter space using a fast estimate of the emulation error and a greedy local optimization algorithm. The reduced basis emulation is found to accurately describe the weak-to-strong pairing cross-over in small grains, the third-order topological phase transition of the interacting Richardson-Kitaev chain, and the complex charge stability diagram of a hybrid quantum dot - superconductor device. RBMs are thus confirmed to be cheap and accurate emulators for the widely encountered superconducting phenomena. Capable of providing orders of magnitude computational speed-up with respect to approaches based only on traditional many-body solvers, they open new possibilities in building and solving models of interacting many-body systems and in better interfacing them with experimental design and data analysis.

16 pages, 17 figures

References in corpus (9)

• The density-matrix renormalization group in the age of matrix product states
• Global sensitivity analysis of bulk properties of an atomic nucleus
• Many-body characterization of topological superconductivity: The Richardson-Gaudin-Kitaev chain
• Model reduction methods for nuclear emulators
• Fast & accurate emulation of two-body scattering observables without wave functions
• Training and Projecting: A Reduced Basis Method Emulator for Many-Body Physics
• Generalizing the calculable $R$-matrix theory and eigenvector continuation to the incoming wave boundary condition
• Applications of reduced basis methods to the nuclear single particle spectrum
• Strongly correlated electrons in superconducting islands with fluctuating Cooper pairs

Cited by in corpus (7)

• Eigenvector Continuation and Projection-Based Emulators
• Surrogate model solver for impurity-induced superconducting subgap states
• Majorana modes in quantum dots coupled via a floating superconducting island
• Reduced basis surrogates for quantum spin systems based on tensor networks
• Emulating generator coordinate method with extended eigenvector continuation: Lipkin-Meshkov-Glick model
• Reduced Basis Method for Driven-Dissipative Quantum Systems
• Optimal Majoranas in Mesoscopic Kitaev Chains