Researcher profile

Andrey V. Chubukov

Andrey V. Chubukov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2026arXiv

Spatially resolved collective modes in d-wave superconductors

We analyze the dispersion of collective modes in a superconductor with $d-$wave symmetry of the order parameter in the presence of long-range Coulomb interaction. We use diagrammatic technique and quasiclassical theory in Keldysh-Nambu formalism to compute longitudinal and transverse pair susceptibilities and extract from them the dispersion of the longitudinal and transverse collective mode. We show that at T=0, the dispersion of the transverse (plasma) mode is the same as in an s-wave superconductor, but at a finite temperature it is softer and has a much larger decay rate due to the partial screening of the Coulomb potential by nodal quasiparticles. We show that the dispersion of the longitudinal mode depends on the direction of momentum with respect to the positions of the nodes of the d-wave gap, while the decay rate of this mode does not depend on momentum. We discuss experimental implications of our results.

preprint2022arXiv

Monte Carlo study of the pseudogap and superconductivity emerging from quantum magnetic fluctuations

The origin of the pseudogap behavior, found in many high-$T_c$ superconductors, remains one of the greatest puzzles in condensed matter physics. One possible mechanism is fermionic incoherence, which near a quantum critical point allows pair formation but suppresses superconductivity. Employing quantum Monte Carlo simulations of a model of itinerant fermions coupled to ferromagnetic spin fluctuations, represented by a quantum rotor, we report numerical evidence of pseudogap behavior, emerging from pairing fluctuations in a quantum-critical non-Fermi liquid. Specifically, we observe enhanced pairing fluctuations and a partial gap opening in the fermionic spectrum. However, the system remains non-superconducting until reaching a much lower temperature. In the pseudogap regime the system displays a "gap-filling" rather than "gap-closing" behavior, consistent with experimental observations. Our results provide the first unambiguous lattice model realization of a pseudogap state in a strongly correlated system, driven by superconducting fluctuations.

preprint2022arXiv

Quantum phase transition in a clean superconductor with repulsive dynamical interaction

We consider a model of electrons at zero temperature, with a repulsive interaction which is a function of the energy transfer. Such an interaction can arise from the combination of electron-electron repulsion at high energies and the weaker electron-phonon attraction at low energies. As shown in previous works, superconductivity can develop despite the overall repulsion due to the energy dependence of the interaction, but the gap $Δ(ω)$ must change sign at some (imaginary) frequency $ω_0$ to counteract the repulsion. However, when the constant repulsive part of the interaction is increased, a quantum phase transition towards the normal state occurs. We show that, as the phase transition is approached, $Δ$ and $ω_0$ must vanish in a correlated way such that $1/|\log[Δ(0)]| \sim ω_0^2$. We discuss the behavior of phase fluctuations near this transition and show that the correlation between $Δ(0)$ and $ω_0$ locks the phase stiffness to a non-zero value.

preprint2022arXiv

Quasiparticle scattering in a superconductor near a nematic critical point: resonance mode and multiple attractive channels

We analyze the scattering rate for 2D fermions interacting via soft nematic fluctuations. The ground state is an s-wave superconductor, but other pairing channels are almost equally attractive. This strongly alters the scattering rate: At energies beyond the pairing gap $Δ$, it is renormalized by contributions from all pairing channels. At energies of order $Δ$, it is determined by the competition between scattering into a gapped continuum and dispersing nematic resonance. The outcome is a "peak-peak-dip-hump" spectrum, similar, but not identical, to the "peak-dip-hump" structure in the cuprates.

preprint2022arXiv

The dynamical exponent of a quantum critical itinerant ferromagnet: a Monte Carlo study

We consider the effect of the coupling between 2D quantum rotors near an XY ferromagnetic quantum critical point and spins of itinerant fermions. We analyze how this coupling affects the dynamics of rotors and the self-energy of fermions.A common belief is that near a $q=0$ ferromagnetic transition, fermions induce an $Ω/q$ Landau damping of rotors (i.e., the dynamical critical exponent is $z=3$) and Landau overdamped rotors give rise to non-Fermi liquid fermionic self-energy $Σ\propto ω^{2/3}$. This behavior has been confirmed in previous quantum Monte Carlo (QMC) studies.Here we show that for the XY case the behavior is different.We report the results of large scale quantum Monte Carlo simulations,which show that at small frequencies $z=2$ and $Σ\propto ω^{1/2}$. We argue that the new behavior is associated with the fact that a fermionic spin is by itself not a conserved quantity due to spin-spin coupling to rotors, and a combination of self-energy and vertex corrections replaces $1/q$ in the Landau damping by a constant. We discuss the implication of these results to experiments.

