## Thermodynamic nature of the 0 -π quantum transition in superconductor-ferromagnet-superconductor trilayers

N. Pompeo, 1 K. Torokhtii, 1 C. Cirillo, 2 A. V. Samokhvalov, 3,4 E. A. Ilyina, 2 C. Attanasio, 2 A. I. Buzdin, 5 and E. Silva 1

1 Dipartimento di Ingegneria, Universit` Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy a 2 CNR-SPIN Salerno and Dipartimento di Fisica 'E. R. Caianiello', Universit` a di Salerno, I-84084 Fisciano (SA), Italy 3 Institute for Physics of Microstructures, Russian Academy of Sciences, 603950 Nizhny Novgorod, GSP-105, Russia 4 Lobachevsky State University of Nizhny Novgorod, Nizhny Novgorod 603950, Russia 5 Institut Universitaire de France and University Bordeaux, LOMA UMR-CNRS 5798, F-33405 Talence Cedex, France

(Dated: August 5, 2014)

In structures made up of alternating superconducting and ferromagnet layers (S/F/S heterostructures), it is known that the macroscopic quantum wavefunction of the ground state changes its phase difference across the F-layer from 0 to π under certain temperature and geometrical conditions, whence the name '0 -π ' for this crossover. We present here a joint experimental and theoretical demonstration that the 0 -π is a true thermodynamic phase transition: microwave measurements of the temperature dependence of the London penetration depth in Nb / Pd 0 84 . Ni 0 16 . / Nb trilayers reveal a sudden, unusual decrease of the density of the superconducting condensate (square modulus of the macroscopic quantum wavefunction) with decreasing temperature, which is predicted by the theory here developed as a transition from the 0 -to the π -state. Our result for the jump of the amplitude of the order parameter is the first thermodynamic manifestation of such temperature-driven quantum transition.

PACS numbers: 74.25.nn, 74.45.+c, 74.78.Fk

## I. INTRODUCTION

Superconductor (S)/Ferromagnet (F) structures are the ideal playground for the search for novel phase transitions: two different, competing orderings come into play in a controllable fashion. In fact, S/F heterostructures can be grown with different S and F thicknesses, d s and d f , respectively, so that the superconducting and ferromagnetic ordering have different effect. Moreover, other external parameters, such as the temperature T , can be used to vary the interactions between F and S. In the recent past, the attention has been focussed on the damped oscillatory behaviour of the Cooper pairs superconducting wavefunction ∆ = | ∆ | e iφ in the ferromagnet: 1,2 the peculiarity of the S/F interaction makes the superconducting wavefunction ∆ oscillate in the F-layer, in addition to the conventional, 'proximity-effect-like', exponential decrease. This unique feature of S/F structures is at the origin of the superconductor-ferromagnetsuperconductor (SFS) π -Josephson junctions (for a review, see Refs. 3,4) characterized by a ground state at phase difference δφ = π between S-layers (one can represent this effect by a sign change of the macroscopic superconducting wavefunction across the F-layer). In this case, one speaks of '0 -π transition'. A sketch of the spatial dependence of the wavefunction in the 0 -and π -states is reported in Fig. 1: in the '0' state, ∆ is only depressed in the F-layer. In the ' π ' state, ∆ has zero value in the middle of the F-layer, and changes sign. The π shift has, among the others, the spectacular consequence of spontaneous supercurrents in ring-shaped structures incorporating π junctions. 5

