Field-induced magnetic instability within a superconducting condensate

See allHide authors and affiliations

Science Advances  19 May 2017:
Vol. 3, no. 5, e1602055
DOI: 10.1126/sciadv.1602055


The application of magnetic fields, chemical substitution, or hydrostatic pressure to strongly correlated electron materials can stabilize electronic phases with different organizational principles. We present evidence for a field-induced quantum phase transition, in superconducting Nd0.05Ce0.95CoIn5, that separates two antiferromagnetic phases with identical magnetic symmetry. At zero field, we find a spin-density wave that is suppressed at the critical field μ0H* = 8 T. For H > H*, a spin-density phase emerges and shares many properties with the Q phase in CeCoIn5. These results suggest that the magnetic instability is not magnetically driven, and we propose that it is driven by a modification of superconducting condensate at H*.

  • Condensed Matter Physics
  • quantum phase transitions
  • superconductivity
  • magnetism
  • neutron scattering
  • strongly correlated electron systems
  • low temperature physics
  • diffraction
  • scaling


Although quantum mechanics was developed to treat nonrelativistic interacting particles, it has also been widely used to describe collective phenomena in condensed matter physics. A macroscopic quantity of interacting particles (~1023) can give rise to novel quantum phases with emergent magnetic and transport properties that are currently hotly debated among theorists and experimentalists alike. Particularly exciting are quantum liquids, including superconducting condensates, that feature coherence phenomena over macroscopic length scales. Pressure, magnetic field, or chemical substitution can induce quantum phase transitions that originate from fluctuations that are governed by Heisenberg’s uncertainty principle. These collective quantum fluctuations can trigger the collapse or emergence of an order parameter, and they can even dominate the finite temperature properties. Quantum fluctuations generally separate phases with different organizing principles (1, 2), leading to magnetic phases with different symmetry as the strength of the fluctuations evolve. This has been demonstrated, for instance, by quantum magnetic insulators (3).

Strong electronic correlations and quantum fluctuations are thought to be at the origin of unconventional superconductivity in many materials. Examples include the layered copper oxides (4), the iron-based pnictides, and the heavy-fermion (HF) superconductors (5). CeCoIn5 is an exceptional material because it features spin-density wave (SDW) order that emerges from the superconducting condensate (6). The material is an ultraclean, ambient-pressure, dx2−y2 superconductor (Tc = 2.3 K) (7) that displays a quasi–two-dimensional Fermi surface (8, 9). Superconductivity is Pauli-limited (10) and coexists at high fields with a SDW that emerges with an ordered moment μ(11 T) = 0.15(5)μB (6, 1114). Magnetism exists only inside the superconducting state and collapses in a first-order transition at the superconducting upper critical field. The ordering wave vector, Q = (0.45, 0.45, 0.5), is pinned along the nodes of the superconducting gap (11) and may result from the condensation of superconducting quasi-particle excitations with increasing magnetic field (15, 16). All this provides evidence for a coupling between superconductivity and magnetic order (11, 12). However, the microscopic interpretation and the nature of the order parameters of this so-called Q phase are still debated (6, 1523).

The Q phase is highly sensitive upon doping and has been suppressed in all previously reported cases (2427). However, recent studies on Nd0.05Ce0.95CoIn5 (28, 29) revealed SDW magnetic order with a similar ordering wave vector as the Q phase for T < TN = 0.8 K at zero field (29). The family of Nd-substituted CeCoIn5, Nd1−xCexCoIn5, displays a rich phase diagram featuring a HF ground state (x > 0.5), superconductivity (x > 0.78), and magnetism (x < 1) (28). The substitution of Nd for Ce in CeCoIn5 introduces localized moments, reduces the hybridization of 4f-bands with the Fermi surface, and weakens the coherence of these heavy bands. This causes the destruction of itinerant magnetism and the emergence of localized moment magnetism at x < 0.5. Small Nd concentrations gradually weaken the superconducting pairing strength, as observed in the reduction in the specific-heat jump (28). The implantation of a small amount of Nd impurities in CeCoIn5 has a similar effect on superconductivity as the substitution with Kondo holes (29) and acts predominantly as random disorder (see the Supplementary Materials). However, the nature of the magnetic order for small Nd substitutions, such as in Nd0.05Ce0.95CoIn5, and its relationship to the SDW in the Q phase is currently not understood.


Evidence for two distinct magnetic phases inside the superconducting condensate

