Robustness of Floquet topological phase at room temperature: a first-principles dynamics study

Ruiyi Zhou; Yosuke Kanai

REVIEW 4 major objections 5 minor 3 references

Robustness of Floquet topological phase at room temperature: a first-principles dynamics study

T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Floquet topological phase survives lattice motion and 300 K, but the drive window narrows.

desk verdict Useful extension of the authors' own RT-TDDFT pumping work, but the 'intact' conclusion rests on labeling Q≈0.92 as W=1 when the Floquet condition is already violated; needs a sharper invariant check. read the letter →

arxiv 2506.12005 v1 pith:SUSTTLMA submitted 2025-06-13 cond-mat.mtrl-sci

Ruiyi Zhou , Yosuke Kanai This is my paper · ORCID

classification cond-mat.mtrl-sci

keywords FloquettopologicalphaseThoulesspumpingtrans-polyacetylenereal-timeTDDFTEhrenfestdynamicsbondlengthalternationPeierlsdistortionwindingnumber

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper asks whether a Floquet topological phase—here, quantized Thouless pumping of electrons in trans-polyacetylene under a sinusoidal electric field—survives when the atomic lattice is allowed to move and when the sample is at room temperature. Using real-time time-dependent density functional theory, with Ehrenfest dynamics for the coupled electron-nucleus motion and snapshot geometries from a 300 K molecular dynamics run, the authors find that the topological phase remains largely intact: the pumped charge per cycle is close to the integer winding number $W=1$ under many driving conditions. However, the range of driving amplitudes and periods that produces the topological phase narrows, and the outcome becomes sensitive to the instantaneous bond length alternation (BLA) of the polymer chain. The practical conclusion is that observing the phase at room temperature is possible in principle but requires careful choice of the driving field, because thermally accessible geometries can push the same drive from topological to trivial response.

What carries the argument

The argument is carried by three linked objects. The first is the Floquet condition: for the Kohn-Sham Hamiltonian to be time-periodic, the initial and final time-dependent orbitals after one cycle must overlap with $|\det(S)|=1$; when this holds, the integrated pumped charge $Q(T)$, computed from the displacement of maximally localized Wannier centers, equals the winding number $W$, an integer that is 1 in the topological phase and 0 in the trivial phase. The second is Ehrenfest dynamics itself, which lets the nuclei respond to the non-equilibrium electron density and thereby breaks perfect Floquet periodicity. The third is the bond length alternation (BLA) index, the averaged difference between C-C and C=C bond lengths, which quantifies the Peierls distortion; the simulations show that as BLA decreases, the Wannier functions of the double bonds become more delocalized, and the emergence of the topological phase becomes more sensitive to the driving period and amplitude.

What would settle it

Compute or measure the pumped charge per cycle at room temperature as a function of driving amplitude and period: if for trans-polyacetylene, or a molecule with a comparable Peierls distortion, there is no driving condition under which the ensemble-averaged $Q(T)$ is an integer within numerical or experimental error, the paper's central claim that the Floquet topological phase remains intact at room temperature would be falsified. A simpler computational version is to run a full Ehrenfest trajectory from a 300 K equilibrated configuration, include nuclear quantum fluctuations, and check whether $Q(T)=1$ survives for the conditions marked as topological in the paper.

Watch

Extended reading notes

Core claim

The central claim is that nonadiabatic Thouless pumping in trans-polyacetylene is robust to both dynamical electron-nuclear coupling and room-temperature thermal disorder, but not unconditionally. When atoms are allowed to respond to the electronic current via Ehrenfest dynamics at 0 K, the Floquet topological phase ($W=1$) is still found in most of the driving-field parameter space, although satisfying the Floquet condition is harder and one case, period $T=150$ a.u. and amplitude $|A|=5\times10^{-3}$ a.u., shows $Q(T)=0.92$ instead of the ideal integer. At 300 K, structures sampled from a first-principles molecular dynamics trajectory have different bond length alternation values, and the same driving field can yield either the topological phase or the trivial phase depending on the instantaneous BLA. Because an experiment at room temperature samples all such geometries, the ensemble-average pumped charge would not be a clean integer even under favorable driving conditions; the paper therefore concludes that the Floquet topological phase remains a real response of the driven system, but the field conditions needed to observe it become more restrictive.

Load-bearing premise

The results assume the mean-field Ehrenfest treatment of the atomic motion, together with the neglect of nuclear quantum fluctuations and the use of static room-temperature snapshots rather than fully coupled hot trajectories, captures the physics that decides whether the pumped charge stays integer within one pump cycle.

