{"id":"9e840296-56f1-43cc-bb40-10aa7605b31e","arxiv_id":"2607.25896","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"β'-Mn3(PO4)2 irreversibly loses long-range order above ~14 GPa; DFT reproduces its anisotropic compression and shows a mechanical instability at the same pressure.","lead":"Researchers compressed the phosphate crystal β'-Mn3(PO4)2 beyond 14 GPa and found it permanently loses crystalline order, with computer calculations showing the crystal becomes mechanically unstable at that pressure. The study provides new data on how complex phosphate frameworks fail under pressure.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The mechanical-instability claim at 14.3 GPa is undercut by an internally inconsistent elastic tensor: Table 2 lists symmetry-forbidden C45, Table 3 replaces it with C35, and the text misstates the ordering of C11/C22/C33. The negative eigenvalue may be an artifact.","rationale":"The reader's weakest assumption—fixed internal coordinates in Rietveld refinement—is a legitimate limitation of the quantitative EoS and compressibility analysis, but it does not threaten the primary qualitative observation of irreversible peak broadening. The elastic-instability claim is different: it is the paper's mechanistic explanation for why the disorder occurs, and it is explicitly highlighted in the abstract and conclusions. The elastic tensors in Tables 2 and 3 contain a symmetry inconsistency (C45 vs C35) and the text contradicts its own C11/C22/C33 ordering. Since the negative eigenvalue at 14.3 GPa is the sole computational evidence for the instability, an error in assembling the tensor could invalidate that central supporting claim. A targeted recalculation is straightforward and would settle the issue. I do not see this as requiring rejection of the experimental result, so the reader's CONDITIONAL verdict remains appropriate.","tokens_in":15705,"tokens_out":12247,"duration_ms":105794,"concrete_test":"Recompute the full elastic stiffness tensor for β'-Mn3(PO4)2 at 14.3 GPa using the Le Page method as implemented in VASP with the same DFT settings, explicitly verifying the monoclinic symmetry (nonzero C15, C25, C35, C46; C45 = 0). Apply the same pressure correction (C_ii → C_ii - P, C_12/13/23 → C_ij + P) and recalculate the eigenvalues of the resulting 6x6 matrix. If the minimum eigenvalue is negative with C35 correctly included, the mechanical-instability claim stands; if not, Table 3 and the Born-criteria statement must be revised. Additionally, check the ambient-pressure tensor with the same procedure to resolve the C45/C35 discrepancy.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's secondary claim—that the crystal becomes mechanically unstable near 14.3 GPa because one eigenvalue of the elastic tensor becomes negative—is the computational link between the observed irreversible disorder and an intrinsic elastic instability. This claim rests on the quality of the reported 6x6 elastic stiffness tensors, and those tensors contain internal inconsistencies.\n\nIn Table 2 (ambient pressure), the header lists 'C45' as an independent constant with value 0.4 GPa. For a monoclinic P2_1/c crystal, the elastic matrix has the symmetry-allowed components C15, C25, C35, C46 nonzero and C45 ≡ 0 by symmetry; Table 3 instead lists C35 = 6.8 GPa and no C45. So either Table 2's C45 is a typo for C35, or the symmetry of the tensor was mishandled. If the elastic matrix was assembled with C45 instead of C35, the eigenvalue spectrum and the generalized Born criteria computed from it would be wrong. The text further states 'C33 > C11 > C22', but the table shows C22 = 159.7 GPa, C11 = 151.5 GPa, C33 = 140.6 GPa, i.e., C22 > C11 > C33. The text's assertion is inconsistent with its own data and would, if true, imply the opposite anisotropy of what the compressibility tensor shows.\n\nThese inconsistencies undermine confidence in the negative eigenvalue at 14.3 GPa. A single spurious entry (C35 vs C45) can change the sign of a near-zero eigenvalue, especially in the pressure-corrected matrix where diagonal terms are reduced by P. If the instability disappears after correcting the tensor, the claimed connection between the observed disorder and elastic instability is unsupported, leaving only a phenomenological observation of irreversible broadening.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a combined high-pressure synchrotron XRD and DFT study of β′-Mn₃(PO₄)₂ up to 20 GPa. The experiments show that the monoclinic β′ phase survives to about 13 GPa, but at 14.1 GPa the diffraction peaks broaden and weaken abruptly, and the degraded pattern persists after decompression to 0.7 GPa. The authors fit a third-order Birch–Murnaghan equation of state (K₀ = 81(2) GPa) and derive the compressibility tensor. DFT reproduces the lattice-parameter evolution, shows that compression is accommodated by Mn–O polyhedral distortions while PO₄ tetrahedra remain rigid, and predicts that three penta-coordinated Mn sites become octahedral before the transition. Elastic-constant calculations are used to claim that the crystal becomes mechanically unstable near 14.3 GPa, correlating with the observed loss of crystallinity. The central experimental finding is that compression produces an irreversible, disordered state rather than a crystalline high-pressure phase.","tokens_in":16174,"tokens_out":4726,"duration_ms":45100,"significance":"If the conclusions hold, the paper provides a useful data point for the high-pressure behavior of complex A₃(PO₄)₂ frameworks: structurally complex, distorted frameworks may undergo irreversible loss of long-range order instead of a cooperative crystalline transition. The work combines synchrotron XRD with DFT, and several computational choices are well grounded: U_eff = 4.3 eV is a literature value with a sensitivity test, the magnetic ordering is selected by energy minimization among several configurations, and the experimental/computational comparison is not fitted to the transition pressure. These strengths make the experimental observation credible. However, the computational support for elastic instability rests on elastic tensors that contain internal inconsistencies, so the central mechanistic link requires verification before the conclusions can be accepted in their present form.","major_comments":[{"comment":"The elastic-constant data are internally inconsistent. For monoclinic P2₁/c, the symmetry-allowed constants include C15, C25, C35 and C46, while C45 is symmetry-forbidden. Table 2 lists C45 = 0.4 GPa and omits C35; Table 3 lists C35 = 6.8 GPa and omits C45. Moreover, the text states 'C33 > C11 > C22', but Table 2 gives C22 = 159.7, C11 = 151.5, C33 = 140.6 GPa. These errors directly affect the eigenvalue analysis and the generalized Born criterion, including the expression displayed in the text, which uses C35. The claim that one eigenvalue becomes negative at 14.3 GPa is therefore not currently supported. Please provide the full 6×6 tensors at both pressures, correct the symmetry-allowed entries, and re-evaluate the Born criteria and eigenvalues.","section":"III.e, Tables 2 and 3"},{"comment":"All quantitative experimental results — lattice parameters, EoS parameters, and the compressibility tensor — come from Rietveld refinements in which all 117 internal atomic coordinates were fixed at their ambient-pressure literature values. The paper acknowledges that the DAC data cannot refine these coordinates. Since the DFT calculations show substantial internal rearrangements (e.g., coordination changes at several Mn sites), the fixed-coordinate model may bias the derived unit-cell parameters and hence K₀, K₀′, and the compressibility axes. The qualitative irreversible-broadening result is unaffected, but the claimed quantitative agreement between experiment and DFT (Section III.b, Figure 5) should be either reassessed with a sensitivity test or stated with this caveat.","section":"III.a and III.b"},{"comment":"The manuscript’s central claim is phrased as a 'loss of long-range crystallographic order.' The diffraction data demonstrate irreversible peak broadening and weakening, but the authors themselves note that the data do not uniquely distinguish partial amorphization, nanocrystallization, severe microstrain, defect accumulation, or unresolved phase coexistence. The abstract and conclusions nevertheless state the stronger interpretation. Please temper the wording to match the evidence, or provide additional characterization (e.g., recovered-sample TEM, diffuse scattering analysis, or a quantitative strain/size analysis) that discriminates among these possibilities.","section":"III.a, Abstract, and Conclusions"}],"minor_comments":[{"comment":"Typo: 'Precipitation was induced by adding by adding ammonium hydroxide' — remove the duplicated 'by adding'.","section":"II.a"},{"comment":"Reference [25] is cited for the bulk modulus of Co₃(VO₄)₂ (122 GPa), but [25] is the PBEsol functional paper. Please check and correct the citation.","section":"III.b"},{"comment":"G. Garbarino is listed in the author-contribution statement but does not appear in the author list. Add the author or correct the contribution statement.","section":"Author contributions"},{"comment":"Please state explicitly whether the reported constants are the raw second-order elastic constants or the pressure-corrected constants used in the Born analysis. The text describes a pressure correction, but it is not clear how it was applied to the values in Tables 2 and 3.","section":"III.e, Table 2 caption"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the experimental observation of irreversible pressure-induced degradation near 14 GPa appears solid and interesting. My main concern is the elastic-tensor inconsistency, which is central to the computational claim of a mechanical instability near 14.3 GPa. This is fixable by re-running the stability analysis, but it must be addressed. I would not reject solely on this basis, because the experimental result and the DFT structural-evolution analysis can stand without the elastic-instability claim if necessary."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The one thing to know: this is the first high-pressure XRD study of β'-Mn3(PO4)2, and the main experimental result — irreversible loss of long-range order above 14.1 GPa — is probably right. The second: the companion computational claim, elastic instability at 14.3 GPa, rests on elastic tensors that are internally inconsistent, so don't repeat that claim until it's been re-derived.\n\nCredit first. The experiment is clean. Peak widths are flat up to 13 GPa, jump at 14.1, and stay broad after full decompression; the irreversibility is the persuasive part. The authors do the right check on non-hydrostaticity — methanol–ethanol turns non-hydrostatic around 10 GPa, yet nothing moves until 14.1 — and they are appropriately cagey about the transformed state, listing partial amorphization, nanocrystallization, severe microstrain, and phase coexistence as all possible. The DFT is solid: full relaxation of the 156-atom cell, Ueff sensitivity test, spin ordering chosen by energy minimization, and good lattice-parameter agreement with experiment. The coordination-evolution result (three MnO5 sites → octahedral by ~5–8 GPa) is useful microscopic detail.\n\nThe biggest soft spot is the elastic tensor reporting; the stress-test note is right. Table 2 lists C45 = 0.4 GPa, symmetry-forbidden for P21/c; Table 3 swaps in C35 = 6.8 GPa. The text then states C33 > C11 > C22 while its own table shows C22 > C11 > C33 — and the compressibility tensor agrees with the table, since the b-axis is least compressible. Any one could be a transcription slip, but the negative eigenvalue at 14.3 GPa is the entire basis for the mechanical-instability claim, and a single wrong entry can flip a near-zero eigenvalue in the pressure-corrected matrix. The authors need to recompute the spectrum from a corrected tensor and say which components actually drive the instability.\n\nSecond, the experimental EoS and compressibility axes come from refinements with all 117 internal coordinates fixed; the paper admits this in Sec. III.a. So K0 = 81(2) GPa and the principal compressibility values carry unquantified systematic bias, and the 'statistically indistinguishable' DFT-vs-experiment EoS comparison overstates the agreement. This doesn't touch the qualitative disorder observation, but it limits how much weight the quantitative section can bear.\n\nMinor: raw data are not deposited, and the disordered product is not identified beyond 'disordered.'\n\nWho gets value: anyone tracking A3(PO4)2 systematics — the parallel with β-Zn3(PO4)2 and the contrast with Co3(PO4)2 is useful. It deserves a serious referee. Send to review, with a firm request for corrected elastic tensors and a re-derived instability claim, plus disclosure of the fixed-coordinate refinement.","headline":"First HP-XRD study of β'-Mn3(PO4)2: the irreversible disorder near 14 GPa looks real, but the elastic-instability claim needs a corrected tensor before it can be used.","tokens_in":16709,"tokens_out":9359,"would_cite":true,"duration_ms":74582,"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":"β′-Mn3(PO4)2 loses crystalline order irreversibly above 14.1 GPa, and density-functional calculations place a mechanical instability at nearly the same pressure, indicating the collapse is intrinsic to the crystal framework.","keywords":["high-pressure X-ray diffraction","manganese phosphate","irreversible structural disorder","elastic instability","compressibility tensor","density-functional theory","polyhedral compression","phosphate frameworks"],"falsifier":"Compress a fresh sample to 20 GPa in a quasi-hydrostatic medium such as helium or neon, then decompress to ambient pressure: if the recovered X-ray diffraction pattern is sharp and matches the starting β′ phase, the claimed irreversibility and intrinsic disorder are wrong. A second check would be measuring single-crystal elastic constants near 14 GPa: if all eigenvalues remain positive, the predicted mechanical instability is absent.","tokens_in":15642,"feed_emoji":"💎","tokens_out":6326,"duration_ms":60420,"temperature":0.7,"pith_summary":"This paper tries to establish that β′-Mn3(PO4)2, the stable ambient form of manganese phosphate, does not survive compression to 20 GPa as a crystalline material. Instead, near 14.1 GPa it develops irreversible structural disorder: X-ray diffraction peaks broaden and weaken, no known Mn3(PO4)2 polymorph can index the pattern, and the broadened pattern persists after pressure is released. The paper also shows through density-functional calculations that the crystal becomes mechanically unstable near the same pressure, with one elastic eigenvalue turning negative at 14.3 GPa. If true, this means the high-pressure limit of this phosphate is a permanent, partially disordered state rather than a new crystal phase, and it connects that behavior to a specific elastic softening of the polyhedral framework.","feed_headline":"Manganese phosphate goes permanently disordered at 14 GPa","feed_subtitle":"New diffraction and theory results tie the collapse to an elastic instability near 14 GPa.","key_machinery":"The load-bearing computational object is the pressure-dependent elastic stiffness tensor of the monoclinic P21/c cell, evaluated at each optimized volume. Mechanical stability is judged through the generalized Born stability criteria; at 14.3 GPa one eigenvalue becomes negative, meaning certain strain modes lower the elastic energy and the crystal is unstable to small perturbations. The complementary structural mechanism is the polyhedral-unit analysis: nearly rigid PO4 tetrahedra versus soft, distorting MnO5/MnO6 polyhedra, with three Mn sites switching from five- to six-fold coordination and one Mn–O bond anomalously stretching. This combination is what the paper uses to tie the observed i","core_discovery":"The central discovery is that β′-Mn3(PO4)2 remains crystalline up to about 14 GPa but then undergoes a loss of long-range crystallographic order that is irreversible upon decompression. The authors report a third-order Birch–Murnaghan equation of state with bulk modulus K0 = 81(2) GPa and strongly anisotropic compression, with the most compressible principal axis lying in the ac plane. Density-functional calculations reproduce the measured compressional behavior and attribute it to distortions of MnO5/MnO6 polyhedra while PO4 tetrahedra remain nearly rigid. Three initially penta-coordinated Mn sites become octahedrally coordinated below 8 GPa, and one MnO6 polyhedron shows anomalous bond len","pith_inferences":["The quantitative lattice-parameter, equation-of-state, and compressibility results could be biased if the internal atomic coordinates shift substantially under pressure, as the DFT calculations indicate they do; a high-pressure single-crystal or better-resolved powder study would test this directly.","The DFT prediction that three Mn sites become octahedrally coordinated between about 4.7 and 7.6 GPa could be probed with Mn K-edge X-ray absorption spectroscopy under pressure, providing an experimental check of the proposed densification mechanism before the 14 GPa collapse.","The anomalous stretching of a single Mn–O bond resembles a soft-mode precursor, so Raman or inelastic X-ray scattering up to 14 GPa might reveal a phonon that softens toward zero frequency, giving a direct dynamical signature of the instability.","A separate compression-decompression cycle that stops just below and just above 14.1 GPa could map the irreversibility boundary and test whether crossing the elastic instability is required to produce the permanent disorder."],"forward_implications":["The crystalline form of β′-Mn3(PO4)2 cannot be used as a stable phase above roughly 14 GPa; any application requiring crystalline integrity would instead produce a permanently disordered material.","The reported bulk modulus and full anisotropic compressibility tensor provide quantitative input for modeling the behavior of manganese phosphates in high-pressure and geophysical settings.","The correspondence between the predicted elastic instability and the experimental disorder onset suggests that other structurally complex A3(PO4)2 frameworks may similarly degrade irreversibly rather than transform to new crystalline phases.","The predicted pressure-driven coordination changes at individual Mn sites offer a microscopic fingerprint that could be observed with element-specific spectroscopies under compression."],"fun_headline_variants":["Mn3(PO4)2 irreversibly disordered above 14 GPa","High pressure makes phosphate permanently disordered","Irreversible pressure collapse in Mn3(PO4)2 at 14 GPa","Permanent disorder in manganese phosphate triggered at 14 GPa","Elastic instability leads to irreversible disorder at 14 GPa"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The quantitative lattice and compressibility results assume that all 117 internal atomic coordinates stay at their ambient-pressure values during compression (Section III.a), because the powder data could not refine them; if those coordinates shift substantially under pressure, as the DFT calculations suggest, the refined lattice parameters and derived elastic quantities could be biased.","fun_headline_variants_meta":{"raw":{"variants":["Mn3(PO4)2 irreversibly disordered above 14 GPa","High pressure makes phosphate permanently disordered","Irreversible pressure collapse in Mn3(PO4)2 at 14 GPa","Permanent disorder in manganese phosphate triggered at 14 GPa","Elastic instability leads to irreversible disorder at 14 GPa"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000227,"raw_usage":{"total_tokens":1372,"prompt_tokens":874,"completion_tokens":498,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":618,"completion_tokens_details":{"reasoning_tokens":408}},"tokens_in":618,"tokens_out":498,"duration_ms":5000,"temperature":1.0,"reasoning_tokens":408,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T01:09:59.227030+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compress a fresh sample to 20 GPa in a quasi-hydrostatic medium such as helium or neon, then decompress to ambient pressure: if the recovered X-ray diffraction pattern is sharp and matches the starting β′ phase, the claimed irreversibility and intrinsic disorder are wrong. A second check would be measuring single-crystal elastic constants near 14 GPa: if all eigenvalues remain positive, the predicted mechanical instability is absent.","supporting_citations":[],"review_version":1}