Advances in Condensed Matter Physics
Volume 2015 (2015), Article ID 835897, 9 pages
http://dx.doi.org/10.1155/2015/835897

Research Article

On New Families of Integrals in Analytical Studies of Superconductors within the Conformal Transformation Method

Ryszard Gonczarek,¹ Mateusz Krzyzosiak,² Adam Gonczarek,³ and Lucjan Jacak¹

¹Faculty of Fundamental Problems of Technology, Wrocław University of Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
²University of Michigan and Shanghai Jiao Tong University Joint Institute, 800 Dongchuan Road, Shanghai 200240, China
³Faculty of Computer Science and Management, Wrocław University of Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland

Received 19 April 2015; Accepted 2 July 2015

Academic Editor: Jaume Gázquez

Copyright © 2015 Ryszard Gonczarek et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We show that, by applying the conformal transformation method, strongly correlated superconducting systems can be discussed in terms of the Fermi liquid with a variable density of states function. Within this approach, it is possible to formulate and carry out purely analytical study based on a set of fundamental equations. After presenting the mathematical structure of the -wave superconducting gap and other quantitative characteristics of superconductors, we evaluate and discuss integrals inherent in fundamental equations describing superconducting systems. The results presented here extend the approach formulated by Abrikosov and Maki, which was restricted to the first-order expansion. A few infinite families of integrals are derived and allow us to express the fundamental equations by means of analytical formulas. They can be then exploited in order to find quantitative characteristics of superconducting systems by the method of successive approximations. We show that the results can be applied in studies of high- superconductors and other superconducting materials of the new generation.

1. Introduction

Over the past few years, research on superconductivity has been gaining new momentum due to the milestone discovery of iron-based superconductors [1, 2]. Along with the ongoing work on high- materials, studies on superconducting systems of the new generation remain one of the most active research areas in superconductivity research. Progress in fabrication and characterization techniques is mutually coupled with efforts to model and predict theoretically the trends and phenomena observed in experiments. An important tool for quantitative theoretical studies on superconducting systems is the gap equation accompanied by the carrier concentration equation. It appears in similar forms in the BCS theory, the Eliashberg equations, the Van Hove Scenario, and the conformal transformation method, with the latter taking fully into account the low-dimensional structure of high- materials. In order to evaluate fundamental thermodynamic functions such as the free energy, the thermodynamic potential, the entropy, the heat capacity, and other quantitative parameters such as the transition temperature, the heat capacity leap, the critical current, the critical magnetic field, and the isotope shift, one should study the gap equation and the carrier concentration equation along with the corresponding formulas for the thermodynamic potential, the critical fields, the concentration of superconducting carriers, and so forth [3–11].

Employing the gap equation with a fixed chemical potential equal to the Fermi energy, in theoretical studies of BCS-type superconducting systems at subcritical temperatures, Abrikosov and Maki have suggested to apply an expansion in the integral equation and restrict it to the first-order terms. The calculated integral allowed them to estimate the solution of the first boundary value problem for the gap of a BCS superconductor as a function of the temperature [12–14].

Contemporary studies focus, however, more on superconducting materials of the new generation. These include high- copper-oxide quasi two-dimensional superconducting systems and various doped superconducting compounds such as spinel- and perovskite-type structures of superconducting compounds of a trivalent rare-earth and a divalent alkali-earth ion. Other classes of superconducting materials that have been actively studied in recent years include novel superconducting compounds of MgB₂ with a C, Al, or Sc substitution, or organic superconductors with a controlled bandwidth and band filling [7, 11, 15–22]. Theoretical description of such superconducting systems requires to take into account spin-fluctuations or strong-correlation effects. These can be included by means of an effective Hamiltonian in the strongly interacting Hubbard model with a given (multiband) one-particle dispersion relation enriched by self-energy corrections, as well as a quite general form of the pairing potential. The latter can be decomposed into an antisymmetric and a symmetric part determining the symmetry of the order parameter [23–30].

In order to address the above problems in a systematic analytical manner, we propose to apply the conformal transformation, which allows us to consider a superconducting system as mapped onto a Fermi liquid in an isotropic momentum space with a naturally introduced function of the density of states [23]. In further discussion, a series expansion of fundamental integrals emerging in the formalism proves to be a powerful tool. In the present paper, we derive formulas for some infinite families of integrals, which allow us to construct such series up to an arbitrary order. We emphasize that until now only few of such integrals were given [12–14, 31–33].