Editorial extensions

If this is right

At 0 K, lattice dynamics induced by the electronic current does not destroy the Floquet topological phase over most of the studied amplitude-period grid, so the topological response is compatible with moving nuclei.
At room temperature, the pumped charge per cycle is no longer guaranteed to be an integer: structure-to-structure variations in BLA can switch the same driving field between $W=1$ and $W=0$.
For any fixed driving condition, the ensemble-averaged pumped charge at room temperature will be less than one electron per cycle, so experimental observation requires selecting a drive whose topological window is wide enough to cover the thermally sampled BLA range.
BLA acts as a useful geometric descriptor: low-BLA, more delocalized structures are harder to pump topologically, consistent with the approach to a gapless metallic limit where the winding number is not defined.
The choice of exchange-correlation functional matters: only a functional that reproduces the Peierls distortion gives a reliable prediction of the topological response, because the bond alternation sets the electronic structure.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the BLA sensitivity carries over to other conjugated polymers, bond-length alternation could serve as a design rule for room-temperature Floquet topological materials: stiff chains that preserve a large BLA under thermal motion should show a wider topological window.
The snapshot approach underestimates dynamical feedback at 300 K; a full Ehrenfest trajectory starting from hot geometries would show whether the simultaneous effect of thermal disorder and electron-nuclear coupling narrows the topological window further than either effect alone.
An explicit numerical test would be to compute the ensemble average $\langle Q(T)\rangle$ over many FPMD snapshots at a fixed drive; the paper's four geometries already suggest the average is non-integer, but a quantitative ensemble would map where the last integer plateau disappears.
The paper's timescale argument suggests that longer pump periods, which are closer to adiabatic, may be more vulnerable to lattice-induced decoherence; examining the period dependence of $|\det(S)|$ at 300 K would test whether the topological window closes from the large-$T$ side.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

First thing to know: this is a genuine extension of the authors' earlier work on nonadiabatic Thouless pumping in trans-polyacetylene, not a new framework. Adding Ehrenfest dynamics and thermal snapshots is a natural next step, and the qualitative conclusion — the Floquet topological phase survives but the driving-field window narrows — is plausible and largely supported. No parameters are fitted; Q(T) comes from first-principles dynamics, so there's no circularity.

What it does well: it actually simulates electron-nuclear coupling and thermal disorder rather than relying on a model Hamiltonian, and it introduces BLA as a structural descriptor to explain why some geometries favor the phase. The Ehrenfest vs RT-TDDFT comparison showing a change in classification at one driving condition is a concrete, useful demonstration. The caveat concerns the abstract's 'remains intact.' The paper labels Q(T)=0.92 with |det(S)|=0.93 as W=1, although their own criterion requires |det(S)|=1 for Q(T)=W. That's thresholding, not invariant computation. It also excludes Structure A at T=150, |A|=2×10^-3 with Q=0.73 because the Floquet condition is not satisfied, so the assignment is inconsistent. The honest statement is that for each geometry there are driving conditions where Q(T) is near 1 and |det(S)| is near 1, but the window is narrowed; 'intact' is an overreach. A direct computation of the Floquet winding number from the evolution operator would settle it.

Soft spots: no convergence tests are reported for time step, cutoff, supercell, or k-points; room-temperature conclusions rest on four snapshots (the authors acknowledge this); and Ehrenfest mean-field dynamics is a timescale-separation argument, which is fair for a single pump cycle but worth flagging. These are moderate, not fatal.

This is a paper for people working on Floquet engineering in realistic materials and on how nonadiabatic effects degrade topological invariants. It deserves a serious referee. The central result can be accepted after the label issue is fixed and convergence data are added; at that point it becomes a solid contribution rather than a conditional one.

Referee Report

4 major / 5 minor

Summary. The manuscript reports first-principles RT-TDDFT and Ehrenfest dynamics simulations of nonadiabatic Thouless pumping in trans-polyacetylene, focusing on how atomic lattice dynamics (at 0 K via Ehrenfest) and thermal fluctuations (at 300 K via FPMD snapshots) affect the Floquet topological phase. The integrated pumped charge Q(T) is computed from time-dependent maximally localized Wannier centers, and the Floquet condition is monitored through |det S| between initial and final Kohn-Sham orbitals after one driving cycle. The main claims are that the Floquet topological phase remains largely intact under electron-nuclear coupling and room-temperature structural disorder, but that the driving-field parameter window (amplitude and period) for observing topological pumping becomes more restrictive and depends sensitively on bond-length alternation. The paper also compares exchange-correlation functionals for reproducing the Peierls distortion.

