# Size of stripe domains in a superconducting ferromagnet

###### Abstract

In a superconducting ferromagnet, the superconducting state appears in the ferromagnetic phase where usually a domain structure has already developed. We study the influence of the superconducting screening currents on a stripe structure with out-of-plane magnetization, in a film of arbitrary thickness. We find that superconductivity always induces a shrinkage of the domains, and there is a critical value of penetration depth below which a mono-domain structure is more stable than the periodic one. Furthermore we investigate the possible different effects of singlet and triplet superconductivity on the domain width, as well as the conditions for the existence of vortices in the domains. The obtained results are then discussed in light of the experimental data of superconducting ferromagnets URhGe, UGe2, and UCoGe.

###### pacs:

75.60.Ch, 74.20.De, 74.25.Ha, 74.70.Tx## I Introduction and model

In the discovered Uranium-based superconducting ferromagnets (SFM) UGe saxena (), URhGe aoki2001 () and UCoGe huy () the Curie temperature is much higher than the superconducting critical temperature which means that superconductivity appears in the ferromagnetic state where it is usually a domain structure (DS) which develops. In previous works krey (); bulaevskii (); sonin1 (); faure (); sonin2 () the influence of the superconducting screening currents on the DS has been studied in the case of thick or bulk systems, when the thickness along the easy magnetization direction is much larger than both the transverse domain width (see Fig. 1) and the London penetration depth . The domain period at equilibrium results from the balance between a positive contribution to the energy density due to domain walls and a negative contribution from the magnetic induction, as the magnitude of both increases when is reduced. An exact energy minimization faure () shows that for the condition (where is an effective domain wall width that parametrizes the wall energy) superconductivity decreases the domain size due to partial penetration of the magnetic field near the domain wall. This energy decrease is proportional to and the formation of the domain wall is favorable when this contribution counterbalances the energy of the domain wall itself, which is proportional to . For the system is in a mono-domain state without any domain wall faure (); sonin1 (); sonin2 (). In the present article we extend Fauré and Buzdin’s work faure () by deriving the expression of the energy valid for any thickness and we discuss the DS in all limits, in particular when . Furthermore we investigate the possible effects of singlet and triplet superconductivity on the DS, as well as the conditions for the existence of vortices in the domains.

Model.— We consider a ferromagnetic film of thickness that can also become superconducting (see Fig. 1). The -axis is chosen perpendicular to the film with the surface edges at . Domain walls parallel to the -plane separate the periodic ferromagnetic structure into domains of equal width and with magnetization alternating along the -axis, i.e. with

(1) |

This means that we consider the domain wall thickness very small compared to and . The energy per surface unit is where the magnetic field energy is given by the relation

(2) |

and the superconducting current energy is expressed in the London limit degennes () (i.e. the superconducting coherence length )

(3) |

Here, is the total area the film surface. The explicit dependences with respect to the domain width result from the Ansatz we made for the magnetization . The contribution from the domain walls may be written as faure ()

(4) |

where is a domain wall characteristical length scale. Usually provides an upper limit of the real domain wall width faure (). Hereafter, the temperature will be considered phenomenologically within an explicit dependence of the London penetration depth . The latter is finite in the superconducting low phase, and diverges at the critical temperature .

Outline of the article.— In the second section we analytically calculate the field distribution and the energy of a periodic stripe DS in a SFM of arbitrary thickness. An exact expression of the energy is obtained as an infinite sum. In section III the equilibrium size of the domains is determined in the superconducting state as a function of the penetration depth. Analytical expressions of the domain width are obtained in limit cases. The stability of the domain structure is also discussed. We then investigate the difference between singlet and triplet superconducting states. Specifically, the triplet pairing is described by a domain wall energy which is temperature dependent, in contrast with the singlet pairing. In section IV we draw the condition for vortices appearance at the center of a domain. Finally in section V we apply our results to the DS in Uranium-based superconducting ferromagnets before concluding.

## Ii Method of solution

First, we minimize the energy with respect to the magnetic induction . This yields the London equation in the film and outside. Once is found by solving the London equation, the resulting total energy, , will be minimized with respect to to determine the domain width at equilibrium.

### ii.1 Magnetic field

