{"id":"22953fbb-95d6-49bc-b68d-0c3e85346327","arxiv_id":"2412.04246","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A strain-engineered magnetostrictive nanomagnet can produce the three-state activation function needed for a ternary stochastic neuron, according to Landau-Lifshitz-Gilbert simulations.","lead":"The authors propose a single magnetostrictive nanomagnet, squeezed by strain, that can act as a ternary stochastic neuron with three random output states (-1, 0, +1). If the simulation is correct, this gives a compact hardware building block for ternary neural networks and other stochastic computing systems.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stress-anisotropy energy values are internally inconsistent by a factor of ~8; Eq. (4) with the stated 100 nm x 2 nm FeGa disk and 50 MPa gives ~76 kT, not the quoted 9.85 kT, so the reported stress-plateau relationship is not reproducible as written.","rationale":"The qualitative mechanism is plausible: compressive stress along y creates an x easy axis, spin-polarized current biases the magnetization toward +/-y, and the resulting time-averaged my versus Is can exhibit a plateau. I do not elevate the macrospin approximation to the primary concern, because the disk is a standard low-barrier geometry and the biaxial-to-uniaxial approximation is actually exact for equal and opposite in-plane stresses, just with a factor-of-two amplitude. The load-bearing problem is internal consistency of the stress-energy scale. Eq. (4) with the stated parameters gives a barrier of about 76 kT at 50 MPa and 152 kT at 100 MPa, while Section 4.2 quotes 9.85 and 19.7 kT. The reported plateau onset (no plateau at 20 MPa, plateau at 40 and 80 MPa) matches the quoted kT scale only if the effective stress is about 6.5 MPa for the 50 MPa label, a factor of about 7.7 smaller than stated. This is not a cosmetic typo: the threshold current at which the neuron leaves the zero state is set by the stress-anisotropy energy, so a factor-of-eight error shifts the predicted operating point by an order of magnitude. Because the paper provides no code and no data, the reader cannot resolve the contradiction. The central claim may survive after correction, but the manuscript as written does not support it reproducibly. The CONDITIONAL verdict is therefore retained, with the additional condition that the authors reconcile the stress-energy scale and provide simulation parameters or code.","tokens_in":9017,"tokens_out":18864,"duration_ms":201975,"concrete_test":"Recompute the stress-anisotropy barrier E = (3/2) lambda_s sigma Omega for the stated 100 nm x 2 nm FeGa disk at sigma = 20, 40, 50, 80, and 100 MPa and compare with the 9.85 and 19.7 kT values quoted in Section 4.2. Then rerun the Section 3 LLG at sigma = 20 MPa exactly as written (Ms=1.32e6 A/m, lambda_s=266.6 ppm, alpha=0.017, T=300 K, 1 microsecond). If Eq. (4) is correct, K/kT is about 30 at 20 MPa and a plateau should appear; if no plateau appears, the implementation of Eq. (3) contains a scale error. This distinguishes a typographical factor-of-eight error from a genuine stress-threshold effect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's own stress-energy accounting is internally inconsistent. Eq. (4) gives E_stress = -(3/2) lambda_s sigma Omega cos^2(theta). With the stated disk (d=100 nm, t=2 nm, Omega=1.57e-23 m^3), lambda_s=266.6 ppm=2.666e-4, and sigma=50 MPa, this yields E = 3.14e-19 J, about 76 kT at 300 K; for 100 MPa it is about 152 kT. Section 4.2 instead quotes 9.85 and 19.7 kT for the same 50 and 100 MPa values, which correspond to sigma about 6.5 and 13 MPa, roughly 7.7 times smaller. The reported Fig. 3 behavior (no plateau at 20 MPa, plateau at 40 and 80 MPa) is consistent with the quoted 3.9/7.9/15.8 kT scale, but not with Eq. (4) as written, which would give 30/60/120 kT and should already produce a plateau at 20 MPa because the initial -y state relaxes to the x easy axis and thermal escape over a 30 kT well is negligible on a 1 microsecond time scale. Thus either Eq. (3), Eq. (4), the material/geometry parameters, or the stress labels in Fig. 3 are in error. Since the plateau width is the defining property of the TSN activation function, this factor-of-eight energy-scale ambiguity makes the central simulation result quantitatively untrustworthy.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a design for a ternary stochastic neuron (TSN) based on a single circular magnetostrictive (FeGa) nanomagnet with zero in-plane shape anisotropy, subjected to uniaxial strain and injected with a spin-polarized current. The authors carry out stochastic Landau-Lifshitz-Gilbert (LLG) simulations and find that when the product of magnetostriction and stress (λ_s σ) is negative, the time-averaged y-component of magnetization <my(t)> versus current Is exhibits a plateau around Is=0, giving three stable states -1, 0, +1. For λ_s σ positive, the activation curve is asymmetric and initial-condition dependent. The paper concludes that this is the first nanomagnetic implementation of a TSN. The central claim is therefore that the strain-induced anisotropy creates a potential well along the x-axis, producing the zero-output plateau needed for ternary behavior.","tokens_in":9420,"tokens_out":8253,"duration_ms":73871,"significance":"If the simulation results are quantitatively reliable, the proposed device would be a simple, compact building block for ternary stochastic neural networks and related probabilistic computing architectures, with potential advantages in area and energy efficiency. The study uses a standard LLG model with material parameters taken from the literature (α=0.017, Ms=1.32×106 A/m, λ_s=266.6 ppm) and does not fit any parameter to a target activation function; the three-state behavior emerges from the physics. The paper also provides a clear qualitative explanation of the plateau mechanism in terms of stress-induced anisotropy energy. However, the quantitative inconsistency in the stress-energy accounting (Section 4.2 versus Eq. (4) and Fig. 3) and the unquantified biaxial-to-uniaxial approximation mean that the specific predicted plateau widths are not reproducible as written; these issues must be resolved before the results can be trusted.","major_comments":[{"comment":"The stress-anisotropy energy values reported in Section 4.2 are internally inconsistent with Eq. (4) and with the stated geometry and stress. For the d=100 nm, t=2 nm FeGa disk (volume Ω=1.57×10^-23 m^3), λ_s=266.6 ppm, and σ=50 MPa, Eq. (4) gives E=(3/2)λ_s σ Ω = 3.14×10^-19 J ≈ 76 kT at T=300 K, not the quoted 9.85 kT. For 100 MPa the correct value is about 152 kT, not 19.7 kT. The quoted values correspond to σ ≈ 6.5 and 13 MPa, roughly a factor of 7.7 smaller. Moreover, Section 4.2 refers to '50 MPa and 100 MPa' while Fig. 3 uses 20, 40, and 80 MPa. If Eq. (4) is correct, even 20 MPa gives a barrier of about 30 kT, which should be sufficient to produce a plateau on the 1 µs simulation time scale; the absence of a plateau at 20 MPa in Fig. 3 is then inconsistent with the stated equations. Since the plateau width is the defining property of the TSN activation function, this factor-of-eight energy-scale discrepancy makes the central simulation result quantitatively untrustworthy. The authors must correct the equations, the geometric/material parameters, or the stress labels, and ideally provide a version of Fig. 3 with the computed energy barriers for each stress value.","section":"Section 4.2, Eq. (4), Fig. 3"},{"comment":"The paper approximates the biaxial strain generated by the piezoelectric gate as a uniaxial strain along the y-axis with a 'larger' magnitude, but it never quantifies this replacement. For a biaxial stress state with σ_xx = -σ_yy, the magnetoelastic energy has the form E = -(3/2)λ_s(σ_xx cos²θ + σ_yy sin²θ)Ω, which is not equivalent to a simple uniaxial term E = -(3/2)λ_s σ_eff cos²θ with an unspecified σ_eff. The effective uniaxial constant depends on the ratio of the two stress components and on the assumed energy expression, and the difference is a factor of order 2 in energy. Because the plateau width is set by the barrier height, this unquantified approximation directly affects the central quantitative result. The authors should either specify the effective uniaxial stress value used in Eq. (3) and how it was derived from the biaxial strain, or implement the full biaxial energy in the simulations.","section":"Section 3, biaxial-to-uniaxial approximation"},{"comment":"The paper repeatedly claims to have 'implemented' a TSN and calls this 'the first and only nanomagnetic implementation of a TSN.' However, the manuscript presents only stochastic LLG simulations; no device is fabricated or measured. The word 'implementation' is thus an overstatement that misrepresents the contribution as an experimental realization. The authors should consistently describe their work as a simulation-based proposal or design, and moderate corresponding novelty claims in the abstract, Section 6, and the title if needed.","section":"Section 3 and Section 6"}],"minor_comments":[{"comment":"The headings 'Positive λsσ product or compressive stress' and 'Negative λsσ product or tensile stress' are reversed with respect to the standard relation: for FeGa with positive λ_s, compressive stress (σ<0) gives a negative λ_s σ product and tensile stress (σ>0) gives a positive product. The text within the sections is correct, but the headings should be swapped to avoid confusion.","section":"Section 4.1 and 4.2 headings"},{"comment":"The sentence 'For the two stress values considered here, 50 MPa and 100 MPa' is not consistent with Fig. 3, which uses 20, 40, and 80 MPa. Please correct the stress values cited in the energy-barrier discussion.","section":"Section 4.2"},{"comment":"Reference [14] lists the year as '2027' for 'Trained ternary quantization'; this appears to be a typo (likely 2017 or 2018). Please correct.","section":"References"},{"comment":"The abstract and some sections contain formatting issues such as 'CIF AR-10' (with a space) and '10 7' (instead of 10^7) for the sample count. Please fix these.","section":"Abstract and text formatting"},{"comment":"The sentence 'We point out that that the contribution here is not just with respect to the activation function' contains a doubled 'that'. Please correct.","section":"Section 6"}],"recommendation":"major_revision","confidential_remarks":"The stress-energy inconsistency is a genuine load-bearing error that undermines the reproducibility of the central simulation result. I would ask the authors to provide a detailed parameter table, the exact form of the stress-field implementation, and ideally the simulation code/data so that the reported plateau widths can be verified. The paper is a simulation-based proposal, not an experimental implementation, and the novelty claim should be calibrated accordingly. If these issues are resolved, the conceptual idea may still be publishable, but in its current form I cannot recommend acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe core idea is genuinely new: use strain-induced anisotropy in a magnetostrictive nanomagnet to create a zero-output plateau in the activation function, which is exactly what a ternary stochastic neuron needs. To my knowledge, it's the first nanomagnetic TSN proposal, and the mechanism is explained clearly. The LLG simulations use standard equations and plausible material parameters, and the qualitative physics — compressive stress with λsσ < 0 pins the magnetization along the hard axis, producing a dead zone around zero current — is sound.\n\nThe quantitative core, though, does not hold up as written. Eq. (4) with the stated disk geometry (100 nm diameter, 2 nm thick FeGa), λs = 266.6 ppm, and σ = 50 MPa gives ≈76 kT at 300 K. Section 4.2 quotes 9.85 kT for 50 MPa and 19.7 kT for 100 MPa — off by a factor of about 8. Worse, Section 4.2 refers to 50 and 100 MPa as the two stress values considered, while Fig. 3 shows 20, 40, and 80 MPa. That is not cosmetic. The plateau width is the defining feature of the TSN, and it depends directly on the stress-energy scale. As written, the simulations are not reproducible, and the central quantitative claim is untrustworthy.\n\nAdd to that the usual simulation-only caveats: no experiment, no code or data, no error bars. The macrospin assumption and the biaxial-to-uniaxial strain approximation are asserted, not validated. And calling a simulation an 'implementation' in the conclusion overstates the result.\n\nThe underlying physics is likely salvageable — the plateau concept is sound, and correcting the energy accounting would almost certainly preserve the qualitative behavior, just with different stress values. But the authors need to reconcile Eq. (4), the quoted kT numbers, and the figure labels, and ideally release the simulation code. This paper deserves a serious referee because the idea is worth engaging with, but it needs major revision before it is reliable. I would not cite it in its current form.\n\nReading group: maybe. The inconsistency makes for a good exercise in calculation-checking, but the paper itself needs fixing first.","headline":"A plausible new nanomagnetic TSN mechanism undermined by an unexplained ~8x error in the stress-energy accounting; needs a serious referee and a substantial revision.","tokens_in":9898,"tokens_out":4926,"would_cite":false,"duration_ms":44579,"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":"A single strained nanomagnet can act as a ternary stochastic neuron.","keywords":["ternary stochastic neuron","magnetostrictive nanomagnet","uniaxial strain","activation function","spin-polarized current","zero-energy-barrier magnet","neuromorphic computing"],"falsifier":"A micromagnetic simulation of the same 100 nm diameter, 2 nm thick FeGa disk under -80 MPa uniaxial compressive stress and spin-polarized current, without the macrospin assumption, would settle whether the zero-current plateau survives in realistic nonuniform magnetization dynamics.","tokens_in":8835,"feed_emoji":"🧲","tokens_out":6307,"duration_ms":57905,"temperature":0.7,"pith_summary":"This paper aims to show that a single zero-energy-barrier magnetostrictive nanomagnet, subjected to uniaxial compressive stress and injected with spin-polarized current, can implement the activation function of a ternary stochastic neuron (TSN) — a neuron whose output randomly takes one of three values, -1, 0, or +1. The crux is that a negative product of magnetostriction and stress (compressive stress on FeGa) creates an energy minimum with the magnetization along the x-axis, so the time-averaged y-component <my(t)> stays at zero for small currents and only rises to ±1 once the current is strong enough. This produces the flat 'plateau' around zero current that makes the 0 state stable, which a standard sigmoid activation cannot do. The authors present stochastic Landau-Lifshitz-Gilbert simulations showing the plateau appears for compressive stresses of 40 and 80 MPa, and they argue this is the first nanomagnetic TSN implementation.","feed_headline":"Strained nanomagnet yields three-state stochastic neuron","feed_subtitle":"Compressive stress creates a zero-current plateau in the activation curve, enabling stable -1, 0, +1 states","key_machinery":"The load-bearing object is the stress-anisotropy energy E = -(3/2)λsσΩcos²θ together with the spin-transfer torque from the injected current. When λsσ<0, the energy minimum sits at θ=90° (magnetization along x, my=0), holding the neuron in the zero state until the current's torque pulls it toward ±y; the dynamics are simulated with the stochastic LLG equation including Slonczewski and field-like torques (relative weights A=1, B=0.3) and thermal noise.","core_discovery":"The central claim is that a zero-energy-barrier (shape-isotropic) magnetostrictive nanomagnet under uniaxial compressive stress — the sign that makes λsσ negative for FeGa — gives a three-level activation function <my(t)> versus spin-polarized current Is, with a plateau around Is=0 that constitutes the stable 0 state of a ternary stochastic neuron. The same device under tensile stress (positive λsσ) instead pins the magnetization near ±y and produces an activation curve that depends on the initial state, which the paper shows is not useful for a TSN. The plateau width grows with stress magnitude, and the resulting activation function also acts as the threshold-based ternary function used in ternary neural networks.","pith_inferences":["Relaxing the macrospin assumption in a micromagnetic simulation could blur the plateau, since strain-induced fields vary across the disk; this is a concrete test of whether the TSN works outside the single-domain idealization.","Ensemble stress non-uniformity will spread plateau widths across many neurons; the paper argues TSN function survives, but the effect on network-level training convergence is left open.","The same strain-anisotropy trick could generalize to higher-radix stochastic neurons by engineering multiple energy minima, though the paper does not pursue that extension.","The positive-λsσ activation curve resembles asymmetric neural-network activations such as ReLU-family functions, suggesting a possible separate use for strained nanomagnets as nonlinear transfer elements."],"forward_implications":["A TSN built this way occupies the same chip area as a binary stochastic neuron but encodes three states, increasing information density.","Stress magnitude tunes the plateau width, giving a voltage-controlled window for the stable 0 state.","The activation function also implements threshold-based ternary functions (Eq. 5), the building block for ternary neural networks that minimize distance between full-precision and ternary weights.","Because the piezoelectric gate is a charged capacitor at steady state, holding the strain consumes no standby power, and a lattice-mismatched substrate could supply the strain without any voltage, at the cost of reconfigurability."],"supporting_citations":[{"why":"Supplies the low-barrier nanomagnet design and tanh-like activation curve of a binary stochastic neuron that the TSN work extends.","marker":"[17]"},{"why":"Provides the target TSN activation function f(x)=1.5tanh(x)+0.5tanh(-3x) with a stable zero plateau that motivates the strain mechanism.","marker":"[20]"},{"why":"Defines ternary weight networks and threshold-based ternary functions used in Section 5 to relate the obtained activation to practical ternary neural networks.","marker":"[13]"},{"why":"Demonstrates localized strain generation in a thin-film piezoelectric, the gate mechanism for applying uniaxial stress.","marker":"[23]"},{"why":"Supplies the spin-transfer-torque model with relative Slonczewski and field-like torque factors (A=1, B=0.3) used in the LLG simulations.","marker":"[25]"},{"why":"Gives the effective magnetic field expressions including the stress field term and the stress-anisotropy energy formula.","marker":"[26]"},{"why":"Provides the FeGa material parameters (Ms, λs, α) used in the simulations.","marker":"[24]"}],"fun_headline_variants":["One nanomagnet, three states: a ternary stochastic neuron","Stress gives nanomagnet a stable zero state for ternary logic","Ternary stochastic neuron from a single strained nanomagnet","Zero-current plateau turns nanomagnet into ternary neuron","Magnetostrictive nanomagnet: ternary stochastic neuron with strain"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction depends on the 100 nm disk behaving as a single macrospin and on treating the gate-induced biaxial strain as a stronger uniaxial strain; if either approximation is wrong, the plateau could shift, narrow, or disappear.","fun_headline_variants_meta":{"raw":{"variants":["One nanomagnet, three states: a ternary stochastic neuron","Stress gives nanomagnet a stable zero state for ternary logic","Ternary stochastic neuron from a single strained nanomagnet","Zero-current plateau turns nanomagnet into ternary neuron","Magnetostrictive nanomagnet: ternary stochastic neuron with strain"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000903,"raw_usage":{"total_tokens":3847,"prompt_tokens":871,"completion_tokens":2976,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":487,"completion_tokens_details":{"reasoning_tokens":2892}},"tokens_in":487,"tokens_out":2976,"duration_ms":18901,"temperature":1.0,"reasoning_tokens":2892,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:35:22.786511+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A micromagnetic simulation of the same 100 nm diameter, 2 nm thick FeGa disk under -80 MPa uniaxial compressive stress and spin-polarized current, without the macrospin assumption, would settle whether the zero-current plateau survives in realistic nonuniform magnetization dynamics.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the low-barrier nanomagnet design and tanh-like activation curve of a binary stochastic neuron that the TSN work extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the target TSN activation function f(x)=1.5tanh(x)+0.5tanh(-3x) with a stable zero plateau that motivates the strain mechanism."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates localized strain generation in a thin-film piezoelectric, the gate mechanism for applying uniaxial stress."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the spin-transfer-torque model with relative Slonczewski and field-like torque factors (A=1, B=0.3) used in the LLG simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the effective magnetic field expressions including the stress field term and the stress-anisotropy energy formula."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the FeGa material parameters (Ms, λs, α) used in the simulations."}],"review_version":1}