{"id":"a376386c-c998-41a7-976f-81a6e6e93aff","arxiv_id":"2608.07884","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The transformation barrier and critical thickness for (110)-to-(100) transitions in 2D noble metal nanosheets are described by a Bell-Evans-Polanyi linear relation.","lead":"Using density functional theory, this paper predicts that very thin (110) sheets of some noble metals spontaneously convert to (100) sheets below a critical thickness. It shows these conversion barriers follow the Bell-Evans-Polanyi relationship, giving a simple rule for when high-energy 2D metal facets become unstable.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 6 predicts Rh critical thickness 2.4 ML, but the same paper reports Rh (110) stable at 2 ML; the BEP critical-thickness claim is internally contradicted.","rationale":"The reader's weakest_assumption correctly identifies the straight-line approximation and constant Delta_ES as the fragile foundation of Eq. 5. My concern is a direct, checkable consequence of that assumption: the derived critical thickness for Rh is 2.37 ML, but the paper's own DFT relaxation and barrier calculations show 2 ML Rh (110) is stable with a finite barrier. This is not an untested extrapolation; it is an internal inconsistency between Eq. 6 and the reported data. If the 2 ML Rh barrier is finite, the BEP model is not 'dominant' for Rh, one of the metals in the abstract and title. If the barrier is absent, the stability claim in Fig. 1(c) is wrong. Either way the central claim needs substantial revision. The Pd and Ag comparisons are suggestive, and the model is transparent, but a paper whose headline generalization is contradicted by its own Table 1 parameters should not be accepted in current form. I agree with the reader that the linearization assumption is the root cause, but the concrete failure at Rh 2 ML is a more decisive way to state it. Verdict: reject current version, with the possibility that a revised manuscript restricted to Pd/Ag and excluding Rh/Ir could be reconsidered.","tokens_in":5289,"tokens_out":11612,"duration_ms":133680,"concrete_test":"Compute the full DFT energy profile E(b) for a freestanding 2 ML Rh (110) sheet with the same VASP settings and locate the maximum along the b-compression path. Report the actual E_b at 2 ML. If the maximum is finite, Eq. 6's predicted n*=2.4 ML is refuted for Rh and the BEP critical-thickness claim must be restricted to Pd/Ag; if the maximum is absent, the paper's assertion that Rh (110) is stable at 2 ML is wrong.","verdict_should_be":"REJECT","load_bearing_attack":"The most load-bearing concern is not merely that the energy curves are approximated as straight lines; it is that this approximation produces a quantitative prediction that contradicts the paper's own DFT results for Rh. Using Table 1 values (m1=-0.427, m2=0.414, Eb=0.234 eV/atom, Delta_ES=0.546 eV/surface atom), Eq. 5 gives E_b'(n)=0.234 - [2m1/(m1-m2)]*Delta_ES/n = 0.234 - 1.0155*0.546/n. At n=2 ML this is -0.043 eV/atom, and Eq. 6 gives n*=2.37 ML. The model therefore says a 2 ML Rh (110) sheet is below the critical thickness and should transform spontaneously. The text, however, states that Rh exhibits a finite barrier down to 2 ML and that Rh nanosheets of all thicknesses retain the (110) facet after relaxation. Since 2 ML < 2.37 ML, this is a direct counterexample to the central claim that the critical thicknesses of these metals are described by the BEP principle. Pd and Ag show rough agreement, but Rh is one of the four metals named in the abstract, so the claim as stated fails.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses DFT to study the thickness-dependent transformation of freestanding (110)-faceted noble-metal nanosheets (Rh, Pd, Ag, Ir, Pt, Au) into (100)-faceted sheets. It proposes a Bell-Evans-Polanyi (BEP)-type relation, Eq. (5), in which the nanosheet transformation barrier decreases linearly with the surface-energy difference ΔE_S divided by thickness, and a critical thickness n*, Eq. (6), below which the (110) sheet is predicted to transform spontaneously. The authors claim that this relation is confirmed by their calculated barriers for Rh and Pd and that the critical thicknesses of Rh, Pd, and Ag are described by the BEP principle.","tokens_in":5572,"tokens_out":7312,"duration_ms":73127,"significance":"If correct, the proposed BEP relation would provide a simple, physically transparent criterion for the metastability limit of high-energy-facet 2D metals, and the paper contributes a useful DFT dataset of nanosheet transformation barriers. The explicit formula in Eq. (5) and the approximate agreement for Pd (7.8 vs 7 ML) and Ag (6.5 vs 5 ML) are of interest. However, the central claim is currently overbroad and is contradicted by the paper's own Rh data, which show a finite barrier at 2 ML despite a predicted critical thickness of 2.4 ML. The paper's value would be substantially strengthened if the claims were restricted to the metals for which the relation is actually supported.","major_comments":[{"comment":"Equation (6) with the Table 1 parameters for Rh (m1=-0.427, m2=0.414, E_b=0.234 eV/atom, ΔE_S=0.546 eV/surface atom) gives n*=2.37 ML. Substituting n=2 ML into Eq. (5) yields E_b'=-0.043 eV/atom, i.e., the model predicts that a 2 ML Rh (110) sheet is below the critical thickness and should transform spontaneously. This directly contradicts the DFT result reported in Fig. 1(c) and the text, which state that Rh (110) sheets remain stable down to 2 ML with a finite barrier. Because Rh is one of the four metals named in the abstract, this internal contradiction invalidates the abstract's claim that the critical thicknesses of these metals are described by the BEP principle.","section":"Eq. (6) and Table 1"},{"comment":"The abstract states that 'the critical thicknesses for these metals are also described by the BEP principle,' referring to Rh, Pd, Ag, and Ir, but the conclusion restricts this statement to Pd and Ag ('the critical thicknesses of Pd and Ag with vanishing energy barriers'). No critical thickness or BEP prediction is reported for Ir, and for Rh the prediction contradicts the DFT result as noted above. Additionally, for Ag the predicted n*=6.5 ML differs from the observed 5 ML by 30%, so even for the favorable metals the agreement is only qualitative. The central claim needs to be restricted to the metals for which it holds, and the quantitative discrepancies need to be addressed.","section":"Abstract and Conclusion"},{"comment":"The derivation of Eq. (5) is built on two assumptions: the bulk energy curves are linear (Eqs. (1)-(2), Fig. 3(e)) and the surface-energy difference ΔE_S is constant along the transformation coordinate. Fig. 2(b)-(e) show that these functions are curved and that the surface-energy difference is b-dependent, and the authors themselves note a local minimum in the Pt bulk curve (Fig. 2(c)). The paper does not quantify how well the straight-line approximation represents the actual curves, and the visual agreement in Fig. 3(g)-(h) is not a substitute for a fit-quality metric. The validity range of the BEP approximation should be stated explicitly.","section":"Eqs. (1)-(5) and Fig. 2"}],"minor_comments":[{"comment":"The text 'In Fig. 3 (g) and (h) we plot the calculated energy barriers (Fig. 2 d-f) versus the derived ones' cross-references Fig. 2(d)-(f), which are surface-energy curves, not barriers; this should be corrected to reference Fig. 3(c)-(d). The caption labels (c) and (d) also appear to be swapped relative to the text's description of Fig. 3(e) and (f).","section":"Fig. 3 caption"},{"comment":"The units of the term 2ΔE_S/n should be spelled out: ΔE_S is per surface atom, and since each sheet has two surfaces and n atomic layers, the factor 2/n converts the energy to eV/atom. A brief explanation of this conversion would improve the derivation.","section":"Derivation of Eq. (4)"},{"comment":"The text contains 'This in consistent with a previous DFT study' which should read 'This is consistent with a previous DFT study', and the conclusion uses 'dominate role' where 'dominant role' is intended.","section":"Typographical errors"},{"comment":"Table 1 would be more complete if it stated the units of m1, m2, b1, and b2, and if it included parameters for Ir or an explanation of their absence, given that Ir is named in the abstract.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's abstract and title overstate the dominance of the BEP principle, and the Rh result is an internally visible counterexample to the critical-thickness claim. The authors should be asked to either reconcile the Rh discrepancy (for example, by showing that the BEP prediction at n=2 is within numerical error or by redefining n*) or to revise the abstract and conclusion so that the critical-thickness claim applies only to Pd and Ag. The paper otherwise fits the journal scope, and the DFT barrier dataset and the explicit formula are potentially valuable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper does something genuinely new—it adapts the Bell–Evans–Polanyi idea to thickness-dependent barriers in freestanding (110) metal nanosheets and writes down an explicit formula (Eq. 5) plus a critical-thickness expression (Eq. 6). The comparison with direct DFT barriers for Pd and Ag is plausible, and the notion of a BEP-dominated stability limit is worth testing. But the central claim is undermined by an internal contradiction with the paper's own data, and the analysis misses one of the four metals named in the abstract.\n\nWhat is actually good: the derivation is transparent, the parameters are tabulated, and the predicted thickness scaling is a clean testable hypothesis. For Pd and Ag, the derived critical thicknesses (7.8 and 6.5 ML) bracket the DFT-optimized values (7 and 5 ML) reasonably. The paper also correctly notes that Pt has a local minimum in the bulk energy curve, which breaks the straight-line assumption; that is a fair caveat.\n\nWhere it falls apart: the stress-test is on target. Using Table 1 for Rh (m1 = -0.427, m2 = 0.414, Eb = 0.234 eV/atom, ΔES = 0.546 eV/surface atom), Eq. 6 gives n* = 2.37 ML. The paper explicitly states that Rh (110) sheets retain their facet down to 2 ML after relaxation. So the model predicts a 2 ML Rh sheet should transform spontaneously, but the DFT relaxation says it does not. That is not a minor discrepancy—it is a direct counterexample to the claim that critical thicknesses are described by the BEP principle. The paper mentions this mismatch but then concludes the principle dominates anyway, which does not follow. Also, the abstract includes Ir, but Table 1 and the derived-thickness discussion cover only Rh, Pd, and Ag. Ir's absence is a gap for a paper claiming four metals.\n\nThe straight-line energy approximation is a further limitation, but the Rh contradiction is the load-bearing issue. A referee could reasonably ask whether the BEP fit is just an interpolation of the thickness range where the barrier is positive, with no predictive power at the point where it should vanish.\n\nWho this is for: researchers working on metastable 2D metals or facet-dependent stability, especially those interested in cheap estimates of synthesis limits. The idea is worth knowing, but the current version overclaims. It deserves a serious referee, but not acceptance without major revision and a direct reconciliation of the Rh result.\n\nRecommendation: send to peer review, but make clear that the Rh contradiction must be addressed, and Ir should be either included or removed from the abstract.","headline":"A clean two-line formula for thickness-dependent transformation barriers in 2D noble metals, but the paper's own Rh results contradict the predicted critical thickness.","tokens_in":6107,"tokens_out":2030,"would_cite":false,"duration_ms":24618,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The Bell-Evans-Polanyi principle, expressed as a straight-line relation between transformation barrier and surface-energy gain, fixes the critical thickness below which ultrathin (110) noble-metal nanosheets spontaneously transform to…","keywords":["two-dimensional noble metals","high-energy facets","Bell-Evans-Polanyi principle","structural transformation","critical thickness","density functional theory","nanosheet stability","metastability"],"falsifier":"Compute the full energy profile of a 6-monolayer Pd (110) sheet along the b lattice parameter with the same DFT settings: the paper predicts a vanishing barrier, so finding an energy barrier larger than the thermal energy at synthesis conditions (or a barrier that does not go to zero at the predicted 7–8 ML critical thickness) would falsify the BEP-dominated description.","tokens_in":5067,"feed_emoji":"📏","tokens_out":7610,"duration_ms":74290,"temperature":0.7,"pith_summary":"Using density functional theory, this paper predicts that freestanding (110) nanosheets of Rh, Pd, Ag, and Ir transform into (100) nanosheets through a b-direction lattice contraction, and that the energy barrier for this transformation falls linearly as the sheet gets thinner. The paper identifies this as a Bell–Evans–Polanyi (BEP) relationship: the barrier plays the role of activation energy, and the surface-energy difference between the (110) and (100) facets plays the role of reaction enthalpy. From that analogy it derives a formula for the critical thickness at which the barrier vanishes, giving approximately 7–8 monolayers for Pd and about 5–6 monolayers for Ag, while Rh keeps a finite barrier down to 2 monolayers. If correct, the result turns the synthesis limit of high-energy-facet 2D metals into a quantity that can be estimated from bulk energy curves and facet surface energies alone.","feed_headline":"Straight-line rule predicts collapse of 2D noble metals","feed_subtitle":"Below 7 (Pd) and 5 (Ag) atomic layers, (110) sheets flip to (100), matching the BEP collapse formula.","key_machinery":"The load-bearing construction is the straight-line approximation of the bulk energy curve along the lattice parameter $b$. During the (110)-to-(100) transformation, $b$ shrinks from $\\sqrt{2}a$ to $a$ while the other in-plane lattice parameter stays fixed. The paper fits the energy of the bulk region to two lines, $y_1=m_1 x+b_1$ and $y_2=m_2 x+b_2$, whose intersection gives the bulk barrier $E_b$. For a sheet of $n$ monolayers, the product-side line is lowered by $2\\Delta E_S/n$ because the two surfaces gain $\\Delta E_S$ per surface atom; intersecting the shifted line with the reactant line gives Eq. (5), and setting the barrier to zero gives the critical-thickness formula. Everything follows from this BEP ansatz.","core_discovery":"The central claim is that the metastability limit of pristine (110)-faceted noble-metal nanosheets is governed by the Bell–Evans–Polanyi principle, not by a more complex thickness-dependent mechanism. Concretely, the DFT-computed transformation barrier $E_b'$ follows $E_b' = E_b - \\frac{2 m_1}{m_1-m_2}\\frac{\\Delta E_S}{n}$, where $E_b$ is the bulk transformation barrier, $m_1$ and $m_2$ are the slopes of the bulk energy versus lattice parameter $b$ on the reactant and product sides, $\\Delta E_S$ is the per-surface-atom surface-energy difference between the (110) and (100) facets, and $n$ is the sheet thickness in monolayers. The paper shows that for Rh, Pd, and Ag the barriers computed from full DFT relaxation agree with this linear expression, and that setting $E_b'=0$ yields critical thicknesses ($n^* = \\frac{2 m_1}{m_1-m_2}\\frac{\\Delta E_S}{E_b}$) close to the values found by direct geometry optimization: 7.8 versus 7 ML for Pd, and 6.5 versus 5 ML for Ag. The authors therefore extend the BEP principle from reaction kinetics to the structural stability of low-dimensional materials.","pith_inferences":["One testable extension is that adsorbates or substrates, by changing the effective $\\Delta E_S$, should shift the critical thickness linearly; a monolayer that lowers the (110) surface energy could stabilize sheets that are otherwise too thin to exist.","The same BEP derivation should apply to other facet pairs, such as (001)-to-(111) transitions in fcc metals, provided the bulk energy along the relevant lattice coordinate is approximately two-line; a systematic failure there would map where the straight-line assumption breaks.","If the relation holds beyond the four metals studied, the critical thickness of any fcc metal nanosheet can be screened from bulk equation-of-state data and ideal-facet surface energies, making the synthesis limit a bulk property plus a surface term.","Platinum and gold are the paper's own exception: their bulk energy curves show a local minimum or a surface-energy maximum, so the two-line model would predict a non-monotonic barrier or a residual barrier even at monolayer thickness; that is a place where the BEP description should be tested rather than assumed."],"forward_implications":["For palladium, (110) nanosheets thinner than roughly 7 monolayers are not metastable: they relax spontaneously to (100) facets, so synthesis below that thickness should not yield the high-energy facet.","For silver the corresponding limit is near 5 monolayers, while rhodium sheets preserve the (110) facet down to at least 2 monolayers because of its larger bulk barrier.","The barrier of a (110) sheet decreases linearly with inverse thickness, which means the dominant destabilizing effect is the growing share of the surface-energy gain, not a change in the transformation mechanism.","For sufficiently thick sheets the barrier converges to the bulk value, so the same formula can be used to decide when surface effects can be ignored in computational studies.","The BEP relation provides a shortcut: critical thickness can be predicted from the bulk energy profile and the (110)/(100) surface-energy difference without running a full thickness-by-thickness DFT scan."],"supporting_citations":[{"why":"Supplies the density functional theory methods used for all total-energy, barrier, and surface-energy calculations.","marker":"[8-9]"},{"why":"The earlier molecular-dynamics estimate of the (001)-to-(111) barrier in ultrathin Au films, which the paper's thickness-barrier study extends and puts on a quantitative footing.","marker":"[6]"},{"why":"Establishes the broken-bond ordering of fcc surface energies (110) > (100) > (111) that motivates the choice of high-energy (110) facets.","marker":"[15]"},{"why":"The DFT study that found a stable body-centered tetragonal phase for Pt, which the paper invokes to explain the local energy minimum in Pt's bulk energy curve.","marker":"[18]"}],"fun_headline_variants":["BEP rule sets collapse limit for 2D noble metals","BEP law governs flip of (110) to (100) in 2D metals","Simple linear rule predicts 2D metal collapse","BEP principle dictates thickness limit of 2D metals","Thin noble metal sheets flip facet via BEP rule"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The formula assumes that the bulk energy along the contraction path is two straight lines with constant slopes and that the surface-energy gain between the (110) and (100) facets stays constant for every thickness; platinum, with its local energy minimum, already shows where this linearity can break.","fun_headline_variants_meta":{"raw":{"variants":["BEP rule sets collapse limit for 2D noble metals","BEP law governs flip of (110) to (100) in 2D metals","Simple linear rule predicts 2D metal collapse","BEP principle dictates thickness limit of 2D metals","Thin noble metal sheets flip facet via BEP rule"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000961,"raw_usage":{"total_tokens":4091,"prompt_tokens":939,"completion_tokens":3152,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":3079}},"tokens_in":555,"tokens_out":3152,"duration_ms":21937,"temperature":1.0,"reasoning_tokens":3079,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:43:21.022841+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full energy profile of a 6-monolayer Pd (110) sheet along the b lattice parameter with the same DFT settings: the paper predicts a vanishing barrier, so finding an energy barrier larger than the thermal energy at synthesis conditions (or a barrier that does not go to zero at the predicted 7–8 ML critical thickness) would falsify the BEP-dominated description.","supporting_citations":[{"cited_title":"Physical Review Letters 2002, 88 (9), 096103","cited_arxiv_id":null,"evidence_quote":"The earlier molecular-dynamics estimate of the (001)-to-(111) barrier in ultrathin Au films, which the paper's thickness-barrier study extends and puts on a quantitative footing."},{"cited_title":"H., Broken-bond rule for the surface energies of noble metals","cited_arxiv_id":null,"evidence_quote":"Establishes the broken-bond ordering of fcc surface energies (110) > (100) > (111) that motivates the choice of high-energy (110) facets."},{"cited_title":"I.; Gall, K., Density functional theory investigation of surface -stress-induced phase transformations in fcc metal nanowires","cited_arxiv_id":null,"evidence_quote":"The DFT study that found a stable body-centered tetragonal phase for Pt, which the paper invokes to explain the local energy minimum in Pt's bulk energy curve."}],"review_version":1}