{"id":"fdd4d692-bf07-4196-bb29-ebd67feeca7f","arxiv_id":"2608.04649","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Time-resolved X-ray diffraction on laser-shocked MgO reveals nanosecond plasticity and a shift from {110} to {100} slip between 95 and 175 GPa, inferred via self-consistent crystal plasticity models.","lead":"Shock experiments on polycrystalline MgO at the European XFEL show that the ceramic deforms plastically within nanoseconds, and that the dominant slip system changes as pressure rises above 100 GPa. The result gives a direct view of how ceramics yield under extreme impacts, relevant for planetary interiors and protective materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inferred slip-system transition rests on an EVPSC model whose rate-independent CRSS assumption is acknowledged to be invalid for shock; at 175 GPa the fitted CRSS values are nearly equal, so the reported 45/55 activity split is fragile.","rationale":"The paper's direct evidence for nanosecond-timescale plasticity in shocked polycrystalline MgO is credible: time-resolved diffraction shows lattice-strain relaxation and texture evolution on the shock transit timescale, and the measured flow stresses are consistent with prior dynamic and static data. The more novel secondary claim, a pressure-mediated switch in dominant slip system between 95 and 175 GPa, is entirely mediated by the EVPSC interpretation. That interpretation is the weakest load-bearing pillar because the authors explicitly acknowledge that EVPSC was not developed for shock deformation, and the specific simplifications listed in Supplement C (axial strain path, constant P-T, rate-independent CRSS, 3000 spherical grains) are all questionable in the 1e6-1e7/s laser-shock regime. The fragility is compounded by the near-equality of the fitted CRSS values at 175 GPa and by the 1-sigma activity uncertainties, which make the inferred 55% {100} activity statistically indistinguishable from a 50/50 split. A rate-dependent re-fit is the decisive check: it directly tests whether the inferred crossover survives when the acknowledged assumption is relaxed. The reader's weakest-assumption identification points to the same model limitation, so I agree with that diagnosis. Because the concern is significant but not demonstrated to be fatal, the appropriate verdict remains CONDITIONAL: the paper should be accepted only if the authors either provide the rate-dependent sensitivity test or soften the 'dominant switch' wording to a strain-rate-dependent trend.","tokens_in":20018,"tokens_out":11945,"duration_ms":150206,"concrete_test":"Re-run the EVPSC optimization at 95 and 175 GPa using a rate-dependent viscoplastic law (e.g., power-law strain-rate sensitivity with n=10-20, or CRSS values from the PNG/MD calculations cited in ref. [22]) over the same CRSS grid, and compare the best-fit relative activities and the 95-175 GPa crossover. If the {100} activity at 175 GPa drops below 50% or the crossover shifts outside the 95-175 GPa bracket, the central transition claim is not robust to the acknowledged model assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The 'change in dominant slip system between 95 and 175 GPa' is not a direct diffraction observable; it is produced by the EVPSC inversion described in the main text section 'Elasto-viscoplastic self-consistent (EVPSC) models' and Supplement C. The authors state that 'EVPSC were not developed to model shock deformation' and then impose purely axial deformation, constant P-T, 3000 spherical grains, and a single strain-rate-independent CRSS per slip system. At laser-shock strain rates near 1e6-1e7/s, CRSS values are expected to be rate dependent, and the issue is acute at 175 GPa: the best-fit CRSS values for the two systems are 1.9 and 2.0 GPa, so the reported 45/55 activity split is a near-balance that could be flipped by a modest rate-induced change in the CRSS ratio. The fit also uses only four Q(hkl) values per P-T condition, and the accepted-grid 1-sigma uncertainties at 175 GPa (7%) make the two activities statistically indistinguishable. Thus the load-bearing inference of a completed switch to {100}<110> dominance is not robust even under the model's own assumptions, and the acknowledged model simplifications could shift the inferred crossover pressure.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports laser-driven shock compression experiments on polycrystalline MgO at three Hugoniot states (27, 95, and 175 GPa) with time-resolved X-ray diffraction at the European XFEL. The authors extract lattice strain parameters Q(hkl) and texture from 2D Rietveld analysis, observe an elastic overshoot followed by a lower flow-state differential stress, and use elasto-viscoplastic self-consistent (EVPSC) simulations to interpret the evolution of lattice strains and texture. The central claim is that MgO deforms plastically on nanosecond timescales and that the dominant slip system changes from {110}<110> to {100}<110> between 95 and 175 GPa, as inferred from the EVPSC activity ratios and the observed texture transition from [100] to [110] maxima.","tokens_in":20321,"tokens_out":3888,"duration_ms":37187,"significance":"If the nanosecond-timescale plasticity claim is robust, it is a significant advance: direct time-resolved diffraction evidence that a hard ceramic can flow plastically within a few nanoseconds under shock, with quantitative differential stress measurements. The independent texture evolution in Fig. 3 and Supplement E is a valuable observable that provides a check on the model. The paper also connects dynamic strength measurements to a broader literature across static and dynamic regimes. However, the slip-system transition, which is a headline claim, is more fragile than the plasticity claim itself; it depends on model inversions whose acknowledged simplifications and statistical uncertainties are not fully quantified.","major_comments":[{"comment":"The central claim of a completed switch to {100}<110> dominance is not statistically supported by the reported best-fit values. At 175 GPa/3000 K (Table S2), the best-fit CRSS values are 1.9(4) GPa for {110}<110> and 2.0(1) GPa for {100}<110>, and the corresponding relative activities are 45(7)% and 55(7)%. Because the two CRSS values agree within their 1σ errors and the activity uncertainties overlap the 50/50 line, the 'change in dominant slip system between 95 and 175 GPa' reported in the abstract and Fig. 4(c) is not a statistically discernible crossover. The authors should report a formal confidence interval for the activity ratio, or soften the claim to 'consistent with an increasing {100} activity' until more pressure points or a rate-dependent model are available.","section":"Fig. 4(c) and Table S2"},{"comment":"The slip activities are partly an inversion of the measured lattice strains because the CRSS values are fitted to the same Q(hkl) data that the EVPSC model then uses to compute the relative activities. The texture comparison is an independent check, but Fig. 3 and Supplement E report only visual agreement; the paper should quantify how well texture alone (without the residual-error weighting in Eq. S1) discriminates between the two slip systems, and report the sensitivity of the activities to the choice of the 1.25 E threshold.","section":"Supplement C, Eq. (S1)"},{"comment":"The authors correctly state that 'EVPSC were not developed to model shock deformation,' yet the model assumes purely axial deformation, constant P-T, 3000 spherical grains with affine interaction, and a strain-rate independent CRSS. Under laser-driven shock with strain rates near 1e6-1e7/s, CRSS is expected to be rate-dependent, and at 175 GPa the two fitted CRSS values are nearly equal; a modest rate-induced change in the CRSS ratio could flip the 45/55 activity split. The authors should provide a sensitivity analysis, e.g., varying the CRSS ratio within the accepted grid or including a power-law rate sensitivity, and show how the inferred crossover pressure shifts. Without such a test, the transition pressure range is not robust.","section":"Elasto-viscoplastic self-consistent (EVPSC) models"}],"minor_comments":[{"comment":"The term 'flow state' is used without a definition; please define it in the caption or main text, and clarify that the time axis is 'time before shock breakout' so that negative times are clear.","section":"Fig. 2 caption"},{"comment":"In the 175 GPa row, the P_VISAR column is blank; use an em dash and add a footnote explaining that the pressure was determined only from impedance matching because VISAR fringes disappeared.","section":"Table I"},{"comment":"The analysis procedure is attributed to Ref. [37], but the detailed methodology is described in Ref. [32], which is stated to be under review; please cite Ref. [32] explicitly at the first use of the Rietveld/MAUD procedure so readers can locate the full method.","section":"Main text, section 'The two-dimensional diffraction data...'"},{"comment":"The sentence describing the acceptance criterion is ambiguous: 'within 1.25 of the minimum E' should be written as E <= 1.25 * E_min or an equivalent explicit inequality.","section":"Supplement C, Eq. (S1)"},{"comment":"The text says the XFEL probe timing was varied in 0.5 to 1 ns increments, but the number of distinct time points per condition is not stated; please specify the total number of shots per condition used in the time series in Fig. 2.","section":"Main text, 'For each P–T condition...'"}],"recommendation":"major_revision","confidential_remarks":"The experimental dataset and the nanosecond-timescale plasticity observation are convincing and valuable. The main issue is that the slip-system transition is a headline claim that is statistically fragile under the model's own assumptions, and the authors' own caveat about EVPSC applicability to shock deformation makes this especially important. I would encourage the editor to request a rate-sensitivity or robustness analysis, or a moderated abstract, before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news here is experimental: in-situ X-ray diffraction of polycrystalline MgO under laser-driven shock at three Hugoniot points, resolving lattice strains, differential stress, and texture as functions of time. The observation that MgO reaches a plastic flow state on nanosecond timescales, with an elastic overshoot at shock entry, is directly supported by the diffraction data and is a solid new benchmark for dynamic strength in a ceramic. The texture evolution from a [100] to a [110] maximum with pressure is a clean qualitative fingerprint, and the comparison with static data and the discussion of prior dynamic measurements are fair and useful. Data are made available, and the MAUD/Rietveld treatment follows established practice. That part deserves real credit.\n\nThe soft spot is the EVPSC inversion, as the stress-test note says. The authors openly state that EVPSC was not built for shock loading, then impose purely axial deformation, constant P-T, 3000 spherical grains, and a strain-rate-independent CRSS per slip system. The two CRSS values are fitted to the same Q(hkl) data that the model then interprets as slip activity, so the inference is partly circular, even though the independent comparison with measured texture helps. The problem is sharpest at 175 GPa: the best-fit CRSS values for {110} and {100} are 1.9 and 2.0 GPa, yielding a 45/55 activity split with 1-sigma errors of 7%. Those two values are statistically indistinguishable, so the claim of a completed switch to {100} dominance is fragile. The pressure at this point also lacks VISAR confirmation; it comes from impedance matching alone, and the model simplifications could shift the crossover pressure by tens of GPa. What is solid is the qualitative trend toward more {100} activity with pressure, which matches static experiments.\n\nWho gets value from this: the dynamic compression community, planetary materials people, and anyone testing polycrystal plasticity models against shock data. The central measurement is a strong benchmark; the slip-system transition should be treated as a motivated hypothesis until tested with off-Hugoniot points, rate-dependent CRSS, or a more direct discrimination between the two slip systems.\n\nRecommendation: this deserves serious peer review. Send it. A good referee will ask for a sensitivity analysis on the CRSS grid and a more explicit statement of how the interpretation depends on the rate-independence assumption. The experimental core is valuable enough that these caveats should be weighed in revision, not used to reject.","headline":"Direct time-resolved XRD shows MgO flows plastically on nanosecond timescales; the inferred slip-system crossover is plausible but rests on a model-dependent fit that is not statistically robust at the highest pressure point.","tokens_in":21318,"tokens_out":1758,"would_cite":true,"duration_ms":22071,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["62.50.-p"],"model":"deepseek-v4-flash","headline":"Laser-shock X-ray snapshots show MgO flowing plastically within nanoseconds and switching slip systems above 100 GPa.","keywords":["magnesium oxide","shock compression","plasticity","slip systems","time-resolved X-ray diffraction","EVPSC","nanosecond dynamics","lower mantle rheology"],"falsifier":"Shock a polycrystalline MgO sample to an intermediate state near 135 GPa and measure its flow-state lattice strains and texture; if $\\{110\\}\\langle 110\\rangle$ slip still accounts for the majority of the EVPSC activity, the transition lies above 135 GPa, and if $\\{100\\}\\langle 110\\rangle$ dominates, it lies below, tightening or contradicting the claimed 95-175 GPa window.","tokens_in":19809,"feed_emoji":"💥","tokens_out":12409,"duration_ms":117726,"temperature":0.7,"pith_summary":"Polycrystalline magnesium oxide, a model ceramic and a major lower-mantle mineral, is shown to deform plastically within nanoseconds when driven by laser shock to pressures up to 175 GPa. Time-resolved X-ray diffraction snapshots, combined with elasto-viscoplastic self-consistent simulations, indicate that the material reaches a plastic flow state on nanosecond timescales and that the dominant slip system switches from $\\{110\\}\\langle 110\\rangle$ to $\\{100\\}\\langle 110\\rangle$ between 95 and 175 GPa. If correct, this means hard ceramics can flow fast enough to matter in high-velocity impacts, and it constrains the rheology of planetary interiors where ferropericlase carries much of the deformation.","feed_headline":"MgO turns plastic in nanoseconds under shock","feed_subtitle":"Flow-state X-ray data show the {100} slip system taking over between 95 and 175 GPa.","key_machinery":"The argument is carried by comparing time-resolved X-ray diffraction with elasto-viscoplastic self-consistent (EVPSC) simulations. From the diffraction images, lattice strain parameters $Q$ for the 111, 200, 220 and 311 reflections give the differential stress through $t = 6G\\langle Q\\rangle$, where $G$ is the pressure-temperature-dependent shear modulus, and orientation-dependent intensities give the texture. The EVPSC model, which treats each of 3000 spherical grains as an inclusion in a homogeneous anisotropic effective medium, is run under purely axial deformation to 10% strain at constant pressure and temperature, with a grid search over critical resolved shear stresses for the two slip systems; accepted solutions minimize a residual between modelled and measured $Q$ values and also reproduce the measured texture maximum. That double match is what lets the paper attribute the [100]-to-[110] texture shift to a change in dominant slip plane rather than to starting texture or arbitrary fitting.","core_discovery":"The paper's central claim is that, under laser-driven shock compression along the B1-phase Hugoniot, polycrystalline MgO undergoes plastic relaxation within the roughly 10 ns drive and that the controlling slip system changes with pressure. At 27(7) GPa / 500 K and 95(7) GPa / 1500 K the diffraction data show an elastic overshoot at shock entry, with differential stress rising to at least 7.4 and 6.9 GPa, before settling to flow-state values of 3.2(0.3) and 4.3(0.9) GPa; at 175(15) GPa / 3000 K only the flow state is captured, at 3.8(1.2) GPa. The flow-state texture's inverse-pole-figure maximum moves from [100] at 27 GPa to [110] at 175 GPa. In the EVPSC fits, the relative activity of $\\{110\\}\\langle 110\\rangle$ slip falls from 77% to 45% while $\\{100\\}\\langle 110\\rangle$ rises from 23% to 55%, placing the mechanism change between 95 and 175 GPa and implying that $\\{100\\}$ slip becomes dominant above roughly 100 GPa.","pith_inferences":["If the slip-system transition is strain-rate sensitive, as the authors suggest, extrapolating these shock results to mantle strain rates may shift the crossover pressure; a rate-dependent critical-resolved-shear-stress law would be needed before using the inferred activities directly in geodynamic models.","The fitted critical resolved shear stresses at 27, 95 and 175 GPa could be checked against atomistic simulations at the same pressures and strain rates; agreement would turn fitted parameters into physical quantities, while disagreement would point to non-axial strain or rate effects.","Repeating the measurement with a different drive duration or strain rate would map how the transition pressure depends on strain rate, separating thermally activated from athermal control of slip in MgO.","The flow-state differential stress at 175 GPa, 3.8(1.2) GPa, is higher than the sub-1 GPa strength inferred from interface instability growth; reconciling these two dynamic measures would clarify what each diagnostic actually measures."],"forward_implications":["Shock-compressed regions of polycrystalline MgO above roughly 100 GPa should be modelled with $\\{100\\}\\langle 110\\rangle$ as the active slip system rather than the ambient-pressure $\\{110\\}\\langle 110\\rangle$ system.","Plastic flow is reached within the ~10 ns laser drive, so hard ceramics should not be assumed to respond only elastically or brittly on these timescales.","The flow-state differential stresses of about 3-4 GPa up to 175 GPa provide direct dynamic-strength benchmarks comparable to static diamond-anvil and gas-gun values.","The elastic overshoot at shock entry means the yield stress at the shock front exceeds the steady flow stress, which matters for interpreting Hugoniot elastic limits in ceramics.","The same single-shot X-ray diffraction plus self-consistent modelling pipeline can be applied to other polycrystalline ceramics to map slip-system activity across pressure, temperature and strain rate."],"supporting_citations":[{"why":"Supplies the elasto-viscoplastic self-consistent modelling procedure and its calibration for periclase, used to convert lattice strains and textures into slip-system activities.","marker":"[43]"},{"why":"Provides first-principles pressure-temperature elastic constants and shear moduli used to extract differential stresses and to set EVPSC elasticity.","marker":"[44]"},{"why":"Reports static rotational diamond-anvil experiments mapping the temperature-pressure dependence of MgO slip systems, giving the static transition range the shock result is compared with.","marker":"[26]"},{"why":"Gives Peierls-Nabarro-Galerkin predictions of critical resolved shear stresses and of the {110}-to-{100} transition, the numerical benchmark for the CRSS trends.","marker":"[22]"},{"why":"Documents dominant {110}<110> slip in polycrystalline MgO up to 47 GPa at ambient temperature, the low-pressure baseline for the mechanism change.","marker":"[9]"},{"why":"Demonstrates time-resolved X-ray diffraction of shock-compressed iron including elastic overshoot before flow, the methodological template for interpreting the MgO stress transients.","marker":"[13]"},{"why":"Reports the low-viscosity, sub-1 GPa strength measurement for MgO at 175 GPa from interface instability growth, the dynamic result this work's flow stress is compared with.","marker":"[18]"},{"why":"Provides gas-gun dynamic yield strengths and Hugoniot data for polycrystalline MgO used here for pressure determination and low-pressure strength comparison.","marker":"[29]"},{"why":"Describes the quantitative single-pulse X-ray diffraction analysis that extracts lattice strain parameters, differential stress, and texture from the shock diffraction images.","marker":"[32]"}],"fun_headline_variants":["MgO's slip mechanism flips above 100 GPa in nanoseconds","Shock-compressed MgO changes slip system in nanoseconds","Nanosecond plastic flow in MgO reveals slip system switch","MgO plasticity shifts to {100} slip above 100 GPa","MgO deforms plastically in ns, slip system changes at ~100 GPa"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The interpretation rests on the EVPSC simplification of the shock: purely axial deformation of 3000 spherical grains to 10% strain at constant pressure and temperature, with a single strain-rate-independent critical resolved shear stress for each slip system, and the paper itself notes the model was not developed for shock deformation; if the real strain path is not axial or the critical stresses depend on strain rate, the inferred slip activities and the 95-175 GPa transition window could shift.","fun_headline_variants_meta":{"raw":{"variants":["MgO's slip mechanism flips above 100 GPa in nanoseconds","Shock-compressed MgO changes slip system in nanoseconds","Nanosecond plastic flow in MgO reveals slip system switch","MgO plasticity shifts to {100} slip above 100 GPa","MgO deforms plastically in ns, slip system changes at ~100 GPa"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001292,"raw_usage":{"total_tokens":5294,"prompt_tokens":983,"completion_tokens":4311,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":4218}},"tokens_in":599,"tokens_out":4311,"duration_ms":31668,"temperature":1.0,"reasoning_tokens":4218,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:48:02.396839+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Shock a polycrystalline MgO sample to an intermediate state near 135 GPa and measure its flow-state lattice strains and texture; if $\\{110\\}\\langle 110\\rangle$ slip still accounts for the majority of the EVPSC activity, the transition lies above 135 GPa, and if $\\{100\\}\\langle 110\\rangle$ dominates, it lies below, tightening or contradicting the claimed 95-175 GPa window.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the elasto-viscoplastic self-consistent modelling procedure and its calibration for periclase, used to convert lattice strains and textures into slip-system activities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides first-principles pressure-temperature elastic constants and shear moduli used to extract differential stresses and to set EVPSC elasticity."},{"cited_title":"Ishimori, S","cited_arxiv_id":null,"evidence_quote":"Reports static rotational diamond-anvil experiments mapping the temperature-pressure dependence of MgO slip systems, giving the static transition range the shock result is compared with."},{"cited_title":"Amodeo, P","cited_arxiv_id":null,"evidence_quote":"Gives Peierls-Nabarro-Galerkin predictions of critical resolved shear stresses and of the {110}-to-{100} transition, the numerical benchmark for the CRSS trends."},{"cited_title":"Merkel, H","cited_arxiv_id":null,"evidence_quote":"Documents dominant {110}<110> slip in polycrystalline MgO up to 47 GPa at ambient temperature, the low-pressure baseline for the mechanism change."},{"cited_title":"Merkel, S","cited_arxiv_id":null,"evidence_quote":"Demonstrates time-resolved X-ray diffraction of shock-compressed iron including elastic overshoot before flow, the methodological template for interpreting the MgO stress transients."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the low-viscosity, sub-1 GPa strength measurement for MgO at 175 GPa from interface instability growth, the dynamic result this work's flow stress is compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides gas-gun dynamic yield strengths and Hugoniot data for polycrystalline MgO used here for pressure determination and low-pressure strength comparison."},{"cited_title":"Ginestet, A","cited_arxiv_id":null,"evidence_quote":"Describes the quantitative single-pulse X-ray diffraction analysis that extracts lattice strain parameters, differential stress, and texture from the shock diffraction images."}],"review_version":1}