Nonmonotonic Evolution of the Superconducting Transition Temperature and Robust Multigap Extended s-wave + s-wave Pairing in Zn-Substituted FeSe Single Crystals

This study demonstrates that Zn substitution in FeSe single crystals induces a nonmonotonic evolution of the superconducting transition temperature driven by enhanced scattering rather than simple pair breaking, while preserving robust multigap extended s-wave + s-wave pairing that imposes stringent constraints on candidate pairing mechanisms.

The Gist

Imagine a microscopic city called FeSe (Iron Selenide). In this city, electrons (the citizens) usually move around chaotically. But when the city gets cold enough, these electrons decide to pair up and dance in perfect unison. This synchronized dancing is called superconductivity, a state where electricity flows with absolutely zero resistance, like a car driving on a frictionless highway.

The big mystery scientists have been trying to solve is: What makes the electrons dance this way, and how do we keep the dance going?

This paper is about a team of scientists who decided to play "musical chairs" with the citizens of this city. They replaced some of the Iron (Fe) citizens with Zinc (Zn) citizens. Zinc is a "non-magnetic" neighbor, meaning it doesn't have a strong personality (magnetic field) that usually disrupts the dance.

Here is what they found, explained simply:

1. The Bumpy Ride (The Non-Monotonic Temperature)

Usually, if you add a new type of citizen to a city, you expect things to get worse gradually. If you add too many strangers, the dance floor gets crowded, and the dancing stops.

However, when the scientists added a little bit of Zinc, the "dance temperature" (the point where superconductivity starts) dropped. But then, as they added more Zinc, the temperature went back up! Finally, if they added even more Zinc, the temperature dropped again.

The Analogy: Imagine a band playing music.

• Too few strangers: The band plays well.
• A few strangers: The rhythm gets a bit off (temperature drops).
• More strangers: Surprisingly, the band finds a new, better rhythm and plays even louder (temperature rises).
• Too many strangers: The crowd gets too chaotic, and the music stops (temperature drops again).

This "up-and-down" behavior told the scientists that the story isn't just about "messing up" the dance. Something more complex is happening with how the electrons interact.

2. The Two-Step Dance (Multigap Superconductivity)

To understand how the electrons are dancing, the scientists looked at the heat the material gave off (specific heat). They found that the electrons weren't just doing one simple dance move.

The Analogy: Think of a ballroom with two different groups of dancers.

• Group A is doing a smooth, steady waltz (an isotropic s-wave gap).
• Group B is doing a more complex, twisting dance that changes speed depending on the direction (an anisotropic extended s-wave gap).

Most other theories suggested the dancers were doing a "figure-eight" (d-wave) or a simple single waltz. But the data proved that both groups are dancing at the same time. This is called multigap superconductivity. The Zinc didn't stop them from dancing together; it just tweaked the rhythm slightly.

