Phase Diagrams and Electromagnetic Propertiesof s–wave Superconductivity of the Extended Hubbard Model with the Attractive Pair–Hopping Interaction

The characteristics of the extended Hubbard model with the pair–hopping interaction J, i.e. the Penson–Kolb–Hubbard model, are analysed for the case of attractive J (J > 0). The effects of on–site U interaction on the properties of the s–wave superconducting state and mutual stability of the superconducting s–wave and magnetic orderings are examined. The electromagnetic and thermodynamic characteristics of the s–wave superconducting phase for arbitrary particle concentration n and interactions are analysed. The system is evaluated for the case of nearest–neighbours hopping on d–dimensional hypercubic d = 2 (SQ) and d = 3 (SC) nonfrustrated lattices. For SQ lattice the effects of phase fluctuations on the properties of the superconducting state are discussed within the Kosterlitz–Thouless (KT) scenario. The KT critical temperatures TKT are calculated and compared to the Tc obtained within the BCS–HFA. The Uemura type plots, TKT vs. λ2(0), are obtained and evaluated as a function of pairing strength and concentration. The results are compared with the case of repulsive J (J < 0) pair–hopping interaction.


Introduction
In this work we are interested in the pairing in the real space characterized by small coherence lengths. Superconductivity with short coherence length is appropriate to a description of the high-temperature superconductivity in materials such as the cuprates, doped bismuthates, ironbased systems, Chevrel phases, organic metals, fullerenes, etc. (for a review, see e.g. Refs. [1,2,[4][5][6] and references therein).
One of the models often used for phenomenological description of dynamics of tightly-bound charged pairs is the attractive Hubbard (AH) model. In the AH model the effective on-site attractive interaction (U < 0) between charge carriers is responsible for pairs creation. Other models used to study the physics of pairing in the real W. R. Czart czart@amu.edu.pl 1 Faculty of Physics, Adam Mickiewicz University in Poznań, Umultowska 8561-614, Poznań, Poland space include the pair-hopping interaction J [8,9,11]. This interaction between the fermion pairs located at different sites (intersite charge exchange) can drive pairs forming and condensation.
In the present paper we analyse the Penson-Kolb-Hubbard (PKH) model which includes both the on-site term (U) as well as the pair hopping interaction J. This is conceptually a simple effective model for studying correlations, electron orderings and superconductivity in narrow band systems with almost unretarded short-range pairing [6, 9, 11-14, 16-19, 21, 23]. The model Hamiltonian is as follows: exchange interaction between pairs, the chemical potential is represented by μ, n iσ = c + iσ c iσ , and In the model J term competes with U interaction. When U = 0 Hamiltonian (1) transforms to the Penson-Kolb (PK) model [8,12,14,17,[25][26][27][28][29], whereas for J = 0 and U < 0 one gets the Hamiltonian of the AH model [1,30,31,33]. So far, the PK and PKH models have been investigated only in a few limits [2,9,11,13,[35][36][37]. The main efforts concerned characteristics at the ground state (T = 0) at halffilling (n = 1) in one dimension (d = 1) [9,11,35]. At T = 0, for both signs of J and U, diagrams of the PKH model on d = 1 and d = ∞ hypercubic lattices have been evaluated within the broken symmetry Hartree-Fock approximation (HFA) and by the slave-boson mean field method in Ref. [11]. As demonstrated, the d = 1 diagrams at n = 1 are occupied by at least nine distinct orderings. Among them: charge-density-waves and superconducting states, site-and bond-located antiferromagnetic phases with their mixed states as well. With increasing lattice dimensionality, the bond-type orderings gradually disappear at d = ∞.
The d = 1 PK model T = 0 phase diagram obtained within the HFA at n = 1 [11] is consistent with the respective diagram derived within exact Lanczos diagonalizations [27,39], density matrix renormalization group method [25,26] and the continuum limit field theory calculations [9,39]. For J > 0, these approaches yield the continuous 2ndorder transition to s-wave pairing phase at J = 0 + without any additional transition within the entire range of J > 0 interaction. This behaviour, at least for alternating lattices, is band-filling independent and it is retained for higher dimensions (including the exactly solvable case of d = ∞) [11].
We have already partially studied the superconducting characteristics of the PK and PKH models for the case of s-wave (Cooper-pair centre-of-mass momentum q = 0) [12,13,18,40,41] and eta-phase (q = Q [16][17][18] superconductivity. Preliminary investigation of the nextnearest neighbour hopping t 2 influence on both eta-and s-phase superconductivity in the 2-dimensional PK model has been also evaluated [18]. We have examined in that paper, the evolutions of the J vs. n phase diagrams and determined the T c and T KT transition temperatures with n and J. We have also derived, for t 2 = 0, the Uemura-type plots for the s-phase. We have analysed, for the PKH model on the nonfrustrated d = 2 SQ lattice, T = 0 phase diagrams involving homogeneous magnetic, charge-ordered and superconducting phases [23,42,43]. Moreover, within the variational approach (where the U term is treated exactly), the superconducting characteristics (for both s-wave and eta-phase) in the zero-bandwidth limit have been analysed [2,6,16,19,44].
In this work we extend our studies of the PKH model and analyse in details the evolution of the model's properties for arbitrary concentration (0 < n < 2) for the attractive J (J > 0). Taking into account both the magnetic orderings and superconducting s-wave states, we investigate the impact of the on-site U interaction on the phase diagrams both at T = 0 and T > 0. The ground state diagrams involving only homogeneous phases are compared with those involving phase separated states. In both types of diagrams we mark a crossover to the Bose-Einstein Condensate (BEC) regime. Moreover, we present here an extensive analysis of the U interaction effect on the electromagnetic as well as thermodynamic characteristics of the s-phase. The evolutions of the properties with n and interaction parameters are determined. We compare the theoretical London penetration depth λ(n) plots with experimental measurements reported by Locquet [48] for La 2−x Sr x CuO 4 films.
For the PKH model on the 2-dimensional lattice, within the Kosterlitz-Thouless (KT) scenario, we investigate the phase fluctuation impact on the s-phase. The d = 2 KT transition temperature (T KT ) is calculated from the superfluid stiffness ρ s evolution with temperature T, which is compared to the KT relation between ρ s and T KT . The KT scenario has already given reliable results for the models with attractive intersite density-density interactions [50, 52,53] and the d = 2 AH model [31,33]. In the former case a correct behaviour of T KT vs. -U-/B in agreement with the available Quantum Monte Carlo (QMC) data was found [33]. The T KT is compared with the critical temperature T c obtained within BCS-HFA. The gap to critical temperature ratio is determined and the Uemura-type "universality plots" [55,57]: T KT vs. ρ s (0) derived within the KT theory are examined. The theoretical Uemura-type plots are compared with the experimental data corresponding to Chevrel phases, doped BaBi0 3 (alternating cubic lattices) and cuprates (SQ lattice) [58,59].
The superfluid properties of the PKH model were studied within a linear response theory [1,30,60,62,64], where the electromagnetic kernel was treated within the HFA-Random Phase Approximation (HFA-RPA) scheme. The HFA-RPA approach is known to be reliable in the case of Hubbard model and its various extensions in the area of the ordered states studied at T = 0, within the whole range of interaction interpolating between the weak and strong coupling limit [1,11,30]. Furthermore, for the electronic models with exclusively intersite interactions in the d = ∞ limit, the HFA leads to exact solutions for any For d < ∞ at T > 0, the HFA is much less reliable especially in the strong coupling limit and in low-dimensional systems as it neglects short-range correlations and phase fluctuation effects.
Hence the HFA neglects phase fluctuation effects and short-range correlations, it is much less reliable in d < ∞ at T > 0, especially for low-dimensional systems in the strong coupling limit.
In a separate work we will study the case of the PKH model with J < 0, which can stabilize eta-pairing.
Organization of the paper is as follows. In Section 2 there are the electromagnetic kernels and basic equations calculated to determine the basic properties of the studied system in the normal and s-wave superconducting phases. In Sections 3, 4, and 5 we present the numerical analysis of the system. The phase diagrams at T = 0 involving s-wave superconducting phase as well as nonordered, antiferromagnetic and ferromagnetic states are given in Section 3.1, and the T > 0 phase diagrams are analysed in Section 3.2. The ground state s-wave phase properties as a function of concentration and coupling strength for SQ and SC structures are shown in Section 4. In Section 5 we devote to superconducting properties analysis at the finite temperatures for the SQ lattice. Going beyond the HFA we calculate the KT critical temperatures T KT and compare them with T c obtained from BCS-HFA. In Section 5. we also demonstrate the Uemura-type "universality plots". Finally, Section 6 contains concluding remarks. In addition, Appendix A is devoted to free energies and self-consistent equations for the phases competing with s-wave paring.

General Formulation
Taking into account the electrons coupling to magnetic field via the vector potential, the PKH model Hamiltonian takes the following form [9, 11-13, 29, 35]: where the Peierls factors incorporate electrons coupling to magnetic field via its vector potential A(r) (diamagnetic effect): e is the electron charge, and < ij > limits the sum to nearest neighbours (nn). Other denotations are given after (1). Following our earlier works on this matter we assume that the effective (phenomenological) model parameters t, U and J contain all possible renormalizations and contributions including those coming from the strong electronphonon or from electrons and other electronic subsystems in chemical complexes or solid couplings [1] (e.g. coupling of electrons with intermolecular vibrations via modulation of the hopping integral [65] or the on-site hybridization in a generalized periodic Anderson model [67,68]). In general, their values and signs could be arbitrary. Particularly interesting is the role of J interaction in a multiorbital model, because of its presence in the iron pnictides [70]. In the strong interaction limit, the J term may also be included in the effective models for Fermi gas on an optical lattice [69].
The self-consistent equations for electron thermal averages and free energy of the PKH model (1) can be determined (within the HFA approximation) by the standard equation of motion or Greens functions approaches [11].
The s-wave phase superconducting order parameter the Fock term where γ k = 2 α cosk α , α = x, y, . . . , and μ are determined by the set of equations: Free energy of the superconducting (S) phase is derived as: The explicit form of the self-consistent equations determining x s , p and μ is the following: where E k = ¯2 k + (−U + J 0 ) 2 x 2 s is a quasiparticle energy, k = k − μ, chemical potential: μ = μ − 1 2 Un, N-number of sites in the system, n = N e N (0 < n < 2) denotes number of electrons per site, J 0 = zJ, k = −tγ k ,t = t + 2pJ/z, z is the number of nearest neighbours (z = 6 and z = 4 for SC and SQ structures respectively), β = 1/k B T.
Free energy of the normal (N) phase has the following form: where μ N and p N are determined as follows: and The equations of the other phases considered in this work are given in Appendix A. However, it should be underlined that all results presented in this work are obtained neglecting effects of the Fock term in all phases, i.e., assuming that p = 0 in each phase. Based on the linear response theory calculations in the weak vector potential A(r) case, the expected value of the total current operator Fourier transform in the direction α (α = x, y, z) [12,13,60,71] takes the given form: The total reaction kernel K αβ (q, ω) has been split up into the paramagnetic and diamagnetic parts, which were subsequently evaluated within the RPA-HFA scheme.
In particular, we get the diamagnetic part as follows: where x s , p and μ are given by (9)- (11). The expression for the paramagnetic part of the kernel is calculated to be where n F (E k ) is the Fermi-Dirac distribution function.
Within London limit the magnetic penetration depth λ is defined from the transverse part of the total kernel which in the static limit reads as follows: where K dia x is given by (15) and for the transverse part of the paramagnetic kernel in the static limit for q → 0 we obtain: At T = 0 K the paramagnetic part is significant when calculating λ for non-local superconductors (Pippard superconductors), common among systems with low T c . The materials with short coherence length, including superconductors with high T c , represent the so-called London limit. For q → 0 and at T = 0 the magnetic penetration depth in the London limit is expressed with the diamagnetic component of the total kernel: We restrict our investigation to analyzing the system in the London limit in which we can apply the local approximation. Derived for this case λ is qualitatively well represented in the weak and strong interactions, moreover, for strong couplings, the results obtained in HFA coincide with those obtained from the perturbation theory.
Having the difference between the normal and superconducting phase free energys and the value of penetration depth we can derive formuls for the thermodynamic critical field H c and the Ginzburg-Landau correlation length ξ GL : where a is the lattice constant, 0 = hc 2e , and we get the estimations for the critical fields: where κ is the Ginzburg ratio: . The HFA expressions determining normal and superconducting ground state energies, E N 0 and E S 0 are: where is the mean kinetic energy of the superconducting (normal) phase.
From (9)-(11) and the condition → 0, where gap amplitude parameter: = (−U + J 0 )x s , one can calculate the HFA transition temperature T c which is an estimation of the temperature of pair-formation [12,33,52].
Due to the phase fluctuation effects, for finitedimensional lattice ( d < ∞), transition to the superconducting phase occurs at temperature lower than T c (HFA). We take into account the phase fluctuation effects by determining another characteristic temperature, the Kosterlitz-Thouless critical temperature (T KT ). For two-dimensional lattices (d = 2) the value of T KT can be derived within the KT theory [12,31,33,72,73]. In terms of the theory the T KT temperature is defined by vortex pair unbinding transition determined by universal jump of the superfluid stiffness (helicity modulus) ρ s at T KT : where Q 0.898 (Monte-Carlo estimates for d = 2 XY model [33,72]). The superfluid stiffness ρ s , which is directly related to the London penetration depth λ and calculated within HFA-RPA scheme using (15), (18) and (17) has the following explicit form: Thus, T KT is determined by (9)-(11) and (26) with ρ s given by (27). Since the formula (27) does not take into account the renormalization coming from the topological excitations (vortex-antivortex pairs) [31,33], calculated this way T KT is only an upper bound estimation of the actual KT transition temperature.
Using our approach one is able to analyse the crossover between BCS-like (extended Cooper pairs ("BCS")) and local Cooper pairs superconductivity (composite bosons, Bose-Einstein Condensate (BEC)). The crossover takes place when increasing pairing interaction (J or -U-(U < 0)) one increases coupling from a weak to strong regime. At the ground state, one can approximate boundary between both regimes from the Leggett criterion [30,74], which requires the chemical potential of the superconducting phase to be located at the bottom of the electronic band, i.e.: where B is the bandwidth of k (B = 2zt) and μ s is determined from the self-consistent (9)-(11) solved at T = 0.

Phase Diagrams of the Model at T ≥ 0
We performed extensive analytical and numerical analysis of the electromagnetic and thermodynamic properties of the s-wave phase of the PKH model (3) for d-dimensional hypercubic lattices.
In this Section we discuss the phase diagrams for arbitrary n and interaction parameters, determined within the HFA in the ground state (Sec. A) and T > 0 (Sec. B), whereas the superfluid characteristics of s-wave phase are analysed in Sections 4 and 5.
For d = 2 we used the square lattice (SQ) density of states expressed with: if -/4t-< 1 and zero otherwise, where K is the complete elliptic integral of the first kind. For d = 3 we use the analytical approximation to D() as calculated numerically by Jelitto [75].

Ground State Phase Diagrams
Below we analyse phase diagrams at T = 0 for the SQ lattice involving the homogeneous: nonordered (N), swave pairing (S), ferromagnetic (F) and antiferromagnetic (AF) states and compare them with the phase diagram involving also phase-separated states. We only analyse the simple magnetic phases but consideration of more complex magnetic orderings would not modify significantly the results. In general, in the ground state, the S states are always stable with respect to the N-phase for J > U/z for arbitrary values of single electron hopping, U, J and n (0 < n < 2). In the figures we mark the crossover to the Bose-Einstein Condensate (BEC) which at T = 0 can be located from Leggett's criterion [30,74] (28). For selected interaction values we also plot the evolution of the order parameters of the respective phases as a function of n. Figures 1 and 2 present the T = 0 phase diagrams U versus n (for J/4t = 1) and J versus n (for U/4t = 2 and U/4t = 3) respectively. The transition lines to N phase are of the 2ndorder, while those between the ordered states are of the 1st-order. As a result of competition between repulsive U and J, the s-wave pairing state is limited on the diagram, Fig. 1, for U < U c (with U > 0), whereas the F and AF phases are stable for U above U c (U > U c ). The F and AF phases occur in the diagram in the limited ranges of concentration n: for low and intermediate values of n and close to half-filling (n = 1), respectively. For U < U c (n = 1) states S may occur within entire range of n. Because of the strong interactions the nonordered phase is not present on the phase diagram. For given concentration n, with increasing on-site U interaction, system exhibits several sequences of transitions: (i) a single 1st-order transition S → F, (ii) a sequence of 1st-order transitions S → AF → F (with increasing U, AF → F transition line asymptotically approaches half-filling) and (iii) at n = 1, a single 1st-order transition S → AF.
The J vs. n phase diagrams, for fixed U at T = 0 is shown in Fig. 2. Analogously to the previous diagram ( Fig. 1), the S-phase may occur in the diagram above a critical J c (U, n). Increasing U reduces concentration ranges occupied by the AF and N phases and expands the range of F state stability in the diagram (comp. Fig. 2a and b). In both diagrams, with increasing J, the system can exhibit several different types of behaviour, depending on n: (i) a single 1st-order transition N → S, (ii) a single 1st-order transition F → S and (iii) a single 1st-order transition AF → S. Figure 3b presents the phase diagram U versus n (for J/4t = 1.0) involving the phase-separated states (PS). This diagram is a preliminary result, a part of a broader analysis of phase diagrams involving PS states, which is soon to be published elsewhere.. In the figures presented here, for the sake of clarity, the μ-dependent plot (Fig. 3a) is given, which helps to determine boundaries between the PS states plotted in Fig. 3b. All homogeneous ordered states are separated by PS states (Fig. 3b). The transitions between homogeneous and the PS states are symbolically named "third order" transitions [6]. In the system that undergoes this type of transition, one of the domains in the PS phase gradually vanishes.
In analogy to earlier presented diagrams involving homogeneous states only (cf. Fig. 1), also in this case (Fig. 3b) the states involving s-wave pairing are stable only below certain critical values U c . For U > U c (n = 1) superconducting S-states can be stable only within a limited range of n (n < 1) and appear in the diagram as a homogeneous or phase-separated states involving AF and F orderings. At half-filling, for U < U c the S-phase may occur within the whole range of concentrations. The homogeneous AF phase occurs only for n = 1 when U > U c . Away from n = 1, within a limited range of U, the AF phase may be stable in the diagram exclusively within the PS(S/AF) and PS(AF/F) states close to n = 1.
As the electron density increases towards half-filling, depending on U, the system exhibits a few distinct transition sequences: (i) a single 1st-order transition F → PS(AF/F) → AF, (ii) a sequence of four 1st-order transitions: S → PS(S/F) → F → PS(AF/F) → AF and (iii) a sequence of two 1st-order transitions: S → PS(S/AF) → AF.
In the plots (Figs. 1, 2 and 3) the crossover to BEC regime is denoted with dashed line. An obvious effect of the repulsive (attractive) U is the reduction (expansion) of the composite bosons (BEC) regime. As we see from the plots, with rising concentration n, the boundaries are moved towards greater values of interaction J for any U. Diagram is plotted for SQ lattice Fig. 2 The J and n phase diagrams plotted for U/4t = 2 and U/4t = 3 at T = 0. All the transitions involving exclusively ordered states are of the 1st-order and those to nonordered (N) phase are of the 2nd-order. Diagrams are plotted for SQ lattice. Other denotations as in Fig. 1 Thus, within a definite range of U an J, the crossover from the weakly to the strongly coupled electron pairs can be achieved by lowering the electron density.

Finite Temperature Phase Diagrams
In this subsection we present a two examples of the finite temperature phase diagrams, as evolution with concentration ( Fig. 4) and evolution with U (Fig. 5) for the repulsive on-site interaction U > 0.
For interaction parameters U/4t = 3.0 and J/4t = 1.0 ( Fig. 4) with increasing temperature the system can exhibit the following phase transitions: (i) a single 2nd-order transition S → N (ii) a sequence of three transitions: 2ndorder S → N and reentrant 2nd-order N → F → N, (iii) a sequence of five transitions: 2nd-order S → N, reentrant 2nd-order N → F → N, reentrant 2nd-order N → AF → N, (iv) a sequence of 1st-order transition F → AF and 2nd-order AF → N and (v) a single 2nd-order transition AF → N.
The diagram T versus U for J/4t = 1.0 and n = 0.9 is given in Fig. 5. As one can notice from the plot, with increasing T the system can undergo either (i) a single 2nd-order

Superfluid Characteristics of s-wave Phase at T = 0
Below we present the evolutions of the ground state superconducting properties of the s-wave phase with model parameters (Figs. 6, 7, 8 and 9). The plots were determined for the parameters' ranges where superconducting S-phase is stable in the phase diagrams analysed earlier in the paper. Due to the electron-hole symmetry of the model all the plots are symmetric with respect to the transformation n → 2−n. As one can see from the figures (cf. Figs. 1-3), where the superconducting characteristics are plotted as a function of U/B for J 0 /B > 0, the S states can still be stable even for U > 0 (repulsibe U) when J 0 − U > 0. With increasing J 0 /B the S-phase boundary shifts towards higher values of U/B.
In Fig. 6 we show plots of the inverse square value (1/λ 2 ) of the London penetration depth and of the order parameter as a function of U/B for fixed values of n and J 0 /B at T = 0 for the SQ and SC lattices. These quantities plotted as a function of n for nonfrustrated 2D SQ lattice and selected values of J 0 /B and U/B are presented in Fig. 7. As the pairing interactions J 0 or -U-(for U < 0) increase, the 1/λ 2 , and other superconducting characteristics evolve smoothly from the limit of weakly interacting single particle carriers (with λ − 2 being proportional to the bandwidth B,   λ − 2 continuously decreases with increasing -U-/B for any -U-and λ −2 ∼ t 2 /|U | for -U-/B >> 1. For fixed J 0 /B > 0 following the increase of -U/B, the λ − 2 first grows, and subsequently it passes through a round maximum and drops at high -U/B proportionally to 2t 2 /-U-+ J (cf. Fig. 6). With increasing J 0 /B the maximum value of λ − 2 increases. Notice that the increase in λ − 2 as a function of J 0 /B for any value of U/B (the increase is linear for large J) is significantly distinct form the behaviour observed for the AH model (J 0 = 0, U < 0).
The following figures present the energy gap E g (0) and the behaviour of the square value of thermodynamic critical field H 2 c at T = 0 for SQ lattice. The quantities are plotted as a function of n for a fixed value of J 0 /B and several values of U/B (Fig. 8a) as well as a function of U for a fixed n and several values of J 0 /B (Fig. 8b). Both characteristics vanish at the empty-band limit (n = 0), reach a maximum at half-filling and increase linearly with increasing -U-/B for U/B >> 1, and in the weak coupling limit they are proportional to exp( −2B J 0 −U ). The plots of G-L coherence length ξ GL and the Ginzburg ratio κ versus U/B and n calculated for SQ lattice for selected values of respective interaction parameters are presented in Fig. 9. In the limit of weakly interacting single particle carriers the As concerns the concentration dependencies of various quantities (see Figs. 7, 8 and 9) one finds that in the low density limit for arbitrary coupling strengths λ −2 ∼ n (as for the fermions in the continuum), (2−n) , for any n. Beyond these limits the concentration dependencies are not universal. Their behaviour highly depend on the features of density of states. Notice, the resulting from the van Hove singularity in D() for SQ lattice, substantial enhancement of and H 2 c close to n = 1 (Figs. 7 and 8) observed in the range of small and intermediate values of (J 0 − U)/B. Similar efect is observed for the local maximum of the Ginzburg ratio κ at n = 1 (Fig. 9a). In addition, in Fig. 9a one can see a significant variation of ξ GL as a function of n in the weak and intermediate coupling strength.

Critical Temperatures for s-wave Phase and the Uemura-Type Plots
Apart form the weak coupling as well as from the case of d = ∞, the phase fluctuations exert a strong impact on the superconducting pairing, and the results display a clear energy scales separation for the pair formation ( ∼ k B T c ) and phase coherence. The effects of the phase fluctuations are most clearly seen for d = 2 lattices, whose phase coherence temperature is determined by k B T KT (26), (27).
Between T c and T KT incoherent s-wave phase occurs, where the pseudogap opens up in the quasiparticle energy spectrum and the system exhibits non-Fermi liquid properties. Not at T c , but only at T KT the real transition to a superconducting state takes place, where the transition to a phase with bound vortex-antivortex pairs occurs. Certain pairing correlations are observed at all temperatures thus T c does not represent a rigorous boundary line. It is rather a reliable estimate of the pair formation and the pseudogap emergence [12,30,36].
For SQ lattice the evolutions of T KT and T c with a change in n for fixed interaction parameters are shown in Fig. 10, whereas Fig. 11 present the plots of transition temperatures as a function of interactions. The T KT temperature can be significantly reduced compared to T c , where for U ≤ 0, J 0 > 0 the highest reduction occurs for low values of carriers density n (cf. eg. Fig. 10 and, for U = 0, Ref. [12]). Moreover, the difference between T c and T KT strongly increases with the increase of on-site attraction U < 0 (cf. Fig. 11). Analogously as found for PK [12] and AH [30] models, also for the PKH model one observes a monotonous increase in T c with increasing interactions, (which is obviously a result of increasing pair binding energy) (cf. Fig. 11 [12]). On the other hand, apart from the weak paring strength regime T KT evolution with interaction parameter is significantly distinct for these models. As we can see in Fig. 11, for the PKH model, similarly as we observe for the AH model, with interaction -U-increase the T KT at first increases exponentially, and subsequently it passes a round maximum and then drops as t 2 /-U- at strong -U-. Analogous behaviour was found for the thermodynamic critical field H c (0) [12,30,33]. On the other hand in the PK model, the model with intersite charge exchange, T KT monotonously increases with J, and the growth of T KT is linear at strong values of J (cf. Ref. [11,12]). Due to the difference in the paring mechanisms, the dynamics of electron pairs is qualitatively different in both models. [11,12,30]. The PKH integrates both of the pairing mechanisms and the interplay between U and J interactions results in the differences in the thermodynamic and electrodynamic properties compared to those in the PK and AH models. In contrast to the BCS case, the phase fluctuations enhance the ratio of the ground state gap (E g = minE k − max(−E k )) to the critical temperature. In Fig. 12a we Notice that for T KT this ratio can be much greater than for T c . We have found, the enhancement of the ratio to be the highest at low concentrations, where it can be higher even by a factor 10 or more. Only for J 0 −U t → 0 we recover the BCS ratio E g (0)/2k B T c ∼ 1.76. Let us stress that the gap ratio goes to infinity for the local pair (the strong coupling) regime if n → 0, since E g (0) remains finite in such a case (cf. Fig. 8a) and T KT/c → 0.
Within the KT scenario we also derived, for SQ lattice, the Uemura-type plots for the superconducting phase of the model, and examples of the T KT/c vs. 1/λ 2 figures for controlling parameter n with fixed values of J 0 /B and U are shown in Figs. 13 and 14. Apart from the weak coupling regime the curves plotted on the figures have shapes analogous to the Uemura's plots experimentaly obtained [55,57] for various classes of the "exotic" superconductors (the superconductors with short-coherence length) such as bismuthates, the organic materials and cuprates. For   (26)). It should be emphasized that T c is not related to ρ s and analogous plots with T c cannot account for the scaling [12, 33, 52,    (⊕)-Tl 2 Ba 2 CuO 6 + δ (for overdoped range). Experimental data taken from [58,59] 53]. Let us point out that for t 2 = 0 [18] the densities of states for SQ and SC lattices are no longer symmetric and the function T KT (1/λ 2 ) does not represent the one-valued function, whereas it exhibits substantial hysteresis due to differences in positions of the maxima values of the T KT (n) and 1/λ 2 (n).
In Fig. 14 we compared our results with experimental data corresponding mainly to SQ lattices (cuprates) and alternating cubic lattices (Chevrel phases, doped BaBi0 3 ) [58,59,76,77]. In order to compare the results we have normalized T c and T KT to unity by scaling both temperatures over respective maximal values of the temperatures: T m c and T m KT , and 1/λ 2 we scaled over a value 1/λ 2 m corresponding to the value of the respective maximal temperature. The experimental data were scaled analogously. Except for the overdoped cases the data compare best with the plots obtained for the values of coupling between strong and intermediate. The experimentally found deviations from the universal dependence in the strongly overdoped regime can be explained either by the asymmetry in D(ε) due to further neighbour hopping in the considered lattices or by the effects of intersite Coulomb interactions [79], which can be different for various families of materials. In Fig. 15 we qualitatively compare the theoretical plots λ(n) obtained for the discussed simple effective model with penetration depth measurements reported by Locquet [48] for La 2−x Sr x CuO 4 films. The experimental data extend from the regime of heavily underdoped to heavily overdoped. The theoretical results fit well the experimental measurements.

Conclusions and Remarks
We have explored in this work the ground state and final temperature phase diagrams as well as the superconducting properties of the Penson-Kolb-Hubbard model. The study have been performed for arbitrary electron density and J > 0. In this model, the intersite charge exchange interaction (a non-local interaction) is responsible for pairs formation and condensation.
The PKH model can be treated as an effective model of the s-wave superconductors with short coherence length, which can be of order of interparticle distance or even a lattice constant and can be applied to isotropic d = 3 structures (such as the barium-bismuthate compounds [80,82], Chevrel phases [1] and fullerides [83][84][85][86][87]), as well as those with layered (d = 2) structures (such as heavy fermion systems and some organic superconductors). Moreover, the model may also lead to some qualitative understanding of the d = 2 high-T c cuprates [76,88,91], iron pnictides and chalcogenides [5,93].
For the SQ lattice we have determined the stability of s-wave pairing state with respect to F, AF and N phases. We have presented the ground state diagrams for homogeneous phases only (without considering PS states) and a preliminary phase diagram involving PS states (Figs. 1-3). In the diagrams the phase transitions to nonordered state are of the 2nd-order, and all the transition lines separating considered ordered states are of the 1storder. Increasing U/4t broadened F phase concentrations n range and reduced the n ranges of AF and N phases occurrence. We have shown that at T = 0 for J > U/z, in the general case of arbitrary U/4t and J/4t the S-phase was always stable with respect to the N phase for any value of the concentration n and single electron hopping.
In the diagram involving PS states (Fig. 3) we have shown that the homogeneous ordered phases were separated by the respective PS states and homogeneous AF phase could occur only at half-filling. For n = 1 the AF ordering could appear only within the phase separated states: PS(S/AF) and PS(AF/F) in a limited extent of concentration n near to the half-filling.
At T = 0 the isotropic s-wave pairing, S-phase, stable for the attractiveJ (J > 0) and the attractiveU (U < 0) is found to exhibit a smooth crossover between the BCS-like regime and the tightly bound pairs limit with increasing pairing interactions. This behaviour is qualitatively different from the one observed for repulsive J (J < 0) when the transition into the eta-phase occurs only above a critical value -J cdependent on band filling (n) and lattice structure. In this case the system never manifests standard BCS-like features [16][17][18]. In the figures we have marked the crossover to the BEC which at T = 0 was located after Leggett's criterion. Within homogeneous S-phase, for a given n, the system could undergo BCS-BEC crossover with decreasing U. The effect of the repulsive (attractive) U is the reduction (expansion) of the composite bosons (BEC) regime. We have demonstrated that in a specific range of U and J the crossover from the weakly to strongly coupled electron pairs could be achieved by lowering the n. With decreasing n the crossover regime is shifted towards lower values of pairing interaction. Notice that the crossover always occurs beyond the range of stability of the PS states involving s-wave domains (cf. Fig. 3). We have also presented examples of the finite temperature phase diagrams involving homogeneous states as a function of concentration and interaction U (Figs. 4 and 5).
The evolution of the thermodynamic and electromagnetic superconducting properties with electron concentration and pair interactions has been analysed for d-dimensional hypercubic lattices (primarily for SQ lattice) and the changes in the behaviour of these characteristics at the crossover between the BCS-like limit and the local pair regime superconductivity have been discussed. We have determined the evolution of the critical temperatures T c (Hartree-Fock), T KT (Kosterlitz-Thouless) as well as the critical fields, the penetration depth, the coherence length and Ginzburg ratio at T = 0 as a function of particle concentration n and interactions.
As the pairing interactions J 0 or − U increase the swave phase characteristics evolve smoothly between the regimes of the weakly interacting single particle carriers and of tightly bound pairs. The λ − 2 increases with increasing J 0 /B for any value of U/B, and the increase is linear for large J. This behaviour is qualitatively different from that found for the AH model (J 0 = 0, U < 0), where λ − 2 continuously decreases with increasing -U-/B and λ −2 ∼ t 2 /|U | for -U-/B >> 1. In the strong interaction regime the electromagnetic and thermodynamic characteristics of the S-phase (for J > 0) become similar to the properties of the eta-phase (for J < 0). For the arbitrary coupling strengths in the low density regime λ −2 ∼ n as for the fermions in the continuum.
The PKH model reduces to PK model for U = 0 and for J = 0 reduces to AH model. The evolution of superconducting characteristics with increasing pairing interaction is essentially different in both systems. For the PK model, in the limit of strong coupling, values of λ − 2 , H c , T KT , κ − 1 increase with J while for the AH model the characteristics decrease with |U |.
Our results obtained for the case U = 0, J > 0 are consistent with those obtained for collective excitations performed using a generalized random-phase approximation [28,36]. The collective mode velocity, in the PK model, rise with J [28] while in the AH model it drops with |U | [36]. Moreover, for the studied models, the conclusions are in qualitatively good agreement with the results of perturbational expansions both in the weak and strong interaction regimes [30,94,95].
The S-phase can be stable for repulsive values of U (0 < U < U c ) as well. In the PKH model, repulsive U modifies the S-phase characteristics (cf. also [17]) and shrinks the extent of the superconducting phase stability. It also can induce 1st or 2nd order transitions to various nonsuperconducting phases.
Apart from the weak coupling regime one finds a substantial increase in the gap to critical temperature ratio resulting from the phase fluctuations and the energy scales separation for the pairs formation ( ∼ k B T c ) and the phase coherence ( ∼ k B T KT ). Between T c and T KT a disordered (incoherent) pair state is possible in which the single-electron excitation spectrum has a gap, but the pairs are phase disordered. The T KT can be substantially lower than the T c .
An important result of our work is the derivation of Uemura-type plots for the model considered. Beyond the regime of weak coupling and low concentrations the results show universal scaling: T KT ∼ 1/λ 2 (0), similarly as it was found for the attractive Hubbard model [33], the PK model [12] and the extended Hubbard model with intersite attraction [52,53]. Apart from the weak coupling regime, the resulting plots have a shape consistent with the experimental Uemura's plots for various classes of superconductors with short-coherence length [17]. Our results support the conclusion [12,30,52,53] that in materials with short coherence length, the ground for the Uemura scaling is a separation of the energy scales of the pairs formation and the phase coherence in the underdopped limit or in the strong interaction regime. from which we get sets of self-consistent equations for each considered ordering type. The order parameters for the phases considered: (i) antiferromagnetic (AF): (ii) ferromagnetic (F): The equations for the phases considered derived for the alternated lattice, (i.e. k + Q = − k ) in the (broken symmetry) HFA: the free energy of the ferromagnetic phase F F : the thermodynamic potential of the ferromagnetic phase Ω F : the set of self-consistent equations: (i) the equation for the number of particles of the ferromagnetic phase: (ii) the equation for the order parameter of the ferromagnetic phase: where E ± F (k) = k ± Ux F , k = k − μ, μ = μ − Un/2. The free energy of the antiferromagnetic phase F AF : the thermodynamic potential of the antiferromagnetic phase Ω AF : the set of self-consistent equations: (i) the equation for the number of particles of the antiferromagnetic phase: (ii) the equation for the order parameter of the antiferromagnetic phase: , (40) and the Fock term: where E ± AF (k) = μ ± 2 + U 2 x 2 AF .