Extracting intrinsic superconducting properties in intercalated layered superconductors using an extended 2D Tinkham model

This study resolves the misclassification of certain intercalated layered superconductors as anisotropic 3D systems by developing an extended 2D Tinkham model that accounts for interlayer misalignment, thereby enabling the accurate extraction of intrinsic bulk 2D superconducting properties and BKT transitions in materials like (Li,Fe)OHFeSe and (CTA)0.5SnSe2.

The Gist

Imagine you are trying to understand how a stack of pancakes conducts electricity when it gets super cold. In the world of physics, these "pancakes" are layers of atoms in a special material called a superconductor. When these materials get cold enough, they lose all electrical resistance, allowing electricity to flow forever without energy loss.

Scientists have long been fascinated by 2D superconductivity—where the electricity flows easily within the flat layers but struggles to jump between them. However, there's a problem: sometimes, when scientists look at these materials, they get confused. They think the electricity is flowing in all directions (3D) when it's actually just flowing in flat layers (2D).

This paper solves that confusion by looking at a specific material called (Li,Fe)OHFeSe and introducing a new way to measure it.

The Problem: The "Wobbly Stack"

The researchers compared two versions of this material:

1. Single Crystals: Think of these as a stack of pancakes where the layers are slightly crooked, tilted, and misaligned. It's a bit of a messy stack.
2. Epitaxial Films: These are like a perfectly stacked, aligned tower of pancakes, where every layer is perfectly flat and parallel to the one below it.

When they tested the "messy" single crystals, the data looked like the electricity was flowing in all directions (3D). But when they tested the "perfect" films, the data clearly showed the electricity was stuck in the flat layers (2D).

The Mystery: Both samples actually are 2D superconductors (they both showed a specific signature called a "BKT transition," which is like a fingerprint proving they are 2D). So, why did the messy crystals look like 3D?

The Analogy: The Blurry Camera

The authors realized the issue wasn't the material itself, but how the "camera" was looking at it.

Imagine taking a photo of a sharp, pointed mountain peak (the perfect 2D behavior).

• If you use a steady camera (the perfect film), the photo shows a sharp, distinct peak.
• If you use a shaky, blurry camera (the misaligned crystal), the sharp peak looks rounded and fuzzy.

In the lab, the "shakiness" comes from the interlayer misalignment. Because the layers in the single crystals are tilted at slightly different angles, the measurement gets "blurred." This blur makes the sharp 2D peak look like a rounded 3D hill, tricking scientists into thinking the material is 3D.

The Solution: The "Extended Tinkham Model"

To fix this, the scientists invented a new mathematical tool called the Extended 2D Tinkham Model.

Think of this model as a photo editing filter.

• The standard model (the old way) assumes the camera is perfectly steady. It tries to fit the blurry photo to a sharp peak and fails, giving wrong answers.
• The Extended Model knows the camera is shaky. It takes the "blur" (measured by how tilted the layers are) and mathematically "un-blurs" the data.

By feeding the model the exact amount of "tilt" in the crystals, they could mathematically reconstruct what the perfect 2D behavior should have looked like.

The Results

1. It Works: When they applied this new "filter" to the messy single crystals, the data suddenly matched the perfect films. They realized the crystals were 2D all along; they just looked 3D because of the tilt.
2. It's Universal: They tested this on another material, (CTA)0.5SnSe2, which also had a messy stack. The old model failed to explain it, but the new "Extended Model" fit the data perfectly.

The Big Takeaway

The paper concludes that many layered superconductors made using "soft" chemical methods (which are great for adding special molecules but often result in messy stacks) have been misclassified as 3D.

The authors provide a simple, effective rule: If you have a layered superconductor with a messy stack, don't just look at the raw data. Use this new "Extended Model" to account for the tilt, and you will uncover the true, intrinsic 2D nature of the material.

This helps scientists stop misidentifying these materials and finally understand their true superconducting properties.

Technical Summary

Technical Summary: Extracting Intrinsic Superconducting Properties in Intercalated Layered Superconductors

Problem Statement
Bulk two-dimensional (2D) superconductivity is characterized by the interplay of symmetry breaking, nontrivial topology, and 2D phase fluctuations. In layered superconductors where the $c$-axis lattice parameter exceeds the superconducting coherence length ($c > \xi_{\perp}$), weak Josephson coupling between layers is expected to result in intrinsic 2D behavior. However, certain intercalated layered superconductors, such as (Li,Fe)OHFeSe and (CTA)$_{0.5}$SnSe$_2$, have been misclassified as anisotropic three-dimensional (3D) superconductors in previous literature. This discrepancy is particularly evident in (Li,Fe)OHFeSe, where some studies report anisotropic 3D behavior while others suggest quasi-2D superconductivity. The authors posit that these conflicting results stem from sample-dependent interlayer misalignments, which blur the angular dependence of the upper critical field ($H_{c2}(\theta)$), leading to incorrect dimensional classifications and underestimated anisotropy values.

