{"id":"7dd4849d-d71a-48e2-95d7-f21af50da1ab","arxiv_id":"1908.05186","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"FePS3 and MnPS3 have low, anisotropic thermal conductivity (k across 0.85 to 1.1 W/mK, k in-plane 2.7 to 6.3 W/mK), and UV Raman reveals new phonon features not seen in visible light.","lead":"Researchers measured the thermal conductivity of two antiferromagnetic semiconductors, FePS3 and MnPS3, and found them to be low and strongly anisotropic, with heat flowing about three times better in-plane than through-plane. The values are the first experimental reports for this material family and matter for designing spintronic devices that generate spin currents from heat.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The TR-MOKE anisotropy claim rests on fitting three unknowns (G, k⊥, k∥) to two modulation frequencies, with no identifiability or sensitivity analysis shown; the reported in-plane values may be non-unique.","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing issue: the TR-MOKE inversion from two modulation frequencies to three unknown parameters is not demonstrated to be unique, and no sensitivity analysis is provided. This is a legitimate concern about the central anisotropy claim, not about the authors' integrity or about disagreement with consensus. The paper has real strengths: it reports the first experimental thermal-conductivity data for FePS3 and MnPS3, supports the FePS3 average value with an independent optothermal Raman measurement, and provides DFT calculations for the Raman mode assignment. However, the optothermal Raman value is a weighted average, and the weights α and β are computed using the TR-MOKE values themselves, so it does not independently verify the in-plane component. The proposed synthetic-data test would settle whether the two-frequency experiment can actually constrain k∥; if the test passes, the concern is resolved. Because the reader's verdict of CONDITIONAL already requires exactly this kind of additional evidence, my read does not change the verdict.","tokens_in":19447,"tokens_out":4379,"duration_ms":48854,"concrete_test":"Recompute the TR-MOKE analysis with synthetic data: generate noiseless and noise-added surface-temperature signals from the stated cylindrical heat-diffusion model at f = 1 and 10 MHz using the reported parameters (FePS3: G = 24 MW/m2K, k⊥ = 0.85 W/mK, k∥ = 2.7 W/mK; MnPS3: G = 23 MW/m2K, k⊥ = 1.1 W/mK, k∥ = 6.3 W/mK), then fit G, k⊥, and k∥ with the same least-squares procedure. Compute the Fisher information matrix and the joint 95% confidence region for k∥. If the confidence interval for k∥ spans values below 1 W/mK or above 6 W/mK, or if the covariance matrix is ill-conditioned, the two-frequency dataset cannot support the reported in-plane conductivity. Additionally, repeat the measurement at 0.5 and 5 MHz; if the best-fit k∥ shifts by more than the stated error bars, the reported anisotropy is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the anisotropic thermal conductivity extracted from TR-MOKE. The paper states that 'By performing measurements at multiple pump modulation frequencies, we can independently determine all three unknown thermal properties of the sample stack: G, k⊥, and k∥,' but then reports data at exactly two frequencies, 1 MHz and 10 MHz. With three free parameters in the cylindrical heat-diffusion model, two modulation frequencies do not by themselves guarantee identifiability; the time-domain traces contain many points, but the information available is set by how sensitively the surface-temperature response depends on each parameter. No sensitivity matrix, parameter-correlation analysis, or confidence contours are presented. If k∥ enters mainly through radial heat spreading at the used spot sizes, its influence can be largely degenerate with the interface conductance G and the through-plane conductivity k⊥. The MnPS3 in-plane value of 6.3 ± 1.7 W/mK is especially exposed because there is no optothermal Raman cross-check for MnPS3, and the power-dependent Raman signal for MnPS3 was too weak to yield a thermal value. The FePS3 optothermal result of 1.35 W/mK is consistent with the TR-MOKE values only after choosing the weight coefficients α = 0.73 and β = 0.27, so it does not independently validate the in-plane component. The anisotropy claim, which is the headline result of the paper, is therefore not yet supported by a demonstrated unique fit.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports phonon and thermal property measurements on exfoliated films of the antiferromagnetic semiconductors FePS3 and MnPS3. Raman spectroscopy was performed with 325 nm (UV), 488 nm (blue), and 633 nm (red) excitation, revealing new UV-only spectral features near 430 and 467 cm-1. Density functional theory calculations for FePS3 were used to propose assignments for these modes. Thermal conductivity was measured by two techniques: steady-state Raman optothermal analysis, yielding an orientation-average value of 1.35 ± 0.32 W/mK for a 400-nm FePS3 film, and time-resolved magneto-optical Kerr effect (TR-MOKE), which is claimed to separately determine through-plane and in-plane conductivities. Reported values are k⊥ = 0.85 ± 0.15 W/mK and k∥ = 2.7 ± 0.3 W/mK for FePS3, and k⊥ = 1.1 ± 0.2 W/mK and k∥ = 6.3 ± 1.7 W/mK for MnPS3. The conclusion emphasizes strong thermal anisotropy and phonon-dominated heat conduction.","tokens_in":19749,"tokens_out":5757,"duration_ms":57805,"significance":"If the reported anisotropic thermal conductivities are correct, they constitute the first experimental thermal transport data for FePS3 and MnPS3, two materials of active interest for antiferromagnetic spintronics. The UV-Raman features, if properly assigned, would also add to the understanding of the vibrational spectra of metal phosphorus trichalcogenides. The paper is honest in reporting the parameter sensitivity of the DFT results, which is a strength, and the thermal conductivity values are direct measurements rather than derived from fitted parameters. However, the central anisotropy claim rests on a TR-MOKE analysis whose identifiability is not demonstrated, and the DFT support for the UV peaks is weakened by the absence of a unique calculation matching both observed modes. The significance of the work is therefore conditional on strengthening the uniqueness and sensitivity analysis of the thermal extraction.","major_comments":[{"comment":"The claim that measurements at two modulation frequencies, 1 MHz and 10 MHz, can independently determine three unknowns (G, k⊥, and k∥) is not supported by any identifiability or sensitivity analysis. With three free parameters in the cylindrical heat-diffusion model, two modulation frequencies do not, by themselves, guarantee a unique fit; the time-varying surface temperature signal at each frequency may be sensitive to only a combination of parameters. In particular, the in-plane conductivity k∥ enters through radial heat spreading, which at the used spot sizes could be largely degenerate with the interface conductance G and the through-plane conductivity k⊥. The paper reports no sensitivity matrix, no parameter correlation contours, and no joint confidence regions. The quoted uncertainties appear to reflect propagated input errors rather than fit degeneracy. Without showing that the data can distinguish k∥ from the other parameters, the reported in-plane values (2.7 W/mK for FePS3 and 6.3 W/mK for MnPS3) and the headline anisotropy claim are not demonstrated to be unique. The authors should provide a sensitivity analysis (e.g., ∂(signal)/∂ln(k∥) at both frequencies) or fit all three parameters simultaneously with reported confidence regions.","section":"TR-MOKE measurements, paragraph beginning 'By performing measurements at multiple pump modulation frequencies'"},{"comment":"The consistency check between the optothermal average (1.35 W/mK) and the TR-MOKE values (0.85 and 2.7 W/mK) is achieved by choosing weight coefficients α = 0.73 and β = 0.27 after the fact. The optothermal Raman measurement alone yields only a weighted average k = αk⊥ + βk∥; the weights are not independently measured but are inferred from the TR-MOKE values. Therefore, this agreement does not provide an independent validation of the anisotropy. The abstract's statement that 'the data obtained by both techniques reveal strong thermal anisotropy' overstates the joint evidence; the anisotropy claim rests solely on the TR-MOKE decomposition, which is itself the subject of the identifiability concern above. To support the claim, an independent measurement of the in-plane conductivity, or a quantitative estimate of α and β from the optothermal experiment's geometry (rather than from the TR-MOKE results), is needed.","section":"Conclusions, paragraph computing α and β"},{"comment":"The DFT assignment of the UV-Raman peaks is not conclusive. The three calculations give lower-mode frequencies of 417, 431, and 417 cm-1 and upper-mode frequencies of 483, 521, and 526 cm-1; no single calculation reproduces both experimental features at 431 and 467 cm-1 simultaneously. The upper-mode frequency varies by about 8% depending on supercell size and the U-J parameter, which the paper acknowledges. The conclusion that these modes 'appear to be the possible candidates' is too strong given this sensitivity. Moreover, the calculations are performed for bulk FePS3 only, whereas the UV-Raman peaks are observed in both FePS3 and MnPS3; no DFT evidence is presented for MnPS3 to support a similar assignment. The authors should either provide a systematic convergence study (e.g., a U-J sweep and phonon dispersion at relevant wavevectors) or substantially hedge the conclusion, noting that the 467 cm-1 feature currently lacks a robust theoretical assignment.","section":"Density functional theory calculations, Table 2 and surrounding text"},{"comment":"The DFT calculations assume the zigzag antiferromagnetic (z-AFM) spin texture, which is the low-temperature ground state, to interpret Raman spectra taken at room temperature, well above the Néel temperature of FePS3 (about 118 K). The paper does not justify why a low-temperature magnetic configuration is appropriate for phonon modes measured in the paramagnetic phase. Table 2 shows that the calculated mode frequencies depend strongly on spin texture (the nonmagnetic calculation gives 417/526 cm-1 versus the z-AFM values), so the choice of magnetic state is not a negligible detail. The nonmagnetic calculation does not match the observed 431/467 cm-1 pair either. Since the material is paramagnetic at the measurement temperature, a paramagnetic or disordered-spin calculation would be a more relevant comparator, or the authors should explicitly argue why the z-AFM frequencies should persist above the magnetic transition. As it stands, the phonon assignment is built on an unexamined assumption.","section":"Density functional theory calculations and Raman temperature dependence"}],"minor_comments":[{"comment":"The sentence 'In order to avoid self-heating effects, we used law laser excitation power in order to avoid self-heating effects' contains a typo ('law' for 'low') and a redundant phrase; it should be revised.","section":"Experimental section, Materials and characterization"},{"comment":"The statement that α > β 'based on the lateral dimensions of the flake size (>10 µm) and the laser spot size (~1.5 µm)' is qualitative. A quantitative estimate of the weight coefficients from the COMSOL model would be more compelling and would avoid the appearance that the weights are chosen post hoc.","section":"Results and Discussion, paragraph on optothermal thermal conductivity"},{"comment":"The column headers in Table 3 are ambiguous: the first column labeled 'k' appears to be used for both the orientation-average value (FePS3 row) and for isotropic or unknown-direction values in the comparison rows. Please clarify the column meaning, e.g., by labeling the columns 'k (average)', 'k⊥', and 'k∥' and explaining in the caption how values for the literature materials were categorized.","section":"Table 3"},{"comment":"The text states that error bars in Figure 7 represent uncertainty in thermal conductivity due to laser spot size uncertainty, but the figure as displayed does not show visible error bars. Please ensure the error bars are drawn and described in the caption.","section":"Figure 7 caption"},{"comment":"The paper refers to an analytical solution of the heat-diffusion equation in cylindrical coordinates from references [42-44] but does not state the specific equations, assumptions (e.g., uniform in-plane conductivity, boundary conditions, heat capacity values, film thicknesses used in the fit), or the exact values of modulation frequencies and spot sizes. Including this information would substantially improve reproducibility and would also help the reader assess the identifiability concern raised in the major comments.","section":"Methods, TR-MOKE description"},{"comment":"The text states that 'Full details of the calculations are provided in the Supplementary Materials,' but the supplementary files as provided contain only figures (XRD, optical, AFM, Raman) and no computational parameters such as pseudopotentials, energy cutoff, k-point sampling, or the phonon calculation method. This is an explicit promise that is not fulfilled. Please either add the missing computational details to the SI or correct the statement in the main text.","section":"Supplementary Information"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reports first-of-kind thermal conductivity measurements for FePS3 and MnPS3, which is potentially of interest to the materials community, but the central anisotropy claim is not yet supported because the TR-MOKE analysis lacks an identifiability study. The issue is fixable with additional analysis (sensitivity, additional modulation frequencies, or confidence regions) and is within the scope of a revision. The DFT section should be revised to either add the promised SI details and a systematic parameter study, or to clearly state that the UV peaks remain unassigned. I recommend major revision rather than rejection, as the underlying measurements may well be correct."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new physics here is the first experimental thermal conductivity of FePS3 and MnPS3, plus two UV-only Raman peaks near 430 and 467 cm−1 that visible Raman does not see. The thermal values are not pulled from a model: the optothermal Raman gives a direct weighted average for FePS3, and the TR-MOKE data are real measured signals. The internal consistency between the two techniques for FePS3 is a plus, and the paper is honest that the DFT mode assignment depends on U-J, supercell size, and spin texture. That honesty earns credit.\n\nThe soft spot is the TR-MOKE extraction. The paper says multiple modulation frequencies independently determine G, k⊥, and k∥, but only two frequencies (1 and 10 MHz) are used for three unknowns. The reader's concern about identifiability is legitimate: without a sensitivity analysis or confidence contours, the in-plane values might be trading off against G and k⊥. This is not a fatal flaw, but it is the load-bearing joint of the anisotropy claim, and MnPS3 has no optothermal cross-check, so its in-plane value of 6.3 W/mK is the most exposed number in the paper. The FePS3 consistency check using α=0.73, β=0.27 is reassuring but not an independent validation of the in-plane component.\n\nThe DFT assignment of the UV peaks is weaker. The 467 cm−1 peak is matched only by the z-AFM 4×4×1 supercell at 483 cm−1, while other parameter sets give 521–526 cm−1. That spread is large enough that the assignment is a suggestion, not a conclusion. The comparison of low-temperature z-AFM calculations with room-temperature data also has a phase mismatch that the paper acknowledges only implicitly.\n\nNone of this sinks the paper. The thermal conductivity measurements are a legitimate new dataset for an important family of antiferromagnetic 2D materials, and the Raman results will be useful to the community. What the paper needs before publication is a proper sensitivity analysis for the TR-MOKE fits, ideally including the correlation between G and the two conductivities, and a more guarded claim about the UV peak assignment.\n\nMy recommendation: send it to peer review. It deserves referee time. The authors should be asked for the TR-MOKE fitting details and a parameter study, not for a rewrite of the science. I would not want this desk-rejected; I would want it revised.","headline":"First thermal conductivity data for FePS3/MnPS3 with real new UV-Raman features, but the TR-MOKE anisotropy fit needs a sensitivity analysis before the headline numbers can be fully trusted.","tokens_in":20326,"tokens_out":757,"would_cite":true,"duration_ms":9306,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"First measurements of heat conduction in the antiferromagnetic semiconductors FePS3 and MnPS3 show heat flows 3 to 6 times faster along the layers than across them.","keywords":["antiferromagnetic semiconductors","quasi-2D materials","Raman spectroscopy","thermal conductivity","spintronics","time-resolved magneto-optical Kerr effect","FePS3","MnPS3"],"falsifier":"Repeat the TR-MOKE measurement on the same flakes with a series of pump spot sizes at fixed modulation frequencies; if the model's best-fit in-plane conductivity shifts by more than the quoted uncertainty when the spot size changes, the two-frequency decomposition is not unique. A stronger check is to measure the in-plane conductivity of a suspended exfoliated flake by a contact method such as a microbridge and compare the result with the 2.7 and 6.3 W m−1 K−1 values.","tokens_in":19230,"feed_emoji":"🔥","tokens_out":9472,"duration_ms":89411,"temperature":0.7,"pith_summary":"This paper reports the first measurements of thermal conductivity in thin crystals of the antiferromagnetic semiconductors FePS3 and MnPS3, which are candidate materials for spintronic devices. The measurements find that heat flows 3.2 times faster along the layers than across them in FePS3 (2.7 versus 0.85 W m−1 K−1) and 5.7 times faster in MnPS3 (6.3 versus 1.1 W m−1 K−1). These are low conductivities, and the authors argue the effect matters because the voltage produced by an antiferromagnetic spin-Seebeck device scales inversely with thermal conductivity. The paper also shows that UV-excited Raman spectroscopy reveals phonon features at about 430 and 467 cm−1 that visible-light Raman cannot see, and assigns them with density functional theory to lattice modes whose frequencies depend on the iron spin arrangement.","feed_headline":"Heat runs 3–6x faster along layers of FePS3 and MnPS3","feed_subtitle":"First measured heat flow in these 2D antiferromagnets gives engineers the numbers needed for spintronic device design.","key_machinery":"The thermal claim rests on two measurement techniques. Time-resolved magneto-optical Kerr effect (TR-MOKE) is a pump-probe method: a thin Co/Pt multilayer on the crystal surface is heated by a modulated laser train, and the polar Kerr rotation of a reflected probe beam tracks the surface temperature decay; the authors fit an analytical cylindrical heat-diffusion model to data taken at 1 and 10 MHz modulation to extract the interface conductance $G$, the through-plane conductivity $k_{\\perp}$, and the in-plane conductivity $k_{\\parallel}$. The steady-state optothermal Raman method uses the Raman peak position as a thermometer: separate calibrations of peak shift versus stage temperature and versus laser power give the coefficients that map laser heating into a thermal conductivity through finite-element solution of the steady heat equation. The phonon assignment uses density functional theory with Hubbard corrections on differently sized supercells and two spin textures, producing candidate modes in the 417–526 cm−1 range that bracket the observed 431 and 467 cm−1 features.","core_discovery":"The central discovery is that quasi-two-dimensional FePS3 and MnPS3 are strongly anisotropic and phonon-dominated heat conductors at room temperature, with through-plane conductivities of $0.85 \\pm 0.15$ and $1.1 \\pm 0.2$ W m−1 K−1 and in-plane conductivities of $2.7 \\pm 0.3$ and $6.3 \\pm 1.7$ W m−1 K−1, respectively. The authors state these are the first experimental values for these antiferromagnetic semiconductors and the first decomposition into in-plane and through-plane components. The anisotropy follows the expected bond picture: weak van der Waals interlayer bonds resist cross-plane transport while covalent intralayer bonds carry heat in-plane. The independent Raman optothermal measurement on FePS3 gave $1.35 \\pm 0.32$ W m−1 K−1, which matches the TR-MOKE components combined as $k = \\alpha k_{\\perp} + \\beta k_{\\parallel}$ with $\\alpha = 0.73$ and $\\beta = 0.27$. A second result is the appearance of a UV-only Raman peak near 467 cm−1 and a shoulder near 430 cm−1, which the calculations associate with modes of $B_u$ and $B_g$ symmetry whose frequencies shift with the magnetic order used in the supercell.","pith_inferences":["Because the extraction uses two modulation frequencies to fix three unknowns, a natural check is to repeat the fits with several pump spot sizes; if the fitted in-plane conductivity moves with spot size beyond the quoted uncertainty, the split between in-plane and through-plane terms would be suspect.","The computed dependence of the 467 cm−1 mode on the iron spin texture raises the possibility of using UV-Raman as a local, non-contact probe of antiferromagnetic order in thin flakes, though the paper does not claim this application.","The reported interface conductance of about 24 MW m−2 K−1, equivalent to tens of nanometers of glass, implies that in real devices the metal–MPX3 contact will dominate total through-plane thermal resistance; thermal design would need to focus on interfaces rather than only the crystal.","Extending the same TR-MOKE protocol to other MPX3 members such as NiPS3 and to monolayer flakes would map how the anisotropy and conductivity evolve with magnetic ion and thickness, testing whether phonon-boundary scattering suppresses the in-plane value at the monolayer limit."],"forward_implications":["A spin-Seebeck device using FePS3 or MnPS3 will develop larger temperature gradients for a given heating power than one built on higher-conductivity materials, directly changing the device's voltage output.","The through-plane values near 1 W m−1 K−1 place this material family among the most thermally insulating crystalline semiconductors, suggesting use as a heat barrier in nanoscale stacks.","The measured tensor gives a benchmark against which first-principles phonon calculations for the MPX3 family can be tested.","The consistency between the Raman optothermal average and the TR-MOKE components ($\\alpha = 0.73$, $\\beta = 0.27$) implies the two methods produce compatible values when the flake is much larger than the laser spot.","Because the metal transducer used in TR-MOKE is thin, the measured in-plane values are attributed to the crystal itself rather than to heat shunting through the metal film."],"supporting_citations":[{"why":"provides the time-domain thermoreflectance heat-diffusion model and measurement approach that TR-MOKE adapts","marker":"[42]"},{"why":"tutorial on time-domain thermoreflectance that underlies the pump-probe extraction of thermal properties","marker":"[45]"},{"why":"demonstrates anisotropic thermal conductivity measurement in layered materials and supplies a comparison framework","marker":"[46]"},{"why":"introduced the optothermal Raman method used for the orientation-average thermal conductivity","marker":"[8]"},{"why":"supplies the optical absorption coefficients that explain the different Raman heating response of FePS3 and MnPS3","marker":"[2]"},{"why":"predicts very low lattice thermal conductivity in the related compound ZnPSe3 and motivates thermal study of the MPX3 family","marker":"[19]"},{"why":"provides the first-principles phonon calculations used to assign the newly observed UV-Raman modes","marker":"[32]"}],"fun_headline_variants":["First heat measurements reveal 3–6x anisotropy in 2D antiferromagnets","Heat flows 3–6x faster in-plane in 2D antiferromagnets","UV Raman uncovers new modes in 2D antiferromagnets","Anisotropic heat flow mapped in FePS3 and MnPS3 films","2D antiferromagnets: heat prefers in-plane by 3–6x"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole separation of in-plane from through-plane heat flow rests on the assumption that measurements at just two pump modulation frequencies, 1 and 10 MHz, can uniquely determine three unknowns: the interface conductance, the cross-plane conductivity, and the in-plane conductivity; the paper presents no sensitivity or identifiability analysis for that split.","fun_headline_variants_meta":{"raw":{"variants":["First heat measurements reveal 3–6x anisotropy in 2D antiferromagnets","Heat flows 3–6x faster in-plane in 2D antiferromagnets","UV Raman uncovers new modes in 2D antiferromagnets","Anisotropic heat flow mapped in FePS3 and MnPS3 films","2D antiferromagnets: heat prefers in-plane by 3–6x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001175,"raw_usage":{"total_tokens":4931,"prompt_tokens":1093,"completion_tokens":3838,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":3726}},"tokens_in":709,"tokens_out":3838,"duration_ms":25755,"temperature":1.0,"reasoning_tokens":3726,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:21:40.932930+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the TR-MOKE measurement on the same flakes with a series of pump spot sizes at fixed modulation frequencies; if the model's best-fit in-plane conductivity shifts by more than the quoted uncertainty when the spot size changes, the two-frequency decomposition is not unique. A stronger check is to measure the in-plane conductivity of a suspended exfoliated flake by a contact method such as a microbridge and compare the result with the 2.7 and 6.3 W m−1 K−1 values.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"demonstrates anisotropic thermal conductivity measurement in layered materials and supplies a comparison framework"}],"review_version":1}