Fit-Free Optical Determination of Electronic Thermalization Time in Nematic Iron-Based Superconductors

This paper introduces a fit-free nematic response function model that combines with the two-temperature model to directly extract and quantify the anisotropic electronic thermalization timescales (110–230 fs) in nematic iron-based superconductors, offering a powerful alternative to complex data-fitting procedures.

The Gist

The Big Picture: Catching a "Hot" Electron Cooling Down

Imagine you have a crowded dance floor (the material) where everyone is dancing wildly because a loud DJ just dropped a beat (a laser pulse). Suddenly, the music stops. The dancers (electrons) are still moving fast and hot, but they need to calm down and return to a normal, cool state.

The scientists in this paper wanted to measure exactly how long it takes for these dancers to stop spinning and settle down. In the world of superconductors (materials that conduct electricity with zero resistance), this "cooling down" time is crucial because it tells us how the material works and why it might be a great conductor.

The Problem: The "Fitting" Nightmare

Usually, to measure this cooling time, scientists have to take a blurry photo of the dancers and try to guess the speed by drawing a curve over the data. It's like trying to guess the speed of a race car by looking at a blurry photo and drawing a line through it. You have to make a lot of assumptions, and if your assumptions are slightly wrong, your answer is wrong. This is called "fitting," and it can be messy and unreliable.

The Solution: The "Nematic Response" Trick

The authors developed a new, clever way to measure this time without drawing curves or making guesses. They call it the Nematic Response Function Model (NRFM).

