Nonthermal magnetization pathways in photoexcited semiconductors

1. Problem Statement

The paper addresses the challenge of stabilizing long-range magnetic order in nominally non-magnetic semiconductors using femtosecond laser pulses. While ultrafast demagnetization and all-optical switching are well-documented, the mechanism for inducing transient magnetic instabilities in non-magnetic materials remains elusive.

• The Gap: Previous theoretical studies (often using Time-Dependent Density Functional Theory, TDDFT) suggest that above-gap photoexcitation can lead to quasi-equilibrium states with broken time-reversal symmetry and magnetic order. However, standard TDDFT implementations are non-dissipative (unitary evolution).
• The Core Question: Can non-dissipative TDDFT simulations accurately reproduce the relaxation pathways toward these transient magnetic states, or is an explicit dissipation mechanism required to capture the transition from the excited state to a lower-energy broken-symmetry phase?

2. Methodology

The author employs a multi-scale modeling approach, moving from a microscopic quantum model to a phenomenological macroscopic description.

A. Minimal Real-Time Spin-Orbital Model

To isolate the fundamental mechanisms, the author constructs a toy model consisting of:

• System: A Hilbert space of four localized $s=1/2$ spins, where one spin ($S_1$) is coupled to an $l=1$ orbital angular momentum ($L_1$).
• Hamiltonian ($H(t)$):
  $$H(t) = H_{SOC} + H_{int} + H_{kick}(t)$$
  • $H_{SOC}$: Spin-orbit coupling ($\lambda$) acting on the coupled spin-orbital pair.
  • $H_{int}$: Anisotropic Heisenberg exchange interactions ($J_{ij}$) between the four spins.
  • $H_{kick}(t)$: A pulsed perturbation representing the laser electric field, acting on the orbital components to transfer angular momentum.
• Dynamics & Dissipation:
  • The system is initialized in an excited eigenstate (mimicking the quasi-equilibrium state after thermalization but before recombination).
  • Dissipative Case: The time evolution is governed by a non-unitary Schrödinger equation incorporating a "quantum friction" term:
    $$i \frac{\partial}{\partial t} |\psi(t)\rangle = \left( H(t) - i \eta [H(t) - \langle H(t) \rangle] \right) |\psi(t)\rangle$$
    Here, $\eta$ is a damping parameter that preferentially decays high-energy components, driving the system toward the lowest-energy sector compatible with symmetries.
  • Non-Dissipative Case: Standard unitary evolution ($\eta = 0$) is used for comparison, mimicking standard TDDFT.

B. Time-Dependent Ginzburg-Landau (TDGL) Model

To bridge the microscopic results with macroscopic observables, a 2D isotropic TDGL model is used to describe the order parameter ($\mathbf{M}$) dynamics:

• Free Energy: A "Mexican hat" potential ($F \propto \Gamma^2 |\mathbf{M}|^2 + \alpha |\mathbf{M}|^4$) where $\Gamma^2 < 0$ indicates instability.
• Equation of Motion: Includes inertial terms, damping ($\gamma$), and an orthogonal driving term ($a_o(t)$) to mimic the transverse spin oscillations observed in the microscopic model.

3. Key Results

A. Microscopic Spin-Orbital Dynamics

1. Role of Dissipation:
   • With Dissipation ($\eta > 0$): The laser pulse triggers a chaotic spin dynamics initially. Over time ($\sim 1000$ fs), the system relaxes to a lower energy state ($\approx -1.28$ eV vs. initial excited state). Crucially, the spins become highly correlated ($\langle S_1 \cdot S_3 \rangle \to 1$), signaling a transition to a coherent, broken-symmetry magnetic state.
   • Without Dissipation ($\eta = 0$): While the laser initiates spin dynamics and angular momentum transfer, the system remains "locked" in the excited state with constant energy. It fails to relax into the lower-energy magnetic ground state.
2. Mechanism Requirements: The emergence of spin dynamics requires the combined effect of three factors:
   • Laser-induced angular momentum injection (via $H_{kick}$).
   • Spin-Orbit Coupling (SOC) to transfer orbital momentum to spin.
   • Dissipation to allow energy relaxation into the broken-symmetry state.
3. SOC Contribution: Despite being the trigger for spin dynamics, the SOC contribution to the total energy remains negligible throughout the process, consistent with ultrafast demagnetization observations.

B. Macroscopic Order Parameter Trajectory

The TDGL model confirms that if the initial fluctuation is driven by orbital angular momentum injection (coherent), the magnetization trajectory can be coherent across the sample. Conversely, if driven by thermal fluctuations, the system forms domains with topological defects (Kibble-Zurek mechanism). The model successfully reproduces the rotating motion of the order parameter observed in the microscopic simulations.

4. Key Contributions

• Identification of the Dissipation Bottleneck: The paper provides strong evidence that standard non-dissipative TDDFT simulations are insufficient for predicting the formation of transient magnetic phases in photoexcited semiconductors. Without a mechanism to relax the photoexcited state (e.g., electron-phonon coupling), the system cannot access the lower-energy broken-symmetry states.
• Minimal Mechanistic Model: The author isolates the specific interplay between SOC, angular momentum injection, and anisotropic exchange, proving that all three are necessary for the onset of magnetic order in this context.
• Qualitative Paradigm: The proposed TDGL trajectory offers a framework for interpreting future ultrafast magnetization experiments, distinguishing between coherent (laser-driven) and incoherent (thermal) magnetization pathways.

5. Significance

• Theoretical Implications: This work challenges the reliance on pure unitary TDDFT for studying light-induced phase transitions. It suggests that to accurately model real materials where photoexcited states relax via many-body interactions, effective dissipative terms must be included in the simulation framework.
• Experimental Guidance: The findings suggest that observing long-lived magnetic order in non-magnetic semiconductors requires not just the right laser pulse (to break symmetry) but also a material environment that supports rapid energy relaxation into the new magnetic ground state.
• Device Design: Understanding these nonthermal pathways is crucial for designing ultrafast magnetic devices with unprecedented operational speeds, as it clarifies the timescales and mechanisms required to stabilize transient magnetic states.