preprint2022arXiv

Twists and turns of superconductivity from a repulsive dynamical interaction

We review recent theoretical progress in understanding spatially uniform $s$-wave superconductivity which arises from a fermion-fermion interaction, which is repulsive on the Matsubara axis, where it is real, but does depend on the transferred frequency. Such a situation holds, e.g., for systems with a screened Coulomb and retarded electron-phonon interaction. We show that despite repulsion, superconductivity is possible in a certain range of system parameters. However, at $T=0$ the gap function on the Matsubara axis, $Δ(ω_m)$, must pass through zero and change sign at least once. These zeros of $Δ(ω_m)$ have a topological interpretation in terms of dynamical vortices, and their presence imposes a constraint on the variation of the phase of the gap function along the real frequency axis, which can potentially be extracted from ARPES and other measurements. We discuss how superconductivity vanishes when the repulsion becomes too strong, and obtain a critical line which terminates at $T=0$ at a quantum-critical point for superconductivity. We show that the behavior of the gap function near this point is highly non-trivial. In particular, an infinitesimally small $Δ(ω_m)$ contains a singular $δ-$function piece $ω_m δ(ω_m)$. We argue that near the critical point superconductivity may be a mixed state with even-frequency and odd-frequency gap components.

preprint2021arXiv

SU(4) symmetry in twisted bilayer graphene - an itinerant perspective

We study symmetry-broken phases in twisted bilayer graphene at small filling above charge neutrality and at Van Hove filling. We argue that the Landau functionals for the particle-hole order parameters at these fillings both have an approximate SU(4) symmetry, but differ in the sign of quartic terms. We determine the order parameter manifold of the ground state and analyze its excitations. For small fillings, we find a strong 1st-order transition to an SU(3)$\otimes$U(1) manifold of orders that break spin-valley symmetry and induce a 3-1 splitting of fermionic excitations. For Van Hove filling, we find a weak 1st-order transition to an SO(4)$\otimes$U(1) manifold of orders that preserves the two-fold band degeneracy. We discuss the effect of particle-hole orders on superconductivity and compare with strong-coupling approaches.

preprint2021arXiv

The Effect of Repulsion on Superconductivity at Low Density

We examine the effect of repulsion on superconductivity in a three-dimensional system with a Bardeen-Pines-like interaction in the low-density limit, where the chemical potential $μ$ is much smaller than the phonon frequency $ω_L$. We parameterize the strength of the repulsion by a dimensionless parameter $f$, and find that the superconducting transition temperature $T_c$ approaches a nonzero value in the $μ= 0$ limit as long as $f$ is below a certain threshold $f^*$. In this limit, we find that $T_c$ goes to zero as a power of $f^*-f$, in contrast to the high density limit, where $T_c$ goes to zero exponentially quickly as $f$ approaches $f^*$. For all nonzero $f$, the gap function $Δ(ω_m)$ changes sign along the Matsubara axis, which allows the system to partially overcome the repulsion at high frequencies. We trace the position of the gap node with $f$ and show that it approaches zero frequency as $f$ approaches $f^*$. To investigate the robustness of our conclusions, we then go beyond the Bardeen-Pines model and include full dynamical screening of the interaction, finding that $T_c$ still saturates to a non-zero value at $μ= 0$ when $f < f^*$.

preprint2020arXiv

Identification of non-Fermi liquid fermionic self-energy from quantum Monte Carlo data