We studied the magnetic order in Nd0.05Ce0.95CoIn5 as a function of field for H||[1 Embedded Image 0]. Figure 1 shows neutron diffraction data at several magnetic fields and T = 40 mK along the (q, q, 0.5) reciprocal wave vector given in reciprocal lattice units. A well-defined magnetic Bragg peak is observed at zero field, consistent with previous measurements (29). With increasing field, the peak first increases in intensity but then completely vanishes around μ0H* ≈ 8 T (Fig. 1D). The key observation is that the peak reappears at higher fields, as shown in Fig. 1 (E and F). The magnetic Bragg peaks were fitted with a Gaussian line shape, showing that the width of the Bragg peaks is resolution-limited at all fields. The peak intensity, Ip, is shown in Fig. 2 as a function of field and provides direct evidence for two distinct magnetic phases: a low-field phase (SDW phase) for μ0H < 8 T and a high-field phase for 8 T < μ0H < 11 T, which we call Q phase because it shares many properties with the Q phase of CeCoIn5. The Q phase in Nd0.05Ce0.95CoIn5 appears for fields larger than 8 T, and its magnetic Bragg peak increases in intensity with increasing field. At μ0H = 11.0(2) T, it collapses in a sharp transition together with superconductivity. This shows that magnetic order is coupled to superconductivity at high fields—similar to the Q phase in undoped CeCoIn5.

Fig. 1 Magnetic order at various magnetic fields.

The scans display diffracted neutron intensity in counts per 30 minutes (cts/30 min) along the tetragonal plane (q, q, 0.5) in reciprocal lattice units (r.l.u.). The reflections correspond to magnetic Bragg peaks at μ0H = 0, 4, 7, 8.5, 9.5, and 10.5 T for (A) to (F), respectively. Background was measured at μ0H = 12 T. All scans were taken at T = 40 mK for H||[1 Embedded Image 0]. The integrated intensity grows for increasing fields up to 4 T, decreases for further field increments, and vanishes around μ0H* ≈ 8 T. For higher fields, magnetic order reappears, increases in intensity, and collapses at μ0H = 11 T.

Fig. 2 Discovery of a magnetic instability inside the superconducting phase.

The position-optimized peak intensity (red) shows two distinct magnetic phases inside the superconducting phase. Both phases vanish at μ0H* = 8.0(2) T (solid and dashed lines are guides to the eye). Because of diffuse scattering, the magnetic intensity remains above the background level (blue) in the vicinity of H*. The high-field phase collapses at the onset of the superconducting phase (resistivity along H||[1 0 0] is shown in green), μ0Hc2 =11.0(2) T, providing evidence for a coupling between superconductivity and magnetic order.

The HT phase diagram of the magnetic phases in Nd0.05Ce0.95CoIn5 is shown in Fig. 3 and shows additional similarities of the high-field phases in pure and doped CeCoIn5: The Q phase in both compounds features very similar HT phase boundaries, and magnetic order is suppressed with increasing temperature. In a further similarity, the field-induced transition from the Q phase to the normal state in Nd0.05Ce0.95CoIn5 remains first-order at Hc2(T) for T < 300 mK. Our results clearly establish that the Q phase in CeCoIn5 is stable under 5% Nd doping of the Ce site.

Fig. 3 A magnetic instability separates two magnetic phases with identical symmetry inside the superconducting condensate.

HT phase diagram of Nd0.05Ce0.95CoIn5 for H||[1 Embedded Image 0]. The color scale represents the background-subtracted magnetic intensity. Orange circles display the magnetic phase boundaries obtained from neutron scattering results (orange dashed lines are guides to the eye). White diamonds show the superconducting phase boundary measured via electrical resistivity and heat capacity (the gray dashed line shows the superconducting phase boundary normalized from CeCoIn5). The color map reveals a field-induced instability at μ0H* = 8.0(2) T. The high-field phase, Q phase, increases linearly in intensity for increasing fields and collapses at the superconducting upper critical field μ0Hc2.

Field-induced magnetic instability separates two SDWs with identical symmetry

As shown in Fig. 1D, a diffuse signal replaces the magnetic Bragg peak around μ0H* ≈ 8 T. This demonstrates that there is a magnetic instability around H* where magnetic order is absent and separates the SDW phase and the Q phase. Such diffuse scattering typically arises from short-range static or dynamic magnetic correlations. The field dependence of the magnetic Bragg peak intensity is linear as a function of field in the proximity of H*, as shown in Fig. 2. Because the magnetic neutron intensity is proportional to the square of the magnetic moment, this means that the ordered magnetic moment M behaves as M = M0|H/H* − 1|β, where β ≈ 0.5 and M0 = 0.19(1)μB for H < H* and M0 = 0.28(2)μB for H > H*, respectively with μ0H* = 8.0(2) T. A mean-field exponent β = 0.5 is also observed over a large field range as in the Q phase of the undoped compound (6) and is expected for quantum phase transitions (30). Although we do not have direct evidence of quantum fluctuations, we point out that the HT phase diagram suggests the presence of a quantum critical region around H* that expands with increasing temperature.

