Condensate Fraction Scaling and Berezinskii-Kosterlitz-Thouless Transition of Superconductivity and Superfluidity

Characterizing the superconducting and superfluid transitions in two-dimensional (2D) many-body systems is of broad interest and remains a fundamental issue. In this study, we establish the {\it condensate fraction} scaling as a highly efficient tool to achieve that and accordingly propose efficient schemes to accurately determine the associated Berezinskii-Kosterlitz-Thouless (BKT) transitions. Using the 2D attractive Fermi-Hubbard model as a testbed and applying numerically exact auxiliary-field quantum Monte Carlo simulations, we access unprecedented system sizes (up to $64\times 64 = 4096$ lattice sites) and perform a comprehensive analysis for the temperature dependence and finite-size scaling of {\it condensate fraction} across the BKT transition. We demonstrate that this quantity exhibits algebraic scaling below the transition and exponential scaling above it, with significantly reduced finite-size effects comparing to the extensively studied on-site pairing correlator. This greatly improves the determination of BKT transition using moderate system sizes. We also extract finite-size BKT transition temperature from condensate fraction, and confirm its logarithmic correction on system size. Based on the accurately determined transition, we reveal that the specific heat displays an anomaly, showing a peak at a temperature slightly above BKT transition. Our findings should be generally applicable to 2D fermionic and bosonic systems hosting superconductivity or superfluidity.