{"id":"29dc10ac-c5e1-4ac6-9bf8-91ec2fb5ee42","arxiv_id":"1908.04581","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"First-principles calculations show CdSe nanoplatelets have a significantly larger negative in-plane thermal expansion than bulk CdSe, driven by flexural ZA and TA-derived E modes.","lead":"This paper computes how CdSe nanoplatelets change size with temperature and finds they shrink when heated, an effect much stronger than in bulk CdSe. The study identifies the specific atomic vibrations that cause this negative thermal expansion, which matters for designing optoelectronic devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative QHA result for the thinnest platelets is unvalidated for flexural-dominated 2D systems; the 2 ML all-temperature negative CTE may depend on missing anharmonic renormalization of the ZA mode.","rationale":"Agree with the reader's weakest_assumption. My stress-test found no internal inconsistency or parameter fitting; the QHA implementation is standard and the bulk benchmark is a real check. However, the strongest claim has two parts: (i) the NTE magnitude enhancement over bulk, and (ii) the 2 ML platelet being NTE over the whole 5-1000 K range. Both are quantitative outputs of QHA. In 2D systems with flexural modes, QHA is known to be less reliable than in 3D because the harmonic approximation at each strain neglects the large anharmonic renormalization of ZA modes; this is not a marginal correction for 2 ML. The paper's own statement in Section II acknowledges that rigorous anharmonicity needs MD. Since no such benchmark is provided, the quantitative claim is not yet secured. This does not refute the paper; it means acceptance should be conditional on one explicit MD/TDEP check. The concrete test above isolates the claim and would settle whether the flexural-mode anharmonicity changes the central conclusion.","tokens_in":12561,"tokens_out":5195,"duration_ms":59829,"concrete_test":"Using the same LDA pseudopotentials and supercell geometry, perform finite-temperature ab initio molecular dynamics (or, equivalently, a self-consistent phonon/TDEP calculation) for the 2 ML and 3 ML F-terminated nanoplatelets at 100, 300, and 600 K, and extract the in-plane CTE from the temperature dependence of the averaged in-plane lattice parameter (or from the stress-fluctuation formula). If the MD/TDEP 2 ML CTE is not negative across this range, or if its magnitude at 300 K is less than roughly half the QHA value in Fig. 4, then the central claim overstates the effect.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative claim — that the in-plane CTE of CdSe nanoplatelets 'appreciably exceeds' the bulk value and that the 2 ML platelet remains negatively expanding up to 1000 K (Fig. 4) — is computed entirely in the quasiharmonic approximation (Section II, Eqs. (1)-(3), Eq. (5)). QHA treats each phonon as independent and harmonic at every strain state, encoding anharmonicity only through the strain dependence of harmonic frequencies. For free-standing quasi-2D systems this is precisely the weak point: the flexural ZA mode, which the paper identifies as a main contributor, has a quadratic dispersion and large out-of-plane amplitudes that drive strong anharmonic coupling between bending and stretching. In the 2 ML platelet these amplitudes are largest, and the paper's own mechanism ('the out-of-plane thermal motion of very thin nanoplatelets results in shrinking') is an anharmonic effect that QHA does not renormalize. The bulk CdSe validation (Fig. 2) does not cover this regime, and no MD, self-consistent phonon, or experimental check for nanoplatelets is provided. If anharmonic renormalization changes the balance between the negative ZA/TA/E contributions and the positive modes, the magnitude—or even the sign at high temperature—could shift. This uncertainty is load-bearing because the headline result is the enhancement of NTE and the all-temperature negativity of the 2 ML platelet.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents first-principles quasiharmonic (QHA) calculations of the in-plane coefficient of thermal expansion (CTE) for free-standing zinc-blende CdSe nanoplatelets of 2–5 monolayers, using DFT (LDA) with norm-conserving pseudopotentials. It also computes CTE for bulk zinc-blende and wurtzite CdSe and validates the bulk results against experimental data for wurtzite CdSe. The central findings are that nanoplatelets exhibit a significantly larger negative in-plane CTE than bulk CdSe and that the 2 ML platelet remains negatively expanding over the entire 5–1000 K range. The authors attribute the enhanced NTE mainly to the out-of-plane flexural ZA mode and to in-plane optical E modes and surface E modes originating from folding of the bulk TA branch with negative Grüneisen parameters. They further report that in-plane polarized optical modes with negative γ appear in a wide range of two-layer quasi-2D systems, which they argue is a first observation.","tokens_in":12810,"tokens_out":5359,"duration_ms":57337,"significance":"If accepted, the paper provides a parameter-free first-principles prediction of an enhanced negative thermal expansion regime in quasi-2D CdSe, with a clear mode-resolved mechanism. The bulk CTE comparison in Fig. 2 gives confidence in the computational methodology, and the Appendix's survey of twelve two-layer systems makes the claim about negative-γ in-plane optical modes a general and falsifiable statement. The machine-checked details of the phonon interpolation, strain dependence, and effective-stress handling in Appendix B are strengths. The main limitation is that the quantitative nanoplatelet predictions, particularly the all-temperature negative CTE of the 2 ML platelet, rely entirely on QHA without direct validation in the flexural-dominated few-monolayer regime.","major_comments":[{"comment":"The load-bearing quantitative claim—that the in-plane CTE of CdSe nanoplatelets appreciably exceeds the bulk magnitude and that the 2 ML platelet remains negatively expanding up to 1000 K—is computed entirely within the quasiharmonic approximation. The QHA assumption stated in Section II (modes remain independent and harmonic at each strain, with anharmonicity entering only through the strain dependence of frequencies) is precisely the weak point for the flexural ZA mode in ultrathin platelets: large out-of-plane displacements drive strong bending-stretching anharmonic coupling, and the paper's own explanation of the 2 ML effect ('the out-of-plane thermal motion of very thin nanoplatelets results in shrinking') is an anharmonic mechanism that QHA does not renormalize. The bulk validation in Fig. 2 does not cover this regime, and no molecular dynamics, self-consistent phonon, or experimental check is provided for nanoplatelets. I recommend that the authors either perform an explicit anharmonic benchmark (e.g., MD or anharmonic perturbation theory) for the 2 ML and 3 ML platelets, or provide a quantitative estimate of how ZA-mode frequency renormalization would alter the sum of Grüneisen parameters that governs the high-temperature CTE. Without this, the magnitude—and potentially the sign at high temperature—of the headline result remains unverified.","section":"Section V, Fig. 4 and Eq. (5)"},{"comment":"The derivation of the CTE connects it to the sum of mode Grüneisen parameters through A0, which is stated as 1/(3B0) in Eq. (5). For the nanoplatelets, the paper later replaces A0 by S11+S12, but the dimensional and notational connection between Eq. (5) and the biaxial stress-strain treatment is not explicit. Please clarify the exact prefactor used for the quasi-2D case and validate it against the standard formula α = (S11+S12) * (∂S_vib/∂T) / A, so that the reader can verify that the 'three times larger' strain-derivative issue mentioned in Section III was properly divided by 3 before solving for uxx(T). If the factor of 1/3 was not applied correctly, all bulk CTE values would be off by a factor of three.","section":"Section II, Eq. (5) and Section III"}],"minor_comments":[{"comment":"In Eq. (5), the quantity A0 is defined as 1/(3B0) but it appears multiplied by 1/V0; please state explicitly whether V0 is the unit-cell volume and whether the final CTE expression is intended as the linear CTE in the cubic case.","section":"Section IV, Eq. (5)"},{"comment":"For the bulk column, specify that the values correspond to zinc-blende CdSe and describe how the effective-thickness correction described in Appendix B enters the reported compliance values for the nanoplatelets.","section":"Table I"},{"comment":"The text states that modes with negative γ reach energies up to ~140 cm^{-1} and that this explains the wider NTE temperature range, but no figure or table directly shows this energy range; consider adding a panel in Fig. 6 or a small table listing the negative-γ modes and their energies.","section":"Section V, phonon spectra text"},{"comment":"The footnote about LDA for quasi-2D systems is attached to a reference entry; please move it to the main text or rephrase so that the argument for choosing LDA is part of the computational methodology.","section":"Reference [29]"},{"comment":"For the Cl- and Br-terminated 3 ML platelets, only the final CTE curves are shown in Fig. 4 with no comment on the differences in the underlying phonon modes beyond the Γ-point frequencies in Table S3; a brief discussion in the text would help the reader connect the mode changes to the CTE evolution.","section":"Section V, Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is well-executed for a QHA study and the bulk validation is convincing. My major revision request is driven by the specificity of the 2 ML all-temperature NTE prediction, which is exactly the regime where QHA's independent-mode assumption is least reliable. I would not reject the paper, but I would want to see either a new calculation (MD or self-consistent phonon) or a careful quantitative discussion of anharmonic renormalization before accepting the headline numbers. The broader claim about negative-γ in-plane optical modes across two-layer systems is interesting and appears well supported by the survey."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper does something genuinely new: first-principles CTE calculations for CdSe nanoplatelets with bulk validation and a clear mechanistic story. The identification of negative Grüneisen parameters for in-plane optical E modes (arising from folding of the bulk TA phonon) is a real contribution, and the survey of twelve two-layer systems adds useful context. The computational detail is careful—proper stress handling in supercells with vacuum (Appendix B), dense phonon interpolation, and bulk wurtzite results that match experiment. The citation pattern is fine; reliance on the author's earlier mode-projection work is not circular because that prior work is itself parameter-free and first-principles.\n\nThe main soft spot is the one the stress-test note flags: QHA for the thinnest platelets. The 2 ML result—negative CTE over the entire temperature range—depends on the flexural ZA mode with large out-of-plane amplitudes, and QHA does not renormalize that anharmonicity. No MD or self-consistent phonon check is provided. I think this is a real caveat for the quantitative 2 ML curve, especially near 1000 K, where the sign could in principle shift if anharmonic renormalization is strong. But it is not load-bearing for the paper's core contribution. The mechanism—negative gamma for TA-derived E modes—is robust, and the bulk validation gives indirect confidence. The paper also does not oversell the prediction; it states the QHA assumptions up front.\n\nAppendix A's generalization to twelve two-layer systems is a bonus but slightly overreaching: negative gamma at a few high-symmetry points is weaker evidence than the full CdSe phonon analysis, and the table lacks uncertainty estimates. That is minor.\n\nThis paper is for computational materials scientists working on 2D thermal expansion and on CdSe nanostructures. It deserves a serious referee and publication, ideally with a sentence in the conclusions acknowledging the QHA limitation for the thinnest platelets.\n\nSend it to peer review.","headline":"Solid first-principles QHA study that finds a new mechanism for negative in-plane CTE in CdSe nanoplatelets; the 2 ML all-temperature result is the least certain piece but the overall argument holds.","tokens_in":13342,"tokens_out":2735,"would_cite":true,"duration_ms":27989,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["65.40.De","63.22.-m"],"model":"deepseek-v4-flash","headline":"The paper claims that CdSe nanoplatelets two to five monolayers thick develop a strongly negative in-plane coefficient of thermal expansion that exceeds the bulk value and persists over a wider temperature range.","keywords":["negative thermal expansion","CdSe nanoplatelets","quasiharmonic approximation","Grüneisen parameter","flexural ZA mode","folded TA phonon","in-plane optical E modes","first-principles phonons"],"falsifier":"Measure the in-plane thermal expansion of freestanding two-to-five-monolayer zinc-blende CdSe nanoplatelets between 5 and 300 K using electron diffraction or nanomechanical deflection; if the CTE is not negative with magnitude at least as large as the bulk, the central claim fails. Alternatively, run full anharmonic molecular dynamics on the two-monolayer platelet and check whether the predicted whole-range negative expansion survives beyond the quasiharmonic approximation.","tokens_in":12340,"feed_emoji":"📉","tokens_out":10081,"duration_ms":79323,"temperature":0.7,"pith_summary":"This paper establishes that a zinc-blende CdSe crystal cut down to a two-to-five-monolayer platelet acquires a strongly negative in-plane coefficient of thermal expansion: heating makes the platelet contract in-plane, with a magnitude that exceeds the negative expansion of bulk CdSe and over a wider temperature range. The prediction comes from first-principles phonon calculations in the quasiharmonic approximation, so no experimental thermal data are fitted. In the two-monolayer limit the contraction persists over the entire computed temperature range. The result matters because nanoplatelet optoelectronics depends on how the band gap shifts with temperature, and thermal expansion is one of the competing contributions.","feed_headline":"Heat shrinks ultrathin CdSe flakes far more than bulk CdSe","feed_subtitle":"First-principles phonon calculations trace the effect to flexural and folded optical modes that pull atoms inward as bonds stretch.","key_machinery":"The machinery is the quasiharmonic approximation expressed through the strain-dependent phonon free energy $F(V_0,T)=E_{\\mathrm{tot}}(V_0)+F_{\\mathrm{vib}}(V_0,T)$, combined with mode-resolved Grüneisen parameters $\\gamma_{jq}=-d\\ln\\omega_{jq}/d\\ln a$ for the in-plane lattice parameter. The negative sign of $\\gamma$ for selected modes, weighted by the thermal occupation factor $\\hbar\\omega/kT$, is what converts vibrational entropy into contraction. In the nanoplatelets the load-bearing objects are the ZA flexural mode and the folded acoustic-like E modes whose displacement patterns are in-plane but whose restoring forces come from bending of Cd-Se bonds.","core_discovery":"On the paper's own terms, the central discovery is that the large negative in-plane thermal expansion in CdSe nanoplatelets is produced by two cooperating groups of modes: the out-of-plane flexural acoustic ZA mode, common to all quasi-two-dimensional systems, and a set of in-plane-polarized optical E modes plus surface E modes that arise when the transverse-acoustic branch of bulk CdSe folds into the slab Brillouin zone. These folded modes inherit the negative Grüneisen parameter of the bulk TA phonon, and the paper reports that this is the first time negative Grüneisen parameters are found for in-plane-polarized optical modes in a quasi-2D system. The same folding mechanism is then checked in twelve two-layer systems, and in-plane optical modes with negative gamma appear in nearly all of them, suggesting the effect is generic for multilayer quasi-2D materials with strong interlayer interaction.","pith_inferences":["Because the mechanism is the folding of a bulk TA branch with negative $\\gamma$, the same enhanced negative in-plane expansion should occur in other zinc-blende II-VI and III-V nanoplatelets; a first-principles check on CdS, ZnSe, or ZnS would be a direct test.","The whole-range contraction of the two-monolayer platelet implies its in-plane sound velocities should stiffen anomalously with temperature, an effect measurable by Raman or Brillouin scattering.","The absence of pressure-induced softening in a system with negative thermal expansion suggests that, in quasi-2D materials, the two effects are not coupled the way they are in framework crystals such as $\\mathrm{Zn(CN)}_2$.","A full anharmonic treatment for the two-monolayer case, where flexural displacements are largest, would reveal whether the predicted negative expansion survives beyond the quasiharmonic approximation."],"forward_implications":["Freestanding CdSe nanoplatelets with 2-5 monolayers should contract in-plane on heating across a wide temperature range, most strongly for the thinnest platelets.","Thermal expansion contributes no more than about $+1.2\\times10^{-5}$ eV/K to $dE_g/dT$ at 300 K, so the large negative band-gap shift seen in CdSe nanoparticles must come from electron-phonon coupling.","Negative Grüneisen parameters for in-plane-polarized optical E modes should appear in most two-layer quasi-2D systems with strong interlayer interaction, including two-layer graphene, BN, SiC, silicene, germanene, blue and black phosphorene, SnS, and TiO$_2$.","Heavier halide termination (Cl, Br instead of F) reduces the magnitudes of the negative Grüneisen parameters and shifts the CTE curves, giving a chemical handle on the thermal response."],"supporting_citations":[{"why":"Classifies the phonon modes of CdSe nanoplatelets and identifies the folded acoustic-like E modes and quasi-Lamb modes that the thermal-expansion analysis uses.","marker":"Ref. 37"},{"why":"Shows that transverse acoustic phonons in zinc-blende semiconductors can have negative Grüneisen parameters, the mechanism inherited by the folded E modes.","marker":"Ref. 35"},{"why":"Provides the experimental CTE data for hexagonal CdSe used to validate the bulk calculation.","marker":"Ref. 6"},{"why":"Determines the bridge position of terminating halogen atoms and the structural model for the nanoplatelet supercells.","marker":"Ref. 5"},{"why":"First-principles quasiharmonic calculation attributing graphene's negative thermal expansion to the flexural ZA mode.","marker":"Ref. 14"},{"why":"Reports negative Grüneisen parameters for in-plane acoustic modes in MoTe2 and WTe2 monolayers, the nearest in-plane analogue.","marker":"Ref. 19"},{"why":"Documents the layer-breathing ZO' mode with negative gamma in multilayer graphene, a comparison point for out-of-plane optical modes.","marker":"Ref. 38"}],"fun_headline_variants":["CdSe nanoplatelets shrink on heating more than bulk CdSe","Folded phonons cause large negative expansion in CdSe flakes","Heat makes CdSe nanoplatelets contract more than bulk","CdSe flakes shrink when heated, folded phonons explain why"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that each phonon mode remains harmonic at every strain state, so all anharmonicity is captured by the strain dependence of mode frequencies; for two-monolayer platelets with large out-of-plane flexural displacements this could miss genuine anharmonic renormalization.","fun_headline_variants_meta":{"raw":{"variants":["CdSe nanoplatelets shrink on heating more than bulk CdSe","Folded phonons cause large negative expansion in CdSe flakes","Heat makes CdSe nanoplatelets contract more than bulk","CdSe flakes shrink when heated, folded phonons explain why"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000867,"raw_usage":{"total_tokens":3703,"prompt_tokens":836,"completion_tokens":2867,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":452,"completion_tokens_details":{"reasoning_tokens":2796}},"tokens_in":452,"tokens_out":2867,"duration_ms":20195,"temperature":1.0,"reasoning_tokens":2796,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:37:26.616648+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the in-plane thermal expansion of freestanding two-to-five-monolayer zinc-blende CdSe nanoplatelets between 5 and 300 K using electron diffraction or nanomechanical deflection; if the CTE is not negative with magnitude at least as large as the bulk, the central claim fails. Alternatively, run full anharmonic molecular dynamics on the two-monolayer platelet and check whether the predicted whole-range negative expansion survives beyond the quasiharmonic approximation.","supporting_citations":[],"review_version":1}