{"id":"5789920b-c6ce-4f2f-87de-c5b8fc4585a9","arxiv_id":"2608.08088","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A physics-guided ML pipeline screens cerium compounds, simulates the top candidates as Ising magnets, and nominates CeGaO3 as a candidate ferromagnet.","lead":"This paper builds a computer pipeline that screens known cerium compounds for likely ferromagnets, then simulates the most promising ones with Monte Carlo and neural networks. If the method holds up, it could speed up the search for permanent magnets made from abundant cerium instead of scarce rare earths.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The under-specified modified-GK f→5d coupling (Eq. 1) is the load-bearing assumption: J_Ce sets both the positive-exchange screen and all reported T_C values, yet its parameters and EuO derivation are not disclosed; a sensitivity or independent DFT+U check is needed.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing point: the material-level mean modified-GK coupling J_Ce is both the selection criterion and the uniform Hamiltonian coupling that fixes T_C in kelvin, and it is under-specified. My independent reading of the manuscript confirms this. The paper explicitly states that the f→5d correction is only proportional to a combination of parameters, with no numerical values given, and the concluding remarks acknowledge that the single effective J is a limitation and that DFT-based confirmation is future work. This is more than a missing implementation detail because the signed value of J is what keeps CeGaO3 and CeF3 in the candidate pool, and the magnitude sets every reported transition temperature. The EuO benchmark validates the Monte Carlo and autoencoder machinery at a fixed J, but it does not validate the exchange model itself because the J input is generated by the same undisclosed procedure. I considered alternative concerns, such as the near-chance balanced accuracy of the Random Forest on minority classes and the three-size finite-size scaling, but these are secondary: the RF stage merely produces a broad pool, and the finite-size analysis is honest about its uncertainty. The paper does have independent support: the EuO benchmark is a useful end-to-end check of the simulation pipeline, and reference 26 independently predicts ferromagnetism in CeGaO3, which mitigates concern about the sign of J for that compound. However, the Ising-candidate label additionally requires uniaxial anisotropy, which is explicitly deferred. Therefore the correct verdict remains CONDITIONAL: the workflow is plausible and the limitations are transparently stated, but the central material-specific claim should not be accepted as a solid prediction until the modified-GK parameters are disclosed and the exchange couplings are confirmed by an independent electronic-structure or sensitivity analysis.","tokens_in":14563,"tokens_out":5453,"duration_ms":56655,"concrete_test":"Recompute the modified-GK couplings for EuO, CeF3, and CeGaO3 from a fully specified version of Eq. (1): authors must report b_fd, Δ_fd, the prefactor, and the angular factor for each Ce–anion–Ce pathway. Then run a sensitivity scan varying the f→5d amplitude from 0 to 2 times the nominal value; if J_Ce for CeGaO3 changes sign or if its ratio to J_EuO shifts by more than about 20%, the screening and all T_C values are not robust. Independently, perform DFT+U or hybrid-functional total-energy calculations for FM and AFM spin configurations of CeGaO3 and CeF3, extract nearest-neighbor J via energy mapping, and compute the magnetocrystalline anisotropy energy; this directly tests whether the positive-exchange screen and the 'Ising ferromagnet' label survive outside the modified-GK model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central assertion that CeGaO3 is a candidate Ising ferromagnet hinges on the material-level modified-GK coupling J_Ce = <J_p> entering Eq. (2) as the uniform J. This single number performs three jobs: it is the positive-exchange filter (J_Ce > 0), it sets the absolute temperature scale via T_C in kelvin, and it defines the magnetic graph connectivity. The section 'Goodenough–Kanamori Screening with f-Orbital Correction' specifies the f→5d channel only as a proportionality, J_{f→5d} ∝ 2 b_fd^2 / [(4 S_i S_j) Δ_fd], 'scaled by the pathway angular factor', with no numerical values for b_fd, Δ_fd, S_i, S_j, no prefactor, and no statement of how the three reported J values (0.5351, 3.6631, 4.7280 meV) were obtained. The EuO benchmark cannot validate this input because J_GK = 0.5351 meV is itself an output of the same under-specified screen; the 12.2% T_C deficit could simply reflect a tuned J rather than an independent test of the exchange model. Since T_C is linear in J, an undisclosed factor of two in the f→5d amplitude would move CeGaO3 from about 327 K to roughly 160 or 650 K, and a sign error would remove the compound from the candidate list entirely. The paper's own concluding remarks concede that 'a single effective exchange parameter J_eff obtained from the average exchange interaction predicted by the modified Goodenough–Kanamori screening model' is a limitation and that DFT+U validation is future work. Reference 26 provides independent first-principles support for ferromagnetic CeGaO3, which reduces the risk of a sign error, but it does not supply the exchange magnitude or the uniaxial-anisotropy information needed for the Ising-candidate label.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a multi-stage computational pipeline for discovering candidate Ising ferromagnets among Ce-based compounds. A Random Forest classifier trained on Materials Project magnetic labels predicts ferromagnetic candidates from structural and SOAP descriptors; a modified Goodenough–Kanamori (GK) screen with an added f-to-5d virtual-excitation channel retains compounds with positive mean exchange J_Ce; Glauber Monte Carlo estimates critical temperatures; conserved-magnetization Kawasaki Monte Carlo on crystal-derived magnetic graphs generates phase boundaries, which are reconstructed by a local-affinity feature and a one-dimensional autoencoder; critical exponents are extracted from the binodal. The pipeline is demonstrated on EuO as an external benchmark and applied to two Ce candidates, CeF3 and CeGaO3, with the central claim that CeGaO3 is a candidate Ising ferromagnet.","tokens_in":14927,"tokens_out":4671,"duration_ms":45354,"significance":"If the full pipeline were validated end-to-end, it would be a useful transferable framework for prioritizing lanthanide-based magnetic candidates: the structural RF stage is fast, the MC/autoencoder machinery is generic, and the EuO benchmark is a genuine external test of the phase-mapping workflow. Strengths of the paper include the explicit EuO benchmark reproducing the 3D Ising exponent within uncertainty (β_AE = 0.337 ± 0.084, β_LA = 0.330 ± 0.063), the careful reporting of phase-boundary errors for both AE and LA methods, the acknowledgment that the effective-J approximation is a limitation, and the commitment to release code and data. However, the Ce-specific predictions rest on an under-specified f-to-5d modification of the GK model whose numerical parameters are not given, and the same J_Ce is used both to select candidates and to define the simulated Hamiltonian. This creates a circularity that the EuO benchmark does not resolve, because the EuO coupling is itself an output of the same screening model. The RF stage also provides near-chance balanced accuracy for the minority classes, so the practical filtering power of that stage is limited.","major_comments":[{"comment":"The f-to-5d correction is specified only as a proportionality, J_{f→5d} ∝ 2 b_fd^2 / [(4 S_i S_j) Δ_fd], with no numerical values for b_fd, Δ_fd, S_i, S_j, no angular prefactor, and no derivation of the three reported couplings (0.5351, 3.6631, 4.7280 meV). Because J_Ce serves simultaneously as the positive-exchange filter (J_Ce > 0) and as the uniform coupling J in Eq. (2), all predicted T_C values in Table 2 scale linearly with this undisclosed input. The EuO benchmark cannot independently validate the model because J_GK = 0.5351 meV is an output of the same screen. The authors should disclose the parameters, provide a sensitivity analysis (for example, varying the f-to-5d amplitude by a factor of two), or replace the screening with a first-principles evaluation for the two Ce candidates. The concluding remark that DFT+U validation is future work is not sufficient to support the central CeGaO3 claim.","section":"Goodenough–Kanamori Screening with f-Orbital Correction, Eq. (1)"},{"comment":"The identification of CeGaO3 as a candidate Ising ferromagnet is circular in structure: the same modified-GK coupling J_Ce > 0 selects the compound, and the same J is then inserted as the uniform coupling in the Ising Hamiltonian whose ferromagnetic phase is simulated. The EuO benchmark validates the MC/autoencoder machinery but not the Ce-specific exchange input. A non-circular test would be to compute J from an independent method (DFT+U or the published first-principles study in Ref. 26) and compare the predicted T_C and magnetic order. Without such a test, the CeGaO3 prediction is inherited from the assumed J rather than independently derived from the material.","section":"Results and Discussion, Ce-Based Candidate Phase Diagrams"},{"comment":"The balanced accuracies for the ternary and combined RF models are 0.434 and 0.453, only slightly above the 1/3 chance level for three classes, and the macro-F1 scores are 0.450 and 0.467. Since the out-of-fold predictions from this model define the 838-compound FM candidate pool, the RF stage contributes little filtering power beyond the majority FM label. The subsequent no-TM-overlap, z≥6, and CIF-availability filters are post-hoc, and the paper does not report how many of the 237 positive-exchange candidates survive each filter or why CeF3 and CeGaO3 are representative of the surviving set. The authors should quantify the selectivity of each funnel stage and report the RF confusion matrix, or explicitly frame the RF stage as a baseline descriptor benchmark rather than a load-bearing screening component.","section":"Ferromagnetic Prediction Using Structural Characteristics"},{"comment":"The CeF3 analysis yields β_LA = -0.102 ± 0.111, which the paper itself describes as 'not physically meaningful for an order-parameter exponent.' While this is an honest report, it means the pipeline's phase-diagram characterization fails for one of the two Ce candidates, and the negative exponent is left unexplained. The paper should analyze why the conserved-magnetization binodal fits fail for CeF3 (for example, lattice connectivity or finite-size effects) or restrict the phase-diagram claims to CeGaO3. As written, the abstract's statement that the framework 'characterize[s] their magnetic phase transitions and critical properties' is not supported for CeF3.","section":"Results and Discussion, Ce-Based Candidate Phase Diagrams"}],"minor_comments":[{"comment":"There are typos in the figure captions and axis labels: 'T emperature' should be 'Temperature', and the axis label 'n 1/' should read 'n^{-1/ν}'. In Figure 10 the exponent labels are missing the β symbol.","section":"Figures 4, 8, and 10"},{"comment":"The seven-feature descriptor includes 'space group,' but the encoding scheme is not specified. Please clarify whether it is one-hot, ordinal, or a numerical space-group number, since this affects the interpretation of feature importances.","section":"Model and Methods, Ferromagnetic Prediction Using Structural Characteristics"},{"comment":"The summation index 'c' in Eq. (1) is used before the channel index is defined. Please list the standard GK terms (σ, π, 90°) explicitly and state the sign rules used, or provide a reference to the specific GK formulation adopted.","section":"Model and Methods, Goodenough–Kanamori Screening with f-Orbital Correction, Eq. (1)"},{"comment":"The finite-size scaling uses ν = 0.876 (3D percolation) for off-critical magnetization sectors, taken from the reference protocol. The paper should justify transferring this exponent to material-derived graphs with non-cubic connectivity, or discuss the sensitivity of the extracted T_bin to this choice.","section":"Model and Methods, Monte Carlo Simulation with Conserved Magnetization"},{"comment":"The text and figures refer to 'ten independent autoencoder cases,' but the main text never defines what distinguishes the cases (initialization seeds, data subsampling, or architecture variants). Please specify this in the Methods section.","section":"Results and Discussion, EuO Benchmark Validation"},{"comment":"The data and code are stated to 'will be made publicly available on GitHub' without a URL or version identifier. A permanent repository link or DOI is needed for the availability claim to be verifiable.","section":"Supporting Information, Data, and Code Availability"}],"recommendation":"major_revision","confidential_remarks":"The paper's novelty lies mostly in the combination of existing tools (RF from Broyles et al., MC/autoencoder from Jang and Yethiraj) with a Ce-specific GK correction. The latter is the weakest link: it is under-parameterized, and the EuO benchmark does not validate it because the EuO coupling is an output of the same model. Reference 26 provides independent first-principles support for a ferromagnetic CeGaO3 ground state; the authors could use that study to calibrate or check the modified-GK J instead of relying solely on the proportionality form. Given the strong EuO benchmark and the transparency about limitations, the issues are fixable within the manuscript's scope, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"My read in short: the genuinely new thing here is the assembly. No one has run a Ce-specific screening funnel that goes from Materials Project structures through an RF structural filter, a modified Goodenough-Kanamori exchange screen, and then a full Glauber/Kawasaki Monte Carlo plus autoencoder phase reconstruction. That integration is useful, and the authors are honest that the components come from prior work. I largely agree with the conditional verdict. The EuO benchmark is the strongest part: the same workflow gives Tc = 60.6 K with J = 0.5351 meV, about 12% below the experimental 69 K, and the extracted beta values (0.337 and 0.330) sit on the 3D Ising value. That validates the MC and autoencoder machinery, not the exchange model.\n\nThe soft spot is exactly where the stress-test note puts it. Eq. (1) states the f-to-5d correction only as a proportionality: no b_fd, no Delta_fd, no prefactor, no derivation of the three J values. This matters because that same J does three jobs: it is the positive-exchange filter, the uniform coupling in the Ising Hamiltonian, and the conversion to kelvin. A missing factor of two in the f-to-5d amplitude would move CeGaO3 from about 327 K to roughly 160 or 650 K. The EuO benchmark does not independently validate that input, since J_GK itself comes from the same under-specified model. Ref. 26 does provide first-principles support for ferromagnetic CeGaO3, which lowers the sign-error risk, but it does not supply the exchange magnitude or the uniaxial-anisotropy data needed for the Ising-candidate label.\n\nThe other weaknesses are real but less central. The RF balanced accuracies (0.590, 0.434, 0.453) mean minority AFM/PM classes are close to chance; the 838-compound FM pool is noisy, though the GK screen and the post-hoc z-at-least-6/no-TM-overlap filters compensate somewhat. The three-size ladders for the Ce compounds are thin for finite-size scaling. CeF3's exponents are unphysical, which the paper acknowledges, and the CeGaO3 autoencoder beta = 0.300 +/- 0.209 has large uncertainty; the local-affinity estimate is the more credible one. None of this is fatal, but it means the Ce-specific predictions are a screening hypothesis, not a demonstrated result.\n\nWho is this for? Computational magnetism and materials-informatics readers who want a template for combining structural ML with MC phase mapping. It is not yet a validation of CeGaO3 as a permanent-magnet candidate. I would send it to peer review with a request for major revision: disclose the GK parameters, provide sensitivity analysis, add larger size ladders, and get one DFT+U check of the exchange couplings and anisotropy.","headline":"A useful integration of existing methods with a solid EuO benchmark, but the Ce-specific claims rest on an under-specified Goodenough-Kanamori coupling and need parameter disclosure plus independent DFT+U validation.","tokens_in":15555,"tokens_out":2844,"would_cite":false,"duration_ms":25999,"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":"A four-stage computational pipeline—Random Forest structural screening, cerium-specific Goodenough–Kanamori exchange filtering, Ising Monte Carlo, and autoencoder phase reconstruction—identifies CeGaO3 as a candidate Ising ferromagnet…","keywords":["cerium-based magnetism","Ising ferromagnet candidate discovery","Random Forest structural screening","Goodenough–Kanamori exchange","Monte Carlo phase diagram","autoencoder phase transitions","uniaxial magnetic anisotropy","rare-earth magnetic materials"],"falsifier":"Compute the Ce–anion–Ce exchange couplings of CeGaO3 and CeF3 with a first-principles method such as DFT+U or hybrid functionals; if the calculated $J_{\\mathrm{Ce}}$ is negative, or far below 4.73 and 3.66 meV, then the predicted ferromagnetism and the reported transition temperatures are artifacts of the assumed model. On the experimental side, magnetization and specific-heat measurements of polycrystalline or oriented CeGaO3 should show a spontaneous moment and an order-disorder feature near 327 K for the central claim to hold.","tokens_in":14299,"feed_emoji":"🧲","tokens_out":8698,"duration_ms":78622,"temperature":0.7,"pith_summary":"The paper tries to establish that the search for cerium-based ferromagnets with uniaxial, Ising-like anisotropy can be accelerated by a four-stage computational pipeline: a structural Random Forest screen, a cerium-specific Goodenough–Kanamori exchange filter, Ising Monte Carlo simulations on crystal-derived magnetic lattices, and an autoencoder that reconstructs the phase diagram from spin snapshots. Applying the pipeline to all 1,011 Ce-containing binary and ternary compounds in a public crystal-structure database reduces the candidate pool to 237 positive-exchange ferromagnets, and two of these, CeF3 and CeGaO3, are carried through full finite-size phase-diagram analysis. The central positive result is that CeGaO3 recovers the three-dimensional Ising order-parameter exponent in the vicinity of the 3D Ising value within uncertainty, identifying it as a candidate Ising ferromagnet worth experimental synthesis and easy-axis validation. A reader should care because cerium is abundant and cheap, and the framework claims to give a reusable route from crystal-structure databases to prioritized magnetic candidates.","feed_headline":"Pipeline flags CeGaO3 as a cerium Ising ferromagnet","feed_subtitle":"This four-stage ML workflow narrows 1,011 cerium compounds and gives CeGaO3 a 3D-Ising signature.","key_machinery":"The load-bearing object is the material-level modified-goodenough–Kanamori coupling $J_{\\mathrm{Ce}} = \\langle J_p \\rangle$, the mean over nonzero Ce–anion–Ce pathway sums of Eq. (1), where each pathway sum includes the standard $\\sigma$, $\\pi$, and $90^\\circ$ terms plus a cerium-specific f→5d virtual channel $J_{f\\to 5d} \\propto 2 b_{fd}^2 / [(4 S_i S_j) \\Delta_{fd}]$. This single number does triple duty: it keeps a candidate only if $J_{\\mathrm{Ce}} > 0$, it sets the uniform coupling in the Ising Hamiltonian $H = -J \\sum_{\\langle ij\\rangle} S_i S_j$, and through $k_B$ it fixes the reported transition temperatures in kelvin. Around this sits the rest of the pipeline: a seven-feature Random Forest (unit-cell volume, density, site count, space group, atomic density, Ce and transition-metal SOAP overlaps) that generates the ferromagnetic candidate pool with out-of-fold predictions, Glauber dynamics for Binder-cumulant $T_C$ estimates, constant-magnetization Kawasaki dynamics on CIF-derived magnetic graphs, graph-distance affinity features $v_i$, and a one-dimensional-latent autoencoder whose latent standard deviation locates the binodal.","core_discovery":"On its own terms, the paper's central discovery is a validated workflow rather than a single material. The modified-GK screen, which adds a positive f→5d virtual-excitation channel to the standard Goodenough–Kanamori superexchange terms, assigns material-level couplings $J_{\\mathrm{Ce}} = +4.7280$ meV to CeGaO3 and $+3.6631$ meV to CeF3; Glauber Binder-cumulant crossings then give $T_C = 326.8\\pm 19.7$ K and $340.7\\pm 15.2$ K respectively. The autoencoder and local-affinity binodal reconstructions for CeGaO3 yield order-parameter exponents ($0.300\\pm 0.209$ and $0.363\\pm 0.057$) consistent with the 3D Ising value, while CeF3's recovered exponents are far from Ising, which the authors read as indicating magnetic excitations rather than a sharp order-disorder transition. EuO, run through the same graph and simulation machinery, gives $T_C = 60.6\\pm 0.9$ K versus the experimental roughly 69 K, providing an external check that the workflow is not grossly off-scale. The result is stated carefully: these are computationally predicted ferromagnets whose uniaxial anisotropy still needs experimental confirmation.","pith_inferences":["If the modified-GK coupling for CeGaO3 is later confirmed by first-principles exchange calculations, the predicted ordering temperature near 327 K would make it an unusually high-$T_C$ cerium-only ferromagnet; the paper itself does not compare this against known Ce intermetallics.","The same pipeline applied to neodymium- and samarium-based compounds could yield a family of rare-earth gap-magnet candidates, since the f-to-5d correction generalizes to other lanthanides.","The non-Ising exponents recovered for CeF3 are a built-in falsifier: experimental single-crystal studies should show no sharp uniaxial transition, suggesting the strong positive $J$ from the f-to-5d channel may be an artifact of the model rather than a material property.","A cheap testable extension would be to rerun the autoencoder stage on the EuO configurations with tanh activations to quantify the claimed bias of that activation choice, since the paper's activation-function analysis is qualitative."],"forward_implications":["If the workflow is sound, CeGaO3 moves from an understudied geometry to a concrete synthesis target for confirming Ce-based Ising ferromagnetism, with a predicted ordering temperature near 327 K and a 3D-Ising-like order-parameter exponent.","The pipeline's 237 positive-exchange candidates constitute a shortlist for higher-cost first-principles and anisotropy calculations, concentrating effort on a small fraction of the 1,011 Ce compounds screened.","EuO's benchmark error (about 12% at the Glauber level) sets a realistic expectation: predicted transition temperatures for new Ce candidates should be read as approximate, with comparable or larger uncertainty.","The framework is claimed to transfer beyond Ce: the same stages with analogous f-orbital corrections apply to other rare-earth and actinide magnetic systems.","The CeF3 result shows the screen can flag a positive-exchange ferromagnet whose simulated phase behavior is not Ising-like, so the pipeline distinguishes a general ferromagnetic candidate from an Ising-anisotropy candidate."],"supporting_citations":[{"why":"Supplies the Random Forest structural-plus-SOAP descriptor and the three-class FM/AFM/PM formulation transferred from uranium compounds to cerium compounds.","marker":"10"},{"why":"Supplies the constant-magnetization Ising protocol, local-affinity feature construction, binodal finite-size extrapolation with exponents 0.630 and 0.876, and the autoencoder architecture used throughout.","marker":"12"},{"why":"Motivates the f-to-5d virtual-excitation channel as an additional positive exchange contribution in rare-earth chalcogenides.","marker":"6"},{"why":"Provides first-principles confirmation of exchange coupling in Eu monochalcogenides, the physical basis for the cerium-specific correction's form.","marker":"7"},{"why":"Provides the experimental EuO Curie temperature near 69 K used as the external ferromagnetic benchmark for the Monte Carlo and autoencoder workflow.","marker":"13"},{"why":"Supplies the crystal-structure entries and magnetic labels that define the screened materials space.","marker":"14"},{"why":"Supplies the Smooth Overlap of Atomic Positions descriptor used to encode local coordination in the Random Forest features.","marker":"17"},{"why":"Reports a first-principles prediction of ferromagnetic order in CeGaO3, used to support the candidate status of the paper's main result.","marker":"26"},{"why":"Gives the three-dimensional Ising critical exponents, including the order-parameter exponent near 0.326, which are the reference values for exponent recovery and phase-boundary errors.","marker":"27"},{"why":"Reports no distinct magnetic ordering in CeF3 single crystals, used to interpret the non-Ising exponents recovered for CeF3.","marker":"28"}],"fun_headline_variants":["ML framework flags CeGaO3 as cerium Ising ferromagnet","AI screen finds CeGaO3 for rare-earth magnet potential","Cerium magnet candidate CeGaO3 predicted by ML","CeGaO3: AI-predicted Ising ferromagnet from cerium screen","ML sifts 1,011 cerium compounds to land on CeGaO3"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything downstream—which compounds stay in the candidate pool, the Ising Hamiltonian, and every temperature in kelvin—rests on the assumption that the average modified-goodenough–Kanamori coupling $J_{\\mathrm{Ce}}$ correctly represents the real exchange in the compound, even though the f-to-5d correction is specified only as a proportionality and the quoted values are asserted without derivation.","fun_headline_variants_meta":{"raw":{"variants":["ML framework flags CeGaO3 as cerium Ising ferromagnet","AI screen finds CeGaO3 for rare-earth magnet potential","Cerium magnet candidate CeGaO3 predicted by ML","CeGaO3: AI-predicted Ising ferromagnet from cerium screen","ML sifts 1,011 cerium compounds to land on CeGaO3"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000861,"raw_usage":{"total_tokens":3790,"prompt_tokens":1051,"completion_tokens":2739,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":667,"completion_tokens_details":{"reasoning_tokens":2640}},"tokens_in":667,"tokens_out":2739,"duration_ms":22050,"temperature":1.0,"reasoning_tokens":2640,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:26:39.404751+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Ce–anion–Ce exchange couplings of CeGaO3 and CeF3 with a first-principles method such as DFT+U or hybrid functionals; if the calculated $J_{\\mathrm{Ce}}$ is negative, or far below 4.73 and 3.66 meV, then the predicted ferromagnetism and the reported transition temperatures are artifacts of the assumed model. On the experimental side, magnetization and specific-heat measurements of polycrystalline or oriented CeGaO3 should show a spontaneous moment and an order-disorder feature near 327 K for the central claim to hold.","supporting_citations":[{"cited_title":"Structure-driven prediction of magnetic order in uranium compounds","cited_arxiv_id":null,"evidence_quote":"Supplies the Random Forest structural-plus-SOAP descriptor and the three-class FM/AFM/PM formulation transferred from uranium compounds to cerium compounds."},{"cited_title":"Unsupervised machine learning method for the phase behavior of the constant magnetization Ising model in two and three dimensions","cited_arxiv_id":null,"evidence_quote":"Supplies the constant-magnetization Ising protocol, local-affinity feature construction, binodal finite-size extrapolation with exponents 0.630 and 0.876, and the autoencoder architecture used throughout."},{"cited_title":"Anomalous transport phenomena in Eu-chalcogenide alloys","cited_arxiv_id":null,"evidence_quote":"Motivates the f-to-5d virtual-excitation channel as an additional positive exchange contribution in rare-earth chalcogenides."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides first-principles confirmation of exchange coupling in Eu monochalcogenides, the physical basis for the cerium-specific correction's form."},{"cited_title":"The magnetic, optical, and transport properties of representatives of a class of magnetic semiconductors: The europium chalcogenides","cited_arxiv_id":null,"evidence_quote":"Provides the experimental EuO Curie temperature near 69 K used as the external ferromagnetic benchmark for the Monte Carlo and autoencoder workflow."},{"cited_title":"P.; Hautier, G.; Chen, W.; Richards, W","cited_arxiv_id":null,"evidence_quote":"Supplies the crystal-structure entries and magnetic labels that define the screened materials space."},{"cited_title":"P.; Kondor, R.; Cs\\' a nyi, G","cited_arxiv_id":null,"evidence_quote":"Supplies the Smooth Overlap of Atomic Positions descriptor used to encode local coordination in the Random Forest features."},{"cited_title":"S.; Farout, M.; Salmani, E.; Erum, N.; Kim, S","cited_arxiv_id":null,"evidence_quote":"Reports a first-principles prediction of ferromagnetic order in CeGaO3, used to support the candidate status of the paper's main result."},{"cited_title":"Improved high-temperature expansion and critical equation of state of three-dimensional Ising-like systems","cited_arxiv_id":null,"evidence_quote":"Gives the three-dimensional Ising critical exponents, including the order-parameter exponent near 0.326, which are the reference values for exponent recovery and phase-boundary errors."},{"cited_title":"V.; Korableva, S","cited_arxiv_id":null,"evidence_quote":"Reports no distinct magnetic ordering in CeF3 single crystals, used to interpret the non-Ising exponents recovered for CeF3."}],"review_version":1}