2. Employed Formalism

Employing the Green function formalism, with the effects of an external magnetic field included, one can derive a set of two fundamental equations in the momentum space, consistent with the mean-field approximation. One of these equations is the gap equation [5, 7, 19–21, 32, 34]:

where and is the superflow, with being the phase of the order parameter and is the vector potential. The other one is the carrier concentration equation:

which determines the chemical potential . This equation allows one to express in (1) in dependence on and study some superconducting systems with a partially filled conduction band. The dispersion relation entering both equations is usually defined with respect to the Fermi level.

In the BCS theory, which was greatly successful in explaining various properties of simple low- superconductors by means of several experimental parameters, an electron system is treated as a Fermi liquid. However, high- materials are complicated, and many fundamental condensed matter physics phenomena play some role in one or another part of their phase diagram. The normal state near the composition corresponding to the highest shows a convincing evidence of being a more composed form of a quantum liquid than the Fermi liquid and should be considered as a non-Fermi liquid. In a wide class of non-Fermi liquids presented in [35] as applicable to modeling of superconducting systems, for some of them, a singular (or marginal) Fermi-liquid behavior is observed in the underdoped (low carrier concentration) regime, whereas the other ones exhibit it in the overdoped (high carrier concentration) regime, near to the optimally doped structure. Therefore, low- superconductors and overdoped high- cuprate superconductors, described by an effective Hamiltonian with a nonlinear dispersion relation , can be discussed in terms of the Fermi liquid after applying a conformal transformation of the momentum space.

In general, such an approach can be applied to anisotropic superconducting systems with an arbitrary dispersion relation for spin-singlet -wave or -wave symmetry states, as well as for the spin-triplet -wave symmetry state, with a fixed carrier concentration. It also allows one to include the effects induced by variations of the carrier concentration , and the pairing potential amplitudes in both the isotropic and the anisotropic channel, in order to discuss the stability of these symmetry states [5, 20, 27, 33–36].

Studies of the pairing potential seem to be particularly interesting and valuable, since a definite pairing mechanism in high- superconducting systems has not been clearly identified yet. However, in the literature, two mechanisms are usually discussed [26]. The first one is based on the single-band Hubbard model (or related models, such as the model) whereas in the other one it is suggested that the role of the electron-phonon interaction is crucial. Quantitative analysis of the pairing interaction in superconducting systems could therefore provide useful insights into the mechanism of Cooper pairs formation in superconducting systems of the new generation [7, 33]. In our approach, we have introduced the pairing potential amplitude in the form . Nevertheless, one should remember that is the momentum-dependent part of the full pairing potential, which is an antisymmetric matrix-function defined in the momentum-spin space. This function is obtained as a result of the Fourier transform of an effective two-particle interaction potential including all many-body effects, and ensuring the spin conservation law [36]. Hence, is a symmetric function of and , and in order to generate attraction needed to form a paired state in the system, this function must be negative. In the formulated approaches (cf. [4, 5, 14, 26]), the function becomes negative and represents the introduced pairing potential amplitude only in a small region of the momentum space near to the Fermi surface, and so (for a spin-singlet -wave state) it can be consider as constant.

2.1. Integrals of the Fundamental Equations for

It is worth to emphasize once again that the conformal transformation method [7, 23, 37] proves that an arbitrary anisotropic system, with a given one-particle energy spectrum, can be considered as an isotropic Fermi liquid in the momentum space. However, it needs then to be taken into account that quantum-mechanical states are not distributed homogeneously. The distribution is defined by a scalar field of the density of states, which for a fixed spin-singlet -wave state, with a given carrier concentration and a pairing potential amplitude , reduces to the density of states function. Therefore, by applying the conformal transformation method to an anisotropic superconductor with a given dispersion relation, one can map quasiparticle excitations onto Landau quasiparticles and consider the system as a Fermi liquid with a fluctuating density of states. Then, (1), after taking into account (2) in the manner discussed above, can be transformed to the form:

