⚛️ quantum physics

Semiclassical Langevin dynamics of long-range dissipative time crystals

This paper develops a semiclassical Langevin approach to analyze finite-size effects in dissipative time-crystalline spin systems, demonstrating that the algebraic scaling of oscillation decay rates and deviation times with system size serves as a robust diagnostic for time-crystal behavior in both spin-1/2 and spin-1 models with long-range interactions.

Original authors: Reyhaneh Khasseh, Rosario Fazio, Angelo Russomanno

Published 2026-07-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Reyhaneh Khasseh, Rosario Fazio, Angelo Russomanno

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a group of people trying to dance in perfect unison. In a perfect world (the "thermodynamic limit"), if they all hold hands and move together, they can keep dancing forever without getting tired. This is what physicists call a Time Crystal: a system that keeps oscillating or "dancing" in time, breaking the usual rule that everything eventually settles down and stops moving.

However, in the real world, people get tired, distracted, or bumped by others. This is dissipation (energy loss). Usually, this noise makes the group stop dancing and just stand still. But, if the dancers are connected by a very strong, long-range "invisible string" (long-range interactions), they might be able to keep dancing together for a very long time, even if the group isn't infinitely large.

This paper is like a sophisticated simulation that asks: "How big does the group need to be, and how strong do the invisible strings need to be, before the dancing stops?"

Here is a breakdown of their findings using simple analogies:

1. The Two Types of "Dancers"

The researchers studied two different scenarios, like two different dance styles:

The Spin-1/2 Model (The "Whispering Circle"):
Imagine a circle of people where everyone whispers instructions to everyone else, but the volume of the whisper drops off the further away the person is. If the whispers travel far enough (long-range), the group can keep a rhythm. The researchers found that even if the whispers aren't perfectly loud to everyone (meaning the connection isn't infinite), the group can still dance for a long time.
- The Discovery: They found a "sweet spot." Even when the whispers get a bit weaker (a specific mathematical threshold), the group still manages to keep dancing longer as the group gets bigger. This was surprising because older theories said the dancing should stop much earlier.
The Spin-1 Model (The "Local Noise, Long-Range Music"):
In this scenario, imagine people are standing in a noisy room (local dissipation), but they are all listening to the same loud music playing from a giant speaker system that reaches everyone equally (long-range Hamiltonian interaction).
- The Discovery: Even though everyone is getting bumped by local noise, the loud music keeps them synchronized. The researchers measured how long the group stays in sync before the noise wins. They found that as long as the music reaches far enough (a specific threshold), the group stays synchronized for a very long time.

2. The "Fading Echo" Analogy

In a small group, the dancing eventually fades away. The researchers measured how fast this fading happens.

Think of a bell ringing. In a small room, the sound dies out quickly. In a massive cathedral, the sound lingers for a long time.
They found that for these "Time Crystal" systems, the "sound" (the oscillation) lingers longer and longer as the system gets bigger.
They calculated a specific "decay rate" (how fast the sound fades). They found that this rate drops significantly as the system grows, but only if the "invisible strings" or "music" reach far enough. If the reach is too short, the sound dies out immediately, no matter how big the group is.

3. The "Magic Number" (The Threshold)

The most exciting part of the paper is finding the limit of this magic.

There is a specific number (called $\alpha$ , roughly 1.2) that acts like a cliff edge.
Below the cliff: The group dances forever (in theory) as it gets bigger. The "Time Crystal" is robust.
Above the cliff: The connection is too weak. The group stops dancing and just stands still.
The Surprise: The researchers found that the "Time Crystal" behavior survives past the point where simple, old-school physics (Mean-Field theory) predicted it should die. It's like finding that a bridge can hold more weight than the engineer's basic calculations suggested, because the bridge has some hidden flexibility (quantum fluctuations) that the simple math missed.

4. The Method: A "Semi-Classical" Simulation

How did they figure this out?

Exact simulations of these quantum systems are like trying to count every single grain of sand on a beach; it's impossible for large groups.
Instead, they used a "Semiclassical Langevin" approach. Imagine watching the dance from a drone. You can't see every tiny twitch of every dancer's finger (quantum details), but you can see the overall flow of the crowd and add a little bit of "random wind" (noise) to simulate the chaos.
This method allowed them to simulate huge groups of dancers that were too big for exact computer calculations, giving them a clear picture of how the "Time Crystal" behaves as it grows.

Summary

The paper shows that Time Crystals are more robust than we thought. Even when the connections between particles aren't perfect or infinite, the system can still maintain its rhythmic, time-breaking dance for a very long time, provided the connections reach far enough. They provided a new way to measure exactly how "tough" this dancing is, showing that it survives in a wider range of conditions than previously believed.

Technical Summary: Semiclassical Langevin Dynamics of Long-Range Dissipative Time Crystals

Problem Statement
Dissipative boundary time crystals (BTCs) represent a phase of matter where continuous or discrete time-translation symmetry is spontaneously broken in the thermodynamic limit, resulting in persistent collective oscillations. While idealized models with infinite-range (fully collective) interactions support stable limit-cycle behavior described by mean-field equations, realistic systems often involve finite-range interactions and local dissipation. A critical open question is the robustness of the time-crystalline phase when the ideal collective structure is weakened. Specifically, it is necessary to understand how finite-size effects influence the lifetime of these oscillations and to identify the thresholds of interaction range (parameterized by a power-law exponent $\alpha$ ) beyond which the time-crystal phase dissolves into a stationary state. Exact numerical simulations using the Lindblad master equation are limited to small system sizes, necessitating alternative methods to probe finite-size scaling in the thermodynamic limit.

