We generalize the Bardeen-Cooper-Schrieffer-Bose-Einstein-condensation (BCS-BEC) crossover of two-component fermions, which is realized by tuning the s-wave scattering length a between the fermions, to the case of an arbitrary effective range re. By using the Nozières-Schmitt-Rink (NSR) approach, we show another crossover by changing re and present several similarities and differences between these two crossovers. Furthermore, the region (re > a/2) where the effective range expansion breaks down and the Hamiltonian becomes non-Hermitian is found, being consistent with the Wigner's causality bound. Our results are universal for low-density interacting fermions with low-energy constants a and re and are directly relevant to ultracold Fermi atomic gases as well as dilute neutron matter. PACS numbers: 03.75.Ss, 03.75.-b, 03.70.+kThe BCS-BEC crossover, which is realized by tuning the s-wave scattering length a in cold atom experiments [1][2][3][4][5][6][7], has been widely accepted as an important concept to understand strongly correlated quantum systems [8][9][10][11][12][13][14][15][16]. Indeed, this phenomenon has been discussed in various systems such as superconductors [17][18][19][20][21][22][23][24] and dense quark matter [25][26][27][28]. In this regard, thermodynamic properties of strongly interacting ultracold Fermi gases have been experimentally investigated by changing a near the unitarity limit [29][30][31][32][33][34][35][36]. The observed quantities are universal for homogeneous twocomponent fermions with a dimensionless coupling parameter 1/(k F a) where k F is the Fermi momentum. The equation of state in this atomic system [34][35][36][37][38][39][40] shows an excellent agreement with a variational calculation for dilute pure neutron matter (PNM) [41][42][43] which has a relatively large negative scattering length a = −18.5 fm [44] and the coupling parameter 1/(k F a) ≃ −0.04 at a subnuclear density ρ ≃ 0.08 fm −3 .On the other hand, there are of course various differences between ultracold Fermi gases and pure neutron matter such as non-locality of the interaction. The most important difference is the magnitude of an effective range r e . While r e in ultracold Fermi gases near a broad Feshbach resonance is negligible, that in neutron matter given by r e = 2.8 fm [44] largely affects system's properties even around the subnuclear density [40]. On the other hand, a narrow Feshbach resonance in ultracold atoms gives a large and negative effective range [45].Since r e is directly related to the phase shift δ(p) (where p is the momentum), one can expect that the negative (positive) effective range induces a strong (weak) attraction [40,[46][47][48][49][50]. In this sense, a natural question arises: How does the superfluid transition behave if one arbitrarily changes the effective range?The purpose of this work is to answer this question and show that another crossover of the superfluid phase transition from the BCS pairing to the molecular BEC occurs when changing the effective range r e . We address -5 ...