Significance. If fully supported, the work would be a valuable first-principles characterization of a molecular Floquet topological pump under realistic conditions. The study's strengths include the use of explicit electron-nuclear dynamics with no fitted parameters, a direct MLWF-based calculation of Q(T) and |det S|, and the identification of bond-length alternation as a structure-sensitive parameter for phase stability. However, as detailed below, the central claim is weakened by an inconsistency between the stated Floquet criterion and the phase assignments, the lack of numerical convergence tests, and the very limited ensemble sampling for the room-temperature conclusions. The significance is therefore conditional on addressing these issues.

major comments (4)

[Results and discussion, Fig. 1 and text after Fig. 2] The text states that Q(T)=W only when |det(S)|=1, yet the topological phase W=1 is assigned to the Ehrenfest case with T=150 a.u. and |A|=5e-3 a.u., for which the manuscript reports |det(S)|=0.93 and Q(T)=0.92. The sentence 'Floquet condition is not lifted appreciably' does not justify replacing the topological invariant with an approximate current integral; Q(T)=W is not defined when the Floquet condition is violated. This phase assignment is load-bearing for the abstract's claim that the Floquet topological phase 'remains intact.' The phase-assignment rule needs to be defined quantitatively, or the invariant W should be computed directly from the Floquet evolution operator rather than inferred from a threshold on Q(T).
[Computational Details and all results] No convergence tests are reported for the 0.1 a.u. integration time step, the 40 Ry plane-wave cutoff, the 55-atom supercell, or the Γ-point k-sampling. Because the central conclusions rely on the proximity of Q(T) to integer values (e.g., Q(T)=0.92 versus 1, and Q(T)=0.73) and on sharp phase boundaries in Fig. 1 and Fig. 4, the numerical uncertainty in these quantities must be established. Without such tests, the physical versus numerical origin of deviations from integer Q(T) cannot be cleanly separated.
[Thermal fluctuation of lattice on electron transport, Fig. 4 and Fig. 5] The room-temperature conclusions are based on only four selected FPMD snapshots. The manuscript itself acknowledges that 'one would ideally calculate the ensemble average of Q(t) to model experiments at room temperature' and that 'even with only four representative geometries' the ensemble average would not yield the 0 K integer value at t=T. Yet the abstract asserts that the Floquet topological phase remains intact at room temperature. The evidence is too sparse to support a general claim; at minimum, the abstract should be qualified to the particular structures studied, or an ensemble average over a larger set of snapshots should be provided.
[Computational Details, Ehrenfest dynamics discussion] The central result concerning lattice dynamics relies on the Ehrenfest mean-field approximation, which the paper acknowledges cannot describe energy exchange such as Joule heating (Refs. 40-41). The timescale-separation argument is plausible, but no numerical test is given to show that mean-field errors are negligible on the few-femtosecond scale of one pump cycle. A comparison with an alternative treatment (e.g., a method including decoherence or electronic friction, or an Ehrenfest simulation with a different nuclear thermostatting scheme) for at least one representative driving condition would strengthen the claim that the observed changes in |det S| and Q(T) are physical rather than artifacts of the mean-field approximation.

minor comments (5)

[Eq. (1)] Equation (1) is typeset with a garbled prefactor ('L!' and missing normalization); please provide the correct Resta-formula expression so that the definition of Q(T) is unambiguous.
[Fig. 4 caption] In the caption of Figure 4, the orange box is described with the same parameters as the pink box (T=125 a.u., |A|=2e-3 a.u.); from the main text, the orange box should correspond to T=150 a.u. and |A|=2e-3 a.u.
[Throughout] There are several typographical errors, including 'first-principal' for 'first-principles' and 'as similarly done' for 'as was similarly done'; these should be corrected.
[BLA definition] The definition of the bond-length alternation index in the text is corrupted by typesetting (the displayed formula is unreadable), making it impossible to verify the normalization; please provide a clean equation.
[Fig. 1 caption] The caption for Figure 1 describes the bottom-panel colors only as 'purple and blue colors' without stating which color corresponds to the trivial phase and which to the topological phase; an explicit legend is needed.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Q(T) is computed from first-principles dynamics and the topological phase label is an inference, not an input.

full rationale

