{"id":"9232183d-d22f-437e-a4cc-445467e18077","arxiv_id":"1908.07006","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For small bosonic samples in a driven optical lattice, stroboscopic dynamics supports the ν=1/2 Laughlin state for U/Jx=10 with ω/Jx≥20, and a slow ramp prepares it with high fidelity.","lead":"This paper uses exact numerical simulations to show that a few interacting bosons in a periodically driven optical lattice can form a Laughlin fractional quantum Hall state, and that a slow ramp can prepare it within about 20 tunneling times. The result maps out a parameter window of interaction strength and driving frequency where the state survives, giving cold-atom experiments a concrete target.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The simulated rotating-frame model omits boundary phasor terms that the appendix says would require extra lab engineering; the experimental-realism claim rests on this untested simplification.","rationale":"I read the paper as an exact-diagonalization study of the driven Bose-Hubbard model Eq. (1), with Laughlin identification based on PES gaps and mode counting. The internal numerics for that model are coherent: the static effective model has the known Laughlin ground states, the Floquet eigenstate analysis is done for Np=4, and the slow ramp shows persistent PES gaps with correct counting for Np=5,6. The weakest step is the translation from the lab-frame derivation in Appendix A to Eq. (1): the boundary phasor terms are dropped explicitly, with an admission that they would require additional engineering. Because all subsequent claims about heating, Floquet eigenstates, and ramp preparation use the simplified model, the experimental-realism part of the central claim is conditional on this simplification being innocuous. This is a limitation openly stated in the manuscript, not an internal inconsistency, and it can be tested directly by including the phasor terms in a small exact simulation. I considered whether the main concern should instead be that the Np=5,6 claims rest on ramp PES rather than Floquet eigenstate PES, but the PES evidence during the ramp is a legitimate dynamical diagnostic and the sharp frequency thresholds are presented mainly for Np=4. I therefore keep the reader's CONDITIONAL verdict unchanged pending the boundary-term test.","tokens_in":14708,"tokens_out":12687,"duration_ms":141687,"concrete_test":"For Np=4 on the 4x8 lattice, simulate the full lab-frame model (A.1)-(A.4) with the R1,R2 periodic boundary conditions and the slow ramp of Eq. (17), without dropping the boundary phasor terms: propagate the two initial wire states |psi_+> and |psi_-> under the time-dependent Hamiltonian, and at t=100T compute the particle-entanglement spectrum of the mixture rho(t), the PES gap, and the overlap with the Laughlin torus state. Run this at (U/Jx=10, omega/Jx=20) and (U/Jx=1, omega/Jx=15). If the Laughlin counting and gap survive with high overlap, the boundary simplification is harmless; if the gap closes or the counting fails, the stated frequency thresholds and the 'experimentally realistic' preparation claim need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the driven-lattice protocol supports and prepares a nu=1/2 Laughlin state at experimentally realistic times. Every numerical result, however, uses Eq. (1), obtained from the rotating-frame Hamiltonian (A.9) by dropping the boundary terms e^{-i*omega*t*(Ly-1)} H_{Ly-1} + e^{i*omega*t*(Ly-1)} H_{-Ly+1} and imposing periodic boundary conditions in both directions. The appendix explicitly states that realizing these terms would require engineering additional non-trivial lab-frame terms, so Eq. (1) is not the lab model with the stated geometry; it is a modified model. At stroboscopic times t=NT these phasors evaluate to 1, so they are not removed by averaging; in a Magnus expansion they first enter at second order, with a scale Jy^2/((Ly-1)*omega). For Jy=5Jx and omega/Jx=15 this is about 0.24Jx, not negligible compared with the engineered J'_y ~ 0.88Jx. These boundary terms are also not captured by the momentum-sector construction used to identify the two degenerate Laughlin-like Floquet states. Thus the claimed preparation protocol remains unverified for the actual lab boundary conditions, and this is the weakest load-bearing link in the 'experimentally realistic' part of the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper numerically studies a periodically driven Bose-Hubbard model on a small tilted lattice with a driving scheme designed to generate a synthetic magnetic flux. Using exact diagonalization of the driven model and particle-entanglement-spectrum diagnostics, it identifies a prethermal regime in which the stroboscopic dynamics supports a ν=1/2 bosonic Laughlin state for Np=4, 5, and 6 particles at intermediate driving frequencies. It further proposes a slow ramp of the y-direction tunneling and the driving amplitude, starting from decoupled wires, and reports preparation of Laughlin-like states on time scales of order 20ℏ/Jx. The paper's central technical content is a parameter scan in U and ω with no fitted parameters; its headline claim is that the protocol is experimentally realistic.","tokens_in":14872,"tokens_out":9397,"duration_ms":96089,"significance":"If the claims hold, this is a useful step toward cold-atom fractional quantum Hall physics: it identifies a finite-size prethermal window in which a driven lattice can host a small bosonic Laughlin state, and it provides a concrete dynamical preparation ramp. The exact-diagonalization calculations are standard and appropriate for the system sizes considered, and the particle-entanglement-spectrum counting is taken from established references rather than fitted, which makes the topological identification credible for the model actually simulated. The paper is also unusually explicit about the approximations made in passing from the lab frame to the simulated model. The main weakness is that the simulated model omits boundary phasor terms that the appendix itself acknowledges would require additional lab-frame engineering, so the experimental-realism portion of the central claim is not yet established.","major_comments":[{"comment":"The Hamiltonian actually simulated, Eq. (1), is obtained from the rotating-frame Hamiltonian (A.9) by dropping the boundary phasor terms e^{-iωt(Ly-1)}H_{Ly-1}+e^{iωt(Ly-1)}H_{-Ly+1} and imposing periodic boundary conditions in both directions. The appendix states that realizing these terms 'would require engineering additional non-trivial terms in the lab frame.' This is not a harmless high-frequency suppression: at stroboscopic times the phasors evaluate to unity, and in a Magnus expansion the boundary terms first contribute at second order with a scale set by Jy^2/[(Ly-1)ω]. For the parameters used in the main text (Jy=5Jx, Ly=8, ω/Jx=15) this scale is about 0.24Jx, which is comparable to the engineered effective hopping J'_y≈0.18Jy≈0.88Jx. Because the central claim includes preparation on experimentally realistic time scales, the authors should either include the boundary terms (or a controlled model of the required lab engineering) in at least the smallest-system simulations, or explicitly restrict the conclusions to the modified model defined by Eq. (1) and remove the direct experimental-realism claim for the stated geometry.","section":"Appendix A, Eqs. (A.9)-(A.11)"},{"comment":"The concluding sentence that 'the stroboscopic dynamics of Np=4,5,6 particles supports the topological ν=1/2 Laughlin state down to ω/Jx=20 for U/Jx=10, and down to ω/Jx=15 for U/Jx=1' is stronger than the evidence presented. Full diagonalization of U_F and particle-entanglement-spectrum analysis of its eigenstates are shown only for Np=4 (Fig. 4); for Np=5 and Np=6 the supporting evidence comes from the time-dependent ramp simulations in Section IV (Figs. 5-7), which demonstrate dynamical preparation but not that the eigenstates of U_F themselves are Laughlin-like in those size sectors. The authors should either add the U_F eigenstate PES for Np=5 and Np=6, or rephrase the conclusion to distinguish eigenstate support from dynamical preparation evidence.","section":"Section III B and Conclusions"},{"comment":"The initial state (7) and the energy observable (8) both use the kick operator K(t=0) approximated by its leading-order expression (A.12). At the intermediate frequencies ω/Jx=15-20 that are central to the paper, the size of the neglected higher-order terms in K is neither estimated nor tested. A quantitative check, for example comparing the leading-order K with a numerically computed micromotion operator for Np=4, would strengthen the claim that the prepared state corresponds to the full driven Floquet dynamics rather than to the truncated micromotion approximation.","section":"Eqs. (7), (8), (A.12)"}],"minor_comments":[{"comment":"The sentence 'By analyzing the effective model from Eq. (1)' should refer to Eq. (3), since Eq. (1) is the driven Hamiltonian and Eq. (3) is the effective model.","section":"Section II B, final paragraph"},{"comment":"The phrase 'overlap ... with high fidelity (better than 1%)' is ambiguous; 'fidelity above 0.99' or 'error below 1%' would be clearer.","section":"Section IV B, caption of Fig. 5(c)"},{"comment":"The notation J'_y≡κ/(2ω) sin(φ/2)J_y would be clearer with explicit parentheses; as printed, it can be misread as κ/(2ω sin(φ/2))J_y.","section":"Section II A, text near Eq. (4)"}],"recommendation":"major_revision","confidential_remarks":"The numerical core is sound for the model actually simulated, and the paper is within the journal's scope. The main question is whether the appendix's acknowledged modification of the boundary conditions can be reconciled with the experimental-realism claims; a small-system test with the boundary phasor terms, or a clearly restricted claim, would be sufficient."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here is my read of Hudomal et al. The genuinely new piece is that they simulate the full driven dynamics of a few interacting bosons rather than relying on an effective model. They find a flattish prethermal window — down to ω/Jx=20 for U/Jx=10, and down to ω/Jx=15 for U/Jx=1 — where the stroboscopic Floquet states show Laughlin-type particle-entanglement spectra, and they show that a slow ramp of Jy and κ can reach those states in roughly 20 tunneling times. The PES diagnostics are sound; the mode counting comes from the external theory of Refs. [52,53], so there is no fitted-parameter circularity. For a small-system numerical study this is a solid, useful contribution.\n\nThe soft spot is the one the authors state themselves in Appendix A. They get from the rotating-frame Hamiltonian (A.9) to the simulated model (1) by dropping the boundary phasor terms e^{-iωt(L_y-1)} H_{L_y-1} and its conjugate, imposing periodic boundary conditions in both directions. At stroboscopic times those phasors evaluate to 1, so they do not average away; in a Magnus expansion they first appear at second order, with size roughly Jy^2/((L_y-1)ω). For their parameters (Jy=5Jx, ω/Jx=15, Ly=8) this is ≈0.24 Jx, about a fourth of the engineered J'_y≈0.88Jx. That is not negligible. The authors write that reproducing these terms would require additional lab-frame engineering, but they do not attempt to bound the effect. Consequently the claim about 'experimentally realistic time scales' is conditioned on an untested simplification. The numerics for the model with these terms neglected stand on their own; the step to the lab is less secure.\n\nA second, smaller caveat is that the initial state in the ramp uses only the leading-order kick operator K(0) from Eq. (A.12). At intermediate frequencies that is probably acceptable, but it does add another approximation to the preparation protocol.\n\nThe math and the data are consistent, and the citation pattern is appropriate — in particular the static ramp of He et al. is credited as the starting point. No code or data are shipped, though the methods are standard enough that independent reproduction should be straightforward.\n\nMy recommendation: send it to peer review. A referee should push for at least an explicit estimate, or a small-system simulation, of the boundary terms. As it stands, the paper is an honest, careful numerical study that will be cited, with the experimental caveat kept in mind.","headline":"Careful exact-diagonalization study showing Laughlin-like states survive in the driven Bose-Hubbard model, but the experimental claim rests on dropping boundary phasor terms the authors themselves flag.","tokens_in":15503,"tokens_out":5935,"would_cite":true,"duration_ms":54119,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper numerically shows that a periodically driven optical lattice can support the ν=1/2 bosonic Laughlin state for 4 to 6 atoms, and that a slow ramp of driving and tunneling can prepare it in about 20 ℏ/Jx.","keywords":["fractional quantum Hall","Laughlin state","driven optical lattices","Floquet engineering","Bose-Hubbard model","particle entanglement spectrum","cold atoms","prethermalization"],"falsifier":"Keep the boundary phasor terms $e^{{-iωt(Ly-1)}}$ and $e^{{iωt(Ly-1)}}$ that the paper drops when imposing periodic boundary conditions in the rotating frame, and exactly simulate the Floquet dynamics of the driven Bose-Hubbard model; if the particle-entanglement gap and Laughlin counting no longer appear at ω/Jx = 15–20 over several hundred driving periods, the central claim would be refuted.","tokens_in":14402,"feed_emoji":"⚛️","tokens_out":13862,"duration_ms":118688,"temperature":0.7,"pith_summary":"The paper asks whether the heating induced by periodic driving in an interacting Bose-Hubbard model necessarily destroys fractional quantum Hall physics. Using exact stroboscopic simulations of four to six bosons in a driven square lattice, it shows the ν=1/2 Laughlin state survives down to a driving frequency ω/Jx = 20 at interaction U/Jx = 10, and down to ω/Jx = 15 at U/Jx = 1, as identified by the gap and mode counting in the particle-entanglement spectrum (the entanglement between subsets of particles). The authors also show that a slow ramp that turns on the driving and the inter-chain tunneling, starting from decoupled wires, prepares the Laughlin state in about 20 ℏ/Jx, an experimentally realistic time scale. If correct, this puts a cold-atom fractional quantum Hall state within reach using the same driving scheme already used to realize the Harper-Hofstadter model.","feed_headline":"Driven lattices can host a bosonic Laughlin state","feed_subtitle":"For 4–6 atoms, the ν=1/2 state survives to ω/Jx=20; a simple ramp reaches it in ~20 ℏ/Jx.","key_machinery":"The central object is the driven Bose-Hubbard model with a sinusoidal density modulation, Eq. (1), whose Floquet high-frequency expansion yields the effective Harper-Hofstadter Hamiltonian of Eq. (3) with a renormalized y-hopping J'_y = κ/(2ω sin(φ/2)) Jy. The argument proceeds by constructing the stroboscopic time-evolution operator U_F (one driving period) and examining the particle-entanglement spectrum (PES)—the negative logarithm of the eigenvalues of the reduced density matrix after tracing out a subset of particles—of the low-lying Floquet eigenstates. The PES is the load-bearing diagnostic: a gap and the Laughlin counting of modes below the gap identify the topological state. For state preparation, the key mechanism is the slow ramp in Eq. (17), which adiabatically connects decoupled wires to the driven lattice.","core_discovery":"The paper claims that the stroboscopic dynamics of Np = 4, 5, and 6 bosons in the driven Bose-Hubbard model with flux density α = 1/4 and driving amplitude κ/ω = 0.5 supports the topological ν = 1/2 Laughlin state for driving frequencies down to ω/Jx = 20 at U/Jx = 10 and down to ω/Jx = 15 at U/Jx = 1. The evidence is the particle-entanglement spectrum of the low-energy eigenstates of the one-period time-evolution operator: a clear entanglement gap with the Laughlin counting of low-lying modes (for example, 10 and 9 modes in the two momentum sectors at Np = 6). The paper further claims that a slow ramp of the y-direction tunneling and the driving amplitude, starting from decoupled wires with atoms on every second wire, prepares this Laughlin state on time scales of order 20 ℏ/Jx, with an overlap to the target Floquet eigenstates of more than 99% (error below 1%) for Np = 4.","pith_inferences":["Beyond the paper: the neglected boundary phasor terms grow with the system circumference; for the small lattices studied (Ly = 8–12) they are presumably minor, but for larger samples or different boundary conditions the stability window could shift.","Beyond the paper: a practical experimental probe could be the particle-entanglement entropy or local density correlations after the ramp; recent detection schemes for fractional excitations in small bosonic systems could be adapted to verify the topological state without measuring a Chern number.","Beyond the paper: applying the same Floquet-plus-entanglement-spectrum analysis to other fillings (for example, ν = 1) or to fermionic atoms could reveal whether the prethermal stability window found here is specific to the ν = 1/2 Laughlin state."],"forward_implications":["The same driving scheme already used to create synthetic magnetic fields in optical lattices (κ/ω = 0.5, φ = π/2) suffices to stabilize a bosonic Laughlin state; no additional magnetic-field engineering is needed.","A concrete prethermal parameter window exists: driving frequencies around ω/Jx = 15–20 (for U/Jx = 1–10) avoid both fast heating and higher-band coupling, giving experiments a target operating regime.","The slow ramp prepares the topological state in about 20 ℏ/Jx, a time scale compatible with current cold-atom experiments; ramp rates up to η/Jx = 0.1 work in the optimal interaction range U/Jx ≈ 5.","Topological features persist on intermediate time scales even in a regime that eventually heats to infinite temperature (for example, ω/Jx = 15 at U/Jx = 10), so the state may be observable before heating dominates."],"supporting_citations":[{"why":"Supplies the experimentally demonstrated driving scheme and parameter values (κ/ω = 0.5, φ = π/2) that the paper adopts.","marker":"[17]"},{"why":"Introduces the driven Bose-Hubbard model with sinusoidal density modulation used as the starting Hamiltonian.","marker":"[43]"},{"why":"Provides the decoupled-wires slow-ramp preparation protocol that the paper extends to the driven interacting model.","marker":"[11]"},{"why":"Gives the particle-entanglement spectrum counting for the Laughlin state on the torus that the paper uses for identification.","marker":"[52]"},{"why":"Establishes the particle-entanglement spectrum as a topological fingerprint for Laughlin states.","marker":"[53]"},{"why":"The heating diagnostic (normalized energy) that defines the prethermal window follows the procedure of this reference.","marker":"[58]"}],"fun_headline_variants":["Driven lattice stabilizes Laughlin state for few bosons","Fast ramp to Laughlin state in driven optical lattice","Stroboscopic drive supports ν=1/2 Laughlin state","Bosons in driven lattice reach topological Laughlin state","Laughlin state prepared in ~20 ℏ/Jx with driven lattice"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The experimental realism of the protocol rests on neglecting the time-oscillating boundary terms $e^{{-iωt(Ly-1)}}$ and $e^{{iωt(Ly-1)}}$ in the rotating-frame Hamiltonian, which the paper admits would require engineering additional non-trivial terms in the lab frame.","fun_headline_variants_meta":{"raw":{"variants":["Driven lattice stabilizes Laughlin state for few bosons","Fast ramp to Laughlin state in driven optical lattice","Stroboscopic drive supports ν=1/2 Laughlin state","Bosons in driven lattice reach topological Laughlin state","Laughlin state prepared in ~20 ℏ/Jx with driven lattice"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000319,"raw_usage":{"total_tokens":1780,"prompt_tokens":908,"completion_tokens":872,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":787}},"tokens_in":524,"tokens_out":872,"duration_ms":9172,"temperature":1.0,"reasoning_tokens":787,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:28:34.963432+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Keep the boundary phasor terms $e^{{-iωt(Ly-1)}}$ and $e^{{iωt(Ly-1)}}$ that the paper drops when imposing periodic boundary conditions in the rotating frame, and exactly simulate the Floquet dynamics of the driven Bose-Hubbard model; if the particle-entanglement gap and Laughlin counting no longer appear at ω/Jx = 15–20 over several hundred driving periods, the central claim would be refuted.","supporting_citations":[{"cited_title":"Aidelsburger, M","cited_arxiv_id":null,"evidence_quote":"Supplies the experimentally demonstrated driving scheme and parameter values (κ/ω = 0.5, φ = π/2) that the paper adopts."},{"cited_title":"Bukov and A","cited_arxiv_id":null,"evidence_quote":"Introduces the driven Bose-Hubbard model with sinusoidal density modulation used as the starting Hamiltonian."}],"review_version":1}