critical temperature @xmath0 and the second critical magnetic field @xmath1 are some of the key characteristics of a superconductor .
these parameters depend on an effective coupling constant with some collective excitations @xmath2 ( @xmath3 is a density of states at fermi level , @xmath4 is an interaction constant ) , on frequency of the collective excitations @xmath5 and on ginzburg - landau ( gl ) parameter @xmath6 ( @xmath7 )
. the larger coupling constant , the larger the critical temperature . for large
@xmath8 we have @xmath9 @xcite ( or @xmath10 in bcs theory ) .
formally the critical temperature can be made arbitrarily large with increasing of the electron - phonon coupling constant . however in order to reach room temperature such values of the coupling constant are necessary which are not possible in real materials .
on the other hand we can increase the frequency @xmath5 due nonphonon pairing mechanisms as proposed in @xcite .
however with increasing of the frequency the coupling constant decreases as @xmath11 . the second critical magnetic field can be enlarged in a `` dirty '' limit , where the mean free path @xmath12 is much less than a cohernce length @xmath13 @xcite . in this case
the gl parameter can be @xmath14 , but still the critical field is low near the critical temperature : @xmath15 .
many different types of superconducting materials with a wide variety of electron pairing mechanisms exist , however all they have the critical temperature limited by values @xmath16 , despite the fact that the highly exotic mechanisms have been proposed .
the coupling constant @xmath8 determines a value of the energy gap at zero temperature @xmath17 , in the same time a ratio between the gap and the critical temperature exists : @xmath18 for presently known materials . thus increasing the zero temperature gap @xmath19 we enlarge @xmath0 .
however the critical temperature can be enlarged by violation the ratio ( [ 0.1 ] ) as @xmath20 due some addition influence on condensate of cooper pairs as shown in fig.([fig1 ] ) .
the simplest examples of such systems are systems where the proximity effect takes place : two superconductors are placed in contact @xcite , or in multi - band superconductors where the interband mixing of order parameters from different bands occurs @xcite . in this cases the cooper pairs are injected into the superconductor ( band ) from another superconductor ( band ) with higher critical temperature ( stronger interaction ) which plays a role of a source of cooper pairs . as a result in a band with lower critical temperature the ratio ( [ 0.1 ] ) can be violated as @xmath21 @xcite as illustrated in fig.([fig1 ] ) .
however the source is function of temperature with own critical temperature . in a case of the contact of two metal the critical temperature of the source drops , however in a case of the multi - band superconductor the proximity effect supports superconducting state in all bands .
and critical temperature @xmath22 are as in bcs theory : @xmath23 .
the gap @xmath24 and the temperature @xmath25 are such that @xmath26 as it can be in two - band superconductors @xcite .
the gap @xmath27 asymptotically aspires to zero as temperature rises , so that @xmath28 .
thus the critical temperature @xmath29 is equal to infinity in this hypothetical situation.,width=302 ]    the hypothetical system where @xmath30 is first proposed in @xcite . in this work
the pairing of massless dirac electrons and holes located on opposite surfaces of thin film of three - dimensional topological insulator is considered .
in such a system both electron - hole coulomb attraction , which leads to pairing , and tunnelling ( or hybridization ) of the electrons between surfaces , which leads to appearance in hamiltonian the `` sources '' of cooper pairs , act jointly . at @xmath31
the gap gradually tends to some value @xmath32 .
thus the hybridization leads to the smearing of a phase transition into the paired state in close analogy with the behavior of a ferromagnetic in an external magnetic field .
this solution means the nonzero energy gap @xmath33 if the electron - hole attraction is switched off or temperature is infinite , that is not correct physically . in its turn
the systems with `` preformed '' cooper pairs have been considered in a work @xcite . in this regime ,
the electrons are paired , but they lack the phase coherence necessary for superconductivity .
the existence of preformed pairs implies the existence of a characteristic energy scale associated to a pseudogap . in a work @xcite a model of hypothetical superconductivity
has been proposed , which demonstrates the principal differences from results of bcs and gl theory due presence of the source term in the model hamiltonian .
unlike the previous works in this model the paired state is absent if the electron - electron coupling is absent @xmath34 .
this means the electron - electron coupling is the cause of the transition to superconducting state only but not the source term . in this theory
the source term determines energy of a cooper pair in some external field which can be called as the _ external pair potential _ by analogy with the terminology in @xcite . in a case of decreasing of cooper pair s energy by the potential the energy gap tends to zero asymptotically as temperature increases .
thus the ratio between the gap and the critical temperature is @xmath30 instead of a finite value in bcs theory . in this work
the free energy functional in the limit @xmath31 has been obtained .
it is shown that the energy gap tends to zero asymptotically as magnetic field increases .
thus the second critical magnetic field is equal to infinity as well as the critical temperature , that is not correct physically , because distance between vortexes can not be less than diameter of the vortex s core .    in present work
we generalize bcs theory in the sense that the invariant under @xmath35 transformation source of cooper pairs ( the external pair potential ) @xmath36 is added to bcs hamiltonian .
we demonstrate that such a hamiltonian describes a hypothetical substance , where an interaction energy between ( within ) structural elements of condensed matter ( molecules , nanoparticles , clusters , layers , wires etc . )
depends on state of cooper pairs : an additional work @xmath37 must be made against this interaction to break a pair . in this model
the potential @xmath37 essentially renormalizes the order parameter so that the ratio ( [ 0.1 ] ) changes as @xmath20 if the breaking of a cooper pair increases energy of the molecular structure ( or creation of the pair lowers the energy ) . in another regime (
the breaking of a cooper pair lowers the energy of the structure ) the renormalization of the order parameter by the potential @xmath37 suppresses a superconducting state and the first order phase transition occurs .
we obtain the free energy of such a system which generalizes landau free energy for the presence of the external pairing potential . for the case @xmath20 in the limit @xmath31
we formulate the effective gl theory which is obtained in high - temperature limit from the general expression for free energy ( unlike the work @xcite where the free energy functional was obtained modifying the ordinary gl expansion ) .
results of our effective gl theory radically differ from results of the ordinary gl theory : the coherence length decreases as temperature rises , the gl parameter and the second critical field are increasing functions of temperature .
according to bcs theory an electron - electron attraction leads to the appearance of nonzero anomalous averages @xmath38 and @xmath39 , which are the order parameter ( pair potential ) of the superconducting state .
the order parameter is determined with some self - consistency equation @xmath40 , which reflects the fact , that superconductivity is a many - particle cooperative effect . in this regime
the charge is carried by pairs of electrons ( current carriers are the pairs with charge @xmath41 ) . to break a pair with transfer of its constituents in free quasiparticle states
the energy @xmath42 is needed .
we can consider the quantity @xmath42 as a work against the effective electron - electron attraction given the fact that the quantity @xmath19 is a collective effect .
the superconducting state of a metal and the breaking of a pair are illustrated in fig.[fig2]a .    in our model
we consider a hypothetical substance , where an interaction energy between ( within ) structural elements of condensed matter ( molecules , nanoparticles , clusters , layers , wires etc . )
depends on state of cooper pairs : if the pair is broken , then energy of the molecular system is changed by quantity @xmath43 , where @xmath44 and @xmath45 are energies of the system after- and before the breaking of the pair accordingly . thus to break the cooper pair
we must make the work against the effective electron - electron attraction and must change the energy of the structural elements : @xmath46 we will call the parameter @xmath37 as the _ external pair potential _
, since it is imposed on the electron subsystem by the structural elements of a substance , unlike the pair potential @xmath19 , which is result of electron - electron interaction and determined with the self - consistency equation @xmath40 . the parameter @xmath37 can be either @xmath47 or @xmath48 and in the simplest case it is not function of the energy gap @xmath49 , @xmath50 is a trivial case corresponding to bcs theory .
moreover we suppose that @xmath37 does not depend on temperature essentially like parameters of electron - phonon interaction .
the condition @xmath51 ensures stability of the cooper pairs ( bound state of the electrons is energetically favorable ) , otherwise transformation ( [ 1.1 ] ) has no sense and such superconducting state can not exist .
if @xmath48 then the breaking of a cooper pair lowers energy of the molecular structure ( or creation of the pair raises the energy ) . in this case
the pairs become less stable .
if @xmath47 then the breaking of the pair increases the energy ( or creation of the pair lowers the energy ) . in this case
the pairs become more stable .
a possible variant of breaking of a cooper pair in these cases is illustrated in fig.[fig2]b and fig.[fig2]c .    without going into the details of interaction of the structural elements
we can write an effective hamiltonian which takes into account the effect of the structure on cooper pairs as some effective external field , like the bcs hamiltonian is an effective hamiltonian describing a system of interacting electrons independently of nature of this interaction .
the order parameter @xmath19 is a complex quantity @xmath52 , where @xmath53 is a phase , and it is the result of a many - particle self - consistent coherent effect . in the same time
the field @xmath37 is an additional parameter imposed on the electron subsystem by the structural elements .
it is easy to see that the following transformations of the order parameter correspond to the transformation ( [ 1.1 ] ) : @xmath54 really the work to break the pair is @xmath55 however the transformations ( [ 1.2 ] ) correspond to the transformation ( [ 1.1 ] ) if @xmath51 only . if @xmath56 then the work is @xmath57 . thus the model based on the transformations ( [ 1.2 ] ) can give `` parasitic '' solutions where the the condition ( [ 1.1 ] ) is not satisfied .
such solutions must be omitted .
the hamiltonian corresponding to the transformations ( [ 1.2 ] ) is @xmath58\\      & \approx&\sum_{\textbf{k},\sigma}\varepsilon(k)a_{\textbf{k},\sigma}^{+}a_{\textbf{k},\sigma }      + \sum_{\textbf{k}}\left[\delta^{+}\left(1+\frac{\upsilon}{2|\delta|}\right)a_{\textbf{k}\uparrow}a_{-\textbf{k}\downarrow }      + \delta\left(1+\frac{\upsilon}{2|\delta|}\right )      a_{-\textbf{k}\downarrow}^{+}a_{\textbf{k}\uparrow}^{+}\right]+\frac{1}{\lambda}v|\delta|^{2},\nonumber\end{aligned}\ ] ] where @xmath59 is bcs hamiltonian : kinetic energy + pairing interaction , energy @xmath60 is counted from fermy surface .
the combinations @xmath61 and @xmath62 are creation and annihilation of cooper pairs operators and we have used anomalous averages : @xmath63 which are the complex order parameter @xmath64 . due the multipliers
@xmath65 and @xmath66 in @xmath36 the energy does not depend on the phase @xmath53 ( @xmath67 ) .
thus both @xmath59 and @xmath36 are invariant under the @xmath35 transformation .
the term @xmath36 is similar to `` source term '' in @xcite , where it means the injection of cooper pairs into the system . on the other hand
, @xmath36 has a form of an external field acting on a cooper pairs only , and @xmath37 is energy of a cooper pair in this field .
eigenfunctions of the operator @xmath68 may not be eigenfunctions of a particle number operator @xmath69 , that reflects the presence of cooper pairs condensate .    in our model
we make the following assumptions .
a change of the chemical potential at transition to superconducting state can be neglected , so that @xmath70 , because in the model we suppose @xmath71 , that occurs for systems with week electron - electron attraction and high density of carriers unlike the systems with the strong attraction and low particle density in the bec regime , so that @xmath72 , where the change of the chemical potential plays important role in formation of the superconducting state @xcite .
in addition we should notice that the coupling constant @xmath4 is a sum of phonon term and coulomb term : @xmath73 . in metals , as a rule , @xmath74 that corresponds to repulsive electron - electron interaction ,
however in such systems the pairing is possible as result of the second order processes , which lead to effective attraction regardless of the sign of interaction @xcite : @xmath75 , where @xmath76 is a coulomb pseudopotential .
such mechanism is not considered in our model . in the model
we suppose a stronger condition @xmath77 , which can occur in nonmetallic superconductors ( for example , in alkali - doped fullerides @xmath78 , where competition between the jahn - teller coupling and hund s coupling takes place @xcite ) .     is needed .
( b ) - the process of braking of a cooper pair if the braking causes the increase of energy of the molecular structure by @xmath37 ( the initial energy level shifts ) .
( c ) - it is the same , but if the braking causes the decrease of the energy.,width=302 ]    in order to find equilibrium value of the order parameter @xmath19 we should obtain system s energy . using the bcs wave function @xmath79 , where @xmath80 , we obtain the energy in a form : @xmath81 where @xmath19 can be supposed real in the absence of magnetic field and current .
the first term is internal kinetic energy , the second term is energy of the electron - electron interaction , the third term is energy of cooper pairs condensate in the external pair field .
the functions @xmath82 and @xmath83 have to minimize the energy , that is @xmath84 : @xmath85=0.\ ] ] there are two different solutions of this equation , which correspond two different physical situations :    1 .
we can suppose @xmath86 then we obtain the @xmath87 functions in a form @xmath88 as in bcs theory with @xmath89 from eq.([1.6 ] ) we can see that if the coupling constant is @xmath90 then the energy gap is nonzero @xmath91 ( if @xmath47 ) that corresponds to results obtained in works @xcite .
thus the operator @xmath36 describes the external source of cooper pairs injecting the pairs into the system . in this case
the operator @xmath36 can have a noninvariant form @xmath92 $ ] . in a work @xcite this operator and its analogues are used to break down the initial symmetry of a noninteracting system and to produce the desired structure ( superconducting , ferromagnetic , solid etc ) - quasiaverages bogolyubov method . in the other case
the operator can have a @xmath93 symmetrical form @xmath94 $ ] , where @xmath37 is the order parameter of another superconductor , for example , the boundary ( proximity ) effect , when a superconductor is placed in contact with a normal metal @xcite or interband mixing of two order parameters belonging to different bands in a multi - band superconductor @xcite occurs .
2 .   on the other hand , in order to obtain the functions ( [ 1.6a ] ) we can suppose @xmath95 then @xmath96 the integration domain in the self - consistency equation ( [ 1.8a ] ) must be cut off at some characteristic phonon energy @xmath97 ( @xmath98 if @xmath99 , outside @xmath90 ) , then we have @xmath100 we can see from eq.([1.8 ] ) , that if the coupling constant is @xmath101 then @xmath102 .
this means , that _ only electron - electron coupling is the cause of superconductivity but not the potential _ @xmath37 . in this case
the operator @xmath36 is an external field acting on the cooper pairs only , and @xmath37 is energy of the pair in this field .
as stated above , the potential @xmath37 is called as _ external pair potential _ since it is imposed on the electron subsystem by the structural elements of a substance ( molecular , clusters etc . if their energy depends on state cooper pairs ) , unlike the pair potential @xmath19 , which is result of electron - electron interaction and determined with the self - consistency equation . if @xmath50 we have usual self - consistency equation @xmath40 in bcs theory . in presence of the external pair potential @xmath103
the self - consistency condition has a form @xmath104 - eq([1.8 ] ) .
thus the potential @xmath37 _ renormalizes _ the order parameter @xmath19 .
+ it should be noticed that @xmath102 is not solution of the self - consistency equation ( [ 1.8 ] ) if @xmath105 ( due the term @xmath106 an uncertainty @xmath107 occurs ) unlike bcs theory ( when @xmath50 ) .
the absence of this solution is consequence of the fact that the derivative @xmath108 does not exist at @xmath102 because @xmath109 and the function @xmath49 is non - differentiable in this point ( the derivative has a jump discontinuity ) , unlike bcs theory where @xmath110 .
thus the regions , where eq.([1.8 ] ) does not have solutions , should be considered as the normal state of a superconductor .
however we can extend the definition of the normal state by the following way .
the order parameter is a complex quantity : @xmath64 .
hence the ratio @xmath65 can be zero due the strong fluctuations of the phase @xmath111 so that average value @xmath112 as demonstrated in @xcite . in this sense @xmath102
can be considered as a solution of the self - consistent equation ( [ 1.8 ] ) .     as function of the potential @xmath37 ( in kelvin ) calculated with eq.([1.8 ] ) , where the parameters @xmath113 , @xmath114 have been used .
if @xmath48 the critical value @xmath115 exists , below which the order parameter is completely suppressed .
the shaded region corresponds to the metastable state , where the energy @xmath116 in fig.([fig5 ] ) is positive .
, width=302 ]    ) as function of the energy gap @xmath19 at different value of the potential @xmath117 in kelvin .
the bold line is energy in bcs theory ( the external pair potential is absent @xmath50).,width=302 ]     shown in fig.[fig4 ] corresponding to different values of the potential @xmath37 .
a region , where @xmath118 , corresponds to metastable states ( it is shaded in fig.[fig3]).,width=302 ]    thus a system described with the hamiltonian ( [ 1.3 ] ) has the energy with two minimums ( [ 1.6 ] ) and ( [ 1.8 ] ) , which correspond to two physically different situations : injection of the cooper pairs into the system from external source - option 1 , and action of the external pair potential caused by dependence of energy of the molecular structure on state of a cooper pair - option 2 . in this work
we will consider the model with the external pair potential only . if @xmath47 then the pairing lowers the energy of the molecular structure , that supports superconductivity @xmath119 . if @xmath48 then the pairing of electrons increases the energy of the molecular structure , that suppresses superconductivity @xmath120 , moreover such a critical value of the potential @xmath121 exists , that if @xmath122 then eq.([1.8 ] ) has not any solutions .
as discussed above this region has been supposed as the normal state of a superconductor .
thus the first order phase transition ( on the parameter @xmath37 ) occurs - fig.[fig3 ] . from the figure
we can see @xmath123 for all @xmath124 , hence the pairs are stable . in the same time eq.([1.8 ] ) has solutions at @xmath125 , however it can shown that in this case the condition ( [ 1.1 ] ) is not satisfied , that is the `` parasitic '' solution occurs in this case .
superconducting state is energetically favorable if the energy ( [ 1.5 ] ) of the superconducting state is less than energy of normal state i.e @xmath126 .
above we have extended the definition of the normal state as @xmath102 due strong fluctuations of the phase @xmath111 so that @xmath127
. then the energy @xmath128 can be determined as @xmath129 . using eq.([1.5 ] ) we have : @xmath130\nonumber\\ & -&v\nu_{f}\left[\int_{-\hbar\omega}^{\hbar\omega}2\varepsilon \textsl{v}_{0}^{2}(\varepsilon)d\varepsilon - g\left(\int_{-\hbar\omega}^{\hbar\omega}\textsl{v}_{0}(\varepsilon)u_{0}(\varepsilon)d\varepsilon \right)^{2}-\upsilon\int_{-\hbar\omega}^{\hbar\omega}\textsl{v}_{0}(\varepsilon)u_{0}(\varepsilon)d\varepsilon\right],\end{aligned}\ ] ] where @xmath131 this renormalization of the normal state is consequence of the extrapolation in the point @xmath102 where the self - consistent equation ( [ 1.8 ] ) does not have solutions .
the expressions for the energies @xmath132 and @xmath128 have not any sense in themselves , only their difference @xmath116 is a physical quantity .
however for a case @xmath47 interpretation for @xmath128 and @xmath133 is possible .
let the electron - electron interaction is switched off : @xmath34 .
in addition it can be shown that @xmath134 .
then charge is still carried by the pairs of electrons ( current carriers are the pairs ) because for their breaking it must be made the work @xmath47 .
this means that the quasiparticle spectrum has a gap @xmath135 . but
this state is not superconducting because the ordering @xmath136 is absent . in our opinion
such state can be interpreted as state with a pseudogap .
for the case @xmath48 such interpretation is impossible , because at @xmath102 the inequality ( [ 1.1 ] ) is not satisfied .
the energy ( [ 1.9 ] ) is shown in fig.[fig4 ] . if @xmath50 we have the energy in bcs theory : the energy has symmetric minimums and a local maximum in a point @xmath102 .
if @xmath47 then the minimums becomes deeper and are located at larger values of @xmath49 .
if @xmath48 then the minimums becomes shallower and are located at smaller values of @xmath49 . at some value of @xmath37
the superconducting state becomes metastable , because the extremum in @xmath102 becomes a local minimum with lower energy . at lager value of @xmath137
the side minimums disappear and we have only one minimum in @xmath102 , that means the superconducting state is suppressed .
the minimums of the function @xmath116 corresponding to different quantities of the potential @xmath37 ( that is values of the function in points determined with eq.([1.8 ] ) ) are shown in fig.[fig5 ] .
we can see the domain of metastable states , where @xmath118 but eq.([1.8 ] ) has solutions .
the region of the solutions @xmath138 corresponding to the metastable states is shaded in fig.[fig3 ] . the minimum of the energy ( [ 1.9 ] )
can be approximated with the expression : @xmath139,\end{aligned}\ ] ] where @xmath140 is a coefficient of the order of one and @xmath141 .
if @xmath50 then we have the energy as in bcs theory .
if @xmath48 then it can be @xmath118 at @xmath142 , that corresponds to the metastable states shown in fig.[fig3 ] and fig.[fig5 ] .
at nonzero temperatures thermal activated quasiparticles appears with some a distribution function @xmath143 , hence the probability that a pair @xmath144 can be involved in formation of superconducting state is @xmath145 @xcite . then generalizing eq.([1.5 ] )
we obtain the free energy @xmath146 in a form : @xmath147-\frac{\lambda}{v}\sum_{\textbf{k},\textbf{p } } u(k)\textsl{v}(k)u(p)\textsl{v}(p)(1 - 2f(k))(1 - 2f(p))\nonumber\\ & + & 2k_{b}t\sum_{\textbf{k}}\left[f(k)\ln f(k)+(1-f(k))\ln ( 1-f(k))\right]-\upsilon\frac{\delta}{|\delta|}\sum_{\textbf{k}}u(k)\textsl{v}(k)(1 - 2f(k)),\end{aligned}\ ] ] where the first term is internal kinetic energy , the second term is energy of the electron - electron interaction , the third term is contribution of entropy , the last term is energy of cooper pairs condensate in the external pair field . minimizing the free energy @xmath148 and @xmath149 we obtain @xmath150 , @xmath151 and @xmath152 where @xmath153 a solutions @xmath154 of the self - consistency equation ( [ 2.3 ] ) at @xmath48 are shown in fig.[fig6 ] , where we can see suppression of superconductivity @xmath155 and the first order phase transition .
the critical temperature ( which is a function of the parameter @xmath37 ) is less than the critical temperature of the superconductor without the external pair potential @xmath156 . at @xmath157 the equation ( [ 2.3 ] )
has no solutions for reasons discussed in section.[zero ] . from fig.[fig6 ]
we can see @xmath123 for all temperatures @xmath158 , that is the work to break a cooper pair ( [ 1.1 ] ) is positive , hence the pairs are stable .    ) if @xmath50 - standard result of bcs theory with the second order phase transition .
squares are solutions of eq.([2.3 ] ) if @xmath159 .
the solutions demonstrate the first order phase transition .
the shaded region corresponds to the metastable states , where the free energy @xmath160 is positive.,width=302 ]    ) if @xmath50 - standard result of bcs theory with the second order phase transition .
squares are solutions of eq.([2.3 ] ) if @xmath161 .
the solutions demonstrate , that the energy gap tends to zero asymptotically as temperature rises .
bold line is asymptotical solution ( [ 2.4 ] ) at large temperature.,width=302 ]    if @xmath47 then the pairing lowers the energy of the molecular structure , that supports superconductivity @xmath162 . in this case
a solution of eq.([2.3 ] ) is such that the order parameter is not zero at any temperature - fig.[fig7 ] . at large temperature @xmath163 , where @xmath164 , the gap is @xmath165 formally the critical temperature of such system is equal to infinity ( in reality it limited by the melting of the substance ) .
it should be noted that if @xmath34 , then for any @xmath37 a superconducting state does not exist ( @xmath102 always ) .
this means the electron - electron coupling is the cause of the transition to superconducting state only but not the external pair potential .
superconducting state is energetically favorable if the free energy ( [ 2.1 ] ) is less than the free energy of normal state i.e @xmath166 . as in section [ zero ]
we can extend the definition of the normal state as @xmath102 and sew together @xmath167 with @xmath168 in this point , namely the normal state s free energy must be such that @xmath169 as in bcs theory . then using eq.([2.1 ] ) and the method of section [ zero ] - eqs.([1.9],[1.10 ] ) , we can write the free energy in a form : @xmath170d\varepsilon - g\left(\int_{-\hbar\omega}^{\hbar\omega}\textsl{v}u(1 - 2f)d\varepsilon \right)^{2}\nonumber\\ & + & 2k_{b}t\int_{-\hbar\omega}^{\hbar\omega}\left[f\ln f+(1-f)\ln ( 1-f)\right]d\varepsilon-\upsilon\frac{\delta}{|\delta|}\int_{-\hbar\omega}^{\hbar\omega}\textsl{v}u(1 - 2f)d\varepsilon]\nonumber\\ % \nonumber\\ & -&v\nu_{f}[\int_{-\hbar\omega}^{\hbar\omega}2\varepsilon \left[u^{2}_{0}f_{0}+\textsl{v}^{2}_{0}(1-f_{0})\right]d\varepsilon - g\left(\int_{-\hbar\omega}^{\hbar\omega}\textsl{v}_{0}u_{0}(1 - 2f_{0})d\varepsilon \right)^{2}\nonumber\\ & + & 2k_{b}t\int_{-\hbar\omega}^{\hbar\omega}\left[f_{0}\ln f_{0}+(1-f_{0})\ln ( 1-f_{0})\right]d\varepsilon-\upsilon\int_{-\hbar\omega}^{\hbar\omega}\textsl{v}_{0}u_{0}(1 - 2f_{0})d\varepsilon],\end{aligned}\ ] ] where @xmath171 if @xmath50 we have the free energy in bcs theory : the energy has symmetric minimums and a local maximum in a point @xmath102 - fig.[fig9]a . at @xmath172
the free energy has only one minimum in @xmath102 - superconducting state is absent .
the free energy ( [ 2.5 ] ) at different temperatures and @xmath173 is shown in fig.[fig8 ] . with increasing of temperature
the minimums becomes shallower and are located in smaller values of @xmath49 . at some temperature
the superconducting states in the minimums becomes metastable , because the extremum in @xmath102 becomes a minimum with lower energy . in the metastable phase the free energy @xmath174 but eq.([2.3 ] )
has solutions . at higher temperature
these minimums disappear , and we have only one minimum in @xmath102 , that means the superconducting state is absent .
the region of the solutions @xmath154 corresponding to the metastable states is shaded in fig.[fig6 ] .
if @xmath47 then the minimums becomes deeper and they are located in larger values of @xmath49 than in the case @xmath50 - fig.[fig9]b .
however the main difference between these cases is that the free energy at @xmath47 has absolute minimums where @xmath175 at all temperatures .
this means the system is in superconducting state at any temperature
. the minimums of the free energy ( [ 2.5 ] ) at @xmath31 are determined with eq.([2.4 ] ) .    ) as function of the energy gap @xmath19 at @xmath173 at different temperatures @xmath176 in kelvin.,width=302 ]    ) as function of the energy gap @xmath19 at absence of the external pair potential @xmath50 at different temperatures @xmath177 in kelvin .
( b ) - the free energy ( [ 2.5 ] ) at @xmath178 at the same temperatures @xmath177.,width=566 ]    let us consider several important limit cases . let the external pair potential is absent @xmath50 and temperature is slightly less than the critical temperature @xmath179 , hence the gap is small @xmath180
. then we can expand the free energy ( [ 2.5 ] ) in powers of the order parameter : @xmath181,\ ] ] where @xmath182 and @xmath183 are some coefficient .
thus we have ginzburg - landau expansion .
the case @xmath47 in a limit @xmath31 is the most interesting .
the energy gap @xmath19 is small in this region
. then we can expand the free energy ( [ 2.5 ] ) in powers of @xmath184 and in powers of @xmath49 : @xmath185.\ ] ] from condition @xmath186 we obtain @xmath187 that coincides with eq.([2.4 ] ) .
it should be noted that in this limit the contribution of the kinetic energy ( the first term in eq.([2.1 ] ) ) is proportional to @xmath188 , in the same time the contribution of the entropy ( the third term in eq.([2.1 ] ) ) is proportional to @xmath184 .
thus main factor , which resists the pairing at @xmath31 , is the entropy , and the contribution of the kinetic energy can be omitted unlike gl expansion . substituting the energy gap ( [ 2.4 ] ) in the function ( [ 2.8 ] ) we have a gain in the free energy for the superconducting state : @xmath189 that determines the thermodynamical magnetic field @xmath190 , which tends to zero asymptotically as temperature rises and is proportional to the potential @xmath37 .
at presence of gradient of the order parameter @xmath191 ( for example , when a current @xmath192 occurs ) electrons are paired with momentums @xmath193 and @xmath194 accordingly .
the nonzero momentum of the pairs causes an additional kinetic energy of the condensate , hence the expansion of superconductor s free energy has a form @xcite : @xmath195 where @xmath196 is the condensate s density . in this case
the density is a function of the momentum , and a critical momentum @xmath197 exists such that @xmath198 . as a rule @xmath199 .
thus we have to consider an operator of the kinetic energy in a form @xmath200,\end{aligned}\ ] ] where @xmath201 using the free energy ( [ 2.5 ] ) with the kinetic energy calculated on the operator ( [ 3.2 ] ) and passing to the limit @xmath202 at @xmath47 we can obtain a free energy : @xmath203,\ ] ] where the coefficients are @xmath204 then we can obtain an equilibrium value of the gap , value of the free energy in this point and the critical momentum of the condensate accordingly : @xmath205    in the general case the vector @xmath206 can run through the whole momentum space ( when the pairing potential is spatially inhomogeneous @xmath207 ) , then we should sum the free energy ( [ 3.4 ] ) over all values of the momentum @xmath206 .
moreover if the system is in a magnetic field , then the momentum @xmath206 must be replaced by @xmath208 , where @xmath209 is a fourier component of a magnetic vector potential @xmath210 , so that @xmath211 is a magnetic induction vector .
then the free energy functional takes a form @xmath212,\ ] ] where the last term is energy of the magnetic field @xmath213 . minimizing the free energy functional ( [ 3.8 ] ) by the magnetic field
@xmath214 we obtain a current in a form @xmath215 where we can suppose @xmath216 in the local ( london ) limit . in a case of a singly connected superconductor using gauge transformations @xmath217 eq.([3.9 ] ) can be reduced to the equation @xmath218 describing meissner s effect . if we consider a superconductor with an inner cavity , then along a closed path lying within the superconductor around the cavity at a distance from the cavity s surface larger than magnetic penetration depth we have @xmath219 . integrating along the path
we obtain magnetic flux quantization @xmath220 , where @xmath221 is a magnetic flux quantum .
thus the free energy functional ( [ 3.8 ] ) describes basic properties of superconductivity .    ) at @xmath222 - line ( 1 ) , the effective ginzburg - landau free energy ( [ 3.10 ] ) at @xmath223 - line ( 2).,width=283 ]    using the transformations @xmath224 and the inverse transformation @xmath225 we can write the functional eq.([3.8 ] ) in real space , however the functional will have a complicated and inconvenient form due to terms @xmath226 and @xmath227 .
nevertheless we have seen above that the functional ( [ 3.8 ] ) describes basic properties of superconductivity , hence the functional can be replaced by an effective gl functional , which has the same symmetry , the same extremes and the same values in these extremes - fig.[fig10 ] .
the effective gl functional has a form @xmath228.\ ] ] the extremum @xmath229 , the value of the free energy @xmath230 and the critical momentum @xmath231 ( if there are only one @xmath206 in the sum and @xmath232 ) are @xmath233 comparing eq.([3.11 ] ) with eqs.([3.6],[3.7 ] ) we can see coincidence of extremes @xmath229 at @xmath234 , coincidence of the values @xmath230 and coincidence of the critical momentums @xmath231 . some different dependence @xmath235 has no basic value .
the current are @xmath236 that coincides with the current ( [ 3.9 ] ) in the local limit .
this means that the magnetic responses will be identical .    in real space
the functional ( [ 3.10 ] ) takes a form @xmath237d\textbf{r},\ ] ] for convenience we can use a dimensionless variable @xmath238 , then @xmath239d\textbf{r},\ ] ] where @xmath240 is thermodynamical magnetic field @xmath190 , @xmath241 is a coherence length , which is determined with properties of a superconductor only ( at @xmath31 it does not depend on @xmath37 ) . the coherence length in bcs ( and in gl )
theory depends on temperature as @xmath242 .
that is at @xmath243 it increases as temperature rises , at @xmath244 it diverges . at @xmath172
the coherence length has physical sense of a correlation radius of fluctuations @xmath245 which decreases as @xmath246 at large @xmath247 .
the superconducting phase at @xmath172 arises fluctuationally in a form of bubble of size @xmath248 within a normal conductor , that is accompanied by an increasing of free energy ( the first term in eq.([2.8 ] ) ) .
the external pair potential @xmath47 withholds the fluctuationally arisen superconducting bubbles with such @xmath19 that the increasing in the free energy at formation of the babble is compensated by the potential @xmath37 ( the second term in eq.([2.8 ] ) ) .
this continues until the superconducting phase does not fill the entire volume of the metal .
thus the coherence length ( [ 3.16 ] ) corresponds to the correlation radius of fluctuations in gl theory .
it should be noted , that due to the transition to superconducting state in a region of size @xmath249 a change of the free energy is @xmath250 that is much less a thermal energy @xmath251 at @xmath31 . in gl theory
the inequality ( [ 3.16a ] ) means strong fluctuations . in a case
@xmath48 the situation is opposite : the normal phase arises fluctuationally in the superconductor , the external pair potential withholds the fluctuationally arisen normal bubbles until the normal phase does not fill the entire volume of the superconductor , thus superconducting phase is suppressed - fig.([fig6 ] ) .
this is physical mechanism of the renormalization of the order parameter ( [ 1.8],[2.3 ] ) by the external pair potential .
small coherence length and strong fluctuations in this model shows similarity with superfluid he ii @xcite . apparently due the strong fluctuations
such regime of superconductivity is possible in 3d systems only .
moreover josephson effect can not be observed in such regime .
in the same time we should notice that the fluctuations at @xmath252 do not interact between themselves , since the free energy ( [ 3.8 ] ) does not contain terms @xmath253 , unlike gl functional where the interaction of order parameter s fluctuations plays principal role near the critical temperature : @xmath254 @xcite . thus fluctuations in our model are fundamentally different from fluctuations within the critical domain in the standard theory .
it should be noted that phonons with lower energies than temperature of the electron gas @xmath255 are perceived by the electrons as static impurities @xcite . as well known the static impurities do not influence on value of @xmath0 in @xmath256-wave superconductors ( the anderson s theorem ) .
however for @xmath257-wave pairing the nonmagnetic impurities destroy superconductivity like magnetic impurities , hence suppression of superconductivity by thermal phonons is observed in such systems and a change to the first order transition occurs @xcite .
thus above - described regime of superconductivity is possible in @xmath256-wave superconductors only .
varying the functional ( [ 3.14 ] ) by @xmath258 and @xmath259 we obtain gl equations @xmath260 where a magnetic field penetration depth is @xmath261 the penetration depth increases with temperature , however , unlike the ordinary gl theory , it is finite quantity at all temperatures
. then gl parameter is @xmath262 the parameter increases with temperature unlike the ordinary gl theory , where the parameter is constant .
this means that at large temperature all superconductors under the external pair potential become type ii superconductors . knowing the thermodynamical magnetic field and the magnetic penetration depth we obtain a pair - breaking current as @xcite : @xmath263 besides the critical temperature an important characteristic of a superconductor are the first @xmath264 and the second @xmath1 critical fields .
the second critical field is @xmath265 @xmath1 is proportional to temperature that differs radically from gl theory , where @xmath266 . in gl theory
the transition to normal state in magnetic field occurs when average distance between vortexes becomes the order of the coherence length @xmath248 .
however in our case the length decreases with temperature @xmath267 .
this decrease compensates the rapprochement of the vortexes .
the first critical field is @xmath268 we can see the relatively strong fall with temperature , that indicates the above - mentioned tendency of transition to type ii superconductivity at large temperature .
the field of melting of the vortex lattice @xcite is @xmath269 obviously at @xmath31 the field decreases with temperature faster than the first critical field @xmath270 as temperature increases .
thus in this regime the vortex lattice is not formed , that is result of above - mentioned strong fluctuations .
in this work we have considered a hypothetical substance , where an interaction energy between ( within ) structural elements of condensed matter ( molecules , nanoparticles , clusters , layers , wires etc . )
depends on state of cooper pairs : an additional work @xmath37 must be made against this interaction to break a pair . in this case
the invariant under @xmath35 transformation source of the pairs @xmath36 has to be added to bcs hamiltonian .
we have obtained a free energy for this model , which has two minimums .
one of them corresponds to injection of the pairs into the system from an external source .
even if the electron - electron interaction is absent the energy gap is nonzero @xmath91 ( if @xmath47 ) that corresponds to results obtained in works @xcite . in another minimum ( [ 1.8],[2.3 ] )
only the electron - electron coupling is the cause of superconductivity , but not the source @xmath36 , however _ the source essentially renormalizes the order parameter _ :
in presence of the potential @xmath103 the self - consistency condition has a form @xmath104 . in this case
we call the source as the _ external pair potential _ , since the potential is imposed on the electron subsystem by the structural elements of matter , unlike the pair potential @xmath19 , which is result of electron - electron interaction and determined with the self - consistency equation . if @xmath50 we have an usual self - consistency equation in bcs theory . if @xmath48 then the pairing of electrons increases energy of the molecular structure , that suppresses superconductivity .
the phase transition conductor - superconductor becomes the first order phase transition with lower critical temperature than at @xmath50 and a region of metastable states occurs . if @xmath47 then the pairing lowers the energy of the molecular structure , that supports superconductivity . in this case at large temperatures @xmath163 the energy gap tends to zero asymptotically as @xmath184 .
thus , formally , the critical temperature is equal to infinity , however the energy gap remains finite quantity . hence the ratio between the gap and the critical temperature
is @xmath30 instead of a finite value @xmath271 for all known materials .
possible realization of this model has been proposed in @xcite .
it should be noted that the normal state @xmath102 in this model ( when @xmath47 ) can be considered as state with a pseudogap : since the quasiparticle spectrum @xmath133 has a gap @xmath135 , then charge is carried by the pairs of electrons , but this state is not superconducting because the ordering @xmath136 is absent .    for the case @xmath47 and
@xmath31 the effective gl free energy functional has been obtained .
the functional describes a superconductor without critical temperature : the order parameter depends on temperature as @xmath184 .
the thermodynamical magnetic field @xmath190 , the first critical field @xmath264 and the pair - breaking current @xmath272 tend to zero asymptotically as temperature increases and they are functions of the external pair potential .
the magnetic penetration depth @xmath4 increases with temperature as @xmath273 and it is a finite quantity at all temperatures .
the coherence length @xmath248 decreases at large temperature as @xmath246 .
the length corresponds to a size of fluctuationally arising bubbles of superconducting phase in a normal conductor ( the correlation radius of fluctuations ) . at @xmath172
the superconducting phase arises fluctuationally in a normal conductor .
the external pair potential @xmath47 withholds the fluctuationally arisen superconducting bubbles .
this continues until the superconducting phase does not fill the entire volume of the metal . in this case
a change of the free energy due to transition to superconducting state in the region of size @xmath249 is much less a thermal energy @xmath251 , that means strong fluctuations . in a case @xmath48 the external pair potential withholds the fluctuationally arisen normal bubbles in a superconductor until the normal phase does not fill the entire volume of the superconductor , thus superconducting phase is suppressed .
this is physical mechanism of above - mentioned renormalization of the order parameter by the external pair potential .
small coherence length and strong fluctuations in this model shows similarity with superfluid he ii . apparently due the strong fluctuations such superconductivity is possible in 3d systems only . in the same time the fluctuations at @xmath252 do not interact between themselves .
thus fluctuations in our model are fundamentally different from fluctuations within the critical domain in the standard gl theory .
moreover such superconductivity is possible in @xmath256-wave superconductors only due large numbers of thermal phonons which are perceived by the electrons as static impurities .    due to
the coherence length is a decreasing function of temperature @xmath267 the gl parameter @xmath274 is an increasing function of temperature @xmath275 unlike the ordinary gl theory , where the parameter is constant .
hence at large temperature all superconductors under the external pair potential become type ii superconductors with the second critical field @xmath1 proportional to temperature . in gl theory the transition to normal state in magnetic field occurs when average distance between vortexes becomes order of the coherence length @xmath248 .
however in our model the length decreases with temperature , that compensates the rapprochement of the vortexes .
this result differs from results of a work @xcite , where the second critical field is infinity like the critical temperature .
the free energy functional ( [ 3.14 ] ) has been obtained in high - temperature limit from the general expression ( [ 2.5 ] ) .
in contrast , in @xcite the free energy functional was obtained by modifying of gl expansion .
another feature of this model is that the field of melting of the vortex lattice is @xmath276 at @xmath31 .
thus in this regime the vortex lattice is not formed , that is a result of the strong fluctuations in such a system .
thus we propose a fundamentally different approach to the problem of increasing of critical temperature .
this approach is not associated with enlargement of the coupling constant or with change of the frequency , but it allows to reformulate the problem in the sense that we change the ratio between the energy gap and the critical temperature as @xmath20 at finite value of the gap . since in this model the second critical field is an increasing function of temperature this gives opportunity to have a superconducting state at very large magnetic fields .