The paper's derivation chain is self-contained. The central observable Q(T) is computed by real-time TDDFT/Ehrenfest dynamics via time-dependent MLWFs (Eq. 1); no parameter is fitted to the target 'topological phase' outcome, and the result is not benchmarked against or forced by its own output. Self-citations (Refs. 22-24, 39, 44-46, 52) supply prior demonstrations, the Wannier-center current formula, and code components, but the formula is the standard Resta/Wannier-center transport expression (Refs. 53-55), and the new lattice-dynamics/thermal-fluctuation results are independent first-principles simulations rather than restatements of those citations. The paper explicitly flags its own limitations: Ehrenfest mean-field dynamics cannot describe energy exchange such as Joule heating (Computational Details, Refs. 40-41), nuclear quantum fluctuations were not included (Conclusion), and only four FPMD snapshots were used instead of an ensemble average (Results). These limitations reduce confidence in the room-temperature conclusion, and the classification of some non-integer Q(T) values (e.g., Q(T)=0.92 with |det S|=0.93 labeled W=1) is a verification concern rather than circularity, because the Floquet-condition theorem Q(T)=W is used as an external inference rule, not as a fitted input. No equation reduces to another by construction, and no fitted parameter is renamed as a prediction.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The central claim rests on standard DFT/RT-TDDFT machinery and on specific modeling choices (SCAN0, Ehrenfest dynamics, four snapshots, threshold on |det(S)|). No free parameters are fitted to data, and no new physical entities are introduced.

assumptions (6)

domain assumption The adiabatic approximation for the exchange-correlation potential in RT-TDDFT is accurate enough for the driven electron dynamics.
Invoked in Computational Details (Refs 56-57); the KS Hamiltonian depends on time-dependent density through the adiabatic XC potential. Standard but not benchmarked against a more accurate many-body method.
domain assumption Ehrenfest mean-field dynamics sufficiently captures electron-ion coupling on the few-femtosecond pump timescale.
Stated in Computational Details; the authors argue the pump timescale is much shorter than the detailed-balance timescale, but acknowledge Ehrenfest cannot describe certain energy exchange such as Joule heating (Refs 40-41). This is load-bearing for the 0K lattice dynamics results.
domain assumption The four selected FPMD snapshots with different BLA values represent the range of room-temperature thermal structures relevant to the topological pump.
Results section (thermal fluctuation subsection); no ensemble average is computed, and the authors state an ideal calculation would ensemble-average Q(t). The conclusion about room-temperature behavior depends on these four structures.
domain assumption The Floquet condition is adequately measured by |det(S)| with a 0.95 threshold, and when it holds Q(T) equals the winding number W.
Results section and Figure 1; the 0.95 threshold is a coloring heuristic, not a derived criterion. The relation Q(T)=W is the standard Thouless result.
domain assumption The Gamma-point approximation with a 55-atom supercell is sufficient to describe the Floquet topological phase in trans-polyacetylene.
Computational Details; adopted from prior work (Refs 22, 24), no convergence study is shown here.
domain assumption Nuclear quantum fluctuations can be neglected; classical nuclei are sufficient for the room-temperature snapshots and Ehrenfest dynamics.
Conclusion states nuclear quantum effects were not considered and could introduce additional challenges; this is an acknowledged limitation.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Robustness of Floquet topological phase at room temperature: a first-principles dynamics study." pith.science (2026). https://pith.science/paper/SUSTTLMA

@misc{pith2026250612005,
  author       = {Pith},
  title        = {Pith review of: Robustness of Floquet topological phase at room temperature: a first-principles dynamics study},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SUSTTLMA}},
  note         = {Machine review of arXiv:2506.12005}
}

read the original abstract

Nonadiabatic Thouless pumping of electrons is studied within the framework of topological Floquet engineering, particularly focused on how atomic lattice dynamics affect the emergent Floquet topological phase in trans-polyacetylene under the driving electric field. As similarly done in the earlier work [Zhou and Kanai, J. Phys. Chem. Lett., 12, 4496 (2021)], the real-time time-dependent density functional theory and Ehrenfest dynamics simulations were used to investigate the extent to which the number of pumped charges remains equal to the topological invariant, the winding number, when the temperature effect of ions and the dynamical coupling of electrons and ions are taken into account. Our theoretical work shows that the Floquet topological phase remains intact but the condition on the driving field necessary for observing the topological phase becomes more limiting.

Figures

Figures reproduced from arXiv: 2506.12005 by the authors.

**Figure 1.** (Top) The determinant of the overlap matrix, S, between the initial and final TD-KS states in a single driving cycle and (Bottom) the time-integrated current over a single driving cycle 𝑄(𝑇) plotted as a function of the driving electric field amplitude |A| and the period 𝑇 for (a) RT-TDDFT and (b) Ehrenfest dynamics. the colors indicate the value of |det(S)| in Top figures; red (>0.95), orange (0.90~0.95), and yello… view at source ↗