We find the magnetic induction by Fourier expansion. Using Maxwell-Thomson equation and symmetry relations and , the field can be written as

(5) |

where we remind that . When solving the London equation in the film and Maxwell’s equation outside, we use the symmetry condition and the continuity condition at the surfaces. This yields, for ,

(6) |

and for ,

(7) |

with

(8) |

Figures 2 and 3 show the distribution of the magnetic field in the normal state and the superconducting state. It is plotted between the domain wall and the domain center. In the normal state, the magnetic field distribution in a wide domain () is concentrated around the domain wall (see Fig. 2(a)) while in a narrow domain () it is nearly uniform and equal to (see Fig. 3(a)). The supercurrent screens the field around the domain wall on a length scale (compare for example the contour in (a) and (b) of Fig. 3). Note that the supercurrent is responsible for a kink that the lines of constant show at the film surface (see Fig. 2(b) and Fig. 3(b)). This is because is discontinuous at the interface when . Flux lines also have a kink at the domain wall since the magnetization is discontinuous there.

### ii.2 Energy of the domain structure

To simplify the expressions of the energy it is convenient to introduce the normalized lengths

(9) |

where is Pearl’s penetration depth which, in the limit , takes the place of as the effective magnetic length scalepearl (). Note that in these notations . Furthermore we will consider the normalized energy and omit the bar henceforth. So the domain wall contribution is and the volume contribution, is

(10) |

where and . For a thick film, i.e. so that , Fauré and Buzdin’s result faure () is recovered as expected. Since the limit has been previously investigated faure () we discuss below the limit , that is , with more details.

#### ii.2.1 Regime

This parameter regime is realized in the normal phase where , as well as inside the superconducting phase, in the vicinity of the critical temperature . When (i.e., when ), the sum in (10) may be expanded in powers of before summing. Hence for ,

(11) |

while for ,

(12) |

#### ii.2.2 Regime and

This parameter regime may be realized in the superconducting phase, for either thin films or large domains. When the sum in (10) can be approximated by an integral. Euler-MacLaurin’s approximation yields

(13) |

where the is defined by Eq. (42) in the Appendix. Then using asymptotic expression (43) of in the limit , one finds

(14) |

and for , asymptotic expression (44) yields

(15) |

## Iii Domain width at equilibrium

### iii.1 General results

The equilibrium size is obtained by minimization of the total energy . Figures 4-6 show as a function of the penetration depth for different values of the normalized wall thickness . In the superconducting state decreases with decreasing temperature from infinity at to a finite value at . As discussed below, just below the domain width always decreases from the normal state value realized for . One can distinguish two regimes. When the domains in the normal state are wide () the decrease is negligible, while when the domains are narrow () their width can drastically shrink down to . Then, in both regimes, reaches a minimum value before diverging at a critical lower bound . This limit corresponds to the situation where the ferromagnetic induction is completely screened by the supercurrent faure (); sonin1 (); sonin2 () and there is no stray field that need to be accommodated by a DS.

The -dependence of the width is similar to the prediction for DS in superconductor-ferromagnet (S/F) bilayers bilayer (); stankiewicz () except for two aspects. First the shrinkage factor in a S/F hybrid is limited and bilayer (). Experiments performed on heterostructures made of a ferromagnetic garnet layer combined with a superconducting layer of Pb tamegai () or Nb vlasko-vlasov () have observed the shrinkage of the domain size due to superconductivity. The smallest observed shrinkage factor was substantially smaller than the theoretical prediction . This circumstance is a non-equilibrium effect and is related to the special dynamics of the vortex structure created in the superconducting layer and coupled to ferromagnetic domains vlasko-vlasov (). The second point is that although the DS can also be unstable in S/F bilayers when decreases, it exists a range of parameters (when the thickness of the F-layer is approximately larger than half the domain width in the normal state) in which remains finite even in the limit stankiewicz (), in contrast with SFM where the DS is never stable for a vanishing .

Hereafter, complementary to the exact numerical solution, we present the analytical results obtained in various asymptotic regimes.

#### iii.1.1 Normal state and vicinity of the transition

In the limit , using expression (11) valid for the vicinity of the superconducting transition where , one finds that the minimum of the energy is realized at