Quantum Monte Carlo (QMC) simulations of correlated electron systems provide unbiased information about system behavior at a quantum critical point (QCP) and can verify or disprove the existing theories of non-Fermi liquid (NFL) behavior at a QCP. However, simulations are carried out at a finite temperature, where quantum-critical features are masked by finite temperature effects. Here we present a theoretical framework within which it is possible to separate thermal and quantum effects and extract the information about NFL physics at $T=0$. We demonstrate our method for a specific example of 2D fermions near a Ising-ferromagnetic QCP. We show that one can extract from QMC data the zero-temperature form of fermionic self-energy $Σ(ω)$ even though the leading contribution to the self-energy comes from thermal effects. We find that the frequency dependence of $Σ(ω)$ agrees well with the analytic form obtained within the Eliashberg theory of dynamical quantum criticality, and obeys $ω^{2/3}$ scaling at low frequencies. Our results open up an avenue for QMC studies of quantum-critical metals.

preprint2020arXiv

Interplay between superconductivity and non-Fermi liquid at a quantum critical point in a metal. II. The $γ$-model at a finite $T$ for $0<γ<1$

In this paper we continue the analysis of the interplay between non-Fermi liquid and superconductivity for quantum-critical systems, the low-energy physics of which is described by an effective model with dynamical electron-electron interaction $V(Ω_m) \propto 1/|Ω_m|^γ$ (the $γ$ model). In paper I [A. Abanov and A. V. Chubukov, Phys Rev B. 102, 024524 (2020)] two of us analyzed the $γ$ model at $T=0$ for $0<γ<1$ and argued that there exist a discrete, infinite set of topologically distinct solutions for the superconducting gap, all with the same spatial symmetry. The gap function $Δ_n (ω_m)$ for the $n$th solution changes sign $n$ times as the function of Matsubara frequency. In this paper we analyze the linearized gap equation at a finite $T$. We show that there exist an infinite set of pairing instability temperatures, $T_{p,n}$, and the eigenfunction $Δ_n (ω_{m})$ changes sign $n$ times as a function of a Matsubara number $m$. We argue that $Δ_n (ω_{m})$ retains its functional form below $T_{p,n}$ and at $T=0$ coincides with the $n$th solution of the nonlinear gap equation. Like in paper I, we extend the model to the case when the interaction in the pairing channel has an additional factor $1/N$ compared to that in the particle-hole channel. We show that $T_{p,0}$ remains finite at large $N$ due to special properties of fermions with Matsubara frequencies $\pm πT$, but all other $T_{p,n}$ terminate at $N_{cr} = O(1)$. The gap function vanishes at $T \to 0$ for $N > N_{cr}$ and remains finite for $N < N_{cr}$. This is consistent with the $T =0$ analysis.

preprint2020arXiv

Interplay between superconductivity and non-Fermi liquid at a quantum-critical point in a metal. I: The $γ$-model and its phase diagram at $T=0$. The case $0 < γ<1$

We analyze a class of quantum-critical models, in which momentum integration and the selection of a particular pairing symmetry can be done explicitly, and the competition between non-Fermi liquid and pairing can be analyzed within an effective model with dynamical electron-electron interaction $V(Ω_m)\sim 1/|Ω_m|^γ$ (the $γ$-model). In this paper, the first in the series, we consider the case $T=0$ and $0<γ<1$. We argue that tendency to pairing is stronger, and the ground state is a superconductor. We argue, however, that superconducting state is highly non-trivial as there exists a discrete set of topologically distinct solutions for the pairing gap $Δ_n (ω_m)$ ($n = 0, 1, 2..., \infty$). All solutions have the same spatial pairing symmetry, but differ in the time domain: $Δ_n (ω_m)$ changes sign $n$ times as a function of Matsubara frequency $ω_m$. The $n =0$ solution $Δ_0 (ω_m)$ is sign-preserving and tends to a finite value at $ω_m =0$, like in BCS theory. The $n = \infty$ solution corresponds to an infinitesimally small $Δ(ω_m)$. As a proof, we obtain the exact solution of the linearized gap equation at $T=0$ on the entire frequency axis for all $0<γ<1$, and an approximate solution of the non-linear gap equation.We argue that the presence of an infinite set of solutions opens up a new channel of gap fluctuations. We extend the analysis to the case where the pairing component of the interaction has additional factor $1/N$ and show that there exists a critical $N_{cr} >1$, above which superconductivity disappears, and the ground state becomes a non-Fermi liquid.We show that all solutions develop simultaneously once $N$ gets smaller than $N_{cr}$.

preprint2020arXiv

