{"id":"33fa3cdd-be95-4a56-a4a9-29fb1f3ad5fe","arxiv_id":"1908.07358","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The Rabi-Stark model produces selective k-photon interactions whose resonance frequencies depend on the bosonic state, with a trapped-ion scheme proposed to observe them.","lead":"This paper shows that adding a Stark coupling term to the quantum Rabi model makes multi-photon transitions selective, so each resonance involves only one chosen number of bosons. It also proposes a trapped-ion setup that could simulate this model and observe the effect.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"General odd-k formula Eq. (5)-(7) is asserted without derivation and omits the second-order resonance shifts already needed at k=3; five-photon peaks deviate ~4%, so extrapolation to arbitrary k is not yet supported.","rationale":"The reader's conditional verdict and weakest-assumption analysis already identify the same soft spot: the perturbation-theory resonance condition is not corrected beyond second order, and the five-photon peaks deviate from the approximate analytic formula. My reading sharpens this into a specific technical objection. The k=3 case is handled carefully: the paper shows that the second-order Hamiltonian (3) induces a Stark shift that moves the resonance, and Fig. 2 confirms that the corrected tilde-delta condition is needed. The generalization to odd k in Eqs. (5)-(7) is not accompanied by a derivation in Appendix A or elsewhere, and Eq. (6) gives only the uncorrected detuning. The five-photon data in Fig. 3 are therefore compared with a formula that the paper itself has implicitly shown to be insufficient at k=3. The observed ~4% systematic shift may well be explainable by the analogous second-order correction, but the paper does not compute it. This does not undermine the core demonstration for k=1 and k=3, nor the trapped-ion implementation, so the appropriate verdict remains conditional rather than reject. The requested check is a finite, well-defined calculation: deriving the fifth-order effective Hamiltonian with the second-order shift and comparing the corrected resonance condition to the three observed peaks. This would settle whether Eq. (5)-(7) is the right general expression or merely a leading-order approximation.","tokens_in":21042,"tokens_out":6196,"duration_ms":65750,"concrete_test":"Compute the fifth-order effective Hamiltonian for the Rabi-Stark Hamiltonian using the same Dyson-series method as in Appendix A (or an equivalent high-order Schrieffer-Wolff transformation), this time including the second-order Stark-shift correction as was done for k=3. Then solve the corrected resonance condition tilde-delta(5)_{N0-} = 0 for omega0/omega at N0 = 2, 3, 4 with g/omega = 0.1 and gamma/omega = 0.9, and compare the roots with the observed peaks -3.227, -5.072, -6.918. If the corrected condition matches the peaks, Eqs. (5)-(7) need to be amended to include the shift; if it does not match, quantify the remaining discrepancy and test whether a fifth-order perturbative treatment truncated at leading order is valid at g/omega = 0.1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of the paper is that the Rabi-Stark model produces selective k-photon interactions whose effective Hamiltonian is given by Eqs. (5)-(7). At k=3, however, the naive resonance condition delta(3)=0 does not match the numerical peak; the authors must go to an interaction picture with respect to the second-order Hamiltonian (3), shifting the detuning to tilde-delta = delta + Delta_e - Delta_g. This correction is essential for the agreement shown in Fig. 2. The general odd-k formula in Eq. (6) contains no such correction, and it is asserted by analogy ('following the same procedure') rather than derived. The text itself acknowledges that the analytic calculation of exact resonance frequencies for higher orders rapidly becomes challenging. Consequently, the five-photon resonances are predicted with the uncorrected formula omega0^c/omega = 5 - gamma(2N0+5), while the numerically observed peaks in Fig. 3 are omega0/omega = -3.227, -5.072, -6.918 versus -3.1, -4.9, -6.7: a systematic ~4% shift that grows with N0. If the second-order Stark shift is responsible, then Eqs. (5)-(7) are not the complete effective description and the claimed resonance condition and (g/omega)^k scaling for arbitrary k are unverified. If it is not solely responsible, some other neglected higher-order or off-resonant process contributes. Either way, the general k-photon prediction is demonstrated convincingly only for k=1 and k=3, and the extrapolation to k=5 and beyond rests on an unproven formula.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the quantum Rabi model with a Stark coupling term (the Rabi-Stark model), H = (ω0/2)σz + ω a†a + γ a†a σz + g(σ+ + σ−)(a + a†), and claims that the interplay of rotating and counter-rotating terms produces selective k-photon interactions whose resonance frequencies depend on the bosonic Fock state through the Stark shift. The authors derive an interaction-picture Hamiltonian, obtain second- and third-order effective Hamiltonians via time-dependent perturbation theory, and assert a general odd-k formula (Eqs. (5)–(7)). They validate the three-photon dynamics against numerical integration of the full Hamiltonian after including a second-order Stark-shift correction, and report numerical five-photon resonances. The second half proposes a trapped-ion implementation using three laser drivings, derives an effective Rabi-Stark Hamiltonian that includes higher-order carrier corrections, and compares the full ion Hamiltonian with the effective model for one- and three-photon processes.","tokens_in":21366,"tokens_out":9688,"duration_ms":99797,"significance":"The selective multiphoton interactions and the trapped-ion simulation, if fully established, would be a useful contribution to quantum control and quantum simulation in the strong and ultrastrong coupling regimes. The paper has clear strengths: the Dyson-series derivations for the second- and third-order Hamiltonians are explicit and detailed; the numerical validations are performed on the full Hamiltonian (1) with no fitted parameters; the trapped-ion mapping in Appendix C keeps important second-order corrections; and dissipative effects are analyzed in Appendix B. The main weakness is that the general odd-k analytical claim, which is the central result, is asserted rather than derived and is incomplete without the second-order resonance shifts that are essential already at k=3; the five-photon numerical check shows a systematic resonance shift that Eqs. (5)–(7) do not predict.","major_comments":[{"comment":"The general odd-k formulas are not supported by the evidence presented. For k=3, the raw resonance condition δ^(3)_5+ = 0 does not locate the resonance in Fig. 2(a); agreement requires shifting to δ̃^(3)_5+ = δ^(3)_5+ + Δ^e_8 − Δ^g_5 from the second-order Hamiltonian (3). Yet Eqs. (5)–(7) contain no analogous shift, and no derivation of the corrected condition for general k is given. For k=5, the predicted qubit frequencies ω0^c/ω = 5 − γ(2N0 + 5) give −3.1, −4.9, and −6.7, while the numerical peaks in Fig. 3 are −3.227, −5.072, and −6.918, a systematic 3–4% shift that grows with N0. This shows that either Eq. (5) is not the complete effective description, or higher-order corrections are required. The claim that Eqs. (5)–(7) describe selective k-photon interactions for arbitrary odd k is therefore not established. The authors should derive the general corrected resonance condition, or explicitly restrict the analytical claim to k=3 and present Fig. 3 as numerical evidence only.","section":"Section II.B, Eqs. (5)–(7) and Fig. 3"},{"comment":"The paper does not quantify the accuracy or selectivity of the effective Hamiltonian at k=5. In Fig. 3 the population transfer is partial, and the authors themselves note that the remaining population goes to states |g,N0+1⟩ and |g,N0−1⟩. No comparison between the full dynamics of Eq. (1) and the effective Hamiltonian (5) at k=5 is provided, and no fidelity or leakage estimate is given. Since the central claim is selectivity, it is necessary to show quantitatively how well Eq. (5) reproduces the full dynamics and over what parameter range, rather than only locating approximate resonance peaks.","section":"Section II.B, Eqs. (5)–(7) and Fig. 3"}],"minor_comments":[{"comment":"The displayed summation expression for δ^(k)_n± is malformed: it contains the undefined symbol δ±_k and the equality to (k−1)ω + δ±_{n+(k−1)/2} does not follow from the sum as printed. The simplified version is understandable, but the sum should be rewritten correctly.","section":"Eq. (6)"},{"comment":"The product in Ω^(k)_n± is written as s=1,3,... without an explicit upper limit or step size; it should specify s=1,3,...,k−2, and the notation δ^(s)_n± should be defined.","section":"Eq. (7)"},{"comment":"The statement that even-k transitions 'will average out as a consequence of the RWA' is misleading: the absence of even-k transitions follows from the exact parity symmetry of Hamiltonian (1), as the authors note later. The parity argument should be stated explicitly.","section":"Section II.B, even-k statement"},{"comment":"After the RWA, the term involving a†a should be σx rather than σ+, since Eq. (60) and the basis {|+⟩,|−⟩} require σx. Please check this notation.","section":"Appendix C, Eq. (59)"},{"comment":"The caption says the dark curve represents 'lower values of log10|δ̃|', but since the log diverges away from zero, the curve should be described as the locus of minima (or approximate zeros) of δ̃.","section":"Fig. 2(a) caption"}],"recommendation":"major_revision","confidential_remarks":"The k=1 and k=3 selective interactions, and the trapped-ion proposal, appear sound and are well supported by independent numerics. The main risk is the general odd-k formula: it is asserted without derivation, omits the second-order shift known to be required at k=3, and the k=5 numerics show a systematic shift. This is fixable by either deriving the corrected general resonance condition or restricting the analytical claim; hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick summary: the central phenomenon—Stark-induced selective multi-photon transitions in the Rabi-Stark model, with Fock-state-dependent resonance frequencies—is convincingly demonstrated for the one- and three-photon cases, and the trapped-ion implementation is carefully validated. I would send this to referees. But the general odd-k formula is asserted, not derived, and that distinction matters here.\n\nThe strengths first. The k=3 effective Hamiltonian is genuinely derived in Appendix A via Dyson series, and the agreement with exact numerics in Fig. 2 is real. Note that the agreement only works after including the second-order Stark shifts, by moving to an interaction picture with respect to Eq. (3); that is an honest, checkable piece of perturbation theory. The one-photon selectivity in Fig. 1 cleanly shows the Stark mechanism. The trapped-ion proposal is a cut above the usual 'in principle' simulator: Appendix C includes second-order corrections from the carrier interaction, and the numerics compare Eq. (10) against the full ion Hamiltonian with realistic parameters. Dissipative effects are addressed in Appendix B, including the honest statement that a five-photon process at g/omega = 0.1 needs kappa/omega < 1e-5. The citations to prior multiphoton-resonance work and to earlier selective interactions are appropriate; Stark-induced selectivity of multi-photon transitions is genuinely new, even if the perturbative machinery is standard.\n\nThe soft spots, in proportion. First, Eqs. (5)-(7) for general odd k are asserted 'following the same procedure' with no derivation. This is not a formality: at k=3 the naive resonance condition delta(3)=0 does not match the numerics, and the paper has to add a second-order level-shift correction. The general formula contains no such correction. Second, the five-photon validation uses the uncorrected resonance condition, and the observed peaks are systematically ~4% off, the deviation growing with N0. The paper calls them 'close,' which is fair, but it does not test whether the k=3-style correction closes the gap. Until that is done, the (g/omega)^k scaling and the resonance condition for arbitrary odd k are plausible extrapolations, not established results. The paper itself concedes that analytic resonance frequencies for higher orders 'rapidly become challenging,' which reads as an honest admission that Eqs. (5)-(7) are conjectural.\n\nWho this is for: readers working on Rabi-Stark models, USC dynamics, and trapped-ion simulation will get real value from the k=3 mechanics and the ion implementation. The general-k claim should not be cited as proven. Recommendation: send to peer review, and have the referee press on the status of Eqs. (5)-(7). With either a derivation or an explicit narrowing of the claim, this becomes a solid PRA-type paper.","headline":"The selective k-photon mechanism is convincingly demonstrated for k=1 and k=3 with an honest, well-validated trapped-ion proposal; the general odd-k formula is an unproven ansatz and should be treated as such.","tokens_in":21907,"tokens_out":6765,"would_cite":true,"duration_ms":63654,"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":"Adding a Stark term to the quantum Rabi model makes odd-order k-photon transitions selective, with resonance frequencies controlled by the boson number; the paper derives this from perturbation theory and shows a trapped ion can simulate…","keywords":["quantum Rabi model","Rabi–Stark model","selective interactions","multi-photon transitions","ultrastrong coupling","time-dependent perturbation theory","trapped-ion quantum simulation","Fock states"],"falsifier":"Scan $\\omega_0$ in a numerical simulation of the full Rabi–Stark Hamiltonian at $g/\\omega=0.1$, $\\gamma/\\omega=0.9$, with initial state $|g,7\\rangle$ and evolution time $t=\\pi/(2\\Omega^{(5)}_{2-})$. If the population of $|e,2\\rangle$ does not peak near $\\omega_0/\\omega=-3.227$ while $|g,6\\rangle$ and $|e,1\\rangle$ stay out of resonance, the effective five-photon Hamiltonian is not capturing the dynamics. The paper itself reports that at $g/\\omega\\approx 0.3$ the simple Jaynes–Cummings-type transfer degrades and population leaks to $|g,N_0\\pm1\\rangle$, so observing that leakage at modestly larger coupling would mark the regime where the perturbative selective-interaction claim breaks down.","tokens_in":20845,"feed_emoji":"⚛️","tokens_out":13376,"duration_ms":122641,"temperature":0.7,"pith_summary":"The paper tries to establish that the quantum Rabi model with a Stark coupling term—the Rabi–Stark model—produces selective $k$-photon interactions in the strong and ultrastrong coupling regimes. The central claim is that the rotating and counter-rotating parts of the linear Rabi coupling jointly generate odd-order multi-photon transitions, and that the Stark term makes each transition resonant for a chosen boson Fock state while neighboring Fock states stay out of resonance. The authors derive this from time-dependent perturbation theory, obtaining an effective Hamiltonian whose $k$-photon strength scales as $(g/\\omega)^k$ and whose resonance frequency depends on the Fock index through the Stark shift. They check the prediction numerically for three- and five-photon transitions and propose a single trapped ion as a simulator of the full model. If correct, this turns the Rabi–Stark model into a way to address one Fock state at a time, with direct use for bosonic state preparation and reconstruction.","feed_headline":"Stark term makes Rabi-model photons interact one Fock state at a time","feed_subtitle":"In strong and ultrastrong coupling, each k-photon transition fires only for one chosen boson number, and a trapped ion can simulate it.","key_machinery":"The central machinery is the interaction-picture Hamiltonian $H_I(t)=\\sum_n\\Omega_n(\\sigma_+ e^{i\\delta^+_n t}+\\sigma_- e^{i\\delta^-_n t})|n+1\\rangle\\langle n|+\\mathrm{H.c.}$, with $\\Omega_n=g\\sqrt{n+1}$ and Stark-shifted detunings $\\delta^\\pm_n=\\omega\\pm[\\omega_0+\\gamma(2n+1)]$. The nonzero $\\gamma$ makes these detunings Fock-state-dependent, which is what converts a degenerate multiphoton ladder into a selective one. A Dyson-series expansion of this Hamiltonian produces, at odd order $k$, the effective $k$-photon Hamiltonian of Eq. (5) with amplitudes $\\Omega^{(k)}_{n\\pm}\\propto (g/\\omega)^k$ and resonance frequencies $\\delta^{(k)}_{n\\pm}=(k-1)\\omega+\\delta^{\\pm}_{n+(k-1)/2}$; even-order terms are either diagonal (inducing a small Stark-like shift of the resonance) or average away under the rotating-wave approximation by parity symmetry. This effective-Hamiltonian construction is what carries the argument, because it converts the dynamical phenomenon of selective $k$-photon exchange into a concrete resonance condition and a calculable interaction strength.","core_discovery":"Adding a diagonal Stark term $\\gamma a^\\dagger a\\sigma_z$ to the quantum Rabi model changes the physical content of the model: for odd $k$, the combined rotating and counter-rotating terms generate effective $k$-photon Jaynes–Cummings/anti-Jaynes–Cummings interactions of the form $H_I^{(k)}(t)=\\sum_n(\\Omega^{(k)}_{n+}e^{i\\delta^{(k)}_{n+}t}\\sigma_+ + \\Omega^{(k)}_{n-}e^{i\\delta^{(k)}_{n-}t}\\sigma_-)|n+k\\rangle\\langle n|+\\mathrm{H.c.}$, with $\\Omega^{(k)}_{n\\pm}\\propto (g/\\omega)^k$ and with resonance condition $\\delta^{(k)}_{n\\pm}=(k-1)\\omega+\\delta^{\\pm}_{n+(k-1)/2}=0$, where the index-dependent piece contains $\\omega_0+\\gamma(2(n+(k-1)/2)+1)$. Consequently a specified initial Fock state $|g,N_0\\rangle$ or $|e,N_0\\rangle$ can be brought into resonant $k$-photon exchange with $|e,N_0+k\\rangle$ or $|g,N_0+k\\rangle$ while other Fock states remain far from resonance. Even-$k$ transitions are absent by parity symmetry, and second-order terms act as a Stark-like energy shift that moves the resonance slightly away from its naive value. Numerical evolution of the full Hamiltonian confirms selective three-photon and five-photon population transfer in the strong and ultrastrong regimes, and the trapped-ion derivation shows how this model can be realized with one ion and three laser drivings.","pith_inferences":["If the same selectivity survives with a nonperturbatively corrected resonance condition, the largest practical payoff would be Fock-state-dependent operations: one could prepare, rotate, or read out a specific boson number state by choosing $\\gamma$ and $\\omega_0$, something the paper motivates through state preparation but does not develop into gates.","The reported five-photon peaks deviate from the second-order analytic values by roughly $0.1$–$0.2\\,\\omega$, suggesting the perturbative formulas are a starting point; a resummed or exact-spectrum version of the resonance condition would be needed to predict peak positions as $g/\\omega$ grows past $0.1$.","A clean experimental test suggested by the parity argument is to scan for even-$k$ resonances: observing a two- or four-photon transition at finite $g$ would contradict the parity-based cancellation, while observing only odd orders would support the effective-Hamiltonian picture.","Since the trapped-ion implementation can choose $\\gamma$ of either sign and can reach $g/\\omega$ beyond $0.1$, it could also probe where the selective Jaynes–Cummings-type description breaks down, for instance by checking the paper's own expectation that at $g/\\omega\\approx 0.3$ population leaks to neighboring Fock states."],"forward_implications":["For each odd $k$, the Rabi–Stark model has selective $k$-photon resonances whose frequencies are shifted by the Stark term, so a chosen Fock state can be addressed without disturbing its neighbors.","The one-photon case inherits the same selectivity: setting $\\omega-\\omega_0=\\gamma(2N_0+1)$ makes only the doublet $\\{|e,N_0\\rangle,|g,N_0+1\\rangle\\}$ resonant, recovering earlier selective Jaynes–Cummings physics.","The strength of $k$-photon processes decreases as $(g/\\omega)^k$, so higher-order transitions require longer times; the numerical results show three-photon exchange at $g/\\omega=0.1$ and five-photon exchange with partial transfer at the same coupling.","Even-$k$ multiphoton transitions are suppressed by parity, so the observable multiphoton spectrum of the Rabi–Stark model is organized by odd orders.","A single trapped ion driven by two sideband fields plus a carrier field reproduces the Rabi–Stark Hamiltonian with independently tunable $\\omega_0$, $g$, and $\\gamma$, allowing the selective interactions to be observed in the strong and ultrastrong coupling regimes."],"supporting_citations":[{"why":"It introduces the Rabi–Stark model and provides the physical context for the Stark coupling.","marker":"[16, 17]"},{"why":"It supplies the spectral and symmetry structure of the Rabi–Stark model used to identify allowed odd-order transitions and spectral-collapse points.","marker":"[15]"},{"why":"It gives the energy-spectrum solutions of the Rabi–Stark model that support the claim that two-photon terms matter only near spectral collapse.","marker":"[19]"},{"why":"They demonstrate multiphoton resonances in the linear quantum Rabi model, the phenomenon this paper generalizes with Stark selectivity.","marker":"[35, 36]"},{"why":"It supplies the Dyson-series integration formula used to derive the second- and third-order effective Hamiltonians and their generalization.","marker":"[55]"},{"why":"They establish the selective one-photon Stark-interaction scheme in cavity QED that this work extends to the multiphoton Rabi–Stark setting.","marker":"[29, 31]"},{"why":"They provide the trapped-ion realizations and proposals of the quantum Rabi model that the simulation scheme builds on.","marker":"[49–51]"},{"why":"It provides the trapped-ion carrier and sideband Hamiltonian used as the starting point for the derivation.","marker":"[2]"}],"fun_headline_variants":["Stark term makes Rabi photons pick one Fock state per jump","Selective k-photon transitions from Stark-shifted Rabi model","Stark term enables Fock-state-selective photon exchange","Rabi model with Stark term: each resonance targets one n","k-photon interactions that fire for a single boson number"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the expansion in successive interactions, keeping only low orders and discarding fast-oscillating terms, stays accurate at the ultrastrong couplings studied here, so the analytic resonance condition locates the true multiphoton peak.","fun_headline_variants_meta":{"raw":{"variants":["Stark term makes Rabi photons pick one Fock state per jump","Selective k-photon transitions from Stark-shifted Rabi model","Stark term enables Fock-state-selective photon exchange","Rabi model with Stark term: each resonance targets one n","k-photon interactions that fire for a single boson number"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000201,"raw_usage":{"total_tokens":1418,"prompt_tokens":1025,"completion_tokens":393,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":641,"completion_tokens_details":{"reasoning_tokens":305}},"tokens_in":641,"tokens_out":393,"duration_ms":4501,"temperature":1.0,"reasoning_tokens":305,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:20:34.293733+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Scan $\\omega_0$ in a numerical simulation of the full Rabi–Stark Hamiltonian at $g/\\omega=0.1$, $\\gamma/\\omega=0.9$, with initial state $|g,7\\rangle$ and evolution time $t=\\pi/(2\\Omega^{(5)}_{2-})$. If the population of $|e,2\\rangle$ does not peak near $\\omega_0/\\omega=-3.227$ while $|g,6\\rangle$ and $|e,1\\rangle$ stay out of resonance, the effective five-photon Hamiltonian is not capturing the dynamics. The paper itself reports that at $g/\\omega\\approx 0.3$ the simple Jaynes–Cummings-type transfer degrades and population leaks to $|g,N_0\\pm1\\rangle$, so observing that leakage at modestly larger coupling would mark the regime where the perturbative selective-interaction claim breaks down.","supporting_citations":[{"cited_title":"Eckle and H","cited_arxiv_id":null,"evidence_quote":"It supplies the spectral and symmetry structure of the Rabi–Stark model used to identify allowed odd-order transitions and spectral-collapse points."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the energy-spectrum solutions of the Rabi–Stark model that support the claim that two-photon terms matter only near spectral collapse."}],"review_version":1}