where denotes averaging over the angles in the scalar product and and . Here, is the dimensionless density of states function, and it is found from the given dispersion relation inherent to the studied anisotropic superconducting system. More precisely, it appears as the final product in the conformal transformation approach and can also include effects of the particle-hole asymmetry and an extra factor, as, for example, the spectral function defined in the Eliashberg equations. Since for the BCS model , we take it in the form in further discussion. According to previous studies [38], the function can substantially change various parameters of a superconducting system. The effects are particularly pronounced if the function has one narrow fluctuation (peak) in the vicinity of the Fermi surface. Moreover, peaks of various shapes yield similar results. Therefore, at this point, we do not need to consider any particular analytical properties of the function , treating it as a locally constant function.

In the next step, introducing the symbols , , , and , (3) can be rewritten as

where the function is redefined into , which can be again modeled in various manners. The factor can be eliminated by employing (4) in the limit cases and .

In order to obtain the full representation of the gap equation of the Abrikosov-Maki type when , we note that

where the differential operator should be read as

Hence, for further calculation, we define the integrals

where for we can replace by . Substituting , differentiating by , and putting , we find that

Let us now define an extended form of which will become useful in our discussion later on:

so . After applying the differential operator in the integral once, then substituting , differentiating by , putting , and completing some algebra, we arrive at the following recursive formula:

Note that substituting into the obtained formula (10), we get

and taking into account that , we find

In order to calculate , after integrating by parts, it can be transformed into the form:

from which in turn, after some algebra, we obtain

Here, the symbol denotes the integral

calculated in [39] and reported in the appendix.

Eventually, employing the relation (A.19), formula (14) reduces to the form:

from which it is easy to find

The obtained formulas (10), (16), and (A.20) allow us to evaluate , for , and eventually find

where

for and . Hence, in particular, for , we have

and, consistently, for , and , their explicit forms read

Furthermore, for ,

which for and yields

Note also that can be written in the form:

and hence the integrals are expressed by means of a linear combination of , where , and .

Since for the formulated problem, we have to find the integrals and ; below we focus on , solely. Employing (18), (A.10), and (A.19), after some algebra, we find

and hence

In particular,

Now, restricting the discussion to the case of , we can rewrite (4) in the form (cf. [14]):

where we include that, for , the gap function and we also assume that can be replaced by in this limit. The obtained equation allows one to find the superconducting gap as a function of the temperature using the method of successive approximations. So, neglecting all terms with in the sum, one gets the solution of the first boundary value problem for (28). Taking it as the first approximation and including the next terms in (28), one can obtain solutions in the subsequent approximations. For example, in the third step, we find

where it is always taken into account that is close to and should be expressed as . Moreover, using relations established in [40], we can derive the reduced specific heat difference in the form:

This allows us to find that, at , the reduced specific heat jump is equal to , whereas the slope of the tangent line to the reduced specific heat difference curve is

which corresponds to the angle of .

Note that the sum (28) is an oscillating power series, whose coefficients for tend to as . Thus, the expression (28) converges if , and then the resulting sequence converges, yielding a solution to the problem.

2.2. Integrals of the Fundamental Equations for

In developed models of superconductors, some additional parameters, such as a supercurrent or an external magnetic field, are included. In that case, one should also consider other expressions describing the superconducting current and the current of normal particles, the free energy difference, the entropy, and the magnetic susceptibility. These expressions allow one to derive a number of quantitative thermodynamic characteristics of the system such as the transition temperature, the isotope shift, the specific heat jump, the critical current, the critical superflow, the London penetration depth, the critical magnetic field, and other quantities important in characterization of superconducting systems [7, 11, 20]. Note that the integrand in the gap equation (4)—as well as in other expressions mentioned above—can be expanded in a double series, with respect to and the additional parameter . For the integrand in the gap equation, it takes the form:

Hence, there appear some new specific integrals, which can be represented as

and rewritten in the equivalent form

However, in order to evaluate such integrals, we yet define their extended forms:

so that . Now, use the following relations for even-order derivatives of , that is,

where

for and , and

for and . Note that , and

and including (9), we eventually find that

and hence putting , we obtain

where the integrals are expressed by means of formula (18), and, for example,

for and . Note also that

and such averages vanish for odd exponents.

3. Conclusions

In the present paper, we have found and discussed some infinite families of specific integrals inherent in models of superconducting systems. The results presented here can be easily applied in analytical studies on superconductors of the new generation, since the gap equation and the carrier concentration equation along with the other particular equations mentioned in the paper are important tools in research of thermodynamic properties of superconductors. They are commonly used to derive fundamental quantitative characteristics for various classes of superconducting systems.