(16) |

which yields the domain width

(17) |

in the normal state (see Fig. 7). The domain width just below the superconducting transition is then related to the normal-state width by

(18) |

#### iii.1.2 Stability of the periodic domain structure

The periodic domain structure is always stable in the normal state where the energy minimum is obtained for a finite width . This is not the case in the superconducting state. In the limit and , the saddle-point equation for the energy is

(21) |

where the function is defined by Eq. (42) in Appendix. There is an energy minimum at a finite only for . Since increases with , this condition defines a lower bound below which the periodic structure is unstable. In the thought experiment where decreases from infinity (normal state) to zero, the period first decreases from but then it increases back before diverging as (see Fig. 4). Note that the limit of vanishingly small has been considered in Ref. sonin1 () and the conclusion of the absence of DS was made. When in the normal state, the approximations of yield

(22) |

while when the normal-state width , the periodic structure exists down to

(23) |

Figures 8 and 9 show that these analytical relations fit well the numerical results in the limit regimes.

One can now understand the difference of shrinkage factor between a narrow domain and a wide domain. For a narrow domain with and , the width in the superconducting state decreases with the penetration depth as . The shrinkage goes on until the DS is destabilized by an excessive superconducting screening. The minimum value that the width can have is then , which is much smaller than . In the opposite limit when and , the actual magnetic length scale is , not . And the domain width behaves as when is smaller than . Thus the minimum , so that a wide domain does not significantly shrink while the penetration depth decreases to the instability value .

### iii.2 The case of triplet superconductivity

As we will see in a next section, the temperature dependence of the domain equilibrium size can provide a phenomenological way of distinguishing triplet from singlet pairing. Here, we describe the general formalism; numerical applications to some specific compounds will be shown later as examples, in the framework of Uranium-based superconductors.

#### iii.2.1 Modification of the wall energy due to the depletion of condensation energy

In the case of a triplet superconductor with Cooper pairs fully spin polarized by the intrinsic magnetization of the ferromagnetic domains, a depletion of the order parameter occurs on a width at both sides of the interface separating two domains of opposite magnetizations. This inhomogeneity results in a loss of condensation energy at a domain wall. At a first approximation the domain wall energy for is modified by replacing the effective domain wall width by the temperature dependent defined by

(24) |

where , and

(25) |

is proportional to the ratio of the superconducting energy at to the magnetization energy. Here, we have assumed , , and . Using the asymptotic expressions relating to the normal state width , one can estimate from experimental data in the limit with

(26) |

and in the limit with

(27) |

#### iii.2.2 Behavior of the domain width below the transition

A triplet order parameter yields a vanishing correction to the domain wall energy at the transition. The temperature dependence of the domain size is in the limit ,

(28) |

and in the limit ,

(29) |

So, as in the singlet case, always decreases with decreasing temperature in the vicinity of . However at a lower temperature and for a large parameter , the behavior for the triplet case can be significantly different from the singlet case. Since for a given there is an universal curve for the singlet pairing, it is possible to detect an experimental signature of triplet superconductivity as a deviation from this curve. This is schematically illustrated in Fig. 10 which shows the temperature dependence of in the singlet and in the triplet case. When temperature is reduced from (on the right side where ) to (on the left side where ), the -curve in the triplet case interpolates between the curve (a) obtained with a constant and the curve (b) for . The triplet-case curve follows the singlet-case curve with at high temperature but at a lower temperature a deviation appears which is maximum at .

In the expression (29) of for the triplet pairing, the first-order term is proportional to while the second-order term is proportional to so one can expect that the singlet/triplet discrepancy can be amplified when reducing the thickness of the film. As shown in Fig. 10 this discrepancy is enhanced when the ratios and become larger. For small enough thickness the difference is not only quantitative but is also qualitative as the triplet-pairing curve can increase with decreasing temperature while the singlet-pairing curve decreases.

## Iv Condition for vortex appearance

In this section we estimate the energy for the creation of a vortex-antivortex pair in the middle of adjacent domains. This additional energy is the sum of (i) the energy decrease when (anti)vortices are driven by screening currents away from domain boundaries to domain centers and (ii) the concomitant increase of the interaction energy between one vortex and its antivortex as they move apart from each other. For simplicity we neglect the vortex-core energy assuming the limit of large .