Here is the analogy:
Imagine the dance floor isn't just a flat room; it's shaped like a long hallway. The dancers move differently depending on whether they are moving down the hallway (let's call this the "Parallel" direction) or across the hallway (the "Perpendicular" direction).

1. The Setup: The scientists shine a laser pulse on the material. This heats up the electrons.
2. The Observation: They measure how the material reflects light in two directions: down the hallway and across it.
3. The Difference: Because the material is "nematic" (meaning it has a preferred direction, like a hallway), the electrons cool down at slightly different speeds in these two directions.
   • Direction A might cool down in 100 femtoseconds (a femtosecond is a quadrillionth of a second—faster than a blink).
   • Direction B might cool down in 110 femtoseconds.

The Magic Moment: The "Valley"

Instead of trying to calculate the speed of each direction separately, the scientists subtract the two measurements from each other.

• Imagine you have two runners. One is slightly faster than the other.
• If you plot the difference between their positions over time, the line goes up, hits a peak, and then goes down.
• In this experiment, the "difference line" creates a deep valley (a minimum point).

The Breakthrough: The scientists discovered that the exact time this valley hits its lowest point is almost exactly the average time it takes for the electrons to cool down.

• Old Way: "Let's guess the curve, adjust the knobs, and hope we get the right answer."
• New Way: "Look at the graph. Where is the bottom of the valley? That's the answer."

It's like finding the exact moment a pendulum stops swinging by looking at the lowest point of its arc, rather than trying to calculate the physics of the swing from scratch.

Why This Matters

1. No Guessing: This method is "fit-free." You don't need to tweak parameters to make the math work. You just look at the data, find the bottom of the valley, and read the time. It's much more reliable.
2. Speed: They tested this on iron-based superconductors (materials that might one day power our cities with loss-free electricity). They found the electrons cool down in about 110 to 230 femtoseconds.
3. Consistency: They compared their "valley method" with the old "curve-fitting method" (called the Two-Temperature Model). The results matched perfectly! This proves their new trick works.
4. Universal Tool: This trick works for any material that has this "hallway" shape (electronic nematicity). It gives scientists a powerful new tool to map out how fast energy moves in these materials without getting bogged down in complex math.

Summary in One Sentence

The scientists invented a clever trick where they subtract two measurements from each other to find a "valley" in the data; the bottom of that valley instantly reveals exactly how fast electrons cool down, removing the need for messy math and guesswork.

Technical Summary

1. Problem Statement

Iron-based superconductors (FBSs) exhibit complex electronic orders, including electronic nematicity (a phase that breaks rotational symmetry without magnetic order). Understanding the dynamics of these systems requires measuring the electronic thermalization time ($\tau_e$), which characterizes how quickly excited electrons redistribute energy and reach a quasi-thermal equilibrium.

• Current Challenges:
  • Traditional methods rely on the Two-Temperature Model (TTM) or multi-exponential fitting of transient optical signals. These approaches require complex, global fitting procedures with multiple free parameters, which can introduce ambiguity and bias.
  • While Angle-Resolved Photoemission Spectroscopy (ARPES) can measure electron-electron scattering rates, it provides momentum-dependent data that differs from the macroscopic, k-space averaged $\tau_e$ obtained via optical methods.
  • There is a need for a robust, fit-free method to directly extract $\tau_e$ and its anisotropy in nematic materials without relying on complex modeling assumptions.

2. Methodology

The authors propose a new framework called the Nematic Response Function Model (NRFM). This method leverages polarization-resolved ultrafast pump-probe spectroscopy to extract thermalization times directly from the data's temporal features.

A. Theoretical Framework

1. Nematic Response Function ($\eta$):

   • The model defines the nematic response as the difference between transient reflectivity signals measured along two orthogonal crystallographic directions (parallel $\parallel$ and perpendicular $\perp$ to the nematic axis):
     $$\eta(t) = \left(\frac{\Delta R}{R}\right)_\parallel - \left(\frac{\Delta R}{R}\right)_\perp$$
   • Assuming the intrinsic response follows exponential decays with characteristic times $\tau_\parallel$ and $\tau_\perp$, and normalizing the amplitudes ($A_\parallel = A_\perp$), the function $\eta(t)$ exhibits a single, pronounced extremum (minimum or maximum).

2. Direct Extraction of Time Constants:

   • Infinitesimal Pulse Limit: For an ideal delta-function pulse, the time position of the extremum ($t_{min}$) corresponds directly to the average thermalization time ($\tau_{avg} = (\tau_\parallel + \tau_\perp)/2$).
   • Finite Pulse Correction: Real experiments use finite-duration laser pulses (Gaussian profile). The authors derive analytical corrections to account for the Instrument Response Function (IRF).
     • The corrected average time is given by: $\tau_{avg} \approx t_{min} - \frac{\tau_{IRF}^2}{\kappa^2 t_{min}}$ (where $\kappa \approx 2.355$).
     • The anisotropy ($\Delta\tau = \tau_\perp - \tau_\parallel$) is derived from the amplitude of the extremum ($\eta_{min}$): $\Delta\tau \approx -e \cdot \eta_{min} \cdot \tau_{avg} \cdot \exp(\dots)$.

3. Experimental Setup:

   • Samples: Epitaxial thin films of $FeSe_{1-x}Te_x$ (specifically $x=0.2$), pure $FeSe$, and optimally doped $Ba(Fe_{0.92}Co_{0.08})_2As_2$.
   • Technique: Polarization-resolved pump-probe measurements using a 35 fs, 800 nm Ti:Sapphire laser. The pump fluence was varied up to $8 mJ/cm^2$.
   • Data Processing: The two polarization channels were renormalized to ensure equal amplitudes ($r=1$), eliminating amplitude mismatch errors that could shift $t_{min}$.

3. Key Contributions

• Fit-Free Extraction: The primary contribution is the development of a methodology that extracts $\tau_e$ by identifying a specific time marker ($t_{min}$) on the nematic response curve, eliminating the need for multi-parameter global fitting.
• Analytical Correction for IRF: The paper provides a rigorous derivation of how finite pulse widths shift the extremum position and offers a closed-form analytical correction to recover the true $\tau_{avg}$.
• Anisotropy Quantification: The model allows for the direct calculation of the anisotropy in relaxation times ($\Delta\tau$) between orthogonal directions, a feature often obscured in standard TTM fits.
• Validation: The NRFM results are cross-validated against independent TTM fits, demonstrating consistency while offering superior robustness against fitting parameter choices.

4. Results

The study was conducted on three samples at temperatures below their superconducting transition ($T_c$):

• $FeSe_{0.8}Te_{0.2}$ (8 K):

  • NRFM yielded $\tau_{avg}$ values of 110–230 fs, consistent with TTM fits.
  • The method successfully captured the fluence dependence (a slight decrease in $\tau_e$ with increasing pump fluence), aligning with Gurzhi's theory for correlated electron systems.
  • NRFM and TTM showed excellent agreement in both mean times and anisotropy trends.

• $FeSe$ (3 K):

  • Similar dynamics were observed, though TTM fits were slightly less accurate ($R^2 \geq 96\%$) compared to the Te-doped sample, likely due to a more pronounced nematic phase.
  • NRFM extracted $\tau_{avg}$ values within one pulse width of the TTM results. Notably, $\tau_e$ showed an unexpected increase with pump fluence, suggesting complex photoinduced dynamics (e.g., Cooper pair breakdown).

• $Ba(Fe_{0.92}Co_{0.08})_2As_2$ (8 K):

  • NRFM and TTM yielded comparable $\tau_e$ values with no clear fluence dependence.
  • The extracted anisotropy ($\Delta\tau$) differed slightly between models but remained within the experimental time resolution (pulse duration).

• General Findings:

  • The IRF correction was found to be negligible ($\sim 3$ fs) for the experimental conditions ($\tau_{IRF} \approx 50$ fs, $\tau_{avg} \approx 100-200$ fs), validating the use of the simpler approximation $\tau_{avg} \approx t_{min}$.
  • The NRFM approach proved robust across a wide range of pump fluences and sample types.

5. Significance

• High-Throughput Analysis: By removing the need for complex fitting, NRFM enables rapid, reliable mapping of relaxation times as a function of temperature, doping, or fluence. This is crucial for studying phase transitions and critical fluctuations in nematic materials.
• Direct Physical Insight: The method isolates the electronic thermalization time from the nematic signal itself, providing a direct probe of anisotropic electron dynamics without the ambiguity of separating electron-phonon coupling constants in TTM.
• Broad Applicability: While demonstrated on iron-based superconductors, the NRFM is applicable to any material exhibiting electronic nematicity, offering a powerful tool for the broader condensed matter community to study ultrafast quasiparticle dynamics.
• Resolution of Discrepancies: The study clarifies that optical $\tau_e$ represents an integral measure of the electron subsystem's relaxation to a quasi-thermal distribution, distinct from the microscopic scattering times measured by ARPES, thereby reconciling different experimental perspectives.

In conclusion, the paper establishes a robust, fit-free optical technique that simplifies the extraction of electronic thermalization times in nematic superconductors, providing a new standard for analyzing ultrafast anisotropic dynamics.