Methodology
The authors develop a semiclassical Langevin approach rooted in the truncated-Wigner approximation (TWA) and discrete TWA (DTWA) frameworks. Instead of solving the full Lindblad master equation, the dynamics are formulated via quantum Langevin equations for local spin operators.

Formalism: The method replaces operator products with products of their expectation values (mean-field approximation) while retaining stochastic noise terms to account for leading-order quantum fluctuations generated by dissipative channels.
Spin-1/2 Model: For a spin-1/2 chain with power-law Lindblad operators (dissipation is long-range, Hamiltonian is local), the dynamics are mapped to stochastic differential equations for local spin components $s_i^\mu$ . The noise is represented by real Gaussian white noises.
Spin-1 Model: For a spin-1 system with local dissipation and long-range Hamiltonian interactions, the dynamics are formulated using a basis of Gell-Mann matrices (the $SU(3)$ generators). The local variables are the expectation values of these matrices, $\lambda_i^\mu$ .
Numerical Implementation: The stochastic equations are solved using an Itô interpretation with a Runge-Kutta scheme for deterministic drift and Gaussian Wiener increments for noise. Crucially, the authors find that for the models studied, sampling over different initial conditions (standard in DTWA) is unnecessary; a single initial condition suffices because memory of the initial state is lost rapidly compared to the oscillation period.
Diagnosis: Time-crystalline order is diagnosed by analyzing the finite-size scaling of oscillation stability. For spin-1/2, the focus is on the decay rate of the oscillation envelope. For spin-1, two probes are used: the "deviation time" (how long the finite-size trajectory stays close to the mean-field limit cycle) and the scaling of the dominant Fourier peak power.

Key Results

Spin-1/2 Model (Power-Law Dissipation):
- In the fully collective limit ( $\alpha=0$ ), the system exhibits robust oscillations where finite-size damping is suppressed as system size $L$ increases.
- For finite power-law exponents $\alpha$ , the oscillation envelope decays exponentially in time, but the decay rate $D$ scales algebraically with system size: $D(L) \sim L^{-\beta}$ .
- A positive scaling exponent $\beta$ indicates that the decay rate vanishes in the thermodynamic limit, signifying a stable time crystal.
- The exponent $\beta$ decreases monotonically as $\alpha$ increases. The time-crystalline regime persists for $\alpha \lesssim 1.2$ . Notably, this threshold is slightly above $\alpha=1$ . While mean-field theory predicts relaxation to a fixed point for $\alpha > 1$ (where mean-field is no longer exact in the thermodynamic limit), the semiclassical Langevin approach predicts persistent oscillations in the range $1 < \alpha \lesssim 1.2$ . This suggests that quantum fluctuations, captured by the method, stabilize the time crystal beyond the mean-field regime.
Spin-1 Model (Local Dissipation, Long-Range Hamiltonian):
- The system displays time-crystalline behavior when the coherent interaction is sufficiently long-range.
- Deviation Time: For $\alpha \leq 1$ , the time $t^*$ at which the finite-size trajectory deviates from the mean-field solution scales as $t^*(L) \sim L^b$ . The exponent $b$ is positive for small $\alpha$ (approx. 0.37 for $\alpha=0$ ) and decreases as $\alpha$ approaches 1, vanishing near $\alpha \approx 1.2$ . This analysis is strictly valid only for $\alpha \leq 1$ , where mean-field theory describes the thermodynamic limit.
- Fourier Peak Scaling: The power of the dominant Fourier peak scales as $P(L) \sim L^s$ . The exponent $s$ is positive for $\alpha \lesssim 1.2$ , indicating time-crystal behavior. Unlike the deviation time analysis, this probe does not rely on the assumption that the thermodynamic limit is mean-field, making it robust for $\alpha > 1$ .
- Both probes identify a critical threshold for the disappearance of the time crystal in the interval $\alpha \in [1, 1.5]$ , consistent with previous cumulant-expansion findings.

Significance and Claims
The paper claims that the semiclassical Langevin dynamics provides a useful finite-size probe for dissipative time crystals in regimes where exact Lindblad simulations are computationally intractable.

Beyond Mean-Field: The approach successfully captures the robustness of time-crystalline order in regimes where mean-field theory fails (specifically for $\alpha > 1$ in the spin-1/2 model). It demonstrates that the time-crystal phase extends beyond the strict mean-field limit, provided the dissipation or interaction remains sufficiently long-ranged.
Diagnostic Tools: The work establishes finite-size scaling exponents (decay rate $\beta$ , deviation time exponent $b$ , and Fourier peak exponent $s$ ) as direct diagnostics for the robustness of the time-crystalline phase.
Methodological Validation: The results for the spin-1 model align with previous findings from cumulant expansions, validating the semiclassical Langevin approach as a reliable tool for studying open quantum systems with long-range interactions.
Physical Insight: The study highlights that while the thermodynamic limit for $\alpha \leq 1$ is described by mean-field dynamics, the finite-size scaling towards this limit is non-trivial and depends on $\alpha$ . Furthermore, the algebraic scaling of the deviation time (rather than logarithmic) suggests regular, non-chaotic dynamics in these time-crystalline phases.

1. The Two Types of "Dancers"

2. The "Fading Echo" Analogy

3. The "Magic Number" (The Threshold)

4. The Method: A "Semi-Classical" Simulation

Summary

More like this