The field dependence of the in-plane component of the wave vector transfer, q, is shown in Fig. 4. The data display a linear decrease in the wave number from q ≈ 0.448 at zero field to q ≈ 0.441 at 10.5 T and appear to be unaffected by the disappearance of magnetic order around μ0H* ≈ 8 T, which is probably related to robust features of the Fermi surface topology. The weak field dependence of the wave vector transfer is in contrast to the Q phase in CeCoIn5 where the wave vector is independent of the magnetic field (11, 12). The symmetry of the magnetic order is identical for the low-field SDW phase and high-field Q phase: At both zero field and μ0H = 10.5 T, an identical amplitude-modulated SDW with field-dependent magnetic moments aligned along the c axis (see Fig. 4, inset) and magnetic moments of μ(0 T) = 0.10(5)μB and μ(10.5 T) = 0.15(5)μB are found, respectively. For increasing fields, the magnetic peak intensity increases and reaches a broad maximum at μ0H ≈ 4 T before it gradually weakens toward H*. At μ0H = 4 T, the magnetic structure is also an amplitude-modulated SDW with a somewhat larger magnetic moment. This initial increase of intensity may result from a nonuniform magnetic domain population under a magnetic field. However, the magnetic structure determination clearly demonstrates that the symmetry of the magnetic structure does not change within the SDW phase.

Fig. 4 Field-dependent wave vector and magnetic structure.

The in-plane component of the propagation vector decreases linearly with increasing fields. Inset: Magnetic structure refined at μ0H = 0 and 10.5 T. An amplitude-modulated structure with a moment configuration perpendicular to the basal plane is found, similar to the Q phase of CeCoIn5.


The identical magnetic symmetry of the ground states in the SDW phase and the Q phase demonstrates that the magnetic instability at μ0H* = 8 T is not purely driven by magnetic fluctuations and that it must emerge from other types of fluctuations. Because the Kondo breakdown in Nd1−xCexCoIn5 occurs around x = 0.5 (28), we exclude the possibility that the observed magnetic instability for x = 0.95 is driven by charge valence fluctuations. On the other hand, it is well known that Pauli-limited superconductors feature instabilities toward complex superconducting phases at high fields, including modulated superconductivity and phases of coexisting magnetic and superconducting order (20, 31, 32). Our results may be explained by a modification of the superconducting condensate that induces a change in the magnetic properties. This is consistent with a recently proposed scenario in undoped CeCoIn5, where SDW order only exists for fields higher than μ0H = 9.8 T and is induced by the emergence of p-wave superconductivity in the d-wave condensate (6, 12, 33, 34). The scenario of novel superconductivity for high fields in the tetragonal plane of CeCoIn5 is expected from a number of microscopic theories (1622). We note that the HT phase diagram of Nd0.05Ce0.95CoIn5 shows some qualitative similarities with the xT phase diagram of La-based cuprate La2−xBaxCuO4, where two superconducting phases exist in the presence of charge and stripe order and are separated by one or more critical points (35). For Nd0.05Ce0.95CoIn5, it needs to be determined whether the magnetic instability is caused by a single or several quantum phase transitions. In this respect, it will be important to study the HT phase diagrams of Nd1−xCexCoIn5 for x > 0.95.

In summary, we report the discovery of a magnetic instability region inside the superconducting phase of Nd0.05Ce0.95CoIn5 that separates two antiferromagnetic phases with identical symmetry. Whereas a low-field SDW is suppressed at μ0H* = 8 T, another SDW emerges from fields μ0H* > 8 T and shares many properties with the Q phase of CeCoIn5, such as the sudden collapse of the SDW at the upper critical field. Our experiment shows that the Q phase in CeCoIn5 is stable under 5% Nd doping of the Ce site. The identical magnetic symmetry of the two magnetically ordered phases suggests a quantum phase transition driven by a modification of the superconducting condensate.