Interplay between superconductivity and non-Fermi liquid at a quantum-critical point in a metal. IV: The $γ$ model and its phase diagram at $1<γ<2$

In this paper we continue our analysis of the interplay between the pairing and the non-Fermi liquid behavior in a metal for a set of quantum-critical (QC) systems with an effective dynamical electron-electron interaction $V(Ω_m) \propto 1/|Ω_m|^γ$ (the $γ$-model). In previous papers we studied the cases $0<γ<1$ and $γ\approx 1$. We argued that the pairing by a gapless boson is fundamentally different from BCS/Eliashberg pairing by a massive boson as for the former there exists an infinite number of topologically distinct solutions for the gap function $Δ_n (ω_m)$ at $T=0$ ($n=0,1,2...$), each with its own condensation energy $E_{c,n}$. Here we extend the analysis to larger $1< γ<2$. We argue that the discrete set of solutions survives, and the spectrum of $E_{c,n}$ gets progressively denser as $γ$ approaches $2$ and eventually becomes continuous at $γ\to 2$. This increases the strength of &#34;longitudinal&#34; gap fluctuations, which tend to reduce the actual superconducting $T_c$ and give rise to a pseudogap region of preformed pairs. We also detect two features on the real axis for $γ>1$ which become critical at $γ\to 2$. First, the density of states evolves towards a set of discrete $δ-$functions. Second, an array of dynamical vortices emerges in the upper frequency half-plane. These two features come about because on a real axis, the real part of the interaction, $V&#39; (Ω) \propto \cos(πγ/2)/|Ω|^γ$, becomes repulsive for $γ>1$, and the imaginary $V^{&#39;&#39;} (Ω) \propto \sin(πγ/2)/|Ω|^γ$, gets progressively smaller at $γ\to 2$. The features on the real axis are consistent with the development of a continuum spectrum of $E_{c,n}$ obtained using $Δ_n (ω_m)$ on the Matsubara axis. We consider the case $γ=2$ separately in the next paper.

preprint2020arXiv

Nematic superconductivity in twisted bilayer graphene

Twisted bilayer graphene displays insulating and superconducting phases caused by exceptional flattening of its lowest energy bands. Superconductivity with highest $T_c$ appears at hole and electron dopings, near half-filling for valence or conduction bands. In the hole-doped case, the data show that three-fold lattice rotation symmetry is broken in the superconducting phase, i.e., a superconductor is also a nematic. We present the mechanism for nematic superconductivity. We take as an input the fact that at relevant dopings the Fermi energy lies in the vicinity of twist-induced Van Hove singularities in the density of states and argue that the low-energy physics can be properly described by patch models with six Van Hove points for electron doping and twelve Van Hove points for hole doping. We obtain pairing interactions for the patch models in terms of parameters of the microscopic model for the flat bands, which contains both local and twist-induced non-local interactions and show that the latter gives rise to attraction in different superconducting channels. For hole-doping, we find two attractive channels, $g$ and $i$-waves, with almost equal coupling constants. We show that in the co-existence state, where both order parameters are non-zero, the three-fold lattice rotation symmetry is broken, i.e., a superconductor is also a nematic. We find two possible nematic states, one is time-reversal symmetric, the other additionally breaks time-reversal symmetry. Our scenario for nematic superconductivity is based on generic symmetry considerations, and we expect it to be applicable also to other systems with two (or more) attractive channels with similar couplings.

preprint2020arXiv

Normal state properties of quantum critical metals at finite temperature

We study the effects of finite temperature on normal state properties of a metal near a quantum critical point to an antiferromagnetic or Ising-nematic state. At $T = 0$ bosonic and fermionic self-energies are traditionally computed within Eliashberg theory and obey scaling relations with characteristic power-laws. Quantum Monte Carlo (QMC) simulations have shown strong systematic deviations from these predictions, casting doubt on the validity of the theoretical analysis. We extend Eliashberg theory to finite $T$ and argue that for the $T$ range accessible in the QMC simulations, the scaling forms for both fermionic and bosonic self energies are quite different from those at $T = 0$. We compare finite $T$ results with QMC data and find good agreement for both systems. This, we argue, resolves the key apparent contradiction between the theory and the QMC simulations.

preprint2020arXiv