Up to now, experiments were directed towards transport observations of the π -shift of the wavefunction: the efforts were concentrated in measurements of the critical current of Josephson junctions with F barrier, SFS Josephson junctions. In order to avoid excessive depression of the superconducting wavefunction and in order to have well-established superconducting electrodes, the typical structure used to detect the π -shift was made up of relatively thick S-layers (on the scale of the superconducting coherence length ξ s ), and F-layers of thickness of the order of the ferromagnetic coherence length, ξ f = √ /planckover2pi1 D /E f ex = √ D /h f (here, D f is the diffusion coefficient in the ferromagnet, E ex is the exchange energy and h represents the exchange field - in appropriate units - acting on the electron's spins 6,7 ). Since ξ f is very small (few nanometers) in strong ferromagnets, weak ferromagnets have been often employed for the ease of controlling the geometrical conditions for the 0 -π transition. The sign change of ∆ has been experimentally observed in SFS Josephson junctions, so that evidence for the π -shift is nowadays robust: measurements of the critical current of SFS junctions with CuNi alloys as F-barrier 6,7 showed 0 -π crossover as a function of d f . Moreover, the subtle role of the temperature T has been revealed: only in the vicinity of the so-called critical thickness condition ( d f /similarequal ξ f ) can the temperature trigger the 0 -π crossover. 6,7 Also, the π -shift was observed in SFS junctions with PdNi F-barrier 8 with d f ∼ 6 nm, even if the transition could not be driven by the temperature in the latter case. So, SFS Josephson junctions with CuNi al-

loys as F-layer remain the only example of temperature mediated 0 -π crossover.

Aim of this work is to assess the thermodynamic nature of the 0 -π transition: while the crossover clearly exists in the relative phase between two superconducting electrodes , no evidence for the change of the order parameter | ∆ | 2 of the macroscopic wavefunction has been brought forward: the finding of a change in | ∆ | 2 at the 0 -π transition would be a direct evidence for the thermodynamic nature of such a transition. In particular, the temperature dependence of the amplitude of the order parameter at the critical thickness has not been studied neither theoretically nor experimentally: measurements of the London penetration depth in Nb/Ni bilayers (hybrids based on a strong ferromagnet) as a function of the Ni layer thickness 9 showed only the d f dependence of the extrapolated zero-temperature superfluid density, similarly to singleT/T c estimates in Nb/PdNi/Nb structures, 10 with no reported temperature dependence.

We are interested in measuring the temperature dependence of the amplitude of the overall macroscopic wavefunction. This aim imposes constraints on the geometrical structure under investigation: while in previous experiments 6-8,11 the superconducting electrodes were thick enough so that the influence of the 0 -π transition on the S-layers superconducting wavefunction ∆ was negligible, here thin S-layer will be more suited.

In this Article we present joint experimental and theoretical results, showing the thermodynamic nature of the 0 -π transition. We perform measurements of the temperature dependence of the effective London penetration depth, ∆ λ eff ( T ), in Nb Pd / 0 84 . Ni 0 16 . / Nb trilayers with relatively small d s ∼ ξ s . In the sample with d f = 2nm, close to the critical thickness, we observe a reentrant jump in the ∆ λ eff ( T ) curve: λ eff increases with lowering temperature, i.e. the superfluid density decreases, a strong experimental evidence of temperature mediated 0 -π transition. As expected, this phenomenon is absent in trilayers with different values of d f . In agreement with the experimental data, the theoretical analysis predicts the upward jump of ∆ λ eff as the temperature decreases providing a qualitative description of the observed phenomenon.

## II. SAMPLES GROWTH AND EXPERIMENTAL SETUP

Nb/Pd 0 84 . Ni 0 16 . /Nb trilayers were grown on Al O 2 3 substrates by ultrahigh vacuum dc diode magnetron sputtering at an Ar pressure of 3 · 10 -6 Torr after obtaining a base pressure of 2 · 10 -8 Torr following the procedure described in Ref. 12: the samples were all prepared in the same deposition run thanks to the presence of a movable shutter in the deposition chamber which selectively covers the substrates, glued by silver paste on the holder, which were kept at room temperature during the deposition process. The typical deposition rates

FIG. 1: (Color online) Schematic spatial dependence of the superconducting order parameter in S/F/S trilayers: the curve represents a sketchy behaviour of the pair wave function. (a) In the even mode (0-phase), due to the symmetry, the derivative of the pair wave function vanishes at the centre of the F-layer. (b) In the odd mode ( π -phase), the pair wave function itself vanishes at the centre of the F-layer, and it has a π -shift in its phase in diametrically opposite points.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0665v1-figure-7f7d17e2-922a-45c6-9f79-f7170be9ef4d.jpg)

were 0.28 nm/s for Nb and 0.40 nm/s for Pd 0 84 . Ni 0 16 . measured by a quartz crystal monitor previously calibrated by low-angle X-ray reflectivity measurements on deliberately deposited thin films of each material. The Nb layers have the same nominal thickness in all the samples of the series, d s = 15 nm, while the thickness of the Pd 0 84 . Ni 0 16 . layer changes from d f = 2 nm to 9 nm (the range where we expect d f ∼ ξ f ). The total thickness of the trilayers is then d = 2 d s + d f . A pure Nb film with d = 30 nm was also grown for comparison (labelled in the following with d f = 0 nm). It has a superconducting critical temperature T c 0 = 7 5 K and . ξ s = √ D / πT s 2 c 0 = 6 nm (here D s is the diffusion coefficient in the superconductor). 13 Estimates of E ex /similarequal 14

meV and ξ f /similarequal 3 nm for the Pd 0 84 . Ni 0 16 . alloy have been reported elsewhere 12 . The complete electrical characterization of S/F/S trilayers has been previously reported in Refs. 14,15. High-resolution transmission electron microscopy showed excellent crystallinity of the Nb layers, roughness at the SF interfaces less than 1 nm, leading to good interface transparency, and suggested some interdiffusion of Nb into the PdNi layer. 16 In a previous study, 17 the local atomic structure was investigated to assess whether the F-layer could induce significant disorder in the S-layers. To this aim, extended X-ray absorption spectroscopy at the Nb K-edge was performed. The results did not show correlations between the structural disorder in the Nb layer and the superconducting properties, such as T c , H c 2 and the microwave results described in the following.