### iv.1 Energy decrease due to the screening current

The force of the screening current acting on a vortex degennes () is , where the quantum of flux T.m. When the vortex moves from the domain wall to the domain center, the current then produces a work equal to

(30) |

where is the coupling energy of the flux carried by one vortex with the magnetization of a ferromagnetic domain.

The work has simple analytical expressions in limit cases. In the regime ,

(31) |

where the Lerch transcendent . If furthermore , expression (31) can be simplified as

(32) |

Thus increases with the width of the domain and, when the latter becomes large compared to the magnetic radius of a vortex, the work is simply equal to the coupling energy .

In the opposite regime , the screening current is vanishingly small so its work is only equal to a fraction of . For wide domains (), the dependence of on the domain width is linear as

(33) |

where Catalan’s constant , while for narrow domains (), the dependence is quadratic as

(34) |

### iv.2 Interaction energy of a vortex-antivortex chain

The interaction energy of a chain composed by one vortex (antivortex) located at the center of every domain of positive (negative) magnetization is the magnetic energy required to separate by a distance vortices and their paired antivortices created at domain walls. In the case or , the energy per one vortex buzdin () has the asymptotic expression

(35) |

We have assumed . In a large- superconductor with , there is always a temperature regime close to in which .

In the opposite limit, for the interaction energy per vortex is

(36) |

while for , it is

(37) |

### iv.3 Condition for the appearance

From the above results we can derive the condition of the vortex stability in a superconducting domain structure. The latter is assured when the total energy per vortex is negative. This condition defines a critical value of magnetization above which vortices appear. One can distinguish four regimes in which has a simple analytical expression. In a film with narrow domains (i.e. ),

(38) |

and

(39) |

Note that Eq. (38) recovers the expected result that vortices appear when the magnetization exceeds the lower critical field of a non-magnetic superconductor defined as . However, when the penetration depth is larger than the domain width, the coupling of the vortex flux with the DS magnetization is not energetically optimum so the critical magnetization (39) is higher than the mono-domain value .

## V Application to Uranium-based superconducting ferromagnets

(J/m) | (nm) | (m) | (mm) | (m) | (nm) | (J/m) | ||
---|---|---|---|---|---|---|---|---|

UGe | 14200 | 15 | 1 | 2 | 4 | 13.6 | 40 | 0.039 |

UCoGe | 51 | 15 | 1.2 | 0.15 | 2 | 45 | 29 | 2.36 |

URhGe | 3070 | 18 | 0.9 | 0.2 | 20 | 3450 | 30 | 6.5 |

According to experimental data (see Table 1), the Uranium-based superconductors UGe, UCoGe and URhGe are in the limit . The effective domain wall thickness can then be easily calculated in this limit from the experimental value of the domain width in the normal state. We roughly estimate the zero-temperature condensation energy with the BCS formula where is the volumic specific heat jump at . Among the three compounds UCoGe has the smallest magnetization resulting into the largest (see Table 1). It is then the most promising candidate where to look for possible observation of the non-monotonic temperature dependence of the domain width in the superconducting phase. Using its estimated parameters we plotted the temperature dependence of the domain width for the singlet and the triplet scenario in a film of half-thickness mm (Fig. 11) and m (Fig. 12). For mm (Fig. 11) the domain width decreases below in both the singlet and the triplet case but there is a significant difference of size at between the two. By reducing the film thickness the discrepancy is amplified. For instance with m (Fig. 12) one obtains a qualitative difference: for the triplet case the width increases by 32% from its value in the normal state while for the singlet case it shrinks by 21%.

## Vi Conclusion

We have done the complete analysis, within the London approximation, of the domain width at equilibrium in a SFM film of arbitrary thickness . We have shown that the ratio follows a universal dependence on the normalized penetration depth and the normalized effective wall thickness (see Figs. 4-6). In addition to the exact numerical dependence, we have derived analytical expressions of this relation in limit cases. In particular, we have recovered the previously published results by Fauré and Buzdin faure () which had been established for and , and we have complemented them with the analysis in the opposite limit and .