Neutron diffraction measurements were carried out on the thermal-neutron lifting-counter two-axis spectrometer D23, the single-crystal four-circle diffractometer D10, the cold-neutron three-axis spectrometer IN12 at the Institut Laue-Langevin (Grenoble, France), and the cold-neutron triple-axis spectrometer RITA-II at the Swiss Spallation Neutron Source SINQ, Paul Scherrer Institut (Villigen, Switzerland). For all measurements, the same single crystal of mass m = 64 mg was used (see also the Supplementary Materials). The experiments required temperatures down to T = 40 mK and maximal fields of μ0H = 12 T, which were achieved using vertical-field magnets with dilution inserts. On D23, an incident neutron wavelength λ = 1.28 Å was used. The experiment on D10 was conducted with a neutron wavelength λ = 2.36 Å. For the experiment, we used the 3He single detector in a four-circle setup and optimized the signal-to-noise ratio via the analyzer option of the instrument. A clean wavelength λ = 4.83 Å was derived from a velocity selector and double-focusing pyrolitic graphite monochromator on IN12. The neutrons were collimated by α = 80′ before scattering at the sample. On RITA-II, we used a wavelength λ = 4.217 Å. The high-energy background was minimized with a collimator (α = 80′), followed by a pyrolytic graphite filter before and a beryllium filter after the sample. The signal was recorded from the central analyzer blade of the RITA-II nine-bladed multianalyzer, and the background in Fig. 2 was obtained from the two neighboring blades.

The magnetic structure at zero field was obtained from the data measured on D10. Six independent magnetic Bragg reflections were analyzed using the three possible irreducible representations. The amplitude-modulated SDW along the c axis with μ(0 T) = 0.10(5)μB results in Rf = 13.6%. Refinements in other representations yield Rf > 43% constantly. For the magnetic structure analysis at μ0H =10.5 T at RITA-II, six independent reflections were used. For the best refinement, in the amplitude-modulated representation with μ(0 T) = 0.15(5)μB, Rf equals 12.6%, whereas models in other representations yield Rf > 30%. In each model, equally populated domains, q = (q, ±q, 0.5), were assumed, and the inset of Fig. 4 shows one of the two domains. For the refinements, the FullProf suite was used (36). The prefactors M0 above and below the H* were obtained by a fit of the field-dependent ordered magnetic moment. This field dependence was recovered from refinements of the integrated intensity in the amplitude-modulated representation that were calculated from the peak intensities shown in Fig. 2. The superconducting upper critical field, Hc2(T), was measured by means of electrical resistivity and heat capacity at Commissariat à l’Energie Atomique et aux Energies Alternatives in Grenoble, France. These experiments were performed on a small part of the single crystal used for neutron diffraction and on small crystals of the same batch. Experimental data are shown in the Supplementary Materials.

The magnetic phase diagram in Fig. 3 was realized via neutron diffraction data measured at the peak-optimized position of the magnetic Bragg peak on the instruments RITA-II and IN12. The phase boundary was obtained from temperature scans at μ0H = 1, 2, 3, 4, 4.5, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, and 10 T, as well as field scans at T = 40, 150, 250, and 300 mK. The field scan at T = 40 mK is shown in Fig. 2, and the temperature scans at μ0H = 7 to 9.5 T are displayed in the Supplementary Materials. The color plot represents the monitor-normalized background-subtracted peak intensity. Points among the measured data points were interpolated from the nearest neighbor-measured peak intensities.


Supplementary material for this article is available at

Supplementary Materials and Methods

fig. S1. Macroscopic results.

fig. S2. Temperature dependence of magnetic intensity in the critical region.

Reference (37)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We acknowledge the Institut Laue-Langevin and the Paul Scherrer Institut for the allocated beam time on D23, D10, IN12, and RITA-II. We thank P. Fouilloux, P. Decarpentrie, B. Vettard, K. Mony, X. Tonon, B. Grenier, M. Zolliker, and M. Bartkowiak for technical assistance. Discussions with M. Sigrist and S. Gerber are acknowledged. Funding: We thank the Swiss National Science Foundation (grant nos. 200021_147071 and 200021_138018). This work was additionally supported by the Swiss State Secretariat for Education, Research, and Innovation through a Collaboration Research Group grant. The research also received funding from the European Community’s Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 290605 (PSI-FELLOW/COFUND). We also acknowledge the Danish neutron society community Danscatt. Author contributions: D.G.M., S.R., J.L.G., and M.K. planned and led the project. The neutron experiments were performed by D.G.M., S.R., J.L.G., E.R., C.N., J.O.B., B.O., and M.K. The resistivity and heat capacity measurements were performed by G.B., G.K., and D.A. The sample was grown by G.L., and the data were analyzed by D.G.M. The manuscript was written by D.G.M., S.R., J.L.G., and M.K. with the input of all coauthors. Competing interests: All authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors. Correspondence and requests for materials should be addressed either to M.K. or to S.R. (emails: michel.kenzelmann{at} and raymond{at}

Stay Connected to Science Advances

Navigate This Article