Microwave measurements were performed using the dielectric resonator technique. 15,18 The quantity experimentally measured is the resonant frequency ν 0 of the resonator incorporating the sample. The sample was placed as an end wall in a cylindrical dielectric-loaded resonator. The resonant frequency of the resonator ν 0 depends on the energy of the electromagnetic field stored in the volumes of the resonator and of the sample portion where the screening of the field, occurring on a length scale λ eff , is not complete. Since the empty resonator gave no additional temperature dependencies in the small temperature range here explored, the resonant frequency of the resonator ν 0 ( T ) changed only as a consequence of the change of the sample screening. Thus, the experimentally accessible quantity is the temperature variation of the effective penetration depth with respect to a given temperature T ref , ∆ λ eff ( T ) = λ eff ( T ) -λ eff ( T ref ), which is obtained from ν 0 through the relation:

$$\Delta \lambda _ { e f f } ( T ) = - \frac { G } { \pi \mu _ { 0 } } \frac { \nu _ { 0 } ( T ) - \nu _ { 0 } ( T _ { r e f } ) } { \nu _ { 0 } ^ { 2 } ( T _ { r e f } ) }, \quad ( 1 ) \, \stackrel { \widehat { \text{pera} } } { \text{ref} } } \\ \text{for } \, i \, \real.$$

where µ 0 = 4 π · 10 -7 H / m and G is a calculated geometrical factor.

The microwave assembly was placed in a 4 He cryostat where temperatures down to 2.8 K were reachable. The cylindrical resonator was loaded with a Rutile (TiO ) 2 cylinder, with negligible temperature dependence of the complex permittivity below 10 K. The resonant mode chosen was the TE 011 , with circular induced currents on the sample. A magnetic field up to µ H 0 = 0.7 T ould be applied perpendicular to the sample plane. The setup has been extensively described previously. 19

Since the total thickness d of the S/F/S structure is smaller than the penetration depth, the London penetration length and the losses (see below) are averaged over the whole sample. In this full-penetration regime one has for the effective penetration depth: 20

$$\lambda _ { e f f } ( T ) = \lambda ( T ) \coth \frac { d } { \lambda } \approx \frac { \lambda ^ { 2 } ( T ) } { d }. \quad \ \ ( 2 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The measured ∆ λ eff directly compares to the square of the superconducting parameter through 1 /λ eff ∝ | ∆ . | 2

Note that due to the small thickness of the F-layer and to the fact that its conductivity is smaller than that of the Nb layers, 12 its contribution to the averaged London penetration depth is negligible. It is also important to stress that the resonant mode used (TE 011 ) induces only in-plane microwave currents: no current across the SF boundaries are involved and, as a consequence, ∆ λ eff is related only to the superconducting order parameter | ∆ , without contributions from tunnelling between lay-| 2 ers. In the mixed state, the microwave currents set in motion the quantized flux lines which then contribute to the field attenuation (losses) and screening. The simultaneous measurement of the resonator quality factor and resonant frequency allows to determine the vortex resistivity ρ v , related to the forces acting on the quantized flux lines and ultimately to the vortex pinning constant k p (see below).

## III. EXPERIMENTAL RESULTS

Raw data for the resonant frequency are reported in Fig. 2a, for the samples with d f = 0, 2, and 8 nm. The flattening of ν 0 close to T c is due to the loss of sensitivity when the samples become nearly electromagnetically transparent, and it is not relevant for the present purposes. The raw data show (i) that for d f = 8 nm the temperature dependence of the screening is smoothened with respect to pure Nb, and (ii) that the sample with d f = 2 nm exhibits a nonmonotonous temperature dependence of the screening. This last effect is the main experimental result of this paper, and we will concentrate on it in the following. The data of Fig. 2a were converted to ∆ λ eff , and plotted as a function of the reduced temperature, t = T/T c in Fig. 2b. In order to have the same reference reduced temperature, t ref = T ref /T c = 0 72, . for all the samples, different T ref had to be set. It is evident that Nb ( d f = 0) and the sample d f = 8 nm do not show any significant feature. Focussing on the sample with d f = 2 nm, the temperature dependence of ∆ λ eff can be divided in two regimes. Close to T c it is more similar to the one of the Nb film, while as T is lowered it shows a jump upward, followed by a smoother behaviour, similar to what shown by the sample with d f = 8 nm. The data have not been scaled vertically: the crossover between the two different regimes is unambiguous. It is worth to stress that the upward jump of the screening with decreasing T implies the decrease of the superfluid density, that is of | ∆ , with decreasing | 2 temperature. This nonmonotonous, reentrant behaviour is clearly the most striking result, and thus it must be thoroughly checked. The check for reproducibility was performed by repeating the measurements after disassembling and reassembling the resonator and the sample holder. While the absolute ν 0 changed (as it is expected), the temperature shift yielded identical results. The results are reported in Fig. 2c, where we show the full temperature range accessible, and they show that the two

measurement sets experimental artifacts.

FIG. 2: (Color online) Temperature dependence of the screening exhibited by the trilayers: (a) Raw data for the resonant frequency of the resonator, ν 0 ( T ), incorporating different samples with d f = 0, 2, and 8 nm. The screening vanishes when the critical temperature is approached. A peak is evident in the sample with d f = 2 nm, showing nonmonotonous screening. The data for the sample with d f = 8 nm have been shifted downward by 4 MHz for clarity. The flat part at the superconducting transition is not a saturation, but the region where the resonator loses sensitivity. (b) Effective penetration depth ∆ λ eff ( ) = t λ eff ( ) t -λ eff ( t ref ) for the same samples, obtained with t ref = 0 72. . With d f = 2 nm, the crossover at t ∼ 0 9 between two different regimes, corresponding to a . decrease of the superfluid fraction with decreasing temperature, is evident. (c) Reproducibility test of the peak in ∆ λ eff ( T ) after disassembling and reassembling the resonator (empty symbols). Full symbols, original data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0665v1-figure-bcb0cab2-ac67-493b-b368-6b1b79629e59.jpg)

