{"id":"6adc72f1-c092-4606-9933-690be6eeb439","arxiv_id":"2505.24444","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Simulating the repeated-application steady state of pulsed DNP, rather than a single polarization transfer, markedly improves agreement with experiments and reveals that the steady orbit differs from the single-transfer trajectory.","lead":"This paper simulates pulsed dynamic nuclear polarization by directly computing the stroboscopic steady state of repeated microwave pulse sequences, and shows that this matches experiments better than the usual single-transfer analysis. It validates the method against data at 0.34 T, 1.2 T, and new 3.4 T measurements.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"W-band agreement is fit-defined: Rmax is adjusted to the same field profile used as validation; a leave-one-out refit would determine if the model predicts or absorbs the data.","rationale":"No significant technical error is found in the steady-state algorithm itself. The Newton-Raphson formulation of Eq. (8) is well posed, and the stability analysis in Eq. (9) is plausible. Propagator squaring is stable. The paper honestly reports that the effective-Hamiltonian method is unstable. Credit should be given for these contributions. The conditionality arises from the strength of the experimental validation. The central claim is that steady-state simulations reproduce experimental field profiles, optimal repetition times, Rabi-frequency dependencies, and transfer curves. The weakest link is that some of the parameters used in those comparisons are adjusted to the same data. Specifically, Rmax is described both as a convergence threshold and as an adjustable parameter, creating an ambiguity: if it is fully converged, the 'adjustment' statement is contradictory; if it is adjustable, the W-band field profile is not an independent test. Similarly, the Rabi ensemble ranges are hand-picked per sequence and per frequency. A leave-one-out calibration/prediction split would settle whether the agreement is predictive. The reader identified the same weakest assumption; I agree. I do not see a reason to change the verdict from CONDITIONAL, but the concern is significant enough that the paper should be revised to clarify the fitting status of these parameters.","tokens_in":1053,"tokens_out":5412,"duration_ms":129125,"concrete_test":"Perform a leave-one-out refit. Use only X-band NOVEL transfer times to set the electron-proton coupling, and use a radial-convergence scan of Eq. (14) before any experimental comparison to set Rmax; fix the Rabi ensemble bounds from the stated B1 conversion efficiencies rather than from observed profiles. Then predict the W-band field profile (Fig. 4), the W-band pulse-length scan (Fig. 5), and the Q-band repetition-time profiles (Fig. 6). If predicted relative peak intensities at +-110/+-170 MHz and the ordering and shapes in Figs. 5-6 shift by more than about 20% or lose the experimental features, the fitted parameters are absorbing the physics; if they survive, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The algorithmic core (Newton-Raphson steady state) is sound: the Jacobian condition-number bound in Sec. 2.3 is parameter-free, and the Spinach implementation is described with data deposit. The load-bearing concern is the validation of the claim of 'remarkably good' agreement. Section 3.3 states the electron-proton coupling 'was adjusted to reproduce the polarization transfer times observed in NOVEL and pulsed SE experiments at 0.34 T'; Section 4.2 states, for the W-band field profile, 'the largest distance in the ensemble (20 A) was adjusted to reproduce the relative peak intensities in the experimental field profile.' The Rabi-frequency ensemble ranges in Secs. 3.5/4.5 (10-20 MHz, 25-35 MHz, per-frequency ranges) are likewise selected rather than predicted. Sec. 3.5's statement that Rmax = 20 A 'was empirically found to converge the radial part of the average' is ambiguous: if the radial average is converged, Rmax is not a free knob; if it is a free knob, then the W-band comparison is partially circular. This needs disambiguation before 'remarkably good' can be read as predictive.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces three algorithms for computing the stroboscopic steady state of a periodically driven dissipative spin system, tailored to pulsed dynamic nuclear polarization (DNP). The first algorithm, propagator squaring, is stable but expensive; the second, based on matrix-logarithm effective Hamiltonians, is numerically unstable for dissipative DNP; the third, Newton-Raphson iteration on the fixed-point equation of the one-period propagator, is claimed to be stable and efficient. The authors implement these in Spinach, benchmark them on CPU and GPU, and compare steady-state simulations with experimental field profiles, repetition-time scans, polarization-transfer curves, and Rabi-frequency dependences for NOVEL, TOP, XiX, and TPPM DNP at 0.34 T, 1.2 T, and 3.4 T, including new W-band XiX DNP data. The central claim is that steady-state simulations reproduce experimental results remarkably well, substantially better than single-transfer simulations, and that the steady orbit differs from the one-time-transfer trajectory.","tokens_in":22281,"tokens_out":4727,"duration_ms":63926,"significance":"The algorithmic contribution is a genuine step forward: the Newton-Raphson steady-state method is clearly explained, the condition-number bound in Sec. 2.3 is parameter-free and system-size independent, and the comparison among three algorithms is useful for practitioners. The code and data are deposited in Spinach and Zenodo, which aids reproducibility. The new W-band XiX DNP experiments provide valuable reference data. The validation of the method against multiple pulse sequences and magnetic fields is ambitious, and the improvement over single-transfer simulations is demonstrated qualitatively. However, the predictive claim is weakened by several hand-adjusted parameters, notably the effective electron-proton coupling and the distance-ensemble cutoff, so the paper's central conclusion that agreement is 'remarkably good' and predictive needs stronger support.","major_comments":[{"comment":"The W-band validation is partially circular by the paper's own description. Sec. 4.2 states that 'the largest distance in the ensemble (20 A) was adjusted to reproduce the relative peak intensities in the experimental field profile,' while Sec. 3.5 says the same value 'was empirically found to converge the radial part of the average.' These statements are in tension: if 20 A is a converged integration cutoff, it is not a free parameter, but if it is adjusted to the data, the W-band comparison is a fit rather than a prediction. The authors should disambiguate this and ideally provide a leave-one-out or cross-validation analysis, e.g., fitting Rmax on one subset of the field profile and testing on the rest, to establish the predictive content of the model.","section":"Sec. 4.2 (W-band field profile) and Sec. 3.5"},{"comment":"The paper reports that the steady-state simulation of NOVEL with a flip-back pulse predicts nuclear polarization 'above the maximum theoretical enhancement factor of 329.' This is a red flag: the theoretical maximum enhancement is a fundamental bound, and exceeding it indicates either an inconsistency in the simulation model (e.g., in how the flip-back pulse or relaxation is treated) or a misstatement. The manuscript does not comment on this discrepancy, yet it appears in a figure supporting the claim of good agreement. The authors should investigate whether this is a modeling error, a unit/definition issue, or a genuine prediction that requires explanation.","section":"Sec. 4.3, Figure 7 (NOVEL with flip-back)"},{"comment":"The claim in Sec. 5 that 'the different pulsed-DNP experiments were successfully simulated with the same set of parameters' is overstated. Some parameters are fixed across experiments, but others are not: T1n,bulk is 26 s at X band but 52 s at Q and W band; the electron Rabi frequency ensemble ranges differ per experiment (e.g., 14-16 MHz for NOVEL at X band, 10-20 MHz for TOP/XiX at Q band, 25-35 MHz for TPPM); and the effective electron-proton coupling is fitted to 0.34 T transfer times. A precise inventory of which parameters are global constants, which are field-specific, and which are fitted to the data being compared would be necessary to support the 'same set of parameters' claim and to assess the degrees of freedom used in the validation.","section":"Sec. 5 and Secs. 3.3-3.5"},{"comment":"The agreement between simulation and experiment is assessed only visually. Given the number of adjustable parameters and the qualitative nature of the comparisons, the central claim of 'remarkably good' agreement would be much strengthened by a quantitative metric, such as normalized root-mean-square deviation between simulated and experimental field profiles, repetition-time curves, or transfer curves, ideally with a statement of which features are used for fitting and which are tested. Without such a metric, it is difficult to distinguish genuine prediction from overfitting, especially in the W-band field profile where Rmax is explicitly adjusted.","section":"Secs. 4.2-4.5, Figures 3-9"}],"minor_comments":[{"comment":"The abstract contains a typo: 'a greement is good' should be 'agreement is good.'","section":"Abstract and Sec. 4.2"},{"comment":"Equation (9) is garbled: the expression 'min 12i T rTr T e e' is not readable as printed. Please rewrite it with proper exponential notation and define all symbols in the surrounding text.","section":"Eq. (9)"},{"comment":"The solvent composition 'd₈-glycerol:D₂O:H₂O 60:30:10' would be clearer with volume/volume percentages explicitly stated, and the isotope subscript formatting should be consistent throughout.","section":"Sec. 3.1"},{"comment":"The phrase 'peaks above the maximum theoretical enhancement factor of 329' is ambiguous; if the plotted quantity is Iz, the maximum possible value is 0.5, so the reader cannot tell whether the simulation is exceeding a thermal-polarization enhancement bound or a different definition. Please clarify the ordinate and the theoretical bound.","section":"Sec. 4.3, Figure 7"},{"comment":"The benchmarks in Table S1 cover systems up to one electron and nine protons, but the main-text statement that Newton-Raphson is 'up to three times' faster than squaring is based on a single representative point at each dimension; reporting error bars or multiple benchmark points would make the performance claim more robust.","section":"Sec. 2.3 and Sec. 4.1"}],"recommendation":"major_revision","confidential_remarks":"The algorithmic core is sound and the paper is likely to become a useful methods contribution. However, the validation section contains a fit-defined parameter for the W-band profile and a statement about exceeding the theoretical maximum enhancement that should be resolved before acceptance. I would suggest the editor require a careful response on these two points, and a quantitative agreement assessment, before sending back to the authors for revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, the algorithmic core is solid and genuinely new: computing the stroboscopic steady state of a pulsed DNP experiment via Newton-Raphson, with a parameter-free conditioning bound and a Spinach implementation, is a practical tool that single-transfer simulations cannot match. Second, the \"remarkably good\" agreement with experiments is real but not entirely predictive, because several parameters are adjusted to the very data the paper compares against.\n\nWhat is actually new: the Newton-Raphson fixed-point method for the DNP period propagator, the numerical demonstration that the matrix-logarithm effective-Hamiltonian route is unstable in dissipative DNP settings, and the first XiX DNP data at 3.4 T. The performance benchmarks are clean and the GPU scaling is useful.\n\nWhat the paper does well: it validates the method across three magnetic fields (X, Q, W) and four pulse sequences (NOVEL, TOP, XiX, TPPM). The repeated-time profiles and Rabi-frequency dependencies are exactly the quantities that single-transfer simulations cannot touch, and the steady-state results reproduce them qualitatively and often semi-quantitatively. The authors are also transparent about many of their modeling choices—they state which relaxation parameters are guessed, and the Discussion acknowledges the model cannot yet predict absolute enhancements.\n\nThe soft spots are proportional. The Rmax issue flagged in the stress test is legitimate: Section 3.5 says 20 Å was \"empirically found to converge the radial part of the average,\" while Section 4.2 says the same 20 Å was \"adjusted to reproduce the relative peak intensities\" of the W-band field profile. Those two statements are in tension and the paper needs to disambiguate them. The electron-proton coupling is also fitted to 0.34 T transfer times, and the Rabi-ensemble ranges are selected rather than predicted. So the W-band comparison is partly circular, and the strongest claims should be softened. That said, the Q-band results were not fitted to Q-band data, and the parameter scan in Figure 5 is a genuine prediction of shapes and relative intensities. The circularity is a blemish, not a fatal flaw, because the method's value does not depend on perfect first-principles prediction.\n\nThis paper deserves a serious referee. The algorithmic contribution is important and the experimental breadth is impressive. I would ask the authors to clarify Rmax, separate fitted from predicted quantities, and add error estimates to the comparisons. The flaws are correctable in revision.","headline":"The steady-state algorithm is a real contribution, but the experimental validation is partly fit-defined, so treat \"remarkably good\" with caution.","tokens_in":22769,"tokens_out":2876,"would_cite":true,"duration_ms":37783,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65H10","15A16","81V55"],"pacs":["82.56.-b","33.25.+k","76.30.-v"],"model":"deepseek-v4-flash","headline":"Pulsed DNP simulations now match experiment by computing the stroboscopic steady state directly, with the Newton-Raphson method as the recommended algorithm.","keywords":["pulsed dynamic nuclear polarization","stroboscopic steady state","Newton-Raphson method","propagator squaring","matrix logarithm","NOVEL","XiX DNP","spin dynamics simulation"],"falsifier":"A concrete falsifier is an ab initio simulation of the same pulsed DNP experiments with a many-proton spin system (e.g., a realistic trityl radical surrounded by tens of protons at crystallographically or MD-derived positions) and independently measured relaxation parameters, rather than the calibrated single-proton coupling and adjusted ensemble cutoff. If such a simulation fails to reproduce the field profiles, repetition-time optima, or Rabi-frequency dependencies that the calibrated single-proton steady-state model reproduces, then the calibrations in Sections 3.3 and 3.5 were absorbing physics omitted from the model. A simpler quantitative check is to compute the W-band field profile with the maximum ensemble distance set to the physically expected value from the trityl proton distribution rather than the adjusted 20 Å and compare the relative peak intensities; if the peaks at ±60 and ±230 MHz remain too strong, the model's agreement was sensitive to the fitted cutoff.","tokens_in":21781,"feed_emoji":"🧲","tokens_out":3603,"duration_ms":32107,"temperature":0.7,"pith_summary":"This paper argues that the standard practice of simulating a single one-time polarization transfer is insufficient for pulsed dynamic nuclear polarization (DNP), and that the relevant quantity is the stroboscopic steady state reached after many repetitions of the microwave pulse sequence. It introduces and compares three algorithms for computing this steady state from the pulsed-DNP propagator, finding the Newton-Raphson method stable and most efficient. Steady-state simulations reproduce experimental field profiles, optimal repetition times, Rabi-frequency dependencies, and polarization transfer curves at X, Q, and W band better than single-transfer simulations. A key observation is that the spin trajectory during the first transfer differs from the steady orbit, which has implications for pulse-sequence design.","feed_headline":"Pulsed DNP simulations now match experiment by computing the steady state directly","feed_subtitle":"Newton–Raphson converges on the stroboscopic orbit that repeated pulse sequences actually reach, outperforming single-transfer simulations.","key_machinery":"The central object is the pulsed-DNP period propagator $P$ that takes the spin density matrix forward by one repetition time $T$, including the microwave pulse sequence and relaxation. The steady state is defined as the fixed point $\\rho_\\infty = P \\rho_\\infty$, with the trace constraint $\\mathrm{Tr} \\rho_\\infty = 1$. The advocated Newton-Raphson algorithm solves this fixed-point equation directly by setting up the residual $f(\\rho) = P\\rho - \\rho$, forming the Jacobian $J = P_{2:n,2:n} - \\mathbb{1}$, and iterating Newton steps until convergence. Its condition number is bounded by $2/(r_{\\min}T)$ where $r_{\\min}$ is the smallest relaxation rate, and this bound does not grow with system size, making the method numerically stable where the matrix-logarithm effective-Hamiltonian approach fails.","core_discovery":"The central claim is that the stroboscopic steady state of a pulsed DNP experiment can be computed accurately and efficiently with the Newton-Raphson method, and that such steady-state simulations reproduce experimental field profiles, optimal repetition times, Rabi-frequency dependencies, and polarization transfer curves better than the standard single-transfer simulations. The paper states: 'We have introduced the simulation of pulsed DNP in the steady state,' 'we recommend the Newton-Raphson method for its stability and complexity scaling,' and 'Agreement with experiments at X, Q, and W band ... is remarkably good.'","pith_inferences":["A testable extension is to apply the Newton-Raphson steady-state approach to adiabatic-passage sequences such as frequency-swept or ramped-amplitude NOVEL, which the authors expect to be equally useful.","The observation that steady orbits preserve the electron-proton dipolar coupling better than single-transfer analysis suggests that optimized sequences could exploit this mismatch to achieve efficient high-field pulsed DNP at lower peak microwave power.","The model's success with a single electron-proton pair, after calibrating the coupling and ensemble range against experiments, implies that the same procedure could be used to calibrate effective spin-system parameters for other polarizing agents and matrices, provided the calibration data exist.","The documented numerical instability of the matrix-logarithm method in dissipative Liouville space suggests that steady-state optimal-control sequence design will require either a stable effective-generator method or a workaround such as direct gradient-based optimization of the fixed-point equation."],"forward_implications":["The optimal repetition time of pulsed DNP sequences can be predicted, which was impossible with single-transfer simulations.","Steady-state simulations reproduce the shapes and relative intensities of field profiles and parameter scans at X, Q, and W bands, including the first XiX DNP data at 3.4 T.","The dependence of the enhancement factor on the electron Rabi frequency, including the contrast between XiX (roughly constant) and TOP (steadily increasing) DNP, is reproduced.","The steady orbit during repeated pulse application differs from the trajectory of the first transfer, so pulse-sequence design should target the steady orbit rather than a single transfer.","Ensemble averaging over electron-proton distances and electron Rabi frequencies improves agreement with experiments, indicating that these distributions matter for quantitative prediction."],"supporting_citations":[{"why":"Provides the Q-band TOP/XiX/TPPM experimental data and the single-transfer simulation methodology that this paper supersedes.","marker":"[22]"},{"why":"Provides the X-band NOVEL experimental data and the polarization transfer times used to calibrate the electron-proton coupling.","marker":"[19]"},{"why":"Provides the X-band experimental context and the bulk proton relaxation time used at X band.","marker":"[51]"},{"why":"Supplies the distance- and angle-dependent nuclear relaxation model used in the simulations.","marker":"[61]"},{"why":"Is the Spinach software library in which the algorithms are implemented and the simulations are run.","marker":"[43]"},{"why":"Supplies the matrix-logarithm method whose numerical instability is analyzed and compared against.","marker":"[37]"},{"why":"Supplies the dynamical-systems result that periodically driven dissipative linear systems have a unique steady orbit to which all initial conditions converge.","marker":"[36]"},{"why":"Supplies the GMRES iterative linear solver used for the Newton-Raphson inverse-times-vector steps.","marker":"[46]"},{"why":"Describes the W-band pulsed EPR spectrometer on which the new 3.4 T XiX DNP experiments were performed.","marker":"[52]"}],"fun_headline_variants":["Steady-state simulations match pulsed DNP experiments","Newton-Raphson computes pulsed DNP steady state efficiently","Pulsed DNP: Direct steady-state calculation beats single-shot analysis","Simulating pulsed DNP: steady state via Newton-Raphson matches data","Pulsed DNP simulations now reach steady state with Newton-Raphson"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The most fragile load-bearing premise is that a single electron coupled to a single proton, with the electron-proton coupling adjusted to reproduce NOVEL transfer times at 0.34 T, plus a simple literature relaxation model and hand-selected ensemble ranges, is sufficient to predict pulsed DNP behavior at other fields. This modeling choice enters in Section 3.3 ('The electron-proton coupling was adjusted to reproduce the polarization transfer times observed in NOVEL...') and in Section 3.5, where the maximum distance in the ensemble is 'adjusted to reproduce the relative peak intensities' of the W-band profile that is then compared to the simulation. If these adjustments are absorbing the physics that the model omits, the agreement with experiments could be partly coincidental.","fun_headline_variants_meta":{"raw":{"variants":["Steady-state simulations match pulsed DNP experiments","Newton-Raphson computes pulsed DNP steady state efficiently","Pulsed DNP: Direct steady-state calculation beats single-shot analysis","Simulating pulsed DNP: steady state via Newton-Raphson matches data","Pulsed DNP simulations now reach steady state with Newton-Raphson"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000974,"raw_usage":{"total_tokens":4127,"prompt_tokens":923,"completion_tokens":3204,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":3116}},"tokens_in":539,"tokens_out":3204,"duration_ms":27031,"temperature":1.0,"reasoning_tokens":3116,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:22:17.708421+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete falsifier is an ab initio simulation of the same pulsed DNP experiments with a many-proton spin system (e.g., a realistic trityl radical surrounded by tens of protons at crystallographically or MD-derived positions) and independently measured relaxation parameters, rather than the calibrated single-proton coupling and adjusted ensemble cutoff. If such a simulation fails to reproduce the field profiles, repetition-time optima, or Rabi-frequency dependencies that the calibrated single-proton steady-state model reproduces, then the calibrations in Sections 3.3 and 3.5 were absorbing physics omitted from the model. A simpler quantitative check is to compute the W-band field profile with the maximum ensemble distance set to the physically expected value from the trityl proton distribution rather than the adjusted 20 Å and compare the relative peak intensities; if the peaks at ±60 and ±230 MHz remain too strong, the model's agreement was sensitive to the fitted cutoff.","supporting_citations":[],"review_version":1}