{"id":"86754f2c-2f22-410e-aa1d-2e5f3d0692c8","arxiv_id":"2504.18828","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A numerical time-dependent Schrödinger study predicts that the n-photon excitation probability of the 8 eV 229mTh isomer scales as intensity^n times pulse width squared in the perturbative regime, with saturation at higher intensities.","lead":"This paper calculates how many laser photons are needed to excite thorium-229 nuclei into their low-energy isomeric state, and how laser intensity and pulse length affect the yield. It finds simple scaling rules for 2-, 3-, and 4-photon excitation and shows that very intense, very long pulses saturate the excitation probability.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative curves, including the 2-photon saturation plateau near 0.3, are computed from a single magnetic-substate initial condition (Eq. 7) with no averaging over the six ground-state substates; an unpolarized 229Th sample would require an incoherent average that may shift the results.","rationale":"The paper's central claim is quantitative: P_f^(n) ∝ Γ_I^2 I^n in the leading-order regime, and high-order effects cap the 2-photon excitation at ~0.3. These numbers come from solving the TDSE with an initial condition that picks one basis state (Eq. 7). An experiment with a real 229Th sample starts from an unpolarized ensemble, i.e., a density matrix proportional to the identity on the ground manifold. The question is whether the single-state calculation represents that average. The angular-momentum algebra in Appendix A shows that U+ contains tensor components with M=±1 and different reduced matrix elements for ground-ground, isomer-isomer, and ground-isomer couplings. For a two-photon process, the amplitude is a sum over intermediate substates weighted by the appropriate 3j symbols and energy denominators. Because the isomeric manifold has J=3/2, the available intermediate and final M values depend on the initial M; for example, the stretched state |5/2> has fewer two-photon channels than |1/2>. Summing |A_{fi}|^2 over final M and photon helicities does not reduce to a constant independent of initial M by the standard 3j orthogonality, because that orthogonality requires summing over complete sets of magnetic quantum numbers, while here the final manifold is incomplete and the intermediate sums are weighted by different reduced matrix elements. Thus the reported curves are likely initial-state dependent. The paper does not provide an averaging prescription, nor does it state which i was used for the figures. This is directly load-bearing for the 0.3 plateau claim. The reader's weakest assumption identified precisely this issue; my analysis agrees. I also examined the Eq. (13) off-by-one S factor flagged by the reader. That factor is block-diagonal in the photon-number index λ and only multiplies each c^{(λ)} by a phase, so it does not affect probabilities |c^{(λ)}|^2; it is a typographical error rather than a numerical flaw. The initial-state averaging is the more substantive concern. The proposed test—rerunning the TDSE for each of the six initial M and averaging—would settle it.","tokens_in":14195,"tokens_out":20862,"duration_ms":201484,"concrete_test":"Reproduce the 2-photon excitation calculation for each of the six ground-state substates |G,M> with M=-5/2,-3/2,-1/2,+1/2,+3/2,+5/2 using the same 10x10 M1 matrix, the same Gaussian pulse with I=10^24 W/cm^2 and ω=E_I/2, and the same numerical parameters as Fig. 4(a); average the six final probabilities with equal weights. Compare this averaged P_f^(2)(Γ_I) curve with the single-state curve that generated Fig. 4(a). If the averaged saturation plateau differs from ~0.3 by more than 20%, or the pulse duration at which high-order effects set in shifts by more than a factor of 2, the reported quantitative predictions do not apply to an unpolarized 229Th target. If the two curves coincide to within those tolerances, the concern is resolved and the single-state result is representative.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II.B and Eq. (7) fix the initial state to one basis vector |φ_i>, and P_f^(n)=|c^(n)|^2 is reported for that single state. The paper never specifies which M substate is used or whether results are averaged over the 6-fold ground manifold. A physical unpolarized 229Th target is an incoherent mixture of all M substates, so the measured excitation probability should be a weighted average over initial M and sum over final M. The 10x10 M1 matrix (Appendix A) contains within-manifold couplings proportional to μ_G and μ_I as well as ground-isomer couplings, so the two-photon amplitude involves different intermediate-state sums for different initial M. The 3j angular-momentum factors make the total transition strength from a given M_G to all allowed M_I dependent on M_G because the J_I=3/2 manifold truncates the final-state sum (e.g., from M_G=5/2 only M_I=1/2 and 3/2 can be reached with two σ photons, while from M_G=1/2 additional channels open). Therefore the leading-order scaling and the high-order saturation plateau (Figs. 4 and the conclusion's ~0.3 limit) may be specific to the chosen initial substate. If the initial state was chosen to maximize the coupling, the reported curves overestimate the unpolarized yield; if it was a typical state, the plateau position and height may still shift. The statement in Sec. III.A that 'the nucleus is not polarized' sets the beam axis as z but does not define the statistical operator. This gap directly affects the central quantitative claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a theoretical study of degenerate multi-photon excitation of the 8-eV isomeric state 229mTh from the 229Th ground state in the direct laser-nucleus interaction. The authors solve the time-dependent Schrödinger equation with a numerical method that sorts contributions by the net photon absorption number, and they use a 10x10 M1 interaction matrix built from measured B(M1), magnetic moments, and excitation energy. They report that, for short pulses, the n-photon excitation probability scales as the square of the FWHM pulse duration and the n-th power of the peak laser intensity, and that for long pulses high-order effects impose an intensity-independent saturation plateau near 0.3 for the 2-photon case, while the 3- and 4-photon limits depend on intensity. The paper argues that current and near-future high-intensity lasers could produce non-negligible isomer populations.","tokens_in":14626,"tokens_out":15684,"duration_ms":144617,"significance":"If the results hold, the paper provides a falsifiable prediction: degenerate two-photon excitation of 229mTh could reach probabilities on the order of 0.1 at 10^24 W/cm^2 with picosecond pulses, and the n-photon scaling law is a clean, parameter-free consequence of the interaction Hamiltonian. A notable strength is that no fitting to the target result is performed: the inputs, B(M1; I→G), the magnetic moments, and the excitation energy, are taken from earlier experiments. The photon-number-resolved TDSE method is a useful conceptual tool, provided the fast-propagation identity is correct. The central short-pulse scaling is analytically sound, but the quantitative saturation curves currently rest on an unverified numerical identity and an unspecified magnetic-substate averaging procedure, so the paper's quantitative claims are not yet fully substantiated.","major_comments":[{"comment":"The initial condition in Eq. (7) fixes a single basis state |φ_i⟩ and P_f^(n) is computed as |c^(n)|^2 for that state, but the index i is never specified and no average over the sixfold ground-state manifold is performed. The statement in Sec. III.A that \"the nucleus is not polarized\" only sets the beam axis as z; it does not define a statistical operator. Because the 10×10 M1 matrix in Appendix A contains within-manifold couplings and the 3j symbols make the two-photon amplitude depend on M_G, an unpolarized sample requires the incoherent average over initial M_G and a sum over final M_I. The quantitative curves and the ≈0.3 saturation plateau in Fig. 4 and the Conclusions may therefore be specific to a single magnetic substate; please state which state was used and provide the ensemble-averaged results or a quantitative bound on the difference.","section":"Sec. III.A, Eq. (7), Figs. 2–8"},{"comment":"The fast-propagator identity is algebraically incorrect as printed. From the stated recurrence T_{I,i+1}=S T_{I,i} S^{-1} one obtains T_{I,i}=S^{i-1}T_{I,1}S^{-(i-1)}, so ∏_{i=1}^n T_{I,i}=S^{n-1}(T_{I,1}S^{-1})^{n-1}T_{I,1}, which is not equal to S^{n+1}(S^{-1}T_{I,1})^n in general. The printed formula also fails to reduce to T_{I,1} when n=1, contradicting the statement that the approximation vanishes for n=1. Since this identity is the basis of the long-pulse saturation calculations (Figs. 3, 4, 6, 8), the derivation must be supplied and the numerical results must be either recomputed with the corrected expression or shown to be unaffected.","section":"Sec. II.B, Eq. (13)"},{"comment":"The two-step approximation and the photon-number truncation (λ_high=12–14, λ_low=−10) are central to the reported numbers, but no quantitative convergence test is provided; the sentence \"By varying the parameters, we ensure the convergency\" is not a numerical error statement. Please report, for at least one representative case, the dependence of P_f^(n) on δt, N, n, and the photon-number cutoffs.","section":"Secs. II.B and III.C"}],"minor_comments":[{"comment":"The phrase \"isomeric state isomeric state\" is duplicated; also clarify that the \"two-level system\" includes the magnetic substates of the ground and isomeric levels, since the numerical basis is 10-dimensional.","section":"Sec. III.A"},{"comment":"The figure captions contain LaTeX artifacts such as \"I /greaterorequalslant 10^20 W/cm^2\"; these should be typeset correctly.","section":"Figs. 2, 5, 7 captions"},{"comment":"The displayed definition of S is incomplete: the block structure should be written explicitly, e.g., each diagonal block is \\tilde{S} e^{-iλωδt} with λ ranging from λ_high to λ_low.","section":"Eq. (12)"},{"comment":"References [55] and [64] contain corrupted author names (\"P/suppress l´ ociennik\" and \"G/suppress lowacki\"); these should be corrected.","section":"References"},{"comment":"The paper should add an explicit limitations paragraph listing the assumptions of a bare nucleus in vacuum, no atomic electrons or hyperfine coupling, and no decoherence, since these are relevant to the experimental interpretation of the predicted yields.","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The central short-pulse scaling is credible and the paper is not circular, but the quantitative long-pulse results depend on an apparently incorrect fast-propagation identity and on an unspecified initial magnetic substate. Please ask the authors to correct the algebra in Eq. (13), state the initial state, and provide ensemble-averaged results; the revised manuscript can then be evaluated on its numerical and experimental claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper is worth a read if you are tracking 229mTh excitation, but treat the quantitative long-pulse curves as provisional. The core scaling law—P_f^(n) ∝ I^n Γ_I^2 for short pulses—is exactly what lowest-order perturbation theory gives, and the authors correctly say so. The new content is the photon-number-resolved TDSE implementation and the first scan of 2-, 3-, and 4-photon excitation probabilities with realistic laser parameters. The parameter inputs (B(M1), moments, energy) come from prior experiments; there is no fitting to the target result. That part is clean.\n\nThe soft spots are two. First, Eq. (13) is not a minor off-by-one: the claimed identity T_{I,n}...T_{I,1} = S^{n+1}[S^{-1}T_{I,1}]^n does not follow from T_{I,i+1}=S T_{I,i} S^{-1}. Even for n=2 the two sides differ by a factor of S, and the statement that the approximation vanishes for n=1 is plainly false for the printed formula. The long-pulse saturation curves (Figs. 4, 6, 8) rely on this fast-propagation trick, and no code or convergence data are provided. As printed, the numerical method is not checkable.\n\nSecond, the initial state is underspecified. Eq. (7) picks one basis vector |φ_i> among the six ground substates, and P_f^(n) is reported for that single state. The appendix says the nucleus is “not polarized”, but that does not define the statistical operator. A real unpolarized sample requires averaging over initial M_G and summing over final M_I; the 3j factors do make the total two-photon strength from M_G=5/2 differ from M_G=1/2 because the J=3/2 manifold truncates the final-state sum. So the ~0.3 saturation limit in the conclusion may be specific to the chosen substate. If the authors intended a single-substate calculation, they should say so; if they meant an unpolarized target, they need the incoherent average.\n\nThe bare two-level nucleus in vacuum is a stated idealization, so I would not count that against them.\n\nThis deserves a serious referee—the topic is timely and the parameter scan is genuinely new—but the referee should demand a derivation or correction of Eq. (13), a specification of the initial substate with a proper unpolarized average (or an explicit single-substate disclaimer), and convergence tests. I would condition acceptance on those.\n\nBest.","headline":"Useful parameter scan for multiphoton 229mTh excitation, but the long-pulse saturation claims rest on an unproven and apparently incorrect numerical identity, and the initial magnetic-substate averaging is unspecified.","tokens_in":15095,"tokens_out":14160,"would_cite":false,"duration_ms":131788,"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":"Two-photon lasers can excite the thorium-229m nuclear isomer","keywords":["multiphoton excitation","thorium-229 isomer","laser-nucleus interaction","time-dependent Schrödinger equation","M1 nuclear transition","nuclear clock","high-intensity laser pulses"],"falsifier":"Measure the 229mTh yield from a bare-nucleus target irradiated by a short, intense pulse at half the isomer energy; if the two-photon probability does not scale as $\\Gamma_I^2 I^2$ in the low-probability regime, or if the saturation limit is found to exceed about 0.3, the central scaling claim would be contradicted.","tokens_in":13984,"feed_emoji":"⚛️","tokens_out":8962,"duration_ms":77158,"temperature":0.7,"pith_summary":"Multi-photon excitation is usually discussed for atoms and molecules; this paper asks whether a bare nucleus can absorb several laser photons at once and reach an excited isomer. The authors study thorium-229, whose first excited state, the 229m isomer, lies only 8.36 eV above the ground state, and solve the time-dependent Schrödinger equation with a scheme that sorts the wavefunction by how many photons have been absorbed. They find that for short, intense pulses the n-photon excitation probability follows a clean power law: proportional to the square of the pulse duration and to the n-th power of the peak intensity. For long pulses, laser-induced photon emission interrupts the climb, and the two-photon excitation probability saturates near 0.3 no matter how high the intensity is. If these numbers hold, current and near-future high-intensity lasers might create measurable amounts of 229mTh directly in a bare-nucleus target, a new path toward studying laser-nucleus interactions and nuclear clocks.","feed_headline":"Two-photon lasers can excite the thorium-229m nuclear isomer","feed_subtitle":"Near-future lasers could populate the 8-eV nuclear isomer; yield scales as pulse width squared times intensity squared.","key_machinery":"The computational engine is a time-dependent Schrödinger equation solver in which the wavefunction is split into components labelled by the net number of absorbed photons, so that the n-photon excitation amplitude can be followed separately. The transfer matrix for one time step is tridiagonal in this photon-number ladder, and a two-step approximation factors a long Gaussian pulse into intervals, reducing the number of matrix multiplications from linear to logarithmic; this makes nanosecond-scale pulses tractable. The laser-nucleus coupling is treated in the multipole expansion, and for 229Th the M1 transition dominates over E2, giving a 10-by-10 coupling matrix between the six ground and four isomeric magnetic substates. The identity that carries the argument is the simple scaling $P_f^{(n)} \\propto \\Gamma_I^2 I^n$ in the leading-order regime, which the numerics establishes and which the authors present as the key rule for estimating isomer yield.","core_discovery":"The central claim is that degenerate multi-photon absorption can drive the 229mTh isomer in the direct laser-nucleus interaction, with a calculable probability. In the leading-order regime where the process is dominated by the net absorption of n photons and emission is negligible, the final excitation probability satisfies $P_f^{(n)} \\propto \\Gamma_I^2 I^n$, with $\\Gamma_I$ the FWHM pulse duration and $I$ the peak intensity. This scaling is verified numerically for the 2-, 3-, and 4-photon cases under the assumption of an isolated two-level nucleus coupled by the M1 transition between the $J=5/2$ ground and $J=3/2$ isomer manifolds. At long pulse durations, higher-order effects become significant and the excitation probability saturates; for the 2-photon case the saturation value is approximately 0.3 across intensities from $10^{20}$ to $10^{24}\\,\\text{W/cm}^2$, while for 3- and 4-photon cases the maximum grows with intensity. The paper further argues that these results imply that high-intensity short-pulse lasers, rather than low-intensity long-pulse ones, are the promising regime for experimental multiphoton isomer production.","pith_inferences":["Outside the paper: an incoherent density-matrix treatment over the ground magnetic substates is the natural next step; if the average differs from the single-state result, the quoted 0.3 saturation should be reinterpreted as an upper bound for polarized targets only.","Outside the paper: the same photon-number-resolved TDSE method could be applied to other low-lying nuclear isomers or to transitions in highly charged ions where the nucleus and shell interact, possibly revealing enhanced multiphoton rates via electronic bridge mechanisms.","Outside the paper: if the scaling law is confirmed experimentally, the $^{229\\text{m}}$Th yield could serve as a diagnostic for the peak intensity and temporal shape of extreme laser pulses, since the probability depends sensitively on both $\\Gamma_I^2$ and $I^n$."],"forward_implications":["At intensities around $10^{23}\\,\\text{W/cm}^2$ and pulse durations of a few hundred femtoseconds, the two-photon channel could produce a non-negligible $^{229\\text{m}}$Th population in a bare-nucleus target, making direct laser-nucleus multiphoton excitation experimentally testable.","The scaling law $P_f^{(n)} \\propto \\Gamma_I^2 I^n$ gives a simple calibration rule: once one point is measured, the n-photon yield at other intensities and pulse lengths is predicted.","The 3- and 4-photon channels are strongly suppressed ($I^3$ and $I^4$) and will likely remain below detection except at the very highest planned intensities.","The saturation at about 0.3 for two-photon excitation means that simply raising intensity or extending pulse length cannot push the isomer population above that in this model; different mechanisms, such as involving atomic electrons or multiple pulses, would be needed.","For nuclear clock applications, the authors advise low-intensity, narrow-linewidth long pulses, but note that practical issues such as linewidth and control remain unsolved."],"supporting_citations":[{"why":"Supplies the original multiphoton excitation theory that the monochromatic limit of the present method recovers.","marker":"[1]"},{"why":"Provides the adopted isomer excitation energy, lifetime, and M1 reduced transition probability used as input for the calculations.","marker":"[33]"},{"why":"Documents the record laser intensity that frames the short-pulse parameter space considered.","marker":"[39]"},{"why":"Surveys current and planned high-power laser facilities, including long-pulse systems, used to select the parameter ranges.","marker":"[42]"},{"why":"Gives the ground-state magnetic dipole moment of 229Th used to build the M1 coupling matrix within the ground manifold.","marker":"[61–63]"},{"why":"Provides the ratio of the isomer to ground magnetic dipole moments, fixing the isomer magnetic moment in the matrix.","marker":"[64]"}],"fun_headline_variants":["Multi-photon lasers can trigger thorium-229m isomer","High-intensity lasers could excite thorium-229m via multi-photon","Multi-photon absorption drives thorium-229m isomer","Laser pulses can directly populate thorium-229m isomer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation starts from a single, fixed magnetic substate of the 229Th ground state, and the reported probabilities assume that this represents the real excitation probability; an unpolarized sample would need an incoherent average over ground substates, and no such average is given.","fun_headline_variants_meta":{"raw":{"variants":["Multi-photon lasers can trigger thorium-229m isomer","High-intensity lasers could excite thorium-229m via multi-photon","Multi-photon absorption drives thorium-229m isomer","Laser pulses can directly populate thorium-229m isomer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001063,"raw_usage":{"total_tokens":4477,"prompt_tokens":982,"completion_tokens":3495,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":598,"completion_tokens_details":{"reasoning_tokens":3423}},"tokens_in":598,"tokens_out":3495,"duration_ms":23027,"temperature":1.0,"reasoning_tokens":3423,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:10:40.587648+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 229mTh yield from a bare-nucleus target irradiated by a short, intense pulse at half the isomer energy; if the two-photon probability does not scale as $\\Gamma_I^2 I^2$ in the low-probability regime, or if the saturation limit is found to exceed about 0.3, the central scaling claim would be contradicted.","supporting_citations":[{"cited_title":"A new scheme for isomer pumping and depletion with high-power lasers","cited_arxiv_id":"2404.07909","evidence_quote":"Provides the adopted isomer excitation energy, lifetime, and M1 reduced transition probability used as input for the calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Surveys current and planned high-power laser facilities, including long-pulse systems, used to select the parameter ranges."},{"cited_title":"Minkov and A","cited_arxiv_id":null,"evidence_quote":"Provides the ratio of the isomer to ground magnetic dipole moments, fixing the isomer magnetic moment in the matrix."}],"review_version":1}