Since the reentrant superconducting superfluid is a very relevant result, we estimated λ 2 ∝ | ∆ | -2 in a different experiment. We applied a magnetic field H perpendicular to the sample plane, thus entering the mixed state of the superconducting structure. Then, quantized flux lines are pinned by defects, which exert a recalling force (per unit length) that can be assumed to be elastic as a first approximation (small displacements δr ): F p = -k δr p , where k p is the vortex pinning constant.

Since k p ultimately depends on the energy gain of the vortex sitting on a defect, as a crude approximation λ is connected to k p by the approximate equality between elastic and condensation energy: 21 1 2 k ξ p 2 s ≈ c p 1 2 µ H ξ 0 2 c 2 ( c p is a constant of order unity). Using for the thermodynamic critical field H 2 c ≈ H H c 1 c 2 one finds:

$$k _ { p } \approx c _ { p } H _ { c 2 } \Phi _ { 0 } / 4 \pi \lambda ^ { 2 } ( T ) \text{ \quad \ \ } ( 3 )$$

where Φ 0 is the flux quantum. The upper critical field H c 2 is obtained from dc resistivity or from the disappearance of the field-dependent microwave signal. 14 The pinning constant k p is easily obtained from microwave measurements of the magnetic field dependent vortexmotion complex resistivity ρ v ( H ) = ρ v 1 ( H ) + i ρ v 2 ( H ), within a wide class of models for vortex relaxational dynamics. 22-25 All the details, as well as the uncertainty intervals, associated to the determination of k p , have been discussed previously. 22 Thus, we can estimate the temperature dependence of the superfluid density using a conceptually different experiment, with the benefit of making the experiment in situ : we measured ρ v ( H ) at several temperatures with the same resonator described above, and we derived k p ( H ) within the GittlemanRosenblum model. 23 Numerical differences in the estimate of k p according to different models 22 are absorbed in c p , and are not expected to affect the temperature dependence.

Fig. 3a presents sample measurements of the magnetic field dependence of the normalized vortex-state resistivity at T = 3 66 K ( . t = 0 72) for the sample with . d f = 2 nm. The field dependence of the derived k p is reported in Fig. 3b. We have chosen to use the maximum value, k p,max , in Eq.(3) to derive ∆ λ eff . Measurements of k p,max at several temperatures, reported in Fig. 3c, were thus converted into estimates for ∆ λ eff using the measured H c 2 ( T ), reported in Fig. 3d. The datum point for ∆ λ eff at T = 3 66 K exactly matches the zero-field data . by taking c p = 1 03, which is kept fixed. . All the remaining points, plotted in Fig. 3e as large purple circles, show then a very satisfactory agreement with the zero-field data, especially taking into account the crudeness of the model. 26 This is an important check of the experimental finding, and concludes the experimental results of this paper: for the sample with d f = 2 nm, where the condition d f ∼ ξ f is approximately fulfilled, 12 the temperature dependence of the effective penetration depth (and thus of the superfluid density), independently probed both in the Meissner state and in the vortex state, exhibits the same nonmonotonous behaviour, with a step-like increase of the superfluid with increasing temperature at t ∼ 0 9. .

## IV. THEORY

With the support of thorough theoretical calculations, we argue in the following that the observed reentrant jump in the superfluid density is a thermodynamic manifestation of the temperature-induced 0 -π transition.

FIG. 3: (Color online) Confirmation of the nonmonotonous behaviour of the superfluid density from measurements in the mixed state: (a) normalized complex vortex resistivity at T = 3 66 K ( . ρ n = 22 µ Ω · cm, Ref. 15). (b) pinning constant k p at T = 3 66 K, and definition . of k p,max . (c) k p,max as a function of T . (d) Upper perpendicular critical field H c 2 ( T ); the continuous line is a fit with H c 2 (0) = 2.53 T and T c = 5.06 K. (e) Comparison of ∆ λ eff ( T ) = λ eff ( T ) -λ eff ( T ref ) (here T ref = 3.2 K) as measured in the Meissner state (red dots) and as derived from k p,max (purple circles) using Eq.(3).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0665v1-figure-0db7e88c-7054-4c01-8347-8568ee6f1abb.jpg)