Methodology
The study employs a comparative approach using two forms of (Li,Fe)OHFeSe: bulk single crystals and epitaxial films.

1. Sample Synthesis and Characterization: Single crystals were synthesized via a hydrothermal ion-exchange method, while epitaxial films were grown on LaAlO$_3$ substrates using matrix-assisted hydrothermal epitaxy (MHE). X-ray diffraction (XRD) rocking curves were used to quantify interlayer alignment, specifically measuring the full width at half maximum (FWHM) of the (006) diffraction peak.
2. Transport Measurements: Electrical transport properties were measured using a four-probe method on a Quantum Design PPMS-9 system. Key measurements included:
   • Angular dependence of the upper critical field ($H_{c2}(\theta)$) from $\theta = 0^\circ$ ($H \parallel c$) to $\theta = 90^\circ$ ($H \parallel ab$).
   • Current-voltage ($I-V$) characteristics to identify Berezinskii–Kosterlitz–Thouless (BKT) transitions.
   • Resistive transitions to determine $T_c$ and $T_{BKT}$.
3. Modeling: The authors developed an extended 2D Tinkham model. Unlike the standard 2D Tinkham model which assumes ideal crystalline alignment ($\Delta\omega = 0$), this extended model incorporates the experimentally measured FWHM values ($\Delta\omega$) as a Gaussian distribution. The measured resistance is modeled as a convolution of the ideal scaled resistance and the angular misalignment distribution:
   $$R_{measure}(H, \theta) = \int g(\omega) \cdot R\left(\frac{H}{H_{c2}(\theta - \omega)}\right) d\omega$$
   where $g(\omega)$ represents the Gaussian distribution of interlayer misalignments.

Key Results

1. Structural Differences: The (Li,Fe)OHFeSe single crystals exhibited a large FWHM of 1.49° for the (006) peak, indicating significant interlayer misalignment. In contrast, the epitaxial films showed a reduced FWHM of 0.12°, demonstrating superior interlayer alignment.
2. Dimensionality Discrepancy:
   • Epitaxial Films: The $H_{c2}(\theta)$ data displayed a cusp-like peak near $\theta = 90^\circ$, fitting well with the standard 2D Tinkham model and yielding a high anisotropy ratio ($\gamma \approx 151$).
   • Single Crystals: The $H_{c2}(\theta)$ data appeared rounded and was better fitted by the 3D anisotropic Ginzburg-Landau (G-L) model, yielding a much lower anisotropy ($\gamma \approx 37$).
3. Confirmation of Bulk 2D Nature: Despite the differences in $H_{c2}(\theta)$ fitting, both sample types exhibited clear BKT transitions, evidenced by the scaled resistive transitions and the power-law behavior ($V \propto I^\alpha$) with $\alpha \approx 3$ at $T_{BKT}$. This confirms that both samples possess intrinsic bulk 2D superconductivity ($c > \xi_{\perp}$), and the "3D" appearance of the single crystals is an artifact of poor crystalline alignment.
4. Validation of the Extended Model:
   • By incorporating the specific FWHM values (1.49° for crystals, 0.12° for films) into the extended 2D Tinkham model, the authors successfully reproduced the experimental $H_{c2}(\theta)$ data for both sample types using a consistent intrinsic anisotropy ($\gamma \approx 128$).
   • The model was further validated on (CTA)$_{0.5}$SnSe$_2$, a material with a large FWHM (4°) that previously deviated from the standard 2D Tinkham model. The extended model provided a significantly better fit without introducing new free parameters, confirming its bulk 2D nature.

Significance and Claims
The paper claims to resolve the long-standing confusion regarding the superconducting dimensionality of intercalated layered superconductors. The authors assert that:

• Interlayer misalignment is a critical factor: Imperfect crystalline alignment, common in materials synthesized via soft chemical methods (e.g., hydrothermal or electrochemical intercalation), distorts $H_{c2}(\theta)$ measurements, leading to the misclassification of 2D superconductors as 3D.
• The Extended 2D Tinkham Model is a robust tool: This model quantitatively captures the blurring effects of interlayer misalignment, allowing for the extraction of intrinsic superconducting properties (such as the true anisotropy $\gamma$) from bulk samples that do not possess perfect crystalline orientation.
• Broad Applicability: The method is presented as a general approach for investigating bulk 2D superconductivity in intercalated compounds, particularly those where high-temperature annealing is not feasible, and where interlayer misalignments are inevitable.

The work does not propose new applications or future experimental directions beyond the scope of characterizing these materials but emphasizes the necessity of accounting for structural imperfections when interpreting transport data in layered superconductors.