{"id":"5a2f3039-a221-474d-8b11-eac25b517a82","arxiv_id":"2607.25591","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Rocksalt BaS–PbS superlattices, superstructures, wires, and solid solutions are predicted to be ferroelectric or antiferroelectric through a chain-structure TO phonon instability.","lead":"First-principles calculations predict ferroelectric, antiferroelectric, and mixed polar order in BaS–PbS rocksalt structures, caused by unstable vibrations in Pb–S chains stretched by barium. The work suggests a new family of non-perovskite ferroelectrics and a complex multi-state energy landscape relevant to memory and relaxor behavior.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PBEsol strain-sensitive phonons are the key risk; verification should target marginal phases.","rationale":"Reader's weakest_assumption and my concern overlap. I set agreement to partial because the direction of the 0.6% lattice error weakens the reader's specific worry: a larger experimental cell would add tensile strain, the very driver of the chain instability. However, functional errors can still affect the phonon frequencies and relative energies directly, and the paper's strongest claims include an essentially infinite set of near-degenerate minima. The best single test is a functional/volume sensitivity calculation on a few representative supercells, including one marginal case. If the instability is robust there, the central claim is supported; if not, the paper should be read cautiously. The paper's internal checks (positive phonons, elastic stability, L11 absence of instability, positive Pb on-site force constant) are good evidence of self-consistency. The recommendation is to keep the CONDITIONAL verdict: accept the qualitative discovery while requiring the sensitivity check and a more careful statement about nonergodicity and the 'infinite number of minima' (which goes beyond the finite set of supercells calculated).","tokens_in":13053,"tokens_out":7946,"duration_ms":86748,"concrete_test":"Recompute the zone-center TO phonon and the relaxed FE/AFE phase energies for Ba3PbS4 (L12), BaPbS2 (L10), (BaS)3/(PbS)1, and (BaS)1/(PbS)5 using the same ONCVPSP pseudopotentials but with the SCAN functional (or HSE for the zone-center mode). Also repeat the PBEsol calculations with the cell volume uniformly expanded by +0.6% and internal coordinates relaxed to mimic experimental lattice constants. If the first three structures retain the soft mode and FE/AFE energy gains of several meV while only the marginal fourth changes, the central discovery survives; if the instabilities vanish across the board, PBEsol is the critical failure and the verdict should be conditional at best.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The claim depends on PBEsol DFT phonon instabilities in strain-sensitive rocksalt superstructures. Section II reports PBEsol lattice parameters 0.6% smaller than experiment, and Fig. 6 shows the TO-mode frequency varies strongly with lattice parameter; a modest functional error could relocate the stability boundary. This is a real risk, but it is not uniform: Pb-S chains in these structures are stretched by Ba substitution, and an increase in the lattice parameter (toward experiment) would stretch them further, so the robust instabilities (e.g., 29.45 meV in Ba3PbS4, Table I) should survive or deepen. The real vulnerability is the proliferation of near-degenerate phases with tiny energy gains—(BaS)1/(PbS)5 at -0.004 meV, six phases within 0.5 meV in Table IV—where PBEsol accuracy, zero-point motion, and thermal fluctuations decide whether the phase exists. The positive-phonon and elastic checks (§III-B) show internal consistency but do not establish transferability. Hence the load-bearing assumption is that PBEsol is quantitatively predictive for marginal chain-order energies, not merely that it finds a qualitative instability.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports first-principles DFT/DFPT calculations predicting ferroelectric (FE), antiferroelectric (AFE), and mixed FE-AFE instabilities in ordered superstructures, superlattices, quantum wires, and disordered solid solutions of the BaS-PbS system with the rocksalt structure. The instability is attributed to softening of transverse-optical phonons in linear -Pb-S-Pb-S- chains, driven by tensile strain from Ba substitution. The authors identify numerous low-energy phases with near-degenerate energies, propose a multi-minimum potential landscape, and suggest possible low-temperature nonergodicity. They also extend the prediction to selenide and telluride analogs.","tokens_in":13297,"tokens_out":4069,"duration_ms":41946,"significance":"If correct, the paper identifies a new family of ferroelectric/antiferroelectric materials in a simple rocksalt system, expanding the known structural classes beyond perovskites and IV-VI ferroelectrics. The work is methodologically solid: it uses standard DFPT, checks phonon and elastic stability, cross-validates against strained bulk PbS and prior FP-LAPW data, and provides parameter-free predictions (no fitted parameters, machine-checked stability criteria). The qualitative mechanism—chain-structure instability from Ba-induced tensile strain—is plausible and supported by the absence of instability in the CuPt structure where chains are broken. However, the quantitative claims of near-degenerate phases and an effectively infinite number of minima rest on energy differences of tens to hundreds of micro-electronvolts, which are at the limit of DFT accuracy, particularly with PBEsol's known lattice-parameter underestimate. The paper is a significant prediction but requires additional validation of the marginal phase stabilities before the strongest conclusions can be accepted.","major_comments":[{"comment":"The central claim of a multi-minimum potential relies on energy differences of 0.4 meV or less among the six phases listed in Table IV (e.g., Cmc21 at -29.85 meV vs R3m at -29.45 meV). PBEsol underestimates lattice constants by ~0.6% (Sec. II), and Fig. 6 shows the TO instability is strongly strain-dependent. A functional error of this magnitude could easily reorder these phases or eliminate the marginal ones (e.g., (BaS)1/(PbS)5 with -0.004 meV in Table II). The authors should provide an estimate of the numerical and functional uncertainty on these energy differences, and ideally test a subset of phases with a more accurate functional (e.g., HSE) or with experimental lattice constants. Without this, the specific ground-state assignments and the claim of 'virtually infinite' near-degenerate minima are not quantitatively robust.","section":"§III.B, Table IV"},{"comment":"The inference from a dense set of unstable q-points to 'a virtually infinite number of minima in the configuration space, separated by potential barriers' is a logical leap. A continuum of harmonic instabilities does not guarantee that each q vector yields a distinct metastable minimum; anharmonic couplings could merge them into a single soft valley or create only a few basins. The paper demonstrates several distinct minima but does not map the energy surface along different q vectors or compute barriers. The 'infinite minima' claim is central to the nonergodicity conclusion and should either be demonstrated by explicit calculations for several q vectors or clearly reframed as a speculative extrapolation.","section":"§III.B, last paragraph"},{"comment":"The strain sensitivity of the TO mode is used to support the mechanism, but it also highlights a vulnerability. PBEsol's 0.6% lattice-parameter underestimate means the calculations are performed at a slightly compressed state compared to experiment. For strongly unstable structures (e.g., Ba3PbS4, -29.45 meV), the instability is robust because stretching would deepen it. However, for marginal cases such as (BaS)1/(PbS)5 (-0.004 meV) or the 3BaS/7PbS SL (-0.161 meV), the instability could be an artifact of the underestimated lattice constant. The paper should either test these marginal structures at experimental lattice parameters or otherwise quantify the strain error margin. This is load-bearing because the purported coexistence of FE, AFE, and mixed phases at near-zero energy cost underpins the nonergodicity prediction.","section":"§III.C, Fig. 6"}],"minor_comments":[{"comment":"The abstract and conclusion emphasize 'discovered' phenomena; it would be clearer to say 'predicted' since no experimental confirmation is provided.","section":"Abstract/Introduction"},{"comment":"Typo: 'calcualated' should be 'calculated'. Also, provide more detail on the pseudopotentials (e.g., valence-electron configurations) and convergence checks for k-point density and plane-wave cutoff.","section":"Sec. II"},{"comment":"Several rows lack polarization values for some phases (e.g., 7BaS/1PbS Amm2, 3BaS/1PbS Amm2). Either fill these in or explicitly state they were not calculated; otherwise the table appears incomplete.","section":"Table II"},{"comment":"The figure has many data points from different structure types; using distinct symbols or a legend would improve readability. The 'quantum wire' point is mentioned in the text but not clearly identified in the figure.","section":"Fig. 4"},{"comment":"The disordered solid solutions are modeled with SQS8 and SQS16, but no supercell sizes or number of atoms are given. This information is essential for reproducibility.","section":"Sec. II / III.B"},{"comment":"The on-site force constant for Pb is reported as +0.02290 Ha/Bohr^2 for one superstructure with 12.5% PbS. Providing the corresponding values for other compositions or for the selenide/telluride analogs would strengthen the conclusion that Pb is not off-center.","section":"Sec. III.C"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the qualitative prediction is likely correct, given the robust instabilities in several structures and the cross-checks against strained bulk PbS. However, the central quantitative claim—near-degenerate FE/AFE/mixed phases and an effectively infinite multi-minimum landscape—rests on energy differences that are within typical DFT error bars for PBEsol. I recommend requesting a robustness analysis: e.g., testing a few marginal phases with a hybrid functional or at experimental lattice constants, and providing error estimates for the energy differences. The 'infinite minima' inference also needs justification beyond the presence of a dense set of unstable q-vectors."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The bottom line: this is a credible first-principles discovery of a new ferroelectric mechanism in rocksalt chalcogenides, but the paper stretches the evidence when it talks about infinite minima and nonergodicity.\n\nThe genuinely new result is the prediction that BaS–PbS rocksalt superstructures and superlattices have soft TO phonons in linear –Pb–S–Pb–S– chains, leading to ferroelectric, antiferroelectric, and mixed chain orders. That is a real extension of the perovskite chain-instability idea to a different structural family, and the paper supports it with thoughtful checks: phonon spectra of the relaxed low-symmetry phases are positive, elastic constants are positive-definite, the instability tracks strained bulk PbS, the CuPt structure (which cannot host infinite Pb–S chains) is stable, and an isolated quantum wire also shows the instability. The link to Ba-induced tensile strain is well argued, and the on-site Pb force constant calculation reasonably rules out simple lone-pair off-centering. The methodology is standard and the numbers are internally consistent.\n\nWhere I would push back: the energy landscape is not as solid as the qualitative mechanism. The paper lists several phases with energy differences below 0.5 meV, and one superlattice, (BaS)1/(PbS)5, with a stabilization of 0.004 meV per molecule. Those numbers are below the likely accuracy of PBEsol, especially given the reported 0.6% underestimate of the lattice parameters. The authors themselves note that the task of finding the ground states in all studied structures is practically unsolvable, and that a dense mesh of wavevectors implies an infinite number of minima—that is an extrapolation, not a calculation. Zero-point motion and thermal fluctuations could easily wash out the shallowest minima. The claim of nonergodicity is suggestive but not demonstrated. I also miss input files or detailed convergence data; for a prediction of a new ferroelectric family, a functional-sensitivity test (e.g., SCAN or HSE on the key superstructures) would go a long way.\n\nOverall: the central qualitative claim holds up well; the quantitative ground-state hierarchy and the speculative physics do not. A serious referee should engage with it, and the paper is worth publishing if the authors temper the overclaim and provide additional validation.","headline":"Credible first-principles discovery of a new ferroelectric mechanism in rocksalt BaS–PbS, but the infinite-minima and nonergodicity claims outrun the evidence.","tokens_in":13766,"tokens_out":2302,"would_cite":true,"duration_ms":22453,"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":"Barium-induced stretching of Pb–S chains in rock-salt BaS–PbS produces ferroelectric, antiferroelectric, and mixed phases at nearly equal energies, according to first-principles calculations.","keywords":["ferroelectricity","antiferroelectricity","rocksalt structure","chain-structure instability","BaS-PbS","soft phonon mode","first-principles calculations","nonergodicity"],"falsifier":"Synthesize an ordered Ba3PbS4 sample and measure its low-temperature structure and dielectric response; if no polar or antipolar distortion appears near the estimated ~300 K transition and no polarization develops, the central claim fails. A simpler check is computational: recompute the phonon spectrum at the experimental lattice constant; if the zone-center TO mode is stable there, the predicted ferroelectricity is an artifact of the lattice underestimate.","tokens_in":12927,"feed_emoji":"⚡","tokens_out":8391,"duration_ms":77489,"temperature":0.7,"pith_summary":"The paper aims to establish that the BaS–PbS system in the simple rock-salt structure is ferroelectric, even though neither BaS nor PbS is ferroelectric on its own. Using first-principles calculations, it shows that barium's larger ionic radius stretches the lattice and destabilizes transverse optical phonons in infinite –Pb–S–Pb–S– chains, so the chains polarize collectively. The same calculations reveal stable antiferroelectric and mixed ferroelectric–antiferroelectric phases whose energies are nearly identical to the ferroelectric phases, so the ground state is often not purely ferroelectric. A sympathetic reader would care because this would expand ferroelectricity to a common, simple crystal structure and imply a glass-like, multi-well energy landscape with many switchable polarization configurations.","feed_headline":"Calculations show barium strain makes rock-salt BaS-PbS ferroelectric","feed_subtitle":"First-principles prediction opens a new switchable material class beyond perovskites.","key_machinery":"The load-bearing mechanism is the chain-structure instability of the transverse-optical phonon in infinite –Pb–S–Pb–S– chains running along the fourfold axes of the NaCl lattice. An optical phonon is a vibration in which different atoms move in opposite directions; here the unstable eigenvector consists of out-of-phase Pb and S displacements along a chain, so each chain develops a one-dimensional electric polarization. Adjacent chains interact weakly and with a sign that favors antiparallel alignment in BaS–PbS. This single mechanism explains the soft mode at the zone center, the instabilities at zone-boundary points, the existence of FE, AFE, and mixed phases, and the absence of the effect","core_discovery":"The central claim is that ordered superstructures, [001] superlattices, a PbS quantum wire in a BaS matrix, and disordered Ba1−xPbxS solid solutions all show a soft transverse-optical phonon whose eigenvector displaces Pb and S atoms out of phase along linear –Pb–S–Pb–S– chains. The soft mode appears only when the structure is stretched by substituting larger Ba atoms; it is not caused by Pb off-centering, since the on-site force constant on Pb is positive. Condensing this mode yields ferroelectric phases with polarization up to roughly 0.15 C/m2, alongside antiferroelectric and mixed phases. The paper concludes that the near-degeneracy of these phases produces a multi-minimum potential with","pith_inferences":["As an extension of the chain mechanism, other rocksalt systems pairing a polarizable heavy cation with a larger cation that stretches the lattice could show the same instability; strontium or calcium would not, as the paper itself notes.","If the near-degenerate FE/AFE energies are real, compressive or tensile epitaxial strain should be able to flip the ground state from ferroelectric to antiferroelectric; this is a tunability lever the paper does not explicitly develop.","The predicted multi-minimum potential implies a distinctive experimental signature: frequency-dependent dielectric response and slow, history-dependent polarization relaxation at low temperatures, behavior typical of relaxor or dipole-glass states.","The quantum-wire result suggests one-dimensional ferroelectric elements can be embedded in a non-ferroelectric rock-salt matrix, which, if grown, would provide a direct nanoscale test of the chain mechanism."],"forward_implications":["Rock-salt chalcogenides can host ferroelectricity, expanding the known structural families of ferroelectrics beyond perovskites and IV–VI compounds.","In BaS–PbS, the ground state is frequently antiferroelectric or mixed FE+AFE; an applied electric field may switch between these near-degenerate states, producing behavior like the irreversible AFE-to-FE transition seen in NaNbO3.","The roughly 29 meV per molecule ferroelectric ordering energy in Ba3PbS4 suggests a Curie temperature on the order of 300 K, making room-temperature realizations plausible.","The same chain mechanism carries over to BaX–PbX with X = Se, Te and to quantum wires, but not to SrS–PbS or CaS–PbS, so the effect is tied to the larger Ba cation's tensile strain.","The near-degenerate infinite set of polarization patterns implies a multi-minimum potential, experimentally observable as nonergodic, history-dependent dielectric response at low temperatures."],"fun_headline_variants":["Barium stretch creates ferroelectricity in rock-salt BaS-PbS","New ferroelectric class predicted: strained rock-salt BaS-PbS","Rocksalt BaS-PbS hosts competing ferroelectric and antiferroelectric states","Strain turns BaS-PbS rock salt into switchable ferroelectric","First-principles reveal ferroelectric order in BaS-PbS beyond perovskites"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The prediction rests on the approximate density-functional phonon description being accurate at the true lattice constant: calculated lattice parameters are about 0.6% smaller than experiment, and the soft-mode frequency depends strongly on strain, so a modest error could erase the instability.","fun_headline_variants_meta":{"raw":{"variants":["Barium stretch creates ferroelectricity in rock-salt BaS-PbS","New ferroelectric class predicted: strained rock-salt BaS-PbS","Rocksalt BaS-PbS hosts competing ferroelectric and antiferroelectric states","Strain turns BaS-PbS rock salt into switchable ferroelectric","First-principles reveal ferroelectric order in BaS-PbS beyond perovskites"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000248,"raw_usage":{"total_tokens":1390,"prompt_tokens":759,"completion_tokens":631,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":503,"completion_tokens_details":{"reasoning_tokens":520}},"tokens_in":503,"tokens_out":631,"duration_ms":6298,"temperature":1.0,"reasoning_tokens":520,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T01:58:16.802427+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Synthesize an ordered Ba3PbS4 sample and measure its low-temperature structure and dielectric response; if no polar or antipolar distortion appears near the estimated ~300 K transition and no polarization develops, the central claim fails. A simpler check is computational: recompute the phonon spectrum at the experimental lattice constant; if the zone-center TO mode is stable there, the predicted ferroelectricity is an artifact of the lattice underestimate.","supporting_citations":[],"review_version":1}