We have found that the domain width always decreases with temperature when the DS enters the superconducting phase. The screening supercurrent induces a decrease of the domain width if the latter is larger than the penetration depth. However, the paramagnetic screening suppresses the DS below a critical value of penetration depth, which means that the system may become mono-domain at a finite temperature if . The domain shrinkage is relatively small when the domain width in the normal state is much larger than the film thickness. Indeed, in this limit, the width decreases as which is only a fraction smaller than since . With the other limit shape of domains, i.e. , the shrinkage can be important since then the domain width decreases as while the critical penetration depth .

Furthermore, we have investigated the effect of the triplet pairing within an effective temperature-dependent renormalization of the domain wall energy. Actually, the DS magnetization alternatively suppresses the order-parameter components with opposite spin-projections. This inhomogeneity results in an additional positive contribution to the domain wall energy, that is absent for the singlet pairing. The supplementary term vanishes at the transition temperature but can be significant at a lower temperature. The variation of the domain width is then no more described by the universal dependence of on obtained for the singlet case. This gives a phenomenological way of distinguishing triplet from singlet pairing in experiments. And, as we have shown, the discrepancy can be amplified by reducing the film thickness (see Figs. 11 and 12).

We have also established that vortices can appear when the magnetization exceeds a critical value . In narrow domains (i.e. ), the latter is equal to the lower critical field of a non-magnetic superconductor if , but in the opposite limit , it is much larger than . For wide domains (i.e. ), we have found that the critical magnetization behaves similarly, with the difference that is then replaced by Pearl’s effective penetration depth .

The available experimental data (Table 1 and Ref. saxena (); ureview (); uge2 (); sakarya (); huy (); ureview (); ucoge (); hykel (); aoki2001 (); ureview (); urhge (); dolocan ()) for UGe, URhGe and UCoGe show that in these compounds and the conditions for the vortices appearance (section IV) are fulfilled. This means that the effects of domain shrinkage (considered in section III for the case where the vortices are absent) should be weakened. Another difficulty to experimentally observe the evolution in the superconducting state may be related with the pinning of the domain wall and/or vortex pinning. At present there are no convincing experimental data on the change of the domain structure below .

It has been recently observed in S/F bilayers tamegai (); vlasko-vlasov () that the coupling between vortices and magnetic domains leads to a strong shrinkage of domains in the presence of the oscillating field used for equilibration of the domain structure. It would be interesting to perform the similar experiments with superconducting ferromagnets to study these coupling effects.

We thank K. Hasselbach, C. Paulsen, J.-P. Brison, and D. Aoki for helpful discussions. This work has been supported by the French ANR project SINUS.

## Appendix A The function

We define the function as the integral

(42) |

### a.1 Limit

The integral can be rewritten as

so by using the approximations in the numerator and in the denominator when , and when ,

Then for

(43) |

### a.2 Limit

In the limit , then

(44) | |||||

## References