Orbital transmutation and the electronic spectrum of FeSe in the nematic phase

We consider the electronic spectrum near $M=(π,π)$ in the nematic phase of FeSe ($T<T_{\rm nem}$) and make a detailed comparison with recent ARPES and STM experiments. Our main focus is the unexpected temperature dependence of the excitations at the $M$ point. These have been identified as having $xz$ and $yz$ orbital character well below $T_{\rm nem}$, but remain split at $T>T_{\rm nem}$, in apparent contradiction to the fact that in the tetragonal phase the $xz$ and $yz$ orbitals are degenerate. Here we present two scenarios which can describe the data. In both scenarios, hybridization terms present in the tetragonal phase leads to an orbital transmutation, a change in the dominant orbital character of some of the bands, between $T > T_{\rm nem}$ and $T \ll T_{\rm nem}$. The first scenario relies on the spin-orbit coupling at the $M$ point. We show that a finite spin-orbit coupling gives rise to orbital transmutation, in which one of the modes, identified as $xz$ ($yz)$ at $T \ll T_{\rm nem}$, becomes predominantly $xy$ at $T > T_{\rm nem}$ and hence does not merge with the predominantly $yz$ ($xz$) mode. The second scenario, complementary to the first, takes into consideration the fact that both ARPES and STM are surface probes. In the bulk, a direct hybridization between the $xz$ and $yz$ orbitals is not allowed at the $M$ point, however, it is permitted on the surface. In the presence of a direct $xz/yz$ hybridization, the orbital character of the $xz/yz$ modes changes from pure $xz$ and pure $yz$ at $T \ll T_{\rm nem}$ to $xz \pm yz$ at $T > T_{\rm nem}$, i.e., the two modes again have mono-orbital character at low $T$, but do not merge at $T_{\rm nem}$. We discuss how these scenarios can be distinguished in polarized ARPES experiments.

preprint2020arXiv

Quantum Phase Transition in the Yukawa-SYK Model

We study the quantum phase transition upon variation of the fermionic density $ν$ in a solvable model with random Yukawa interactions between $N$ bosons and $M$ fermions, dubbed the Yukawa-SYK model. We show that there are two distinct phases in the model: an incompressible state with gapped excitations and an exotic quantum-critical, non-Fermi liquid state with exponents varying with $ν$. We show analytically and numerically that the quantum phase transition between these two states is first-order, as for some range of $ν$ the NFL state has a negative compressibility. In the limit $N/M\to \infty$ the first-order transition gets weaker and asymptotically becomes second-order, with an exotic quantum-critical behavior. We show that fermions and bosons display highly unconventional spectral behavior in the transition region.

preprint2020arXiv

The pairing glue in cuprate superconductors from the self-energy revealed via machine learning

Recently, machine learning was applied to extract both the normal and the anomalous components of the self-energy from photoemission data at the antinodal points in Bi-based cuprate high-temperature superconductors [Y. Yamaji {\it et al.}, arXiv:1903.08060]. It was argued that both components do show prominent peaks near 50 meV, which hold information about the pairing glue, but the peaks are hidden in the actual data, which measure only the total self-energy. We analyze the self-energy within an effective fermion-boson theory. We show that soft thermal fluctuations give rise to peaks in both components of the self-energy at a frequency comparable to superconducting gap, while they cancel in the total self-energy; all irrespective of the nature of the pairing boson. However, in the quantum limit $T \to 0$ prominent peaks survive only for a very restricted subclass of pairing interactions. We argue that the way to potentially nail down the pairing boson is to determine the thermal evolution of the peaks.

preprint2020arXiv

Thermodynamic signatures of an antiferromagnetic quantum critical point inside a superconducting dome