3. The "Ghost" Effect (Why Zinc Didn't Break the Dance)

In physics, there's a rule called the Anderson Theorem. It says that if you have a simple dance (like a standard waltz), adding non-magnetic strangers (Zinc) shouldn't hurt the dance at all. But if the dance is complex and relies on the dancers being on "opposite sides" of the room (sign-changing), adding strangers breaks the connection immediately.

The Surprise: The Zinc did change things, but it didn't destroy the superconductivity. The two groups of dancers kept their relative sizes the same, even as more Zinc arrived.

The Analogy: Imagine a dance where two groups hold hands. If you put a stranger between them, they usually let go. But here, the Zinc acted like a ghost. It passed through the crowd without breaking the hand-holding between the two main groups. This suggests the dancers aren't holding hands in a way that is easily broken by strangers. It implies the "dance move" (the pairing mechanism) is likely robust and doesn't rely on the fragile "opposite signs" that usually get destroyed by impurities.

4. The Big Takeaway

The scientists concluded that FeSe is a very special, resilient material.

• It doesn't just have one way of conducting electricity; it has multiple lanes (multiband) working together.
• The "dance" is a mix of simple and complex moves.
• Even when you mess with the city by adding Zinc, the superconductivity is tough enough to survive and even adapt.

Why does this matter?
Understanding how these electrons dance helps scientists design better materials for future technology, like ultra-fast computers or lossless power grids. This paper tells us that the "secret sauce" of FeSe isn't a fragile, single trick, but a robust, multi-layered system that can handle a bit of chaos without falling apart.

In short: The scientists added a new ingredient to a superconductor, expecting it to ruin the party. Instead, the party got weird, then better, then worse, proving that the electrons are dancing in a complex, resilient, two-step rhythm that is much harder to break than anyone thought.

Technical Summary

1. Problem Statement

The superconducting pairing mechanism in iron-based superconductors, specifically FeSe, remains a subject of intense debate. While FeSe has a simple crystal structure, its pairing symmetry (e.g., $s_{\pm}$, $d$-wave, or sign-preserving $s$-wave) and gap structure (nodal vs. nodeless) are unresolved.

• The Gap in Knowledge: Previous substitution studies focused heavily on the Se site (e.g., Te, S) or magnetic Fe-site dopants (e.g., Co, which may carry magnetic moments). The effect of nonmagnetic impurity substitution at the Fe site in high-quality single crystals is less explored.
• Theoretical Context: According to Anderson's theorem, nonmagnetic impurities should not suppress $T_c$ in conventional $s$-wave superconductors. However, in sign-reversing states (like $d$-wave or $s_{\pm}$), nonmagnetic impurities act as strong pair breakers via interband scattering, leading to a monotonic suppression of $T_c$.
• Objective: This study aims to probe the pairing mechanism of FeSe by systematically substituting Fe with nonmagnetic Zn in high-quality single crystals ($Fe_{1-x}Zn_xSe$) and analyzing the evolution of $T_c$, transport properties, and thermodynamic gap structure.

2. Methodology

The researchers synthesized and characterized a series of high-quality single crystals with Zn concentrations ranging from $x = 0$ to $x = 0.023$.

• Synthesis: High-quality single crystals were grown using the flux method with $AlCl_3$ and $KCl$ as fluxes in a two-zone tube furnace. This method ensured millimeter-sized, plate-like crystals with high crystallinity.
• Structural Characterization:
  • XRD: Confirmed high crystallinity and pure phase (only $(00l)$ reflections) with no impurity peaks.
  • EDS: Verified homogeneous distribution of Zn and confirmed the substitution of Zn at Fe sites (no interstitial occupation).
• Physical Property Measurements:
  • Magnetization: Zero-field cooled (ZFC) and field-cooled (FC) measurements to determine $T_c$ onset and bulk superconductivity. Hysteresis loops ($M-H$) were used to calculate critical current density ($J_c$) and lower critical field ($H_{c1}$).
  • Transport: Resistivity measurements ($\rho(T)$) under various magnetic fields to determine upper critical fields ($H_{c2}$) and magnetoresistance (MR).
  • Specific Heat: Low-temperature specific heat ($C_p$) measurements were performed to extract the electronic specific heat ($C_e$) and analyze the superconducting gap structure.
• Data Analysis:
  • $H_{c2}$ was fitted using the Werthamer-Helfand-Hohenberg (WHH) model.
  • Specific heat data was analyzed using single-band isotropic $s$-wave, $d$-wave, and two-gap models (combining isotropic $s$-wave and anisotropic extended $s$-wave).

3. Key Contributions

1. Discovery of Nonmonotonic $T_c$ Evolution: The study reveals that $T_c$ does not decrease monotonically with Zn doping. Instead, it initially drops, recovers to a local maximum at $x \approx 0.014$, and then drops again. This behavior contradicts the simple pair-breaking model expected for sign-changing pairing states under nonmagnetic impurity scattering.
2. Confirmation of Robust Multigap Superconductivity: Through specific heat analysis, the authors definitively rule out single-gap and $d$-wave models. They establish that the superconducting state is best described by a two-gap scenario consisting of an isotropic $s$-wave gap and an anisotropic extended $s$-wave gap.
3. Evidence for Weak Interband Scattering: The relative weights of the two gap components remain nearly constant across the entire doping range. This suggests that Zn substitution induces only weak interband scattering, preserving the multiband nature of the superconductivity.
4. Constraint on Pairing Symmetry: The resilience of superconductivity against nonmagnetic impurities, combined with the specific gap structure, strongly disfavors a pure sign-changing ($s_{\pm}$ or $d$-wave) state in the strong scattering limit, pointing instead toward a predominantly sign-preserving pairing state with significant anisotropy.

4. Key Results

• Structural Properties:
  • Zn substitutes Fe sites, causing a non-monotonic change in the $c$-axis lattice parameter, which remains close to 5.50 Å.
  • Crystals maintain high purity and bulk superconductivity up to $x = 0.023$.
• Superconducting Transition ($T_c$):
  • Parent FeSe ($x=0$): $T_c \approx 9.8$ K.
  • Low Doping ($x=0.006$): $T_c$ drops to $\approx 7.86$ K.
  • Intermediate Doping ($x=0.014$): $T_c$ recovers to $\approx 9.05$ K.
  • High Doping ($x=0.023$): $T_c$ suppresses again to $\approx 7.49$ K.
  • This "dome-like" or non-monotonic behavior suggests a competition between impurity scattering and modifications to the Fermi surface pockets.
• Transport Properties:
  • Resistivity: The structural transition temperature ($T^*$) tracks $T_c$. The Residual Resistivity Ratio (RRR) decreases with Zn doping, indicating increased disorder.
  • Magnetoresistance: Shows a crossover from quadratic ($B^2$) to linear dependence at high fields. The effective mobility decreases with Zn doping, confirming enhanced scattering.
  • Upper Critical Field ($H_{c2}$): Well-described by the orbital-limited WHH model, indicating negligible spin-paramagnetic effects.
• Specific Heat and Gap Structure:
  • Gap Symmetry: Single-band isotropic $s$-wave, $d$-wave, and single-gap extended $s$-wave models failed to fit the data.
  • Two-Gap Model: The data is perfectly fitted by a weighted sum of an isotropic $s$-wave gap ($\Delta_s$) and an anisotropic extended $s$-wave gap ($\Delta_{es}$).
  • Parameters: For $x=0$, the gaps are $2\Delta_s/k_BT_c \approx 2.0$ and $2\Delta_{es}/k_BT_c \approx 4.9$ with weights of 45% and 55%, respectively.
  • Stability: The relative weights of these two components remain nearly invariant ($\sim$40-45% vs. 55-60%) despite Zn doping, indicating that the multiband character is robust.
  • Coupling Strength: The specific heat jump ratio $\Delta C/\gamma_n T_c$ ranges from 1.30 to 1.76, exceeding the weak-coupling BCS limit (1.43), suggesting intermediate coupling strength.

5. Significance

This work provides critical constraints on the pairing mechanism in FeSe:

• Refutation of Simple Pair-Breaking: The non-monotonic $T_c$ evolution under nonmagnetic doping rules out a simple $s_{\pm}$ or $d$-wave scenario where nonmagnetic impurities act as strong pair breakers.
• Support for Sign-Preserving Anisotropic $s$-wave: The robustness of the multigap state and the lack of rapid $T_c$ suppression suggest a sign-preserving order parameter (likely $s_{++}$) with strong anisotropy. The system likely exists in a crossover regime where multiple pairing channels compete.
• Multiband Physics: The results highlight that the electronic structure of FeSe is highly sensitive to impurity-induced scattering, which redistributes contributions among Fermi surface pockets, leading to complex $T_c$ behavior rather than simple suppression.
• Material Design: The demonstration that FeSe can host robust multigap superconductivity even with significant nonmagnetic disorder underscores the stability of its superconducting ground state, offering insights for engineering iron-based superconductors with higher $T_c$ or specific gap symmetries.