the interest in nonminimal theories , i.e. , theories that couple the gravitational field to other fields using cross terms containing the curvature tensor , started to rise as alternative theories of gravity a long time ago .
the nonminimal field theories are based on five classes , divided accordingly to the types of fields that couple nonminimally to gravitation . the first class of nonminimal models deals with the coupling of scalar fields with the spacetime curvature .
they arose first within the scherrer
jordan  thiry
brans - dicke theory where a scalar field couples nonminimally to the ricci scalar , and the different authors @xcite had different motivations to set up such a theory , see @xcite for a historical perspective .
in addition a scalar field conformally coupled to gravitation was suggested in @xcite , see @xcite for a review on the subject .
the second class is based on the modeling of nonminimal interactions of the electromagnetic field with curvature , usually called nonminimal einstein - maxwell model , see , e.g. , @xcite .
the third class comprises einstein  yang - mills models with @xmath0 symmetry @xcite .
the fourth class embodies einstein  yang - mills  higgs models @xcite .
and the fifth class covers models in which an axion pseudoscalar field appears , for instance coupling to the electromagnetic and gravitational fields in what may be called nonminimal einstein - maxwell - axion models @xcite .    within these nonminimal theories exact solutions for electric stars @xcite , magnetic stars @xcite ,
wormholes @xcite , electric black holes @xcite ( where star - like solutions can also be found ) , and regular magnetic black holes @xcite , with a @xmath1 wu - yang ansatz ( see @xcite for the original ansatz and work , and @xcite for a full study of its properties ) have been constructed .
various aspects of star and black hole physics were discussed in these papers in all five classes of the mentioned nonminimally coupled theories .
regular black holes are black holes without spacetime curvature singularities and have been constructed with general relativity coupled to some form of matter or within alternative and extended theories of gravity . for a sample of regular black holes
see @xcite .    in this vein , we want to further investigate nonminimal regular magnetic black holes in spacetimes with a cosmological constant , extended thus the work where no cosmological constant is considered @xcite .
we discuss new examples of such exact solutions of a @xmath2 nonminimal @xmath1-symmetric theory of a gauge field of magnetic type satisfying the wu - yang ansatz .
depending on the sign and value of the cosmological constant the solutions are asymptotically de sitter ( ds ) , asymptotically flat , or asymptotically anti - de sitter ( ads ) .
the new features appearing with the @xmath2 term in the structure of the black holes with a regular center are certainly of interest .
if the dark energy is indeed a cosmological constant term and if a fundamental theory is nonminimal then our study displays the role of the dark energy on the regularization of the internal structure of black holes .
the paper is organized as follows . in sec .
[ sec2 ] we set up the general formalism for a nonminimal einstein  yang - mills theory . in sec .
[ specific ] we present a specific choice for the nonminimal parameters and display the important reduced master equation . in sec .
[ generic ] we investigate the properties and the horizons of the exact solution .
a thorough analysis is made of the solutions , in particular of the regular magnetic black hole solutions , for positive cosmological constant in sec .
[ g0 ] , for zero cosmological constant in sec .
[ eq0 ] , and for negative cosmological constant in sec .
[ l0 ] . in sec .
[ disc ] we conclude .
we study a nonminimal einstein  yang - mills theory with @xmath1 symmetry and a wu - yang - type ansatz as proposed in @xcite .
this theory is a generalization of the nonminimal einstein - maxwell theory with @xmath3 symmetry studied in @xcite .
the basic elements of this einstein  yang - mills theory are given now .
the nonminimal einstein  yang - mills theory can be described in terms of the action functional @xmath4 here @xmath5 is the determinant of the metric tensor @xmath6 , @xmath7 is the ricci scalar , @xmath2 is the cosmological constant , and @xmath8 is the coupling constant ( we are putting the speed of light and the gravitational constant equal to one , @xmath9 , @xmath10 , so that @xmath11 ) .
latin indices without parentheses run from 0 to 3 .
group indices , @xmath12 , run from 1 to 3 , and when repeated should be summed with a kronecker delta metric .
the nonminimal susceptibility tensor @xmath13 is defined as @xmath14 where @xmath15 and @xmath16 are the ricci and riemann tensors , respectively , and @xmath17 , @xmath18 , @xmath19 are the phenomenological parameters describing the nonminimal coupling of the yang - mills field with the gravitational field .
the @xmath1 yang - mills field is described by a triplet of vector potentials @xmath20 , where the group index @xmath12 runs over three values .
the yang - mills field components @xmath21 are connected to the potentials @xmath22 by the formulas @xmath23 here @xmath24 is a covariant spacetime derivative , and the symbols @xmath25 denote the real structure constants of the gauge group @xmath1 .    the variation of the action ( [ act ] ) with respect to the yang - mills potential @xmath26 yields @xmath27 with @xmath28 the tensor @xmath29 is the non - abelian analog of the excitation tensor in electrodynamics .
this analogy allows us to consider @xmath13 as a susceptibility tensor .
the gauge covariant derivative @xmath30 of an arbitrary tensor @xmath31 is defined as follows @xmath32 the variation of the action functional eq .
( [ act ] ) with respect to the metric @xmath6 yields @xmath33 the effective stress - energy tensor @xmath34 appearing in eq .
( [ einmaster ] ) can be divided into four parts , @xmath35 the first term @xmath36 is the stress - energy tensor of the pure yang - mills field .
the definitions of the other three tensors are related to the corresponding coupling constants @xmath17 , @xmath18 , @xmath19 . so @xmath37 \left[f^{(a)}_{mn}f^{mn(a ) } \right ] \
, , \label{ti}\end{aligned}\ ] ] @xmath38 \nonumber\\ & { } -\frac{1}{2}g_{ik}\biggl[{{\hat{d}}}_{m } { { \hat{d}}}_{l}\left(f^{mn(a)}f^{l(a)}_{\ n}\right)-r_{lm}f^{mn ( a ) } f^{l(a)}_{\ n } \biggr]\nonumber \\ & { } - f^{ln(a ) } \left(r_{il}f^{(a)}_{kn } + r_{kl}f^{(a)}_{in}\right ) \nonumber\\ & { } - r^{mn}f^{(a)}_{im } f_{kn}^{(a ) } - \frac{1}{2 } { { \hat{d}}}^m{{\hat{d}}}_m \left(f^{(a)}_{in } f_{k}^ { \
n(a)}\right ) , \label{tii}\end{aligned}\ ] ] @xmath39 \ , . \label{tiii}\end{aligned}\
] ] now we consider the master equations for static spherically symmetric spacetimes and objects .
let us consider a static spherically symmetric spacetime with metric @xmath40 here @xmath41 are spacetime spherical coordinates and @xmath42 is a function depending on the radial variable @xmath43 only .
the gauge field is considered to be characterized by the wu - yang ansatz ( see , e.g. , @xcite ) .
we use a position dependent basis @xmath44 , @xmath45 , @xmath46 for the generators of the @xmath1 group , instead of the standard basis @xmath47 , @xmath48 , @xmath49 .
the generators @xmath44 , @xmath45 , @xmath46 are defined in terms of @xmath47 , @xmath48 , @xmath49 as @xmath50 and satisfy the commutation relations @xmath51={\bf t}_{\varphi } \,,\quad \left[{\bf t}_{\theta } \ , , { \bf t}_{\varphi}\right]={\bf t}_{r } \
, , \quad \left[{\bf t}_{\varphi},{\bf t}_{r}\right]={\bf t}_{\theta}\,.\label{des6}\ ] ] the wu - yang - type ansatz is of the form ( see @xcite for details ) @xmath52 the magnetic parameter @xmath53 is a nonvanishing integer .
the field strength tensor has only one nonvanishing component , @xmath54 clearly , it is a magnetic - type solution and it does not depend on the parameters @xmath2 , @xmath17 , @xmath18 and @xmath19 .
then , the nonminimal gravitational field equations reduced to spherical symmetry have the form @xmath55 , \label{g1}\end{gathered}\ ] ] @xmath56 , \label{g2}\end{gathered}\ ] ] @xmath57\ , , \label{g3}\end{gathered}\ ] ] where a prime denotes a derivative with respect to the radial variable @xmath43 .
equation  ( [ g3 ] ) can be deduced from eq .
( [ g2 ] ) upon differentiation .
clearly , eqs .  ( [ g1 ] ) and ( [ g2 ] )
coincide when @xmath58 , i.e. , @xmath59 this is our first choice for a relation between the parameters @xmath17 , @xmath18 , and @xmath19 .
we search for solutions with @xmath60 and @xmath61 .
equation  ( [ g1 ] ) is compatible with these conditions at @xmath62 , when @xmath63 , i.e. , @xmath64 , or @xmath65 this is our second choice for another relation between the parameters @xmath17 , @xmath18 , and @xmath19 .
these two requirements , eqs .
( [ q1231 ] ) and ( [ q1232 ] ) , restrict the number of the three nonminimal coupling constants , @xmath17 , @xmath18 , and @xmath19 , to just one independent coupling constant , @xmath66 , say .
we then put @xmath67 and so @xmath68 @xmath69 we assume that @xmath70 .      with these choices the remaining master equation for the gravitational field ( [ g1 ] )
can be rewritten as , @xmath71 =   - \frac{{q_{m}^2}}{r^2 } -\lambda r^2 \ , , \label{nn}\ ] ] where @xmath72 with @xmath73 being the magnetic charge .
the solution to eq .
( [ nn ] ) is of the form @xmath77 thus , we deal with a four - parameter family of exact solutions .
the first parameter is the nonminimal parameter of the theory @xmath66 , the second parameter is the cosmological constant @xmath2 , inbuilt into the theory , the third parameter , @xmath78 , is the magnetic charge of the wu - yang gauge field , and the fourth parameter , @xmath79 , is the asymptotic mass of the object , that makes its appearance as a constant of integration .
we consider @xmath80 @xmath81 @xmath82 @xmath83 in what follows .
the limiting case , @xmath84 , gives the minimal coupled theory .
the corresponding solution , taken as the limit @xmath84 of eq .
( [ n00 ] ) , is the magnetic reissner - nordstrm solution with a cosmological constant , @xmath85 this magnetic solution has curvature singularities at @xmath86 .
in addition , for @xmath87 , eq .  ( [ n00 ] ) yields spacetime curvature singularities at finite positive @xmath43 . for @xmath70
there are no curvature singularities as shown below .
that is why we choose @xmath70 .
near the center the metric function @xmath88 behaves as follows : @xmath89 thus displaying that @xmath60 , @xmath61 and @xmath90 .
this means that the point @xmath86 is a minimum of the regular function @xmath88 independently on the sign and value of the cosmological constant @xmath2 , and independently of the mass value @xmath79 . since @xmath60 ,
the curvature scalar is regular in the center : @xmath91 .
the quadratic scalar @xmath92 , and other curvature invariants are also finite in the center for @xmath70 as we are considering .
thus the spacetime is truly regular at the center , and this is due to the nonminimality of the model
. the cosmological constant does not contribute to the regularity at the center .
all the objects displayed are regular objects .
we can anticipate some features of the function @xmath88 .
@xmath88 , independently of the parameters @xmath66 , @xmath2 , @xmath78 , and @xmath79 , has a central small cavity . where does the boundary of this small cavity is situated ?
the simplest definition of this boundary can be associated with the nearest maximum of the function @xmath88 .
the position of this maximum , which gives the width of the small cavity from @xmath86 up to this maximum , with the corresponding value for @xmath42 , is predetermined by the values of all four parameters .
however , it is the nonminimal parameter @xmath66 that provides the creation of this central small cavity . for @xmath43
greater than this first maximum , there exist a zone in which the first derivative @xmath93 is negative .
thus , we deal with a repulsion zone here . the size of the repulsion zone is set by the parameter @xmath78 , roughly of the order of @xmath78 .
in addition when the mass @xmath79 is small enough and the cosmological constant is negative , @xmath94 , the maximum can be converted into an inflection point and the repulsion barrier vanishes . for a wide range of values of the parameters @xmath66 , @xmath95 , @xmath96 , and @xmath2 there exists a newtonian - type attraction zone of the metric function @xmath88 , where it behaves as @xmath97 .
when @xmath98 , this zone starts at infinity and finishes at some point .
when @xmath99 , this newtonian attraction zone appears only if the mass of the object exceeds some critical value .
finally , depending on the sign and value of the cosmological constant the solutions are asymptotically ds , asymptotically flat , or asymptotically ads .
the most important feature of our investigation is the analysis of the horizons as functions of the parameters .
cauchy , event , and cosmological horizons can appear , and in some instances they can coincide one with another or the three altogether .
the details are given below .
when @xmath100 , and the roots @xmath101 are real and positive , we deal with horizons at these radii .
since the equation @xmath100 can be reduced to an algebraic equation of order six , see eq .
( [ n00 ] ) , @xmath102 , one can be sure that the number of horizons is not more than six .
we intend to show that in the model under consideration the number of horizons can not be more than three .
moreover , all the existing solutions for the horizons @xmath101 , namely , @xmath103 , @xmath104 and @xmath105 , depend on sign and value of @xmath2 , i.e. , the cosmological constant participates in the formation of the causal structure of the spacetime .
generically the three horizons that appear are the cauchy horizon , the event horizon , and the cosmological horizon . to study the causal structure of the corresponding spacetimes we have , first , to count the horizon number , second , to classify the horizons , and third , to visualize them in figures .
thus in order to find the location of the horizons , i.e. , to find @xmath101 explicitly one should solve @xmath106 i.e. , the sixth - order equation @xmath107      another important property is the number of the horizons . instead of using eqs .
( [ revnb0])-([revnb ] ) and seeing how many horizons there are by brute force , one can use a more efficient method by introducing one auxiliary function @xmath108 in the following context .
indeed , eq .  ( [ revnb ] ) can be rewritten in the form @xmath109 @xmath110 to count the horizon number we have to determine the number of points in which the plot of the function @xmath111 is crossed by the horizontal line @xmath112 , the mass line . from a physical point of view
this procedure shows how many horizons can appear when the mass of the object is equal to a given @xmath79 .
there are only four types for the function @xmath113 .
they are depicted in fig .
[ figc ] .
we consider the cosmological constant @xmath2 to be positive , zero , and negative .
it is reasonable to start with the analysis of the case @xmath114 , i.e. , asymptotically ds , which is the most interesting and complex . then we study the case @xmath115 , i.e. , asymptotically flat , in order to recall the partially known results .
finally , we further analyze the case with @xmath94 , i.e. , asymptotically ads .
let us consider spacetimes with ds asymptotics .
the model with @xmath114 has a variety of structures since a cosmological horizon appears in addition to the usual inner and outer horizons .
all spacetimes are asymptotically ds .
we start with the analysis of the number of horizons .
we focus on eqs .
( [ revn00d2m])-([revn00d ] ) using the form @xmath116 @xmath117 where we have written @xmath118 in order to make explicit that the cosmological constant is positive .
clearly , the auxiliary function @xmath119 takes negative infinite values at @xmath120 , and positive infinite values at @xmath121 .
the extrema of @xmath119 are given through the equation @xmath122 which in turn from eq .
( [ 1revn ] ) yields the equation @xmath123 in terms of the auxiliary variable @xmath124 , eq .
( [ revn010 ] ) gives the cubic equation @xmath125 .
it has three roots @xmath126 , @xmath127 , and @xmath128 , for which the product @xmath129 is negative , and the sum @xmath130 is positive .
there are three cases : ( i ) there are three real roots and two of them are positive , ( ii ) there are two coinciding positive roots and one negative , and ( iii ) one real root is negative and other two roots are complex conjugated .
this means that three different situations are available .
the configuration can have two extrema ( one maximum and one minimum ) , or it can have two coinciding extrema , or still it can have no extrema .
the points of inflection of @xmath119 are given through the equation @xmath131 which in turn from eq .
( [ 1revn ] ) yields the bicubic equation @xmath132 in terms of the auxiliary variable @xmath133 , eq .
( [ 1revn01 ] ) gives the cubic equation @xmath134 .
it has three roots @xmath135 , @xmath136 , and @xmath137 , for which the product @xmath138 is positive , and the sum is equal to zero , @xmath139 .
this means that only two cases are available : ( i ) two real roots are negative and one real root is positive , and ( ii ) the real part of a complex conjugated pair of roots is negative and the only real root is positive . since negative and complex roots are not physical , we conclude that there exists only one inflection point of the graph of the function @xmath119 .
the sketches of the corresponding three plots of the function @xmath119 are presented in fig .
[ figc ] , panels ( a ) , ( b ) and ( c ) . the mass line @xmath140 crosses the corresponding plots at least in one point , meaning that there is at least one horizon .
more precisely , depending on the values of parameters @xmath66 , @xmath79 , @xmath78 , and @xmath141 , we deal with three independent simple horizons , or with one simple and one double horizon , or with one triple horizon , or even with one simple horizon .      the analysis of horizons and extrema shows that there are two critical masses , @xmath142 and @xmath143 , say , for which the corresponding horizons coincides with one of the extrema . for these masses
the horizon is in a fact a double horizon , it is degenerate .
one coincidence takes place when the plot of the function @xmath88 touches the line @xmath144 in its minimum , and another coincidence happens when @xmath88 touches the line @xmath144 in its maximum .
the critical masses can be calculated as @xmath145 @xmath146 where the radii @xmath147 and @xmath148 are obtained as the corresponding solutions of the pair of equations @xmath149 and @xmath150 .
thus , two critical parameters @xmath142 and @xmath143 appear in the course of analysis .
these critical masses can also coincide . in this case
we deal with a triple horizon , when the two extrema , the inflection point and all three zeros of the function @xmath88 coincide .
the family of solutions under discussion is described by four positive parameters with different intrinsic units : @xmath66 , @xmath79 , @xmath78 , and @xmath2 . in order to clarify the fine details of these solutions we use below auxiliary and effective parameters and auxiliary and effective variables .
first , we introduce four auxiliary parameters with units of length @xmath151 second , we construct three effective positive parameters without units , namely , @xmath152 third , since @xmath99 , in fact here @xmath153 , we introduce the variable @xmath154 without units through @xmath155 we thus use the cosmological scale for the description of the internal structure of the regular nonminimal black holes .    in terms of these variables
the metric @xmath88 , given in eq .
( [ n00 ] ) , is now @xmath156 and takes the form @xmath157 its first derivative is @xmath158}{(x^4+\gamma)^2 } \ , , \label{d13}\ ] ] where @xmath159
. clearly , only the three parameters without units , @xmath160 , @xmath161 , and @xmath162 , are the ones that matter in the search for horizons , @xmath163 , and in the search for extrema , @xmath164 .    in this analysis
the masses @xmath142 and @xmath143 given in eqs .
( [ d170])-([d17 ] ) are important .
they give the cases in which there is a double horizon .
this double horizon can be an extreme black hole horizon or an extreme cosmological horizon , each with a corresponding radius @xmath165 and @xmath166 .
the critical masses can coincide giving a triple horizon .
the plots of @xmath88 are presented in fig .
below we describe the features of the solutions for the different characteristic masses .
\(i ) the case @xmath167    when the mass @xmath79 vanishes , we deal with a regular nonminimal object that is not a black hole . at @xmath168
there is a minimum with @xmath60 ; then the function @xmath156 reaches a maximum at @xmath169 , where @xmath170^{\frac14 } \ , .
\label{d19}\ ] ] the zone @xmath171 is a central attraction zone .
the zone @xmath172 is a repulsion zone .
there is one inflection point of curve @xmath156 in this zone .
this spacetime has only one horizon , namely , the cosmological horizon at @xmath173 , where @xmath174 can be found as follows .
define @xmath175 and @xmath176 as @xmath177 -\frac13 \beta\ , , \label{n13u0}\ ] ] @xmath178    when @xmath179 , one obtains that @xmath180 } \,.\label{d21}\end{aligned}\ ] ]    when @xmath181 , one obtains that @xmath182    when @xmath183 , one obtains that @xmath184 this case is illustrated by curve i in fig .
[ fig5 ] .
\(ii ) the case @xmath185    inside the range @xmath186 there exists a specific mass value , @xmath187 which distinguishes the models with no newtonian - type attraction zone from the ones with a newtonian - type attraction zone , i.e. , with a @xmath97 dominating term in the gravitational potential . for masses in the range
@xmath185 there is no newtonian - type attraction zone .
the corresponding radius @xmath188 appears as a solution of the pair of equations , @xmath189 and @xmath190 ; i.e. , this point of the graph of the function @xmath156 has a coincidence of extrema points and an inflection point .
this limiting case , where the mass is @xmath191 is illustrated by curve ii in fig .
[ fig5 ] and is representative of the range @xmath192 .
the nonminimal regular objects within this mass range are not black holes , since the single simple horizon is a cosmological horizon .
\(iii ) the case @xmath193    clearly , within this range the following features can be visualized , see curve iii in fig .
there is , first , a central small cavity ( between the center @xmath86 and the nearest maximum ) , second , a repulsion barrier ( between the nearest maximum and the following minimum ) , third , the zone of newtonian - type attraction , i.e. , similar to the @xmath97 behavior of the newtonian gravitational potential ( between the minimum and the final maximum ) , and fourth , the zone of cosmological repulsion ( between the final maximum and infinity ) .
the nonminimal regular objects within this mass range are not black holes , since the single simple horizon is a cosmological horizon .
\(iv ) the case @xmath194    when the mass @xmath79 has this value , the nonminimal coupled regular black hole is extremal with the inner and the even horizon coinciding .
the minimum coincides with a double horizon .
there is still an outer cosmological horizon .
this is shown in curve iv of fig .
[ fig5 ] .
\(v ) the case @xmath195    when the mass @xmath79 is in this range , we deal with regular nonminimal black holes inside a cosmological horizon .
specific features of these solutions are the following , see the curve v in fig .
each solution of this type displays a small cavity near the regular center , a repulsion barrier , a cauchy horizon , a black hole region , an event horizon , a zone of newtonian - type attraction of @xmath97 plus corrections , an intermediate maximum , and a zone of cosmological repulsion .
\(vi ) the case @xmath196    when the mass @xmath79 has this value , the nonminimal coupled regular black hole is extremal with the even horizon coinciding with the cosmological horizon .
the maximum coincides with a double horizon .
there is still an inner cauchy horizon .
this is shown in curve vi of fig .
[ fig5 ] . in this case
the black hole is a cosmological extremal supermassive regular black hole as the whole visible universe is swallowed by this supermassive object .
\(vii ) the case @xmath197    when the mass exceeds the second critical mass , the cauchy horizon of the regular nonminimal black hole becomes the cosmological horizon .
we deal now with a supermassive black hole with central small cavity , a unified repulsion barrier , and without a standard zone of newtonian - type attraction , see curve vii on fig .
[ fig5 ] . since for @xmath198
the cauchy horizon is also a cosmological horizon , the black hole is ultramassive . in such universe
there is only one horizon , which is a cosmological one , together with a repulsion zone . in this context
this ultramassive black hole is of a new type , for which all three horizons coincide , namely , the cauchy , event and the cosmological horizons .
\(viii ) the case @xmath199    when the critical masses @xmath142 and @xmath143 coincide , two separatrices convert into one curve , depicted in panel ( b ) of fig .
[ fig5 ] , curve viii .
this curve illustrates the situation , when four specific points of the function @xmath156 coincide , namely , the maximum , the minimum , the inflection point and the horizon .
we deal now with a triple horizon .
when @xmath200 , the solutions behave like curve iii ( but without extrema ) . when @xmath201 the curves look like the curve vii ( but without extrema ) .
this case has been studied in @xcite and we briefly review the results . all spacetimes are asymptotically flat .      in this case we do not need to introduce variables without units as we did in the previous case .
a typical sketch of the auxiliary function @xmath202 looks like the curve depicted in panel ( d ) of fig .
there is a critical mass @xmath203 , corresponding to the minimum of that curve .
depending on the value of @xmath79 of the object , the line @xmath112 can cross the curve @xmath204 zero times , one time or two times .
the corresponding spacetimes have no horizons , when the mass @xmath79 is less than the critical mass , @xmath205 , it has one horizon , a double one , when @xmath206 , and it has two horizons , when @xmath207 .    in more detail , for @xmath115 , eqs .
( [ revn00d2m])-([revn00d ] ) turn into @xmath208 @xmath209 when the mass line @xmath112 touches the curve @xmath204 at the point of minimum , the mass is a critical mass , @xmath210 , given by @xmath211 where @xmath212 is the critical radius
. this radius can be found by imposing @xmath213 and @xmath214 , yielding @xmath215 equations  ( [ m10])-([d33 ] ) then give @xmath216 the parameter @xmath212 is convenient for the analyses . using it
we obtain , in particular , @xmath217 and @xmath218 .
now we give the main features of the four distinct types that appear for @xmath115 , see fig .
[ fig3 ] .
\(i ) the case @xmath167    when @xmath167 we deal with a pure nonminimal object . the metric function @xmath88 given in eq .
( [ n00 ] ) is @xmath219 and its derivative has also a simple form , @xmath220 .
the plot of the function @xmath88 , see curve i in fig .
[ fig3 ] , is located above the horizontal line @xmath221 and has a maximum at @xmath222 , and a minimum at the center .
thus , a newtonian - type attraction zone is absent now , and there is a repulsion barrier that prolongs itself into infinity .
there is no equilibrium point and there are no horizons in the spacetime .    .
the plots of the metric function @xmath88 for the nonminimal regular models with @xmath115 . note that @xmath60 for all curves so that the spacetimes are regular .
all the curves contain a small cavity near the center , a repulsion barrier , and a flat asymptotic behavior .
curve i displays the behavior of @xmath88 , when @xmath167 .
curve ii is for the case @xmath223 .
there are no horizons in the spacetime , the object is a regular star .
curve iii , for case @xmath206 , is the critical case .
there is a double horizon at the minimum of the metric function @xmath88 , i.e. , the cauchy and event horizons coincide .
curve iv describes a black hole with mass exceeding the critical mass @xmath207 , it contains a cauchy horizon and an event horizon .
curves ii , iii , and iv , which correspond to massive objects , have zones of newtonian - type attraction.,height=226 ]    \(ii ) the case @xmath223    when @xmath205 there are no horizons in the spacetime , and we deal with regular stars .
when the mass @xmath79 is in this range one can find the point of minimum of @xmath88 , @xmath224 , say , which is an equilibrium point , in which the zones of attraction and repulsion make contact , see curve ii in fig .
since at this point @xmath225 , the gravitational force does not act on any massive particle , i.e. , a particle can be at rest at this point .
any small perturbation puts this particle to oscillate near this equilibrium point .
\(iii ) the case @xmath206    in the critical case @xmath206 the horizon is double and has radius @xmath226 given in eq .
( [ d33 ] ) . the metric function @xmath88 given in eq .
( [ n00 ] ) has for this case the simple form @xmath227 \,.\label{n11}\ ] ] the multiplier @xmath228 shows explicitly that the horizon at @xmath226 is a double horizon . in other words , we deal with an extremal regular nonminimal black hole . as @xmath229
the spacetime is asymptotically flat .
the plot of the function ( [ n11 ] ) is presented in fig .
[ fig3 ] curve iii .
the curve contains a zone of newtonian - type attraction , which starts at infinity and finishes at the horizon .
the interior of this extremal black hole is regular .
it contains a zone near the center which is a nonminimal small cavity , and it also contains a nonminimal repulsion barrier which is locate between the sphere on which @xmath88 has a maximum and the double horizon .
\(iv ) the case @xmath230 when @xmath230 , the interior of the regular black hole contains a nonminimal small cavity and a repulsion barrier that are separated from the exterior by a nonminimal well , see curve iv in fig .
the edge of the well nearest to the center is the cauchy horizon , the distant edge is the event horizon , it separates the interior from the zone of newtonian - type attraction with a @xmath231 potential term .
the case @xmath232 has some features similar to the @xmath233 case , with in addition some more structure .
all spacetimes are asymptotically ads .
let us consider spacetimes with ads asymptotics , @xmath234 .
the case @xmath94 , like the @xmath235 case , is simpler than the case @xmath114 and there is no need to introduce variables and parameters with no units .
a typical sketch of the auxiliary function @xmath236 looks like the curve depicted in panel ( d ) of fig .
[ figc ] . in this case
there is also a critical mass @xmath203 related to the minimum .
depending on the value of @xmath79 of the object , the line @xmath112 can cross the curve @xmath237 zero times , one time or two times .
the corresponding spacetimes have no horizons , when the mass @xmath79 is less than the critical mass , @xmath205 , it has one horizon , a double one , when @xmath206 , and it has two horizons , when @xmath207 . since the spacetime is asymptotically ads for sufficient large radii there is a cosmological attraction .    in more detail , for @xmath234 , eqs .
( [ revn00d2m])-([revn00d ] ) turn into @xmath238 @xmath239 where we have set @xmath240 . since the second derivative of the function @xmath236 is positive everywhere , and @xmath241 , a typical sketch of this function is like panel ( d ) in fig .
this function has features similar to the function @xmath202 analyzed previously .
we can find the critical value of the mass @xmath203 .
it is @xmath242\,,\label{1ren00}\ ] ] where @xmath212 is the real positive solution of the following bicubic equation @xmath243 in this case there is only one real positive solution .
indeed , if we introduce the quantity @xmath244 as @xmath245 , we see that the product of the roots of eq .
( [ 1revn010 ] ) , @xmath246 is positive , and the sum @xmath247 is negative .
this can happen in three instances : first , when there are two complex roots and one real positive root , say , @xmath248 ; second , when there are three real roots , two of them being negative , and one positive , @xmath248 ; third , when @xmath249 , and @xmath248 . in any case
, there is only one real positive solution @xmath250 of eq .
( [ 1revn010 ] ) . technically , in each of the three mentioned above situations , the unique real positive root @xmath126 can be extracted from the three corresponding solutions of the standard cardano formula ,
namely , @xmath251 where @xmath252^{\frac13}\ , , \nonumber\\ \nu \equiv \left[-\frac{{\cal q}}{2 } - \sqrt{\frac{{\cal p}^3}{27 } + \frac{{\cal q}^2}{4 } } \right]^{\frac13}\,,\label{cardano2}\end{gathered}\ ] ] @xmath253 the number of real roots is known to be controlled by the discriminant @xmath254 .
there are three real roots , when @xmath255 , two roots coincide when the discriminant vanishes @xmath256 , and there is one real root when @xmath257 .
let us consider one interesting example exactly solvable .
we put @xmath258 , i.e. , @xmath259 then @xmath260 , and we obtain from ( [ cardano1])([cardano3 ] ) three real roots @xmath261 the root @xmath126 is positive and thus the minimum of the function @xmath236 gives @xmath262 , i.e. , latexmath:[\[r_{\rm c }   = \frac{1}{\sqrt{3|\lambda|}}\sqrt{\sqrt{3 + { 9}{q_{m}^2 }    @xmath203 in eq .
( [ 1ren00 ] ) in this particular case studied .
\(i ) the case @xmath167    when @xmath167 , the metric function is @xmath264 clearly , for @xmath167 and negative cosmological constant @xmath265 everywhere , and thus there are no horizons .
the derivative @xmath266 \label{d5}\ ] ] is equal to zero at the center @xmath86 , and so this point is a local minimum for the metric function @xmath88 .
now , there are still further features depending on whether @xmath267 , @xmath268 , or @xmath269 .
when @xmath267 , there are two additional extrema , the maximum at @xmath270 \right\}^{\frac14}\,,\label{d7}\ ] ] and the minimum at @xmath271 \right\}^{\frac14}\,.\label{d9}\ ] ] in the limiting case @xmath272 the minimum drifts to infinity , and the maximum coincides with the value @xmath273 , thus recovering the case with @xmath115 .
when @xmath268 , the two radii above coincide with an inflection point , i.e. , @xmath274 yielding @xmath275 .    when @xmath269 there are no additional extrema .
this latter case , @xmath269 , is illustrated by curve i in fig .
[ fig4 ] .
\(ii ) the case @xmath276    for @xmath277 there appears another specific mass @xmath278 , with @xmath279 denoting triple , for which @xmath280 and @xmath281 . at this radius ,
the curve of the metric function @xmath88 has two extrema coinciding with the inflection point .
for @xmath282 the objects correspond to regular stars since there are no horizons , see curve ii in fig .
we note that this curve is of the same type as the curve for the case @xmath167 and @xmath268 . in the range
@xmath276 there are no horizons in the spacetime also , and so we also deal with regular stars , of which curve ii is representative .    . the plots of the metric function @xmath88 for the nonminimal regular models with @xmath94 .
note that @xmath60 for all curves so that the spacetimes are regular .
all the curves contain a small cavity near the center , a repulsion barrier , and an ads asymptotic behavior .
curve i displays the behavior of @xmath88 when @xmath167 .
curve ii , for @xmath282 is also representative of the case @xmath283 .
there are no horizons , the objects are nonminimal stars .
curve iii , for @xmath284 also represents nonminimal stars , in addition it has a flat well with a repulsive barrier followed by a newtonian - type attraction zone .
curve iv is the critical case @xmath206 with a double horizon , it describes an extremal black hole .
curve v describes a black hole with @xmath230 , it has a cauchy horizon and an event horizon .
, height=226 ]    the plots of the metric function @xmath88 for the nonminimal regular models with @xmath115 .
note that @xmath60 for all curves so that the spacetimes are regular .
all the curves contain a small cavity near the center , a repulsion barrier , and a flat asymptotic behavior .
curve i displays the behavior of @xmath88 , when @xmath167 .
curve ii is for the case @xmath223 .
there are no horizons in the spacetime , the object is a regular star .
curve iii , for case @xmath206 , is the critical case .
there is a double horizon at the minimum of the metric function @xmath88 , i.e. , the cauchy and event horizons coincide .
curve iv describes a black hole with mass exceeding the critical mass @xmath207 , it contains a cauchy horizon and an event horizon .
curves ii , iii , and iv , which correspond to massive objects , have zones of newtonian - type attraction .
( iii ) the case @xmath284    when @xmath284 there are no horizons in the spacetime .
the plot of @xmath88 in this case has a minimum at @xmath224 , for which @xmath225 , see curve iii in fig .
since at @xmath224 one has @xmath225 , massive particles suffer no force and can be at rest there .
however , any small perturbation puts the particle to oscillate near this equilibrium position .
\(iv ) the case @xmath206    when @xmath206 , we deal with an extremal regular black hole in an asymptotically ads spacetime .
the horizon is a double horizon , see curve iv in fig .
there are no specific points up to the second order that distinguish the zone of newtonian - type attraction and the zone of cosmological attraction .
\(v ) the case @xmath230    when @xmath230 we deal with regular black holes with a cauchy horizon and an event horizon in an asymptotically ads spacetime , see curve v in fig .  [ fig4 ] .
the solution for the metric potential @xmath88 that we have found , namely @xmath285 , presents a new exact spherically symmetric static solution of a nonminimal @xmath1 einstein  yang - mills theory with a cosmological @xmath2 and a wu - yang ansatz . the expression for @xmath88 shows explicitly a four - parameter family of exact solutions .
these parameters are the nonminimal parameter @xmath66 , the cosmological constant @xmath2 , the magnetic charge @xmath78 , and the mass @xmath79 .
the most important feature of the family of exact solutions is that it has solutions with horizons depending on the relative values of the parameters .
we would like to stress the following details in this context .    for @xmath114 , depending on the values of the parameters , the black hole solution can have three horizons , the cauchy horizon , the event horizon and the cosmological horizon .
when @xmath286 , the first critical mass value , there is only one horizon , the cosmological horizon .
this means that a typical profile of the metric function @xmath88 contains a central small cavity , a repulsion barrier , a zone of rest for matter near the point of minimum , a newtonian - type attraction zone , a zone of cosmological acceleration and a zone beyond the cosmological horizon . when @xmath287 , there is a double extremal horizon , formed by the cauchy and event horizons together , in addition to the cosmological one .
when @xmath288 there are then three separate horizons , the cauchy , event and cosmological horizons . when @xmath289 , there is a separate cauchy horizon , and there is another horizon coincidence , the event horizon and the cosmological horizon come together . in this case
the black hole is a cosmological extremal supermassive regular black hole as the whole visible universe is swallowed by this supermassive object . for @xmath290
the cauchy horizon is also a cosmological horizon , the black hole is ultra massive . in such a universe
there is only one horizon , which is a cosmological one , together with a repulsion zone . in this context
this ultramassive black hole is a new type , in which all three horizons coincide - the cauchy , event , and cosmological horizons .
there is also the case for which @xmath291 where the three horizons coincide .    for @xmath233
, there is no cosmological horizon and depending on the values of the parameters , the black hole solution can have two horizons - the cauchy horizon and the event horizon . when the mass is below a critical mass @xmath203 , there are no horizons . for
@xmath292 a double horizon appears , and when @xmath79 exceeds @xmath203 there are the cauchy horizon and event horizons .    for @xmath234 ,
there is also no cosmological horizon and depending on the values of the parameters , the solution can have two horizons - the cauchy horizon and the event horizon .
there is also a critical mass @xmath203 , below which there are stars and above which there are regular black holes with two horizons .
the critical case is a black hole with a double extremal horizon .
geodesics of massive and massless test particles in these regular black hole geometries display a great variety of cases .
test particles can be in a state of rest on the bottom of the well , or oscillate around this minimum , other particles can overcome the repulsion barrier and be lodged in the central small cavity , and several other trajectories are possible .
however , particles can not leave the well into the exterior since there is a black hole horizon , i.e. , there is a one - way membrane .
abb and aez thank financial support from the program of competitive growth of kazan federal university ( kfu ) project no .
0615/06.15.02302.034 and from the russian foundation for basic research grant ( rfbr ) no .  14 - 02 - 00598 .
abb acknowledges financial support provided under the european union s framework program 7 ( fp7 ) european research council ( erc ) starting grant  the dynamics of black holes :  testing the limits of einstein s theory " grant agreement no .
jpsl thanks fundao para a cincia e tecnologia ( fct ) - portugal for financial support through project no .
pest - oe / fis / ui0099/2014 .
h. f. m. goenner , `` on the history of unified field theories .
1930  ca .
1965 ) '' , living rev .
relativity * 17 * , 5 ( 2014 ) . see also , h. f. m. goenner , `` some remarks on the genesis of scalar - tensor theories '' , gen .
. gravit . * 44 * , 2077 ( 2012 ) ; arxiv:1204.3455 [ gr - qc ] .
f. w. hehl and yu .
n. obukhov , `` how does the electromagnetic field couple to gravity , in particular to metric , nonmetricity , torsion , and curvature ? '' , in _ gyros , clocks , interferometers ... : testing relativistic gravity in space _ , edited by c. lmmerzahl , c. w. f. everitt , and f. w. hehl , lect .
notes phys .
* 562 * , 479 ( 2001 ) ; arxiv : gr - qc/0001010 .
a. b. balakin and j. p.
s. lemos , `` non - minimal coupling for the gravitational and electromagnetic fields : a general system of equations '' , classical quantum gravity * 22 * , 1867 ( 2005 ) ; arxiv : gr - qc/0503076 .
a. b. balakin , h. dehnen , and a. e. zayats , `` nonminimal einstein  yang - mills
higgs theory : associated , color and color - acoustic metrics for the wu
yang monopole model '' , , 124011 ( 2007 ) ; arxiv:0710.5070 [ gr - qc ] .
a. b. balakin , v. v. bochkarev , and j. p. s. lemos , `` nonminimal coupling for the gravitational and electromagnetic fields : black hole solutions and solitons '' , , 084013 ( 2008 ) ; arxiv:0712.4066 [ gr - qc ] .
t. t. wu and c. n. yang , `` some solutions of the classical isotopic gauge field equations '' , in _ properties of matter under unusual conditions , in honor of edward teller s 60th birthday _ , edited by h. mark and s. fernbach ( interscience , new york 1969 ) , p.  349
j. m. bardeen , `` non - singular general - relativistic gravitational collapse '' , in _ abstracts of the 5th international conference on gravitation and the theory of relativity _ , edited by v. a. fock et al .
( tbilisi university press , tbilisi , 1968 ) , p. 174 .
d. horvat , s. iliji , z. narani , `` regular and quasi black hole solutions for spherically symmetric charged dust distributions in the einstein - maxwell theory '' , classical quantum gravity * 22 * , 3817 ( 2005 ) ; arxiv : gr - qc/0409103 .