In order to explain the experimental results, we aimed at calculating the effect of the 0 -π crossover on the superfluid density of an F-layer sandwiched between two S-layers. In particular, two main issues were of interest: (i) assessing the existence of a double transition with lowering the temperature, first in the 0 -state and, at lower temperatures, in the π -state and, most important, that (ii) the superfluid density in the π -state was smaller than in the 0 -state, around the transition. The theoretical model adopted for the calculations should be able to incorporate the main experimental features.

The model is set up as follows. Since the F-layer is a weak ferromagnet we use the Usadel equations 27 for the averaged anomalous Green's functions F and F s for the F- and S-regions, respectively (see Ref. 3 for details). The presence of magnetic disorder (always present in magnetic alloys, and responsible for the main mechanism of the temperature induced 0 -π transition 7,28,29 ) was taken into account with the introduction of the magnetic scattering rate τ -1 s . Thin S-layers ( d s /lessorsimilar ξ s ) were introduced and yielded a nearly constant order parameter in the S-layers. We assumed transparent interfaces. We expect that the main correction in case of nonideal interfaces is a reduction of the effective F thickness due to some interdiffusion processes. We first note that in practice the exchange field h acting on the electron's spins in the ferromagnet and the magnetic scattering rate τ -1 s are much larger than the superconducting critical temperature. Second, the assumption of transparent interfaces yields the boundary conditions 30 at x = ± d / f 2:

$$F _ { s } = F, \, \sigma _ { s } \, \partial _ { x } F _ { s } | _ { s } = \sigma _ { f } \, \partial _ { x } F | _ { f }, \, \text{ \quad \ \ } ( 4 )$$

where σ f and σ s are the normal-state conductivities of the F- and S-metals. Applying the method developed in Ref. 31, we find the solution of the non-linear Usadel equation in the F-layer near the superconducting critical temperature obtaining the self-consistency equation for the superconducting order parameter ∆, which in turn yields the expansion of the free energy near the critical temperature. For the 0-state ( π -state) we should choose the even (odd) anomalous Green's functions F . As a result we obtain the expansion of the free energy F 0 ,π ( T ) in the 0- or π -states (indicated by the superscript) near the critical temperature as well as the critical temperatures T 0 ,π c of the two states, and the stable state is then determined. The expansion of the free energy reads:

$$F ^ { 0, \pi } ( T ) = E _ { 0 } \left [ a ^ { 0, \pi } \frac { T - T _ { c } ^ { 0, \pi } } { T _ { c } ^ { 0, \pi } } \Delta ^ { 2 } + \frac { b ^ { 0, \pi } } { 2 } \, \Delta ^ { 4 } \right ], \quad ( 5 ) \\ \prod _ { \lambda \, \dots \, } \, \pi \, \quad \mu / \infty \, \lambda \, \colon = \lambda _ { \lambda \, \dots \, } \dots \, \dots \, \lambda \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \dots \, \lambda \, \dots \, \dots \, \lambda \, \dots \, \dots \, \dots \, \lambda \, \dots \, \dots \, \dots \, \lambda \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \delta \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dists \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dts \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dclots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \ dclots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \,\dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dclots \, \dclots \, \dclots \, \dclots \, \dclots \, \dclots \, \dclots \, \dclots \dclots \dclots \dclots \dclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclipclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclslclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclcl|clclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclcldclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclclcl|\dots}$$

where E 0 = N (0) Ad f is determined by the electron density of states N (0) in S-layer and by the area A of the cross section of the junction, and the superscripts 0 and π label the quantities for the 0 -and π -state, respectively.

The critical temperatures T 0 ,π c of the transitions into 0 -or π -states are given by the expressions: 3,7

