{"id":"c2fcf10e-614e-4cf0-b32f-28c16b896a8e","arxiv_id":"2501.10720","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Optimal control using CMA-ES with cubic-spline pulses prepares Laughlin-type FQH states in optical lattices with higher fidelity and shorter time than prior adiabatic or Bayesian-optimized ramps, in exact-diagonalization simulations.","lead":"This paper simulates optimal-control ramps that prepare small fractional quantum Hall states in optical lattices faster and with higher fidelity than prior protocols. The authors show numerically that optimized tunneling and tilt sequences reach about 99 percent fidelity in about 31 tunneling times for two atoms, and extend the approach to up to four atoms in larger lattices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline fidelities are computed in the effective HHH model, not the Floquet sequence of the experiment; without a Floquet-level validation the experimental relevance of the optimized ramps remains unestablished.","rationale":"The reader's weakest assumption is exactly the one I find most load-bearing: the optimized protocols are evaluated in the effective Harper-Hofstadter-Hubbard model rather than the Floquet-driven Hamiltonian that would actually be implemented. The paper is internally consistent and the optimal-control results are plausible within the simulated model, but the abstract and introduction promise schemes 'specifically built on' the optical-lattice experiment and 'well suited' for realistic preparation. That external claim depends on the effective-Hamiltonian approximation being quantitatively accurate for the fast, non-adiabatic pulses found by CMA-ES. The paper's own concluding remarks acknowledge this gap by suggesting Floquet heating and higher-band effects as future work. A Floquet-level simulation is therefore the decisive check. Since the reader already assigned a conditional verdict on this basis, my assessment does not change that verdict.","tokens_in":17086,"tokens_out":11512,"duration_ms":133571,"concrete_test":"Re-run the two optimized schemes for the 4x4 two-boson system using the full time-periodic Hamiltonian of the experiment in Ref. [17], with the drive parameters mapped so that the stroboscopic lowest-band dynamics matches Eq. (2) at every instant; include at least the first excited band in the simulation. Compute the final overlap with the target FQH state and the excited-band population. If the fidelity stays above ~95% and the excited-band fraction remains small, the effective-Hamiltonian shortcut is validated; otherwise the experimental fidelity claim needs to be downgraded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III opens by stating 'We directly simulate the effective Hamiltonian instead of modeling the Floquet sequence as performed in the experiment [17].' All reported fidelities (99.0% at 54τ, 99% at 31τ, 96.4% at 65τ, 98.1% at 86τ) are therefore results for the time-dependent Harper-Hofstadter-Hubbard model, not for the periodically driven optical lattice. The optimized pulses are non-adiabatic and deliberately pass through excited states (Fig. 2), so there is no adiabatic protection against Floquet heating or interband transitions. The conclusion itself concedes that applying the protocols to the driven system 'would be interesting' and that higher bands 'beyond the tight-binding regime' remain to be studied. Since the central selling point is an experimentally ready preparation scheme for the Leonard et al. setup, the missing Floquet-level simulation is the load-bearing gap. The numerical results are internally consistent, and the GRAPE cross-check in Appendix A supports the optimization methodology, but these checks validate the effective model, not the experimental setting.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes optimal-control protocols for preparing bosonic Laughlin-type fractional Chern insulator states in finite optical lattices, motivated by the experiment of Léonard et al. [17]. The control fields are the tunneling amplitudes and linear gradients of the Harper-Hofstadter-Hubbard Hamiltonian (2); they are parametrized by cubic splines and optimized with CMA-ES. For the 4x4 lattice with two bosons, a four-step and a two-step scheme reach fidelities of 99.0% at 54τ and 99% at 31τ, respectively; for 6x6 (N=3) and 4x8 (N=4) lattices the two-step scheme reaches 96.4% at 65τ and 98.1% at 86τ. The paper also reports robustness to control noise, time-correlated noise, and static disorder, and cross-checks the optimization with GRAPE.","tokens_in":17305,"tokens_out":10018,"duration_ms":100817,"significance":"The computational study is careful within its stated model: exact diagonalization is used, the cost function is a well-defined fidelity to externally fixed ground states, the robustness tests are extensive, and the GRAPE comparison strengthens confidence in the CMA-ES results. If the high fidelities survive a full simulation of the Floquet-engineered experiment, the protocols would be a meaningful step toward faster FQH preparation. At present, however, the central quantitative claim applies only to the effective HHH model; the experimental relevance, which is the paper's main selling point, remains to be established.","major_comments":[{"comment":"The paper states \"We directly simulate the effective Hamiltonian instead of modeling the Floquet sequence as performed in the experiment [17].\" All reported fidelities (99.0% at 54τ, 99% at 31τ, 96.4% at 65τ, 98.1% at 86τ) are therefore results for the time-dependent Harper-Hofstadter-Hubbard model, not for the periodically driven optical lattice. The optimized pulses are non-adiabatic and deliberately pass through excited states (Fig. 2), so there is no adiabatic protection against Floquet heating or interband transitions. The conclusion itself concedes that applying the protocols to the driven system \"would be interesting\" and that higher bands \"beyond the tight-binding regime\" remain to be studied. I ask the authors to provide a Floquet-level simulation for at least the 4x4 system with the same optimized envelopes, or, alternatively, to explicitly reframe the claims as effective-model control results and remove the direct comparison with the measured 43(6)% fidelity of Ref. [17].","section":"III (opening paragraph)"},{"comment":"Scheme II is initialized with Δx = 100ℏ/τ, which is far outside the parameter range -4ℏ/τ ≤ Δx,y ≤ 4ℏ/τ stated in Eq. (4) and used to define the \"fair comparison\" with Refs. [17,22]. The control field plotted in Fig. 3(d) does not appear to extend to such a value, so the text, figure, and parameter constraint are mutually inconsistent. Since the speed advantage of Scheme II relies on this initialization, the authors must clarify the actual initial tilt and justify that it is experimentally accessible within the quoted range, or correct the bound.","section":"III.B and Eq. (4)"}],"minor_comments":[{"comment":"There is a typo: \"paradiagmatic\" should be \"paradigmatic\".","section":"II (first paragraph)"},{"comment":"The boundary conditions of the lattices are never stated explicitly; the figures suggest open boundary conditions, but this should be confirmed because it affects the spectrum, degeneracies, and the definition of the target ground states.","section":"III.A and IV"},{"comment":"The white-noise amplitude σ in Sec. III.C.1 and the OU volatility σ in Appendix B are not directly comparable, because the stationary variance of the OU process is σ²/(2θ); the authors should state the convention used in Fig. 11.","section":"III.C.1 and Appendix B"},{"comment":"The GRAPE appendix reports \"the minimal duration is estimated T2 = 4τ\" without describing how minimality was determined; please specify the discretization, stopping criterion, and search procedure.","section":"Appendix A"},{"comment":"For the 6x6 and 4x8 systems, please specify the precise form of the initial product state (which column and which sites are occupied) and the boundary conditions used in the exact diagonalization.","section":"IV (Figs. 8-9)"},{"comment":"The number of spline knots, initial parameter values, population size, and number of CMA-ES generations are not reported; providing these details, together with a data/code availability statement, would make the optimization reproducible.","section":"General (reproducibility)"},{"comment":"Please clarify whether the Ref. [22] fidelity values at T=10τ and T=20τ are simulation results and whether they include the same disorder model as the present work.","section":"Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The paper is technically solid within the effective model, and the missing Floquet simulation is the main reason I cannot recommend acceptance in the current form. The Δx = 100ℏ/τ discrepancy also needs a definitive resolution. I found no circularity: the target states are fixed ground states defined independently of the optimization. The authors should be given the opportunity to add the Floquet-level check or to re-scope the paper's claims before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a competent numerical optimal-control paper: it does what it says, and it says what it does. The main result is clear — CMA-ES with cubic-spline control fields prepares Laughlin-type states in the Harper-Hofstadter-Hubbard model with ~99% fidelity in about 31τ for the 4×4 two-boson system, and 96–98% fidelity for 3–4 particles in 6×6 and 4×8 lattices. That is a real, if incremental, advance over the Bayesian-optimization results in Ref. [22]: better fidelity in less time, plus the extension to larger systems. The GRAPE cross-check in the appendix is a nice methodological touch, and the robustness tests against control noise and disorder are carefully done and appropriately averaged. The target states are ground states of the same Hamiltonian, so the high fidelity is a legitimately optimized result, not a fitted prediction.\n\nThe soft spots are all in proportion. The big one is that all of these numbers are computed in the effective HHH model, not the Floquet sequence of the Leonard et al. experiment. The paper says this explicitly in Section III, and the conclusion flags the open question of higher bands and Floquet heating. But that means the headline claim of an \"experimentally ready\" protocol is not yet established: these ramps are non-adiabatic and deliberately pass through excited states, so they have no adiabatic protection against the driving-induced effects. The comparison to the experimental 43% fidelity is also apples-to-oranges, since that number includes the full Floquet dynamics. A Floquet-level simulation of the optimized ramps is the obvious missing piece, and I suspect the fidelities will drop when it is done. Still, this is a limitation the authors are honest about, not a hidden flaw.\n\nMinor issues: no code or data is provided, so the numbers are not independently reproducible; and the \"scalability\" claim rests on only 4 particles, which is honest but modest. The self-citation is fine. The paper is clearly written and the methods are standard enough that the results should be trustworthy.\n\nWho gets value from this: people working on cold-atom FQH simulation, optimal control of many-body systems, and adiabatic vs. diabatic state preparation. I'd bring it to a reading group as a good example of what optimal control can do for a genuinely hard many-body problem, and what it can't yet do for the real experimental setting. It deserves a serious referee; the main request in a report should be a Floquet-level validation or an explicit statement that the protocol is only for the effective model.","headline":"Solid numerical optimal control for FQH state preparation in the HHH model, with honest limitations: the headline fidelities are not validated in the actual Floquet-driven experiment, but the paper is transparent about this and the internal results hold up.","tokens_in":17875,"tokens_out":2167,"would_cite":true,"duration_ms":23663,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Optimized control ramps prepare a two-boson Laughlin-type fractional quantum Hall state in an optical lattice at 99% fidelity in 31 tunneling times, using only four experimental knobs.","keywords":["fractional quantum Hall states","Laughlin-type states","optical lattices","ultracold atoms","quantum optimal control","CMA-ES","Harper-Hofstadter-Hubbard model","state preparation"],"falsifier":"Run the two-step optimized ramp in the actual 4×4 two-atom experiment and measure the fraction of atoms that return to the initial state after reversing the ramp; a return probability well below 90% would show that the effective-model simulation misses heating or higher-band losses. A cheaper numerical check is to simulate the same ramp with the full Floquet Hamiltonian and compare the final fidelity with the effective-model value.","tokens_in":16887,"feed_emoji":"⚛️","tokens_out":17259,"duration_ms":150769,"temperature":0.7,"pith_summary":"The paper claims that small Laughlin-type fractional quantum Hall states—topological phases realized by strongly interacting bosons in a magnetic field—can be prepared in optical lattices much faster than adiabatic protocols allow, by numerically optimizing just four experimental knobs: the two tunneling amplitudes and the two tilt gradients. Using cubic-spline parameterized ramps refined by the covariance-matrix-adaptation evolution strategy, the authors reach about 99% fidelity for the two-boson 4×4 Laughlin-type state in 31 tunneling times, while the original experiment reported about 43% fidelity at 100 tunneling times and the previous optimized protocol reached 94.5% at 100 tunneling times. The same two-step scheme prepares three-boson states in a 6×6 lattice at 96.4% fidelity in 65 tunneling times and four hard-core bosons in a 4×8 lattice at 98.1% in 86 tunneling times. If correct, this makes small Laughlin-type states practical to create with existing quantum-gas-microscope setups and opens a route toward larger strongly correlated topological states.","feed_headline":"Control ramps hit 99% fidelity for fractional quantum Hall states","feed_subtitle":"A two-step control ramp reaches near-unit fidelity in 31 tunneling times, beating the previous 100-tunneling-time benchmark.","key_machinery":"The load-bearing object is the ramp protocol itself: the time-dependent control fields $t_x(t)$, $t_y(t)$, $\\Delta_x(t)$, $\\Delta_y(t)$ of Hamiltonian (2), each represented by a cubic spline through evenly spaced interpolation points. CMA-ES, the covariance-matrix-adaptation evolution strategy, is a gradient-free optimizer that proposes candidate spline parameters, computes the fidelity $F = |\\langle\\psi_f|\\psi_{\\rm target}\\rangle|^2$, and iteratively updates its search distribution to maximize $F$. The two-step protocol couples $t_y$ with $\\Delta_y$ in the first stage and $t_x$ with $\\Delta_x$ in the second, letting the particles delocalize along one axis and then the other. The mechanism that makes the preparation fast is that the optimizer is free to exploit excited states as intermediate stages rather than being constrained to the instantaneous ground state. Robustness is handled by adding white-noise control perturbations or static disorder to the dynamics and, in one variant, using the disorder-averaged fidelity as the cost function.","core_discovery":"The paper establishes that optimized smooth ramps of tunneling amplitudes and tilt gradients can prepare the target fractional Chern insulator ground state without following the adiabatic gap. All fidelities are computed within the effective Harper-Hofstadter-Hubbard model, the tight-binding Hamiltonian with interactions, flux $\\phi = 2\\pi\\times 0.26$, and on-site interaction $U = 8\\hbar/\\tau$ for the 4×4 case, whose ground state is the Laughlin-type state at filling $\\nu = 1/2$. In the four-step scheme, where one control is varied at a time, fidelity reaches 99.0% in total time $T = 54\\tau$; in the two-step scheme, where two controls are varied simultaneously, the same state is prepared at 99% fidelity in $T = 31\\tau$. The scalability claims are that three bosons in a 6×6 lattice reach 96.4% fidelity in $T = 65\\tau$ and four hard-core bosons in a 4×8 lattice reach 98.1% in $T = 86\\tau$.","pith_inferences":["Beyond the paper, a full periodic-drive simulation of the same ramps would test whether the effective-Hamiltonian fidelities survive Floquet heating and higher-band coupling.","Beyond the paper, the same few-knob spline-plus-evolutionary-strategy recipe could be transferred to other target states, such as different fillings or non-Abelian states, since the method does not rely on the specific Laughlin gap structure.","Beyond the paper, the reported robustness is average fidelity over noise realizations; a worst-case or certified fidelity would be a stricter and more useful guarantee for applications."],"forward_implications":["The 4×4 two-boson Laughlin-type state can be prepared at 99% fidelity in 31 tunneling times, more than three times faster than the previous 100-tunneling-time optimized protocol and at higher fidelity than the original experiment's 43%.","The same two-step scheme generalizes beyond the minimal system: 96.4% fidelity for three bosons in a 6×6 lattice and 98.1% for four hard-core bosons in a 4×8 lattice.","Because the optimized ramps tolerate control noise and static disorder, and can be made more robust by disorder-aware training, they are compatible with the site-resolved controls of current quantum gas microscopes.","Smooth spline controls are preferable to piecewise-constant numerical pulses, which oscillate rapidly and would be difficult to implement experimentally.","The protocols deliberately use excited states during the ramp, so preparation time is not controlled by the many-body gap alone; this is qualitatively different from standard adiabatic preparation."],"supporting_citations":[{"why":"Provides the experimental platform and benchmark: the two-atom 4×4 optical-lattice realization of a Laughlin-type FQH state with 43(6)% fidelity at T=100τ, which the optimized protocols are designed to outperform.","marker":"[17]"},{"why":"Gives the previous best optimization benchmark (94.5% at 100τ, 78% at 20τ, 53% at 10τ) and the piecewise-linear parameterization that this work replaces with cubic splines.","marker":"[22]"},{"why":"Supplies the covariance-matrix-adaptation evolution strategy algorithm used to optimize the spline control parameters.","marker":"[31]"},{"why":"Software implementation of the evolution strategy used for the numerical optimizations reported in the paper.","marker":"[81]"},{"why":"Provides the gradient ascent pulse engineering algorithm used in the appendix to cross-check that the smooth CMA-ES protocols reach comparable fidelity without oscillatory pulses.","marker":"[23]"},{"why":"Established that the strongly interacting Harper-Hofstadter-Hubbard model hosts bosonic fractional Chern insulator states at filling ν=1/2, the target state class.","marker":"[45]"},{"why":"Supplies the Streda-marker criterion (CStr≈1/2) used to identify the fractional Chern insulator target states in the 6×6 and 4×8 systems.","marker":"[58]"},{"why":"Identifies the parameter regime, flux α≈0.3 in the hard-core limit, in which the four-boson 4×8 system hosts the FCI target state.","marker":"[79]"}],"fun_headline_variants":["Optimal ramps prepare FQH states 3x faster than adiabatic","Two-step control hits 99% fidelity in 31 tunneling times","Ramped gradients and hoppings speed FQH state prep","Control scheme accelerates Laughlin state creation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All results are simulated in a simplified model that ignores the periodic shaking used in the real experiment; if that shaking heats the atoms or drives them into higher bands, the measured fidelity will be lower than the simulated 99%.","fun_headline_variants_meta":{"raw":{"variants":["Optimal ramps prepare FQH states 3x faster than adiabatic","Two-step control hits 99% fidelity in 31 tunneling times","Ramped gradients and hoppings speed FQH state prep","Control scheme accelerates Laughlin state creation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000525,"raw_usage":{"total_tokens":2533,"prompt_tokens":939,"completion_tokens":1594,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":1523}},"tokens_in":555,"tokens_out":1594,"duration_ms":13283,"temperature":1.0,"reasoning_tokens":1523,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:02:04.282598+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the two-step optimized ramp in the actual 4×4 two-atom experiment and measure the fraction of atoms that return to the initial state after reversing the ramp; a return probability well below 90% would show that the effective-model simulation misses heating or higher-band losses. A cheaper numerical check is to simulate the same ramp with the full Floquet Hamiltonian and compare the final fidelity with the effective-model value.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the covariance-matrix-adaptation evolution strategy algorithm used to optimize the spline control parameters."},{"cited_title":"Motruk and F","cited_arxiv_id":null,"evidence_quote":"Software implementation of the evolution strategy used for the numerical optimizations reported in the paper."},{"cited_title":"Gerster, M","cited_arxiv_id":null,"evidence_quote":"Supplies the Streda-marker criterion (CStr≈1/2) used to identify the fractional Chern insulator target states in the 6×6 and 4×8 systems."},{"cited_title":"Barkeshli, N","cited_arxiv_id":null,"evidence_quote":"Identifies the parameter regime, flux α≈0.3 in the hard-core limit, in which the four-boson 4×8 system hosts the FCI target state."}],"review_version":1}