{"id":"6a4aabc0-44b4-433d-ae27-bf4a200e307e","arxiv_id":"1908.05342","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Molecular dynamics simulations show that epitaxial stress from a coherent NiAl second phase changes the martensitic transformation in Ni63Al37 from first-order to continuous, critical-like behavior.","lead":"Simulations show that layering a shape-memory alloy with a non-transforming metal can switch its sudden, jerky phase change into a smooth, continuous one. This suggests a path to hysteresis-free shape-memory materials for actuators and sensors.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Continuous order parameter may reflect disorder-broadened first-order behavior, not a true second-order transition; the paper's evidence does not rule this out.","rationale":"The reader identified coherency as the weakest assumption, and coherency is indeed necessary for the proposed mechanism. I agree with that concern but see a more fundamental uncertainty inside the simulation evidence: the observable used to claim second-order behavior, continuous order parameter plus zero hysteresis, is also the signature of a first-order transition smeared by quenched disorder in a constrained random alloy. This is especially relevant because the transforming layer is a random 63 at.% Ni alloy only a few nanometers thick at 65% NiAl, so local composition and stress variations are unavoidable. The paper's free-energy-landscape argument (Fig. 1C) is the strongest independent support, but it rests on integrating stress-strain curves (Eq. 1), a procedure that is only valid for reversible paths and is not validated for these transforming systems, making it secondary support. The scaling analyses are honestly reported as approximate, but they do not discriminate: the beta range 0.45-0.7 spans multiple universality classes, and the system-size comparison is presented as 'identical' rather than used for finite-size scaling. A targeted specific-heat or finite-size-scaling measurement would settle the ambiguity. I therefore do not move the verdict; the paper should remain conditional pending that check, but the reason is the untested criticality hypothesis rather than only the realizability of coherent interfaces.","tokens_in":10361,"tokens_out":6701,"duration_ms":78962,"concrete_test":"Recompute the enthalpy and specific heat for the 65 at.% NiAl nanolaminate at the three system sizes in Supplemental Table 1, and analyze the specific-heat peak height and position under finite-size scaling. A second-order transition would show a peak height increasing with L^{alpha/nu} and a pseudocritical temperature shifting as L^{-1/nu}; a disorder-broadened first-order transition would show a saturated peak and no critical shift. As a supporting check, repeat cooling at two rates separated by an order of magnitude: if the apparent continuity and hysteresis depend on rate, kinetic broadening rather than a true change of order is indicated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that internal epitaxial stress changes the transformation order depends on interpreting the continuous strain evolution in Fig. 2 and the near-zero hysteresis in Fig. 3 as thermodynamic criticality. That interpretation is not uniquely supported. The transforming Ni63Al37 layer is a random alloy confined between NiAl layers; in such a system, quenched disorder and strain gradients can sequentially trigger small martensite variants at slightly different local transformation temperatures, producing a continuous average order parameter and suppressed hysteresis while each local event remains first-order. This is a disorder-broadening mechanism, not a Landau-barrier-free second-order transition. The paper's checks do not separate the two: beta is quoted as 0.45-0.7 (Sec. IV C), a range too wide to identify a universality class; gamma is obtained from fitting only a few volume fractions (Supplemental Fig. 6); and the larger-system comparison (5M and 20M atoms) is only described as 'nearly identical' rather than subjected to finite-size scaling. For a genuine critical point, the specific-heat peak should grow and its position should shift with system size; for a disorder-broadened first-order transition, the continuous-looking response is essentially size-independent. Because no latent-heat, susceptibility, or finite-size-scaling analysis is reported, the observation most directly supporting the order-change claim is ambiguous.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript uses molecular dynamics simulations of nanolaminates composed of B2 NiAl and off-stoichiometric Ni63Al37 to argue that epitaxial stress from a coherent second phase can change the martensitic transformation in Ni63Al37 from first-order to second-order-like. The authors compute free energy landscapes along strain paths, cool and heat the laminates at high NiAl volume fractions, and report a continuous transformation strain, near-zero thermal hysteresis, power-law martensite cluster-size distributions, and critical exponents beta around 0.45 to 0.7, tau around 1.3 to 1.6, and gamma near 0.85. They interpret these signatures as evidence for a critical point and attribute the behavior to free energy landscape engineering via coherency constraints. The paper also reports stiffness softening near the transformation temperature as a practical consequence.","tokens_in":10570,"tokens_out":9251,"duration_ms":94182,"significance":"If correct, the central claim is significant: it would be one of the few demonstrations that martensitic transformation order can be tuned by internal material design rather than by an external field, with implications for low-hysteresis shape-memory alloys and actuators. The paper has several strengths: direct MD evidence of continuous strain evolution, explicit contrast with bulk Ni63Al37, a physically motivated free-energy-landscape explanation, and a set of supplemental analyses including cluster scaling, system-size comparisons, and stiffness response. The authors are also appropriately cautious in noting that the universality class remains to be established. However, the evidence for true thermodynamic criticality is not yet conclusive; the distinction from disorder-broadened first-order behavior is the central correctness risk of the paper.","major_comments":[{"comment":"The continuously varying order parameter and near-zero hysteresis shown in Figs. 2 and 3 do not uniquely establish a thermodynamic second-order transition: in a random alloy with strain gradients, local first-order regions can transform sequentially at slightly different temperatures, producing a continuous average order parameter and suppressed hysteresis while each local event remains first-order. The manuscript's checks do not separate these scenarios. In particular, beta is quoted as 0.45 to 0.7, a range spanning several universality classes; gamma is fit from only a few volume fractions (Supplemental Fig. 6); and the 5M/20M-atom comparison is described only as nearly identical without a finite-size-scaling analysis of a susceptibility or specific-heat peak. A disorder-broadened transition would show essentially size-independent continuous response, whereas a true critical point would show a specific-heat or strain-susceptibility peak that grows and shifts with system size. I recommend adding a quantitative finite-size-scaling analysis of the strain susceptibility or excess enthalpy, and if possible a latent-heat estimate, and tempering the second-order claim accordingly.","section":"§IV A/C, Figs. 2 and 3, Eq. (2)"},{"comment":"The exponent analysis is too fragile to carry the universality-class conclusion. The reported beta range of 0.45 to 0.7 includes mean-field, 3D Heisenberg, and tricritical values, so it cannot distinguish mechanisms. For gamma, the fit in Supplemental Fig. 6 uses volume fractions from 60 to 75 at. % NiAl, yet the transformation character is discontinuous at 60% and continuous only at higher fractions; including non-critical compositions in a scaling fit can bias gamma. No confidence intervals or fit residuals are reported. The paper's own sentence in §IV C that additional work will be required to definitively assess the universality class is appropriate, but the abstract and conclusion go on to claim scaling-law behavior consistent with the mean-field universality class without the same caveat. Please either provide a more robust exponent analysis, including joint Ms and beta fits and sensitivity to the composition range, or weaken the universality-class statements to clearly indicate that the exponents are illustrative.","section":"§IV C, Supplemental Fig. 6"},{"comment":"The proposed mechanism rests on the assumption that the NiAl/Ni63Al37 interfaces remain coherent and defect-free throughout the transformation, so that the two phases share the same in-plane lattice parameters. The simulations enforce this condition via periodic boundary conditions and do not test its stability against dislocation nucleation, interface sliding, or decohesion under the transformation-induced stresses. If coherency is lost, the composite free energy becomes a weighted average of two first-order landscapes and the barrier reappears. The paper should explicitly state this idealization as a limitation, and ideally test coherency stability through larger or free-standing cells or seeded defects, or restrict the conclusions to coherent nanolaminates under ideal interface conditions. Without such a test, the engineering-level claim that epitaxial stress can change the order of the transformation is not fully established.","section":"§III, Fig. 1D, and Methods"}],"minor_comments":[{"comment":"The strain rate and cooling/heating rate values appear corrupted in the text as 1∙10! ps-1 and 1∙10%% K/s; please replace them with the correct numeric values and state whether the cooling rate is the same for all systems.","section":"§III and §IV"},{"comment":"The differential form includes all six stress and strain components, but the free energy landscapes in Fig. 1 appear to be computed under a uniaxial or biaxial strain path; please specify which components are integrated and how shear components are constrained, and note that integrating stress-strain at a finite rate gives a quasi-static path on a constrained free-energy surface rather than necessarily the equilibrium Helmholtz free energy.","section":"Eq. (1)"},{"comment":"The text N ∝ V(, has a corrupted exponent and should read N ∝ V^{-tau}; also define how the cluster volume V is measured and whether the cluster-size histogram is logarithmically binned.","section":"Eq. (3) and cluster analysis"},{"comment":"The polyhedral template matching RMS cutoff of 0.12 is a free parameter with no sensitivity test; because the cluster size distribution and fractal dimension depend on structural classification, please report at least one alternative cutoff and confirm that the reported exponents are stable.","section":"§II, PTM analysis"},{"comment":"The statement near zero hysteresis should be quantified; define the hysteresis measure, for example the difference between forward and reverse transformation temperatures, and report it for each composition with error bars.","section":"Fig. 3 and §IV A"},{"comment":"The conclusion that the transformation can be described via scaling laws should be qualified by the caveat in §IV C that the universality class is not definitively established; the abstract's claim of power-law exponents comparable to other ferroic transitions should either be supported by confidence intervals or described as a qualitative comparison.","section":"§V, Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The central observation of a continuous, low-hysteresis transformation in the simulated coherent laminates is likely real, but the strongest claim, that this represents a true second-order or critical transformation, is not yet fully supported by the evidence. The disorder-broadening alternative is the main risk, and the requested finite-size and susceptibility analysis would substantially strengthen the paper. The coherency idealization also deserves explicit discussion even if it cannot be fully tested in the present simulations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nThe thing to know: this paper shows, in MD, that putting a coherent NiAl second phase inside Ni63Al37 takes the martensitic transformation from first-order to continuous. The order parameter grows smoothly, hysteresis drops to zero above 62.5% NiAl, and the free energy landscape develops a flat barrierless region. That is a real, new result, and it's the reason to read the paper. The FELE idea is the authors' own, but the step from external-field-tuned criticality (Xiao et al.) to internal, material-level control is not a trivial extension.\n\nWhat it does well: the direct observations in Figs. 2 and 3 are clean, the cluster-scaling analysis (τ and fractal dimension) is a nice independent probe, and the authors try three system sizes. They also state plainly that β is hard to fit and report a wide range instead of hiding it.\n\nWhere it's soft: the quantitative criticality claim. Treating the continuous response as true second-order requires ruling out disorder-broadening of a first-order transition. The paper doesn't quite do that. β = 0.45–0.7 does not pin down a universality class, γ is extracted from a handful of volume fractions, and the three system sizes are only described as \"nearly identical\"—there's no actual finite-size scaling of the order parameter or specific heat. A specific-heat peak that grows and shifts with system size would separate criticality from the quenched-disorder alternative. The free energy landscape in Fig. 1C is the strongest evidence against that alternative, but the landscape calculation itself is path-dependent, and the coherency it relies on is enforced by periodic boundary conditions, never tested against dislocation nucleation or interface sliding. So the mechanism is plausible, and the engineering claim is explicitly contingent on coherent, defect-free interfaces.\n\nAlso: one interatomic potential (Farkas) underlies everything, and no code or data is shared. These are minor in the context of a demonstration paper, but they would be cheap to fix.\n\nBottom line: the qualitative finding is probably right and worth taking seriously. The criticality characterization is not yet definitive. If I were the editor I'd send it to a referee with a request to focus on finite-size scaling and the disorder-broadening confound, and to soften the universality-class language. For a materials science audience this is a thought-provoking paper, well worth discussing in a group. I'd cite it if I worked on martensites. Send it out.","headline":"A solid MD demonstration that internal epitaxial stress can make martensitic transformations continuous, but the criticality evidence doesn't yet rule out disorder-broadened first-order behavior.","tokens_in":11142,"tokens_out":2830,"would_cite":true,"duration_ms":29619,"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":"Epitaxial stress turns a first-order metal transition continuous.","keywords":["martensitic transformation","second-order phase transition","free energy landscape engineering","epitaxial stress","coherent interfaces","Ni-Al alloys","molecular dynamics","critical exponents"],"falsifier":"Cool an actual 65 at. % NiAl NiAl/Ni63Al37 laminate, or run a molecular dynamics cell with interfaces free to decohere or with pre-existing dislocations, and measure the strain-temperature curve: a discontinuous jump in transformation strain, a nonzero cooling-heating hysteresis loop, or a barrier in the computed free energy landscape at 65 at. % NiAl would all show the order has not been changed by the epitaxial stress.","tokens_in":10116,"feed_emoji":"🔁","tokens_out":5151,"duration_ms":51170,"temperature":0.7,"pith_summary":"The paper claims that embedding a non-transforming second phase inside a martensitic alloy, with coherent interfaces that lock the two phases to the same in-plane lattice parameter, reshapes the free energy landscape enough to remove the nucleation barrier. In simulations of Ni63Al37 nickel-aluminum alloy laminated with NiAl, the transformation order parameter grows continuously rather than jumping once the NiAl fraction reaches 65 at. %, and thermal hysteresis drops to near zero. If correct, this gives materials designers a built-in, external-field-free route to second-order, critical martensitic behavior, with scaling laws and continuously variable strain instead of abrupt shape changes. The authors intend the mechanism as general: any alloy with a suitably matched second phase could be tuned this way.","feed_headline":"Epitaxial stress turns abrupt metal transition continuous","feed_subtitle":"Simulations of NiAl/Ni63Al37 show hysteresis vanishes above 62.5% NiAl, pointing to low-fatigue shape-memory alloys.","key_machinery":"The load-bearing object is the coherent nanolaminate: a periodic stack in which a thin layer of ordered NiAl shares a defect-free interface with the transforming Ni63Al37 alloy, forcing both phases to adopt the same in-plane lattice parameter. Because the total free energy is then nearly the weighted sum of the two components as functions of lattice parameter, one component's preference for austenite cancels the other's preference for martensite, flattening the double-well landscape into a single shallow well. This 'free energy landscape engineering' is what converts the abrupt, hysteretic first-order jump into a continuous second-order transition; the scaling laws for strain order parameter, cluster sizes, and susceptibility then follow from that flat landscape.","core_discovery":"The central discovery is that a martensitic transformation, normally universally first order, can be made continuous by internal epitaxial stress. Coherent integration of B2 NiAl (which never wants to be martensitic) with Ni63Al37 (which does) forces the two phases to share in-plane lattice parameters, so the composite free energy becomes a weighted combination of two opposing landscapes; at 65 at. % NiAl this combination is flat and barrierless. The order parameter—transformation strain—then rises continuously with cooling, with near-zero hysteresis for NiAl fractions above 62.5 at. %, martensite domains that interpenetrate with rough fractal interfaces, and cluster-size distributions following a power law. The extracted critical exponents ($\\beta$ between 0.45 and 0.7, $\\tau$ between 1.3 and 1.6, $\\gamma \\approx 0.85$) fall in the range of other critical ferroic transitions and are consistent with mean-field universality.","pith_inferences":["The same landscape-flattening argument should apply to any martensitic pair where one phase is stable and the other is unstable at all temperatures; the natural next test is to scan the lattice-parameter mismatch between layers and map where the barrier vanishes.","A direct experimental test would be thermal cycling of sputter-deposited NiAl/Ni63Al37 superlattices: if hysteresis persists above 62.5 at. % NiAl, coherency was not maintained and the engineering premise of the mechanism failed.","The continuous, barrierless transformation implies the transition temperature itself may become tunable by layer thickness and composition, since the flat landscape is a balance that shifts with volume fraction.","If the mean-field assignment holds, finite-size rounding in small samples should follow Ginzburg-like criteria, with order-parameter fluctuations growing near the critical temperature; this is checkable in the existing simulation trajectories."],"forward_implications":["In Ni-Al, any NiAl fraction above about 62.5 at. % should yield a continuous transformation with negligible thermal hysteresis, while lower fractions remain first order.","The transformation strain can be varied continuously with temperature, so actuation could in principle be proportional rather than switch-like, reducing functional fatigue from abrupt shape changes.","The observed exponents place the transition near the mean-field universality class, meaning predictions from Landau-type descriptions may transfer to this composite.","The effect should be manufacturable by layered deposition (coherent metallic superlattices) or by coherent precipitates in conventional metallurgy, not only by external stress.","Stiffness should show a sharp V-shaped dip near the critical temperature at high NiAl fractions, matching the simulated biaxial modulus and the reported gum-metal behavior."],"supporting_citations":[{"why":"Establishes the NiAl/Ni63Al37 nanolaminate model and volume-fraction methodology on which the present study builds.","marker":"[27]"},{"why":"Shows that interfacial stresses from coherent phases produce ultra-low stiffness, motivating the landscape-engineering picture.","marker":"[25]"},{"why":"Provides the hetero-epitaxial integration approach and tunable martensitic response used here.","marker":"[26]"},{"why":"Supplies the interatomic potential used for all nickel-aluminum simulations.","marker":"[29]"},{"why":"Provides the experimental Fe-Pd critical martensite under external stress that the material-level mechanism extends.","marker":"[17]"},{"why":"Reports elastically confined titanium-alloy behavior with continuous, higher-order martensitic transformation that this mechanism explains.","marker":"[23]"},{"why":"Underlies the stress and strain integration used to compute the free energy landscapes.","marker":"[35]"}],"fun_headline_variants":["Epitaxial stress turns abrupt martensite into continuous","Stress engineering makes martensite transformation second order","Alloy stress flips martensite from first to second order","Epitaxial stress induces continuous martensite in Ni-Al","Metal martensite transition becomes continuous under stress"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The mechanism assumes the two phases stay atomically coherent, sharing the same in-plane lattice parameter across a defect-free interface, for the whole transformation; if dislocations, interface sliding, or decohesion relax that constraint, the free energy reverts to a weighted average with a barrier, and the transformation becomes first order again.","fun_headline_variants_meta":{"raw":{"variants":["Epitaxial stress turns abrupt martensite into continuous","Stress engineering makes martensite transformation second order","Alloy stress flips martensite from first to second order","Epitaxial stress induces continuous martensite in Ni-Al","Metal martensite transition becomes continuous under stress"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000452,"raw_usage":{"total_tokens":2242,"prompt_tokens":880,"completion_tokens":1362,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":1282}},"tokens_in":496,"tokens_out":1362,"duration_ms":10234,"temperature":1.0,"reasoning_tokens":1282,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:16:17.125889+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Cool an actual 65 at. % NiAl NiAl/Ni63Al37 laminate, or run a molecular dynamics cell with interfaces free to decohere or with pre-existing dislocations, and measure the strain-temperature curve: a discontinuous jump in transformation strain, a nonzero cooling-heating hysteresis loop, or a barrier in the computed free energy landscape at 65 at. % NiAl would all show the order has not been changed by the epitaxial stress.","supporting_citations":[],"review_version":1}