$$\int _ { 0 } ^ { \prime \prime } \ln \left ( \frac { T _ { c } ^ { 0, \pi } } { T _ { c 0 } } \right ) & = \Psi \left ( \frac { 1 } { 2 } \right ) - R e \left \{ \Psi \left ( \frac { 1 } { 2 } + \Omega _ { 0, \pi } \right ) \right \}, \quad ( 6 ) \\ \text{d} \quad \text{$\Phi$} \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi`} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \col. \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$} \colon \text{$\Phi$.}$$

where Ψ is the digamma function, and Ω 0 ,π is the depairing parameter:

$$\Omega _ { 0, \pi } = \frac { \varepsilon T _ { c 0 } } { 2 T _ { c } ^ { 0, \pi } } \begin{cases} k \tanh ( k s _ { f } ), & 0 - \text{phase} \\ k \coth ( k s _ { f } ), & \pi - \text{phase} \end{cases},$$

with ε = σ ξ /σ d ξ f 2 s s s f , s f = d / ξ f 2 f and k 2 = 2(i + 1 /τ h s ). The explicit expressions for the coefficients a 0 ,π and b 0 ,π read:

$$\begin{array} { c } \text{$d$} \\ \text{le} \\ \text{is} \\ \text{for} \\ \text{for} \\ \text{$d$} \\ \text{$e$} \\ \text{nt} \\ \text{$a$} \\ \text{$g$} \\ \end{array} \begin{array} { c } \text{and $bv,\pi$ read} \colon \\ \text{is} \\ \text{is} \\ \text{for} \\ \text{$d$} \\ \text{$e$} \\ \text{nt} \\ \text{na$} \\ \text{$g$} \\ \end{array} \begin{array} { c } a ^ { 0, \pi } = 1 - R e \left \{ \Omega _ { 0, \pi } \Psi ^ { ( 1 ) } ( 1 / 2 + \Omega _ { 0, \pi } ) \right \} \,, \\ b ^ { 0, \pi } = \frac { - 1 } { ( 4 \pi T _ { c } ^ { 0, \pi } ) ^ { 2 } } \, R e \left \{ \Psi ^ { ( 2 ) } ( 1 / 2 + \Omega _ { 0, \pi } ) - \frac { \Omega _ { 0, \pi } } { 6 k ^ { 2 } } \times \\ \text{$e$} \\ \text{nt} \\ \text{na$} \\ \text{na$} \\ \text{$g$} \\ \end{array} \begin{array} { c } \left [ i \mp \frac { i + 4 / \tau _ { s } h } { \cosh \gamma \pm 1 } \left ( 1 \pm \frac { \gamma } { \sinh \gamma } \right ) \right ] \Psi ^ { ( 3 ) } ( 1 / 2 + \Omega _ { 0, \pi } ) \right \} \,, \\ \text{and} \\ \text{$a$} \\ \text{$g$} \\ \text{$g$} \\ \text{$g$} \\ \text{$g$} \\ \text{where $\Psi ^ { ( n ) } ( z ) = d ^ { n } \Psi ( z ) / d z ^ { n }, \, \text{and $\gamma = 2 k s _ { f }. } \end{array}$$

where Ψ ( n ) ( z ) = d n Ψ( ) z /dz n , and γ = 2 ks f .

The functional F 0 ,π ( T ) provides the complete description of the S/F/S trilayers near the critical temperature. The equilibrium energy of the system is:

$$\text{es.} \quad F ^ { 0, \pi } ( \Delta _ { 0, \pi } ) = - E _ { 0 } \left [ a ^ { 0, \pi } \left ( T _ { c } ^ { 0, \pi } - T \right ) / T _ { c } ^ { 0, \pi } \right ] ^ { 2 } / 2 b ^ { 0, \pi }. \ \ ( 7 )$$

FIG. 4: (Color online) Critical temperatures T 0 ,π c vs. the thickness of the F-layer, and existence of the 0 -π transition line: the figure reports the typical dependence of T 0 ,π c on Flayer thickness d f for the even mode (0 -phase, solid red line) and for odd mode ( π -phase, dashed blue line). The inset zooms in the shaded region of the ( T , d c f ) phase diagram, to reveal the existence of the temperature-induced 0 -π transition (dashed green vertical arrow). Symbols + show the 0 -π transition line T 0 ( d f ). Here we choose: d s = 2 ξ s ; σ /σ f s = 0 12; . ξ /ξ s f = 3 ( ε = 0 18), and . √ b π /b 0 /similarequal 1 25. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0665v1-figure-e36101bd-2f1b-4e7b-ac56-cf3af45766cd.jpg)

Thus, the first order transition between 0 -and π -states occurs at F 0 (∆ ) = 0 F π (∆ ), thus determining the tranπ sition line T 0 as:

$$\frac { T _ { c } ^ { 0 } - T _ { 0 } } { T _ { c } ^ { \pi } - T _ { 0 } } = \frac { a ^ { \pi } T _ { c } ^ { 0 } } { a ^ { 0 } T _ { c } ^ { \pi } } \sqrt { \frac { b ^ { 0 } } { b ^ { \pi } } } \,. \quad \quad ( 8 )$$

The crossing of the curves T 0 c ( d f ) and T π c ( d f ) occurs at a value d ∗ f : for d f &gt; d ∗ f it is T π c &gt; T 0 c . At d f &lt; d ∗ f , but at thicknesses near d ∗ f , the decrease of the temperature determines first the transition from the normal to the superconducting 0 -state and then, with the further decrease of the temperature, the transition from the 0to the π -state (see inset of Fig. 4).

Thus, we have calculated the transition temperature 3,7 in the 0 -and π -states, as a function of the ratio d /ξ f f , T 0 ,π c ( d /ξ f f ), and we have found and calculated the transition line between the 0 -and the π -states, T 0 ( d /ξ f f ). The typical T 0 ,π c ( d f ) dependence is presented in Fig. 4. In this diagram, one must stress that near the critical thickness, where T 0 c = T π c , the 0 -π transition line T 0 emerges when F 0 (∆ ) = 0 F π (∆ ). π Only when d f is below, but close to, the critical thickness, the decrease of the temperature determines first the transition from the normal to the superconducting 0-state and then, with the further decrease of the temperature, the transition from the 0-to the π -state, as shown in the inset of Fig. 4.

To complete the picture of the temperature induced 0 -π transition, and to find the quantity directly observed in experiments, we calculate the superconducting order

/s32

FIG. 5: (Color online) Theoretical temperature dependence of the Ginzburg-Landau parameters of the S/F/S structure. (a) Schematic temperature dependence of the gap ∆ 2 0 ,π (upper scale, red lines) and the Ginzburg-Landau energy F GL (∆ 2 0 ,π ) (lower scale, blue lines): T π c &lt; T 0 c and ( a π ) 2 /b π &gt; a ( 0 ) 2 /b 0 ; (b) schematic temperature dependence of the penetration depth λ ∼ 1 ∆ / 0 ,π ; (c) dependence of the superconducting gap jump ∆ ( 2 0 T 0 ) / ∆ ( 2 π T 0 ) = ( b π /b 0 ) 1 / 2 on the F-layer thickness d f /lessorsimilar d ∗ f . The parameters are the same as in Fig. 4.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0665v1-figure-b4b298ac-ad2e-45a3-a248-e61c2f6c1ba4.jpg)

parameter as a function of the (lowering) temperature. In fact, at the transition, see Fig. 5a, the superconducting order parameter jumps from ∆ 0 to ∆ π according to:

$$\Delta _ { \pi } ^ { 2 } ( T _ { 0 } ) = \frac { a ^ { \pi } } { b ^ { \pi } } \left ( \frac { T _ { c } ^ { \pi } - T _ { 0 } } { T _ { c } ^ { 0 } } \right ) = \Delta _ { 0 } ^ { 2 } ( T _ { 0 } ) \sqrt { \frac { b ^ { 0 } } { b ^ { \pi } } }. \quad ( 9 )$$

It is an essential result that, since b 0 /b π &lt; 1 in the whole range of reasonable values for the various parameters (see Fig. 5c), at the transition the order parameter in the π -state has a smaller value than in the 0 -state, whence the reentrant behaviour of the superfluid density and the upward jump of the London penetration depth 1 /λ 2 ∼ ∆ , see Fig. 5b. So our theoretical description completely 2 recover the experimental findings of Section III.

## V. CONCLUSION

In summary, we have shown for the first time the thermodynamic nature of the 0 -π transition, by monitoring the temperature dependence of the order parameter close to the critical thickness condition. The experimental demonstration relied on the observation, by two different experiments, of a reentrant jump in the Cooper pair density in the Meissner and in the mixed states. The observation was performed in an S/F/S heterostructure with the F-layer close to the critical thickness, and with thin S-layers. An accurate theoretical treatment allowed us to find and calculate the 0 -π transition line, as well as the existence of the reentrant jump in the superfluid

density at the temperature-induced 0 -π transition. We have concluded that the observed jump was related to the first order transition from 0- to π -state on cooling.

## Acknowledgments

The authors thank A.S. Mel'nikov and C. Meneghini for stimulating discussions, and S. Sarti and R. Loria

- 1 A. I. Buzdin, L. N. Bulaevskii, and S. V. Panyukov, Pis'ma Zh. Eksp. Teor. Fiz. 35 , 147 (1982) [JETP Lett. 35 , 178 (1982)].
- 2 A. I. Buzdin and M. Yu. Kupriyanov, Pis'ma Zh. Eksp. Teor. Fiz. 53 , 308 (1991) [JETP Lett. 53 , 321 (1991)].
- 3 A. I. Buzdin, Rev. Mod. Phys. 77 , 935 (2005).
- 4 A. A. Golubov, M. Yu. Kupriyanov, and E. Ilichev, Rev. Mod. Phys. 76 , 411 (2004).
- 5 S. M. Frolov, M. J. A. Stoutimore, T. A. Crane, D. J. Van Harlingen, V. A. Oboznov, V. V. Ryazanov, A. Ruosi, C. Granata, and M. Russo, Nature Phys. 4 , 32 (2008).
- 6 V. V. Ryazanov, V. A. Oboznov, A. Yu. Rusanov, A. V. Veretennikov, A. A. Golubov, and J. Aarts, Phys. Rev. Lett. 86 , 2427 (2001).
- 7 V. A. Oboznov, V. V. Bolginov, A. K. Feofanov, V. V. Ryazanov, and A. I. Buzdin, Phys. Rev. Lett. 96 , 197003 (2006).
- 8 T. Kontos, M. Aprili, J. Lesueur, F. Genˆt, B. Stephanidis, e and R. Boursier, Phys. Rev. Lett. 89 , 137007 (2002).
- 9 T. R. Lemberger, I. Hetel, A. J. Hauser, and F. Y. Yang, J. Appl. Phys. 103 , 07C701 (2008).
- 10 N. Pompeo, K. Torokhtii, E. Silva, J. Supercond. Nov. Magn., accepted for publication (2014).
- 11 H. Sellier, C. Baraduc, F. Lefloch, and R. Calemczuk, Phys. Rev. B 68 , 054531 (2003).
- 12 C. Cirillo, E. A. Ilyina, and C. Attanasio, Supercond. Sci. Technol. 24 , 024017 (2011).
- 13 C. Cirillo, S. L. Prischepa, M. Salvato, C. Attanasio, M. Hesselberth, and J. Aarts, Phys. Rev. B 72 , 144511 (2005).
- 14 E. Silva, N. Pompeo, and S. Sarti, Supercond. Sci. Technol. 24 , 024018 (2011).
- 15 K. Torokhtii, C. Attanasio, C. Cirillo, E. A. Ilyina, N. Pompeo, S. Sarti, and E. Silva, Physica C 479 , 140 (2012).

for their help in the microwave measurements. This work was supported, in part, by European NanoSC COST Action MP1201, by French ANR grant 'MASH', by the RFBR (grant n.13-02-97126), by the program 'Quantum Mesoscopic and Disordered Structures' of the RAS, and by Ministry of Education and Science of RF and Lobachevsky State University (agreement 02.B.49.21.0003).

- 16 C. Meneghini and T. Neisius, private communication.
- 17 N. Pompeo, K. Torokhtii, C. Meneghini, S. Mobilio, R. Loria, C. Cirillo, E. A. Ilyina, C. Attanasio, S. Sarti, and E. Silva, J. Supercond. Nov. Magn. 26 , 1939 (2012).
- 18 N. Klein, U. Dahne, U. Poppe, N. Tellmann, K. Urban, S. Orbach, S. Hensen, G. Muller, and H. Piel, J. Supercond. 5 , 195 (1992).
- 19 N. Pompeo, K. Torokhtii, E. Silva, Meas. Sci. Rev. 14 , 164 (2014).
- 20 M. Tinkham, Introduction to Superconductivity , Dover Publications, 2004.
- 21 M. Golosovsky, M. Tsindlekht, and D. Davidov, Supercond. Sci. Technol. 9 , 1 (1996).
- 22 N. Pompeo and E. Silva, Phys. Rev. B 78 , 094503 (2008).
- 23 J. I. Gittleman and B. Rosemblum, Phys. Rev. Lett. 16 , 734 (1966).
- 24 M. W. Coffey and J. R. Clem, Phys. Rev. Lett. 67 , 386 (1991).
- 25 E. H. Brandt, Phys. Rev. Lett. 67 , 2219 (1991).
- 26 This is a crude estimate based on core pinning. Magnetic pinning mechanisms are likely to appear in SFS trilayers. However, we do not expect a strong effect on the temerature dependence in a reduced temperature range, such as the one in which we are interested in.
- 27 K. D. Usadel, Phys. Rev. Lett. 25 , 507 (1970).
- 28 V. N. Kushnir, S. L. Prischepa, J. Aarts, C. Bell, C. Cirillo, and C. Attanasio, Eur. Phys. J. B 80 , 445 (2011).
- 29 M. Faur´ e, A. I. Buzdin, A. A. Golubov, and M. Yu. Kupriyanov, Phys. Rev. B 73 , 064505 (2006).
- 30 M. Yu. Kupriyanov and V. F. Lukichev, Zh. Eksp. Teor. Fiz. 94 , 139 (1988) [ Sov. Phys. JETP 149 , 1163 (1988)].
- 31 A. I. Buzdin, Phys. Rev. B 72 , 100501 (2005).
N. Pompeo K. Torokhtii C. Cirillo A. V. Samokhvalov E. A. Ilyina C. Attanasio A. I. Buzdin E. Silva
Thermodynamic nature of the 0-PI quantum transition in superconductor-ferromagnet-superconductor trilayers
In structures made up of alternating superconducting and ferromagnet layers
(S/F/S heterostructures), it is known that the macroscopic quantum wavefunction
of the ground state changes its phase difference across the F--layer from 0 to
PI under certain temperature and geometrical conditions, whence the name "0-PI"
for this crossover. We present here a joint experimental and theoretical
demonstration that the "0-PI" is a true thermodynamic phase transition:
microwave measurements of the temperature dependence of the London penetration
depth in Nb/Pd_0.84Ni_0.16/Nb trilayers reveal a sudden, unusual decrease of
the density of the superconducting condensate (square modulus of the
macroscopic quantum wavefunction) with decreasing temperature, which is
predicted by the theory here developed as a transition from the 0- to the
PI-state. Our result for the jump of the amplitude of the order parameter is
the first thermodynamic manifestation of such temperature-driven quantum
transition.