In order to apply the results presented in the paper to superconductors of the new generation studied within the conformal transformation method [7], it is enough to take into account a model form of the function —which so far has been treated as constant—in a local interval around , where the integrands in the formulas under discussion are strongly concentrated. The model function should possess a narrow fluctuation (peak) in the vicinity of the Fermi surface. Since various particular shapes of the fluctuation produce similar results, we propose to take the function in the following normalized form:

where we also include (12). The fixed parameter is the height of the local fluctuation, whereas denotes its half-width. These should be estimated in accordance with a given form of the dispersion relation (cf. [7, 38]). Subsequently, the form of fluctuation can be substituted into all integrals and . Because of the function there appear some extra integrals and in the equations under consideration. In such a way, the proposed form of allows us to employ the derived integrals to calculate thermodynamic parameters of composed superconductors [3, 4, 6, 8–10, 12, 14, 21].

Let us also emphasize that although the case corresponds to the Van Hove Scenario, the case can be properly discussed in the conformal transformation method only [7, 21–23, 30]. One should then include additional functions , for , where , which appear after averaging the dimensionless kernel of the density of states with over the angles and can be modeled similarly to .

Appendix

The Integral

As we have shown in [39], the formula for the family of integrals

for , can be derived by analytical calculations. In order to calculate , we introduce the integrals:

where

and is a complex number such that . Then, for each , the function is a continuous function of and , where ranges over the region and , for all . Moreover, for every value of , the function is an analytic function of , and finally the integral , for all , is almost uniformly convergent in the domain . Therefore, is an analytic function of and one can differentiate under the integral sign [41]:

In particular, it yields

Let us first derive a formula for . Using a series expansion of , for , we obtain

which, after changing the summation index to and using the formula

transforms into

where the summation can be extended to running from to since for we obtain zero terms. If is sufficiently large, we can integrate the series term-by-term to get

where the change of variables in the integrals has been performed.

The precise argument which allows term-by-term integration goes as follows. First, we note that term-by-term integration of the series (A.8) is justified for the integration from to , where since the series (A.8) is uniformly convergent for and any . This allows us to integrate (A.8) term-by-term also over the interval provided the right-hand side of (A.9) converges absolutely [41]. It is the case in the half-plane , and hence in that half-plane the formula (A.9) is valid.

In order to find the coefficients standing by the same power of in (A.9), let us introduce

for and . Thus, in particular,

For example,

Now using (A.10) and the formula

for , where is the Riemann -function, we arrive at

By analytic continuation, this formula holds for all and in fact yields the analytic continuation of even to the left of the half-plane with poles at .

Recall that , for , , and . Differentiating (A.14) and then putting , we obtain

Next, using the Riemann functional equation [31, 42]

we arrive at the formula

Now, using formula (A.17) and formulas

where is the Euler constant (), we finally obtain

with . Consistently, we can find in particular cases, and for we have

in agreement with other results [3–5, 12, 14, 15, 31, 39, 43, 44]. Moreover, for subsequent ’s, we obtain

where

Moreover,

and , for , which can be exploited in numerical estimations.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This study is supported by MNiSW in 2015/2016.

References

Y. Kamihara, H. Hiramatsu, M. Hirano et al., “Iron-based layered superconductor: LaOFeP,” Journal of the American Chemical Society, vol. 128, no. 31, pp. 10012–10013, 2006.

Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, “Iron-based layered superconductor La[O1−xFx]FeAs (x = 0.05−0.12) with T c = 26 K,” Journal of the American Chemical Society, vol. 130, no. 11, pp. 3296–3297, 2008.

A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill, 1971.

A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski, Methods of Quantum Field Theory in Statistical Physics, Fizmatgiz, Moscow, Russia, 1962.

A. V. Svidzynskii, Space-Inhomogeneous Problems of the Theory of Superconductivity, chapter 1, Nauka, Moscow, Russia, 1982.

M. Cyrot and D. Pavuna, Introduction to Superconductivity and High-T c Materials, chapter 7.1, World Scientific, Singapore, 1992.

M. Krzyzosiak, R. Gonczarek, A. Gonczarek, and L. Jacak, “Conformal transformation method in studies of high- superconductors—beyond the van Hove Scenario,” in Superconductivity and Superconducting Wires, D. Matteri and L. Futino, Eds., chapter 5, Nova Science Publishers, Hauppauge, NY, USA, 2010.