**Figure 2.** Electron dynamics from RT-TDDFT and Ehrenfest dynamics simulation of trans-polyacetylene with the driving field of T=150 a.u./|A|=5 × 10−3 a.u. (as marked by the red box in [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

**Figure 3.** Average spread values of C=C (red) and C-C (green) Maximally-localized Wannier functions (MLWFs) in the electronic ground state for structures with different Bond Length Alternation (BLA) values. The structures are taken from the trajectory of FPMD simulation of trans-polyacetylene at 300K. The vertical dashed line represents the BLA value of the equilibrium (0K) structure. The error bars indicate the standard devia… view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: (Top) Time-integrated current over one driving cycle Q(T), and (Bottom) determinant of the overlap matrix, S, between the initial and final TD-KS states in a single driving cycle. These values are plotted as a function of the driving electric field amplitude |A| and th…

**Figure 5.** Figure 5: Time-integrated current as a function of time, 𝑄(𝑡), of trans-polyacetylene for the structures of various BLA values with the driving field of (a) T=125 a.u./|A|=2 × 10−3 a.u. (as marked by the pink square box in [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

3 extracted references · 2 canonical work pages

[18]

F.; Schoop, L

Khoury, J. F.; Schoop, L. M., Chemical Bonds in Topological Materials. Trends in Chemistry 2021, 3, 700-715. 19. Martín Pendás, A.; Contreras-García, J.; Pinilla, F.; Mella, J. D.; Cardenas, C.; Muñoz, F., A Chemical Theory of Topological Insulators. Chemical Communications 2019, 55, 12281-12287. 20. Castro, A.; De Giovannini, U.; Sato, S. A.; Huebener, H...

arXiv 2021
[35]

G., Canonical Dynamics: Equilibrium Phase-Space Distributions

Hoover, W. G., Canonical Dynamics: Equilibrium Phase-Space Distributions. Physical review A 1985, 31, 1695. 36. Shirakawa, H.; Ikeda, S., Infrared Spectra of Poly (Acetylene). Polymer Journal 1971, 2, 231-244. 37. Runge, E.; Gross, E. K. U., Density-Functional Theory for Time-Dependent Systems. Physical Review Letters 1984, 52, 997-1000. 38. Yabana, K.; B...

work page 1985
[51]

A.; Rubio, A., Propagators for the Time-Dependent Kohn–Sham Equations

Castro, A.; Marques, M. A.; Rubio, A., Propagators for the Time-Dependent Kohn–Sham Equations. The Journal of chemical physics 2004, 121, 3425-3433. 52. Yost, D. C.; Yao, Y.; Kanai, Y., Propagation of Maximally Localized Wannier Functions in Real-Time Tddft. The Journal of chemical physics 2019, 150, 194113. 53. Resta, R., Quantum-Mechanical Position Oper...

work page 2004

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.

[1] [18]

F.; Schoop, L

Khoury, J. F.; Schoop, L. M., Chemical Bonds in Topological Materials. Trends in Chemistry 2021, 3, 700-715. 19. Martín Pendás, A.; Contreras-García, J.; Pinilla, F.; Mella, J. D.; Cardenas, C.; Muñoz, F., A Chemical Theory of Topological Insulators. Chemical Communications 2019, 55, 12281-12287. 20. Castro, A.; De Giovannini, U.; Sato, S. A.; Huebener, H...

arXiv 2021

[2] [35]

G., Canonical Dynamics: Equilibrium Phase-Space Distributions

Hoover, W. G., Canonical Dynamics: Equilibrium Phase-Space Distributions. Physical review A 1985, 31, 1695. 36. Shirakawa, H.; Ikeda, S., Infrared Spectra of Poly (Acetylene). Polymer Journal 1971, 2, 231-244. 37. Runge, E.; Gross, E. K. U., Density-Functional Theory for Time-Dependent Systems. Physical Review Letters 1984, 52, 997-1000. 38. Yabana, K.; B...

work page 1985

[3] [51]

A.; Rubio, A., Propagators for the Time-Dependent Kohn–Sham Equations

Castro, A.; Marques, M. A.; Rubio, A., Propagators for the Time-Dependent Kohn–Sham Equations. The Journal of chemical physics 2004, 121, 3425-3433. 52. Yost, D. C.; Yao, Y.; Kanai, Y., Propagation of Maximally Localized Wannier Functions in Real-Time Tddft. The Journal of chemical physics 2019, 150, 194113. 53. Resta, R., Quantum-Mechanical Position Oper...

work page 2004