Authors: Vladislav Pokorný and Tomáš Novotný
Phys. Rev. Research 3, 023013 – Published 2 April 2021
Abstract: Interacting single-level quantum dots connected to BCS superconducting leads represent a well-controllable system to study the interplay between the correlation effects and the electron pairing that can result in a 0−π (singlet-doublet) quantum phase transition (QPT). The physics of this system can be well described by the impurity Anderson model. We present an analysis of the statistics of the perturbation expansion order of the continuous-time hybridization expansion quantum Monte Carlo algorithm in the vicinity of such a first-order impurity QPT. By calculating the moments of the histograms of the expansion order which deviate from the ideal Gaussian shape, we provide an insight into the thermodynamic mixing of the two phases at low but finite temperatures, which is reflected in the bimodal nature of the histograms.