R. Szczęśniak and D. Szczęśniak, “Superconductivity in the presence of a pseudogap induced by local charge density waves,” Chemistry of Metals and Alloys, vol. 4, pp. 50–57, 2011.

R. Szczęśniak and A. P. Durajski, “The characterization of high-pressure superconducting state in Si₂H₆ compound: the strong-coupling description,” Journal of Physics and Chemistry of Solids, vol. 74, pp. 641–646, 2013.

R. Szczęśniak and A. P. Durajski, “On the critical temperature and the energy gap in dense SiH₄(H₂)₂ at 250 GPa,” Solid State Communications, vol. 153, no. 1, pp. 26–30, 2013.

R. Gonczarek and M. Krzyzosiak, “Model of superconductivity in the singular fermi liquid,” in Progress in Superconductivity Research, O. A. Chang, Ed., chapter 6, Nova Science Publishers, Hauppauge, NY, USA, 2008.

K. Maki, “Effect of Pauli Paramagnetism on Magnetic Properties of High-Field Superconductors,” Physical Review, vol. 148, no. 1, pp. 362–369, 1966.

A. A. Abrikosov, “On the magnetic properties of superconductors of the second group,” Zhurnal Eksperimentalnoi i Teoreticheskoi Fiziki, vol. 32, p. 1442, 1957 (Russian), Reprinted in: Journal of Experimental and Theoretical Physics, vol. 5, p. 1174, 1957

A. A. Abrikosov, Fundamentals of the Theory of Metals, Nauka, Moscow, Russia, 1962.

D. Kasinathan, K.-W. Lee, and W. E. Pickett, “On heavy carbon doping of MgB₂,” Physica C: Superconductivity and Its Applications, vol. 424, no. 3-4, pp. 116–124, 2005.

J. Kortus, O. V. Dolgov, R. K. Kremer, and A. A. Golubov, “Band filling and interband scattering effects in MgB₂: carbon versus aluminum doping,” Physical Review Letters, vol. 94, Article ID 027002, 2005.

W. S. Agrestini, C. Metallo, M. Filippi et al., “Substitution of Sc for Mg in MgB₂: effects on transition temperature and Kohn anomaly,” Physical Review B, vol. 70, no. 13, Article ID 134514, 2004.

H. Mori, T. Okano, M. Kamiya et al., “Bandwidth and band filling control in organic conductors,” Physica C: Superconductivity, vol. 357-360, no. 2, pp. 103–107, 2001.

R. Gonczarek, M. Krzyzosiak, and A. Gonczarek, “Islands of stability of the d-wave order parameter in s-wave anisotropic superconductors,” European Physical Journal B, vol. 61, no. 3, pp. 299–308, 2008.

R. Gonczarek and M. Krzyzosiak, “Conformal transformation method and symmetry aspects of the group in a model of high- superconductors with anisotropic gap,” Physica C, vol. 426, no. 1, pp. 278–282, 2005.

R. Gonczarek, L. Jacak, M. Krzyzosiak, and A. Gonczarek, “Competition mechanism between singlet and triplet superconductivity in the tight-binding model with anisotropic attractive potential,” European Physical Journal B, vol. 49, no. 2, pp. 171–186, 2006.

R. Gonczarek, M. Krzyzosiak, L. Jacak, and A. Gonczarek, “Coexistence of spin-singlet s- and d-wave and spin-triplet p-wave order parameters in anisotropic superconductors,” Physica Status Solidi B: Basic Research, vol. 244, no. 10, pp. 3559–3571, 2007.

R. Gonczarek, M. Gladysiewicz, and M. Mulak, “On posibble formalism of anisoropic Fermi liqued and BCS-TYPE superconductivity,” International Journal of Modern Physics B, vol. 15, no. 5, p. 491, 2001.

F. C. Zhang and T. M. Rice, “Effective Hamiltonian for the superconducting Cu oxides,” Physical Review B, vol. 37, no. 7, pp. 3759–3761, 1988.

D. Y. Xing, M. Liu, Y.-G. Wang, and J. Dong, “Analytic approach to the antiferromagnetic van Hove singularity model for high-T c superconductors,” Physical Review B, vol. 60, no. 13, pp. 9775–9781, 1999.