- (1) S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, R. K. W. Haselwimmer, M. J. Steiner, E. Pugh, I. R. Walker, S. R. Julian, P. Monthoux, G. G. Lonzarich, A. Huxley, I. Sheikin, D. Braithwaite, and J. Flouquet, Nature 406, 587 (2000).
- (2) D. Aoki, A. Huxley, E. Ressouche, D. Braithwaite, J. Flouquet, J.-P. Brison, E. Lhotel, and C. Paulsen, Nature 413, 613 (2001).
- (3) N. T. Huy, A. Gasparini, D. E. de Nijs, Y. Huang, J. C. P. Klaasse, T. Gortenmulder, A. de Visser, A. Hamann, T. Görlach, and H. v. Löhneysen, Phys. Rev. Lett. 99, 067006 (2007).
- (4) U. Krey, Int. J. Magn. 3, 65 (1972).
- (5) L. N. Bulaevskii, A. I. Buzdin, and S. S. Crotov, Solid State Commun. 48, 719 (1983).
- (6) E. B. Sonin, Phys. Rev. B 66, 100504(R) (2002).
- (7) M. Fauré and A. I. Buzdin, Phys. Rev. Lett. 94, 187202 (2005).
- (8) E. B. Sonin, Phys. Rev. Lett. 95, 269701 (2005); M. Fauré and A. I. Buzdin, ibid. 95, 269702 (2005).
- (9) P. G. de Gennes, Superconductivity of metals and alloys, Perseus Books (Reading, Massachusetts, 1999).
- (10) J. Pearl, Appl. Phys. Lett. 5, 65 (1964).
- (11) A. I. Buzdin, A. S. Melnikov, A. V. Samokhvalov, T. Akashi, T. Masui, T. Matsuda, S. Tajima, H. Tadatomo, and A. Tonomura, Phys. Rev. B 79, 094510 (2009).
- (12) L. N. Bulaevskii and E. M. Chudnovsky, Phys. Rev. B 63, 012502 (2000); L. N. Bulaevskii, E. M. Chudnovsky, and M. Daumens, Phys. Rev. B 66, 136502 (2002); E. B. Sonin, Phys. Rev. B 66, 136501 (2002); M. Daumens and Y. Ezzahri, Phys. Lett. A 306, 344 (2003).
- (13) A. Stankiewicz, S. J. Robinson, G. A. Gehring, and V. V. Tarasenko, J. Phys.: Condens. Matter 9, 1019.
- (14) T. Tamegai, Y. Nakao, and Y. Nakajima, J. Phys.: Conf. Ser. 150, 052263 (2009); T. Tamegai, Y. Nakao, S. Mohan, and Y. Nakajima, Supercond. Sci. Technol. 24, 024015 (2011).
- (15) V. Vlasko-Vlasov, U. Welp, W. Kwok, D. Rosenmann, H. Claus, A. A. Buzdin, and A. Melnikov, Phys. Rev. B 82, 100502(R) (2010).
- (16) C. Kittel, Introduction to Solid State Physics (Wiley, Berkeley 1996), 7th ed.
- (17) L. D. Landau and E. M. Lifschitz, Electrodynamics of Continuous Media (Nauka, Moscow, 1982).
- (18) C. Pfleiderer, Rev. Mod. Phys. 81, 1551 (2009); D. Aoki, T. D. Matsuda, F. Hardy, C. Meingast, V. Taufour, E. Hassinger, I. Sheikin, C. Paulsen, G. Knebel, H. Kotegawa, and J. Flouquet, arXiv 1012.1987v1 (2011).
- (19) K. Oikawa, T. Kamiyama, H. Asano, Y. Onuki and M. Kohgi, J. Phys. Soc. Japan 65 3229 (1996); N. Tateiwa, T. C. Kobayashi, K. Hanazono, K. Amaya, Y. Haga, R. Settai, and Y. O. Nuki., J. Phys.: Condens. Matter 15, L17 (2001); E. D. Bauer, R. P. Dickey, V. S. Zapf and M. B. Maple, J. Phys.: Condens. Matter 13, L759 (2001).
- (20) S. Sakarya, N. H. van Dijk, and E. Brück, Phys. Rev. B 71, 174417 (2005).
- (21) N. T. Huy, D. E. de Nijs, Y. K. Huang, and A. de Visser, Phys. Rev. Lett. 100, 077002 (2008); P. de la Mora and O Navarro, J. Phys.:Condens. Matter 20, 285221 (2008); A. Gasparini, Y. K. Huang, J. Hartbaum, H. v.Lohneysen, and A. de Visser, Phys. Rev. B82, 052502 (2010); A. Gasparini, Y. K. Huang, N. T. Huy, J. C. P. Klaasse, T. Naka, E. Slooten, and A. de Visser, J. Low Temp. Phys. 161, 134 (2010).
- (22) D. Hykel, Microscopie à micro-SQUID: étude de la coexistence de la supraconductivité et du ferromagnétisme dans le composé UCoGe, PhD Thesis, University of Grenoble (2011).
- (23) K. Prokes, T. Tahara, Y. Echizen, T. Takabatake, T. Fujita, I. H. Hagmusa, J. C. P. Klaasse, E. Brück, F. R. de Boer, M. Divis, and V. Sechovsky, Physica B 311, 220 (2002).
- (24) V. O. Dolocan, Imagerie magnétique des supraconducteurs non conventionnels, PhD Thesis, University of Grenoble (2005).