{"id":"6dc1d0a9-7bf8-47a8-9c3d-3aa74be80539","arxiv_id":"2506.10260","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"First-principles Wigner transport calculations show Sn2S3 has anisotropic, low lattice thermal conductivity caused by rattling of lone-pair-bearing Sn(II) atoms, with optical phonons carrying most heat along the b-axis.","lead":"Using first-principles simulations, the authors predict that the mixed-valent compound Sn2S3 conducts heat poorly along two of its three crystal axes, with the weakly bonded Sn(II) atoms rattling inside the lattice. The work suggests that mixed-valent compounds with lone-pair electrons are promising targets for low-thermal-conductivity thermoelectric materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"AIMD sampling is too short to converge the Sn(II) rattler's IFCs; the quantitative kappa tensor and T^-0.68 exponent are therefore not yet established.","rationale":"I agree with the reader that the quantitative claims depend critically on the finite-temperature IFCs. The paper's independent evidence for rattling—MSD, frozen-potential curves, pCOHP, and normalized IFC anisotropy—supports the mechanism, so that part of the central claim is robust. What is load-bearing and unsecured is the specific thermal conductivity tensor and its temperature exponent, which are the headline numerical results. The methods section reports a single 3 ps AIMD run per temperature for TDEP/hiphive extraction; the separate 50,000-step trajectory run is used for displacements and MSD, not for IFC fitting, so it does not cure the sampling problem. Because phonon lifetimes and the Wigner coherence term depend sensitively on both the harmonic reference and the linewidths (Eqs. 2-3), a biased IFC set could change kappa substantially and alter the apparent T dependence. The proposed 30 ps test is computationally heavy but feasible for a 320-atom cell and directly settles whether the reported numbers are converged. Since this is exactly the condition that should gate acceptance of the quantitative claims, I do not move the reader's conditional verdict; I would keep it unchanged pending this check.","tokens_in":13323,"tokens_out":5985,"duration_ms":74867,"concrete_test":"Recompute the TDEP IFCs from new AIMD runs at 300 K and 500 K with at least 30 ps of production time (discarding the first 5 ps) in the same 2x4x2 supercell, then rerun the 3ph+4ph+OD Wigner calculation on the same 11x26x7 q-grid. Compare kappa_300 K and the fitted T exponent to the reported values. If the 300 K tensor changes by less than ~10% and T^-0.68 is reproduced within fitted uncertainty, the sampling concern is settled; if not, the quantitative claims require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims—kappa = 1.73/6.55/1.59 W/m/K at 300 K, the 3.86 anisotropy ratio, the 63% optical contribution along b, and the weak T^-0.68 dependence—all inherit from TDEP/hiphive IFCs extracted from only 3 ps of AIMD in a 2x4x2 supercell (Methods, Computational details). At a 2 fs time step, 3 ps yields 1500 snapshots, and for low-frequency Sn(II) rattler modes around 2-11 meV (periods roughly 0.4-2 ps) this covers only a few vibrational periods. No equilibration discard is reported, and consecutive snapshots are strongly correlated, so the effective sample for fitting cubic and quartic IFCs in a 320-atom cell is small. If the slow, large-amplitude Sn(II) motion is undersampled, the harmonic reference and anharmonic linewidths are biased, which would change phonon lifetimes, the Wigner coherence term, and the fitted temperature exponent. The qualitative rattling mechanism is separately supported by MSD, pCOHP, and potential-energy evidence, so the mechanism is not at risk; the specific low-kappa numbers and weak-T scaling, however, are not yet reliable.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports first-principles calculations of lattice thermal conductivity in mixed-valent Sn2S3, a quasi-1D van der Waals material, using temperature-dependent effective potentials (TDEP), with phonon transport solved via the Boltzmann transport equation and the Wigner formulation. The authors predict strongly anisotropic conductivities of 1.73 (a), 6.55 (b), and 1.59 (c) W m-1 K-1 at 300 K, a weak temperature dependence of T^-0.68 in the full Wigner treatment, and a dominant (63%) optical-phonon contribution to conduction along the b-axis. The low and anisotropic transport is attributed to rattling of the lone-pair-bearing Sn(II) atoms, supported by mean-square-displacement analysis, potential-energy barriers, projected COHP, low average phonon frequencies, and large Grüneisen parameters for Sn(II)-dominated modes.","tokens_in":13586,"tokens_out":6414,"duration_ms":73738,"significance":"If the quantitative results are reliable, the paper makes a useful contribution by extending the rattling paradigm to mixed-valent compounds and by demonstrating, with standard and well-established computational tools, that optical phonons can dominate heat conduction along a specific crystallographic direction. The qualitative mechanism is supported by several independent and parameter-free analyses (MSD, pCOHP, potential wells, frozen-phonon potentials), which is a genuine strength. However, the numerical values that carry the central claims depend sensitively on interatomic force constants extracted from extremely short AIMD trajectories, and the lack of any comparison with available experimental thermal-conductivity data leaves the quantitative predictions unvalidated. The paper is therefore of moderate significance in its current form; the mechanism is plausible but the numbers are not yet established.","major_comments":[{"comment":"The TDEP interatomic force constants (IFCs) used for the thermal-conductivity calculations are extracted from AIMD runs of only 3 ps (1500 steps at a 2 fs time step) in a 2x4x2 supercell. For the low-frequency Sn(II) rattling modes with energies around 2-11 meV, whose periods are roughly 0.4-2 ps, this amounts to only a few vibrational periods, and no equilibration discard or convergence test with respect to simulation length is reported. Because the reported conductivity tensor (1.73/6.55/1.59 W m-1 K-1), the T^-0.68 dependence, and the 63% optical contribution along b all derive from these IFCs, the quantitative reliability of the central claims is not established. The authors should increase the AIMD sampling (for example, to at least 20 ps), test convergence of the extracted cubic and quartic IFCs, and preferably demonstrate that the predicted conductivity is stable with respect to simulation length.","section":"Methods, Computational details"},{"comment":"The paper cites Ref. [41] as an experimental measurement of low thermal conductivity of Sn2S3, and the SI describes thermal-diffusivity measurements on a single-crystal sample, but no experimental conductivity values are given anywhere and no comparison is made with the calculated kappa. Without a quantitative benchmark to experiment, the claim of 'low and anisotropic' conductivity is not contextualized, and the accuracy of the computational pipeline cannot be assessed. The authors should provide the experimental values (from Ref. [41] and from their own SI measurements) and compare them with the computed tensor.","section":"Results and Discussion, first paragraph; SI Fig. S6"},{"comment":"Two different AIMD simulations are described: one of 3 ps in a 2x4x2 supercell (320 atoms) for TDEP extraction, and another of 50000 steps in a 2x3x2 supercell (160 atoms) at 300 K for atomic trajectories and MSD analysis. While these are distinct runs serving different purposes, the manuscript does not clearly state that, and the juxtaposition of '3 ps' and '50000 steps' is confusing. Please clarify which simulation was used for which analysis and ensure the reader is not misled about the length of the MD used for the IFCs.","section":"Methods, Computational details and Results, third paragraph"}],"minor_comments":[{"comment":"The text refers to 'COHP analysis as shown in Fig. 1b', but Fig. 1b displays XRD patterns; the actual pCOHP plot appears only in SI Fig. S2. Either move the COHP panel to Fig. 1 or correct the cross-reference.","section":"Fig. 1"},{"comment":"There are several typographical errors, including 'Brillion' (Brillouin), 'stats' (states), 'oose-Einstein' (Bose-Einstein) in the SI, and 'avoid-crossing' should be 'avoided crossing'. A careful proofread is needed.","section":"Throughout"},{"comment":"The equation for the wave-like (coherence) contribution contains a redundant denominator factor; please check it against the standard Wigner transport expression (e.g., Simoncelli et al., Phys. Rev. X 12, 041013 (2022)).","section":"Eq. (3)"},{"comment":"The statement that '63% of the thermal conductivity along the b-axis is contributed by optical phonons' is made without defining how acoustic and optical modes are separated in the cumulative analysis. Specify the criterion (e.g., based on the phonon branch index or on the eigenvector projection) used to partition the conductivity.","section":"Fig. 3c and text"},{"comment":"The average-frequency analysis reports values for the total system (22.57 meV), Sn(II) (11.21 meV), and S (28.91 meV), but omits the average frequency for Sn(IV). Providing this value would make the comparison between Sn(II) and Sn(IV) more complete.","section":"Results, first paragraph of 'Sn(II) rattler' discussion"},{"comment":"The definition of the normalized trace of the IFC is difficult to parse because of missing parentheses and indices; please rewrite the formula in a clean, standard notation.","section":"Eq. (4)"},{"comment":"The manuscript does not report the convergence of the thermal conductivity with respect to the q-point grid beyond stating that 'a 11x26x7 q-point grid was used'. A table in the SI showing kappa as a function of grid size would strengthen the reliability of the reported values.","section":"Methods, Computational details"},{"comment":"The SI states that isotope scattering is included, but the main text never mentions this. Add a sentence in the Methods section acknowledging the inclusion of isotope effects.","section":"SI, first paragraph"}],"recommendation":"major_revision","confidential_remarks":"The paper's central mechanistic claim (rattling of Sn(II) driven by lone pairs) is plausible and supported by multiple qualitative indicators, but the quantitative results are built on an exceptionally short AIMD trajectory for IFC extraction (3 ps). This is a common pitfall in TDEP-based studies, and the authors can likely fix it with longer simulations and convergence checks. The absence of any experimental comparison for kappa is also problematic and should be remedied. The paper is within the scope of the journal; with the required revisions it could become a solid contribution. I would ask the editor to encourage the authors to address the sampling and benchmarking issues thoroughly rather than treating them as minor."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth knowing: this is a solid, standard first-principles phonon transport paper on Sn2S3, and the central qualitative claim — that Sn(II) rattling drives the low thermal conductivity — is probably right. The numbers, however, are not yet established.\n\nWhat is genuinely new: this is the first calculation I know of that combines four-phonon scattering and Wigner coherence for Sn2S3, and the attribution of the low kappa to Sn(II) rattling is supported by several independent indicators: pCOHP shows the lone-pair anti-bonding character, the potential-energy barriers are lowest for Sn(II), the MSD is anisotropic and largest for Sn(II), the IFC traces show weak bonding, and the Gruneisen parameters are large on the low-frequency Sn(II) modes. That is a well-rounded mechanistic case, and it is not circular — nothing is fitted to the target kappa.\n\nThe soft spots are real but not fatal. The main one: the IFCs used for the quantitative kappa tensor (1.73/6.55/1.59 W/m/K at 300 K, the 3.86 anisotropy, the 63% optical contribution, and the T^-0.68 exponent) come from only 3 ps of AIMD in a 320-atom cell. Sn(II) rattler modes have periods of order 0.4–2 ps, so that is only a few vibrational periods. No equilibration discard or convergence tests are reported. That makes the computed lifetimes and the Wigner term provisional. The qualitative mechanism is not at risk, but the specific numbers should be treated with caution.\n\nSecond, the paper never compares its computed kappa to the experimental value from Ref. [41]. That is an obvious and necessary benchmark, and its absence is conspicuous. Third, there is a technical inconsistency: the Methods say 3 ps / 320 atoms for the TDEP extraction, while the Results describe a separate 50000-step / 160-atom run for trajectories. The authors need to clarify which run feeds which quantity.\n\nMinor issues: the DFT-D1 dispersion correction is dated, and the SI mentions experimental thermal diffusivity measurements that are never connected to the calculated values. These are easy fixes.\n\nOverall: a capable study with a plausible mechanism and provisional quantitatives. It deserves peer review, but a referee should push for longer AIMD, convergence checks on the IFC fit, and a direct comparison with experiment before the numbers are accepted.","headline":"A competent anharmonic-phonon study with a plausible Sn(II) rattling mechanism, but the quantitative kappa values rest on very short AIMD and are never benchmarked against experiment.","tokens_in":14102,"tokens_out":1965,"would_cite":false,"duration_ms":26124,"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":"Rattling Sn(II) atoms, pushed by lone pairs, are the paper's explanation for Sn2S3's low and weakly temperature-dependent lattice thermal conductivity.","keywords":["Sn2S3","mixed-valent compound","rattling","lone pair","lattice thermal conductivity","Wigner transport equation","phonon coherence","anisotropic thermal transport"],"falsifier":"Measure the lattice thermal conductivity of a single crystal of Sn2S3 along the a, b, and c axes from 300 to 700 K: if the temperature exponent is close to $T^{-1}$ rather than $T^{-0.68}$, or if the c-axis value is far from about 1.6 W/m/K, the rattling/Wigner explanation would be in trouble. Inelastic neutron or X-ray scattering that fails to find the predicted flat low-frequency optical branches below about 12 meV and their avoided crossing would also sever the proposed mechanism.","tokens_in":13149,"feed_emoji":"🌡️","tokens_out":11617,"duration_ms":128111,"temperature":0.7,"pith_summary":"This paper argues that the mixed-valent compound Sn2S3 owes its intrinsically low and weakly temperature-dependent lattice thermal conductivity to rattling of its Sn(II) atoms, caused by repulsion from lone-pair electrons. Using molecular dynamics simulations and a full Wigner treatment of phonon transport, it computes lattice thermal conductivities of about 1.73, 1.59, and 6.55 W/m/K along the a, c, and b axes at 300 K, with the temperature dependence flattening to $T^{-0.68}$ once wave-like coherence contributions are included. The authors show that Sn(II) rattling generates flat low-frequency optical phonons that scatter acoustic heat carriers, and that along the strongly bonded b-axis roughly 63% of the heat is carried by optical phonons. The result matters because it identifies mixed-valent compounds with lone-pair cations as a design route to low-thermal-conductivity thermoelectrics and shows that optical phonons can dominate heat conduction in quasi-1D materials.","feed_headline":"Rattling tin atoms keep Sn2S3's heat flow low and flat","feed_subtitle":"Lone-pair wobble yields ~1.6 W/m/K along c and pushes 63% of chain heat into optical phonons","key_machinery":"The load-bearing object is the rattling mode of the Sn(II) sublattice in quasi-1D Sn2S3: weakly bonded, lone-pair-bearing atoms whose large anisotropic displacements produce flat, low-frequency optical phonon branches. The quantitative engine is the Wigner transport equation, which adds off-diagonal wave-like (coherence) heat-flux terms to the standard Boltzmann particle picture; this addition is what changes the temperature dependence from about $T^{-1}$ to $T^{-0.68}$. The paper also uses the normalized trace of interatomic force constants to show that Sn(II)-S bonding is weaker and more anisotropic than Sn(IV)-S bonding, and mode-resolved group velocities to show that optical branches along the b-axis are fast enough to carry heat.","core_discovery":"The paper's central claim is that the intrinsically low and weakly temperature-dependent lattice thermal conductivity of Sn2S3 originates from rattling of Sn(II) atoms, not from ordinary acoustic-phonon scattering. Lone-pair electrons on adjacent Sn(II) atoms repel the Sn(II) sublattice, making Sn(II)-S bonds weak and anisotropic; Sn(II) then vibrates with much larger and more anisotropic displacements than Sn(IV) or S, at a low average frequency of 11.21 meV compared with 28.91 meV for S. These rattling vibrations create flat low-frequency optical branches that hybridize with acoustic phonons and show an avoided crossing, producing strong anharmonicity. Including both particle-like and wave-like heat conduction via the Wigner transport equation gives lattice thermal conductivities of 1.73, 1.59, and 6.55 W/m/K along a, c, and b at 300 K, an anisotropy ratio of 3.86, and a weak temperature dependence $T^{-0.68}$ rather than the conventional $T^{-1}$; along b, optical phonons carry about 63% of the total heat.","pith_inferences":["Inferred: if Sn(II) rattling is the cause, chemical substitutions that tighten the Sn(II) cage or quench the lone-pair repulsion should raise the lattice thermal conductivity and steepen its temperature dependence back toward $T^{-1}$; this is directly testable but not computed in the paper.","Inferred: the same lone-pair rattling mechanism should appear in other mixed-valent compounds with stereochemically active lone pairs, such as Sn2Se3 or Pb-based analogues, making the $T^{-0.68}$ signature and the large optical-phonon share a screening target for future materials searches.","Inferred: because the wave-like coherence contribution here is tied to high-frequency optical phonon pairs, isotope substitution or pressure that shifts those optical frequencies should change the coherence term more than the particle-like term, giving an experimental handle on the Wigner channel."],"forward_implications":["If the rattling picture is correct, Sn2S3's lattice thermal conductivity stays low at high temperature, so its thermoelectric figure of merit should deteriorate more slowly than in conventional $T^{-1}$ phonon-gas materials.","The heat-flow anisotropy ratio of 3.86 means crystal orientation matters: along a and c, heat is throttled by van der Waals gaps, while along b, optical phonons carry most of the heat.","Optical phonons cannot be neglected in quasi-1D and mixed-valent compounds; the common assumption that acoustic modes dominate lattice thermal conductivity fails along the b-axis of Sn2S3.","The rattling fingerprint—low average frequency, large anisotropic atomic displacement, and flat low-energy optical branches—offers a screening criterion for finding other low-thermal-conductivity mixed-valent thermoelectrics."],"supporting_citations":[{"why":"supplies the Wigner transport equation used for both particle-like and wave-like (coherence) contributions, including the weak $T^{-0.68}$ dependence","marker":"[51,52]"},{"why":"provides the temperature-dependent effective potential method for extracting finite-temperature interatomic force constants from molecular dynamics trajectories","marker":"[46,47]"},{"why":"used to extract the harmonic, cubic, and quartic force constants that feed the scattering calculations","marker":"[48]"},{"why":"solves the phonon Boltzmann transport equation for the particle-like thermal conductivity","marker":"[49]"},{"why":"adds four-phonon scattering rates to the transport calculation","marker":"[50]"},{"why":"reports the experimentally measured thermal and thermoelectric properties of Sn2S3 that the calculations aim to explain","marker":"[41]"},{"why":"supplies the phonon dispersion and density of states used to identify the flat low-frequency optical branches and the rattling signature","marker":"[53]"},{"why":"defines the normalized trace of interatomic force constants used to show that Sn(II)-S bonding is weak and anisotropic","marker":"[15]"}],"fun_headline_variants":["Rattling Sn(II) atoms give Sn2S3 low, anisotropic heat flow","Lone-pair repulsion rattles Sn(II), slashing Sn2S3 thermal conductivity","Sn2S3's 1.6 W/mK along c from rattling tin atoms","Optical phonons carry 63% of heat in Sn2S3 due to rattling Sn(II)","Weak bonds and lone pairs make Sn2S3 anisotropic heat insulator"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands on the assumption that a 3-picosecond simulation of 320 atoms samples enough of the slow, large-amplitude wobbling of Sn(II) to determine how it scatters phonons; if that wobble is slower or more collective than the simulation window, the computed heat-flow values and their temperature dependence could change.","fun_headline_variants_meta":{"raw":{"variants":["Rattling Sn(II) atoms give Sn2S3 low, anisotropic heat flow","Lone-pair repulsion rattles Sn(II), slashing Sn2S3 thermal conductivity","Sn2S3's 1.6 W/mK along c from rattling tin atoms","Optical phonons carry 63% of heat in Sn2S3 due to rattling Sn(II)","Weak bonds and lone pairs make Sn2S3 anisotropic heat insulator"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000609,"raw_usage":{"total_tokens":2852,"prompt_tokens":975,"completion_tokens":1877,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":1762}},"tokens_in":591,"tokens_out":1877,"duration_ms":15911,"temperature":1.0,"reasoning_tokens":1762,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:31:41.931815+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the lattice thermal conductivity of a single crystal of Sn2S3 along the a, b, and c axes from 300 to 700 K: if the temperature exponent is close to $T^{-1}$ rather than $T^{-0.68}$, or if the c-axis value is far from about 1.6 W/m/K, the rattling/Wigner explanation would be in trouble. Inelastic neutron or X-ray scattering that fails to find the predicted flat low-frequency optical branches below about 12 meV and their avoided crossing would also sever the proposed mechanism.","supporting_citations":[{"cited_title":"Eriksson, E","cited_arxiv_id":null,"evidence_quote":"used to extract the harmonic, cubic, and quartic force constants that feed the scattering calculations"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"solves the phonon Boltzmann transport equation for the particle-like thermal conductivity"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"adds four-phonon scattering rates to the transport calculation"},{"cited_title":"Saito, K","cited_arxiv_id":null,"evidence_quote":"reports the experimentally measured thermal and thermoelectric properties of Sn2S3 that the calculations aim to explain"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the normalized trace of interatomic force constants used to show that Sn(II)-S bonding is weak and anisotropic"}],"review_version":1}