R. Szczęśniak, “Pairing mechanism for the high- superconductivity: symmetries and thermodynamic properties,” PLoS ONE, vol. 7, no. 4, Article ID e31873, 2012.

E. Pavarini, I. Dasgupta, T. Saha-Dasgupta, O. Jepsen, and O. K. Andersen, “Band-structure trend in hole-doped cuprates and correlation with ,” Physical Review Letters, vol. 87, no. 4, Article ID 047003, 2001.

O. K. Andersen, A. I. Liechtenstein, O. Jepsen, and F. Paulsen, “LDA energy bands, low-energy hamiltonians, , and ,” Journal of Physics and Chemistry of Solids, vol. 56, pp. 1573–1591, 1995.

O. K. Andersen, S. Y. Savrasov, O. Jepsen, and A. I. Liechtenstein, “Out-of-plane instability and electron-phonon contribution to s- and d-wave pairing in high-temperature superconductors; LDA linear-response calculation for doped CaCuO₂ and a generic tight-binding model,” Journal of Low Temperature Physics, vol. 105, no. 3-4, pp. 285–304, 1996.

R. Gonczarek and M. Krzyzosiak, “Some universal relations between the gap and thermodynamic functions plausible for various models of superconductors,” Physica Status Solidi B: Basic Research, vol. 238, no. 1, pp. 29–39, 2003.

I. S. Gradstein and I. M. Ryzhik, Tables of Integrals, Sums, Series and Products, chapter 3.73, Nauka, Moscow, Russia, 1987.

R. Gonczarek, “Influence of the Fermi liquid interaction on the behaviour of superconductors in the strong magnetic fields,” Acta Physica Polonica A, vol. 71, p. 59, 1987.

R. Gonczarek, M. Krzyzosiak, A. Gonczarek, and L. Jacak, “New classes of integrals inherent in the mathematical structure of extended equations describing superconducting systems,” International Journal of Modern Physics B, vol. 29, no. 17, Article ID 1550117, 12 pages, 2015.

R. Gonczarek, “Critical current and depairing superflow in the two-dimensional superconducting Fermi system,” Acta Physica Polonica A, vol. 75, p. 675, 1988.

C. M. Varma, Z. Nussinov, and W. van Saarloos, “Singular or non-Fermi liquids,” Physics Reports: A Review Section of Physics Letters, vol. 361, no. 5-6, pp. 267–417, 2002.

R. Gonczarek, M. Gladysiewicz, and M. Mulak, “Equilibrium states and thermodynamical properties of D-wave paired HTSC in the Tight-Binding model,” Physica Status Solidi (B), vol. 233, pp. 351–363, 2002.

R. Gonczarek, M. Krzyzosiak, and M. Mulak, “Valuation of characteristic ratios for high-T c superconductors with anisotropic gap in the conformal transformation method,” Journal of Physics A. Mathematical and General, vol. 37, no. 18, pp. 4899–4911, 2004.

R. Gonczarek and M. Mulak, “Enhancement of critical temperature of superconductors implied by the local fluctuation of EDOS,” Physics Letters, Section A: General, Atomic and Solid State Physics, vol. 251, no. 4, pp. 262–268, 1999.

R. Gonczarek and J. Czerwonko, “On a class of integrals appearing in subcritical thermodynamics of superconductors,” Bulletin de l'Académie Polonaise des Sciences. Série des Sciences Mathématiques, vol. 48, p. 403, 2000.

M. Mulak and R. Gonczarek, “Discontinuous phase transitions in S-paired fermi systems,” Acta Physica Polonica A, vol. 92, no. 6, pp. 1177–1189, 1997.

E. C. Titchmarsh, The Theory of Functions, Oxford University Press, Oxford, UK, 1968.

H. Bateman and A. Erdelyi, Eds., Higher Transcendental Functions, vol. 1, McGraw-Hill, New York, NY, USA, 1953.

R. Gonczarek, “On Abrikosov's and Maki's theories of type-II superconductors in the paramagnetic green function approach,” Acta Physica Polonica A, vol. 73, no. 5, pp. 683–704, 1988.

J. Bardeen, L. N. Cooper, and J. R. Schrieffer, “Microscopic theory of superconductivity,” Physical Review, vol. 106, no. 1, pp. 162–164, 1957.