Recent experiments in unconventional superconductors, and in particular iron-based materials, have reported evidence of an antiferromagnetic quantum critical point (AFM-QCP) emerging inside the superconducting dome of the phase diagram. Fluctuations associated with such an AFM-QCP are expected to promote unusual temperature dependencies of thermodynamic quantities. Here, we compute the $T$ dependence of the specific heat $C(T)$ deep inside a fully gapped $s^{+-}$ superconducting state as the AFM-QCP is approached. We find that, at the AFM-QCP, the specific heat $C(T)$ vanishes quadratically with temperature, as opposed to the typical exponential suppression seen in fully-gapped BCS superconductors. This robust result is due to a non-analytic contribution to the free-energy arising from the general form of the bosonic (AFM) propagator in the SC state. Away from the AFM-QCP, as temperature is lowered, $C(T)$ shows a crossover from a $T^2$ behavior to an exponential behavior, with the crossover temperature scale set by the value of the superconducting gap and the distance to the QCP. We argue that these features in the specific heat can be used to unambiguously determine the existence of AFM-QCPs inside the superconducting domes of iron-based and other fully gapped unconventional superconductors.

preprint2020arXiv

Valley magnetism, nematicity, and density wave orders in twisted bilayer graphene

We analyze density-wave and Pomeranchuk orders in twisted bilayer graphene. This compliments our earlier analysis of the pairing instabilities. We assume that near half-filling of either conduction or valence band, the Fermi level is close to Van Hove points, where the density of states diverges, and study potential instabilities in the particle-hole channel within a patch model with two valley degrees of freedom. The hexagonal symmetry of twisted bilayer graphene allows for either six or twelve Van Hove points. We consider both cases and find the same two leading candidates for particle-hole order. One is an SU(2)-breaking spin state with ferromagnetism within a valley. A subleading inter-valley hopping induces antiferromagnetism between the valleys. The same state has also been obtained in strong coupling approaches, indicating that this order is robust. The other is a mixed state with $120^\circ$ complex spin order and orthogonal complex charge order. In addition, we find a weaker, but still attractive interaction in nematic channels, and discuss the type of a nematic order.

preprint2019arXiv

Raman response in the nematic phase of FeSe

Raman experiments on bulk FeSe revealed that the low-frequency part of $B_{1g}$ Raman response $R_{B_{1g}}$, which probes nematic fluctuations, rapidly decreases below the nematic transition at $T_n \sim 85$K. Such behavior is usually associated with the gap opening and at a first glance is inconsistent with the fact that FeSe remains a metal below $T_n$, with sizable hole and electron pockets. We argue that the drop of $R_{B_{1g}}$ in a nematic metal comes about because the nematic order drastically changes the orbital content of the pockets and makes them nearly mono-orbital. In this situation $B_{1g}$ Raman response gets reduced by the same vertex corrections that enforce charge conservation. The reduction holds at low frequencies and gives rise to gap-like behavior of $R_{B_{1g}}$, in full agreement with the experimental data.

preprint2019arXiv

The interplay between superconductivity and non-Fermi liquid at a quantum-critical point in a metal

Near a quantum-critical point, a metal reveals two competing tendencies: destruction of fermionic coherence and attraction in one or more pairing channels. We analyze the competition within Eliashberg theory for a class of quantum-critical models with an effective dynamical electron-electron interaction $V(Ω_m) \propto 1/|Ω_m|^γ$ (the $γ$-model) for $0 < γ<1$. We argue that the two tendencies are comparable in strength, yet the one towards pairing is stronger, and the ground state is a superconductor. We show, however, that there exist two distinct regimes of system behavior below the onset temperature of the pairing $T_p$. In the range $T_{cr} < T < T_p$ fermions remain incoherent and the density of states $N(ω)$ displays &#34;gap filling&#34; behavior in which the position of the maximum in $N(ω)$ is set by temperature rather than the pairing gap. At lower $T < T_{cr}$, fermions acquire coherence, and $N(ω)$ displays a conventional &#34;gap closing&#34; behavior. We argue that the existence of the two regimes comes about because of special behavior of fermions with frequencies $ω= \pm πT$ along the Matsubara axis. Specifically, for these fermions, the component of the self-energy, which competes with the pairing, vanishes in the normal state. We further argue that the crossover at $T \sim T_{cr}$ comes about because Eliashberg equations allow an infinite number of topologically distinct solutions for the onset temperature of the pairing within the same gap symmetry. Only one solution, with the highest $T_p$, actually emerges, but other solutions are generated and modify the form of the gap function at $T \leq T_{cr}$. Finally, we argue that the actual $T_c$ is comparable to $T_{cr}$, while at $T_{cr} < T < T_{p}$ phase fluctuations destroy superconducting long-range order, and the system displays pseudogap behavior.