{"id":"d385f92b-6685-4029-80cf-4c27bfa2085c","arxiv_id":"1908.00602","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Surface inversion breaking in collinear magnets is predicted to produce a bound ferroelectric polarization through grad(m^2) magnetoelectric coupling and through DMI-driven spin canting.","lead":"Surface-induced electric polarization is shown to arise in magnetic materials whose bulk is collinear, through the gradient of the magnetization and through surface Dzyaloshinskii-Moriya interactions. The paper maps two phenomenological mechanisms and estimates measurable polarization ranges in thin films and multilayers.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Detectability claim rests on an unquantified magnetoelectric coupling; without a material-specific γχE the predicted polarization spans ten orders of magnitude and may fall below detection.","rationale":"The reader's weakest_assumption identifies the unquantified magnetoelectric coupling γχE as the load-bearing premise, and my stress test agrees. The paper is a legitimate phenomenological GL analysis: the symmetry argument for the allowed coupling is sound, and the DMI route is a standard extension. However, the quantitative claim of detectability is conditional on γχE exceeding ~10^-19 sm/A, while the cited literature range spans ten orders of magnitude. This is not an internal inconsistency, but it is a missing material-specific verification that undermines the broad claim in the abstract and Discussion. The reader's CONDITIONAL verdict is appropriate: accept the symmetry-based mechanism as plausible but require material-specific coupling information before the detectability prediction can be endorsed. I found no fatal flaw, so I do not change the verdict. The proposed DFT test would directly settle whether a representative collinear ferromagnet supports the effect at a detectable magnitude.","tokens_in":11818,"tokens_out":7166,"duration_ms":78705,"concrete_test":"Perform a first-principles DFT calculation for a (001) surface slab of a simple collinear ferromagnet such as bcc Fe, with spin-orbit coupling explicitly disabled to isolate the γ∇(m^2) mechanism. Compute the planar-averaged magnetization profile m(z) and the electric polarization profile p(z) using the modern theory of polarization (or a slab Berry-phase calculation). Fit the relation p(z) = -χE γ d(m^2)/dz using a DFT-derived χE for the same slab to extract the effective product γχE. If the fitted value is smaller than 10^-19 sm/A, the paper's stated detection threshold is not met for that material, and the broad detectability claim would require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that collinear magnetism generates a detectable surface polarization depends on the coupling term p·γ∇(m^2) in Eq. (3). The Introduction explicitly states 'the underlying assumption is that we deal with predominantly magnetic materials with non-zero coupling between magnetic and ferroelectric order parameters,' but γ is never computed, measured, or constrained for any specific material. The paper cites a literature range for γχE of 10^-23 to 10^-13 sm/A (Sec. II, Refs. 29-32) and then states that γχE > 10^-19 sm/A would be needed for a measurable polarization. In footnote 34 the authors acknowledge that 'we focus on the necessary value of coupling constant that would provide measurable polarization,' effectively converting the quantitative prediction into a conditional statement. Since the amplitude of p(z) in Eq. (4) is directly proportional to γχE, a material with γχE at the lower end of the cited range yields a polarization orders of magnitude below current experimental sensitivity, while the Discussion asserts 'The polarization we predict is detectable within the current accuracy of the experimental techniques.' Without a material-specific evaluation of γ, the qualitative symmetry argument stands but the detectability claim is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a Ginzburg-Landau theory for electric polarization induced near surfaces of collinear magnetic materials. In Sec. II the authors solve the one-dimensional GL equation for the magnetization profile near a (001) surface and couple it to polarization through a term p·(γ∇(m^2)+...). They obtain a closed-form p(z) peaked a few correlation lengths from the surface. In Sec. III they consider two DMI-based routes in thin films: an in-plane anisotropy case with a one-dimensional twist and an easy-axis anisotropy case with an axisymmetric soliton; they solve the resulting equations numerically and plot polarization profiles. The paper estimates possible polarization magnitudes and asserts that the predicted polarization is detectable with current experimental techniques.","tokens_in":12118,"tokens_out":6374,"duration_ms":59001,"significance":"If the qualitative claim holds, the paper identifies a general mechanism: any magnetic material with a nonzero magnetoelectric coupling γ and a surface-induced magnetization gradient develops a boundary-localized polarization even in the absence of DMI. This is a useful complement to material-specific first-principles studies. The analytic part of Sec. II is clean, the numerical scheme in Appendix B is described in sufficient detail to be reproducible, and no result is fitted to a target polarization. The main limitation is that the coupling γ is not constrained for any material, so the quantitative predictions span many orders of magnitude and the detectability claim is conditional rather than unconditional.","major_comments":[{"comment":"The prefactor in Eq. (4) is off by a factor of two. Using m(z)=m∞ tanh[(z+a0)/(√2 ξ)], one finds d(m^2)/dz = √2 (m∞^2/ξ) tanh u (1−tanh^2 u), so p(z) = −χE γ √2 (m∞^2/ξ) tanh u (1−tanh^2 u), not −χE γ m∞^2/(√2 ξ) tanh u (1−tanh^2 u). Consequently the maximum polarization should be pmax = −2√2/(3√3) χE γ m∞^2/ξ, which is twice the value implied by the displayed Eq. (4). The approximate estimate pmax ≈ −χE γ m∞^2/ξ remains numerically valid, but the printed formula needs correction.","section":"Sec. II, Eq. (4)"},{"comment":"The claim in the Discussion that 'the polarization we predict is detectable within the current accuracy of the experimental techniques' is not supported by the analysis. The polarization amplitude is proportional to γχE, and the paper itself gives the literature range 10^-23 to 10^-13 sm/A for γχE and the threshold γχE > 10^-19 sm/A for measurability. Because no material-specific value or independent constraint on γ is provided, a material at the lower end of this range yields p ~ 10^-10 µC/cm^2, far below current experimental sensitivity. Footnote 34 states that the authors 'focus on the necessary value of coupling constant that would provide measurable polarization,' which is precisely a conditional statement; the unconditional detectability assertion in the Discussion overreaches. The authors should either identify a material or material class where γχE exceeds the threshold, or explicitly reframe the prediction as conditional on the coupling.","section":"Secs. II and IV; footnote 34"},{"comment":"The numerical DMI calculations are made dimensionless and illustrated with A=1, D=1, k=1, γ=1, χE=1.5, but no mapping is given from these parameter choices to physical materials, and no systematic variation of D, k, the profile function f(z), or the anisotropy is presented. The estimates p ~ 10^-10 to 10^4 µC/cm^2 in Secs. III A and III B rest on the same unconstrained product DγχE as in Sec. II, and on the assumed gradients Δφ/Δz ≈ 10^5 m^-1 and |∂θ/∂z| ≈ 10^4 m^-1, whose compatibility with the chosen dimensionless parameters is not demonstrated. The statement that both profiles f(z) lead to the same physics qualitatively is asserted but not shown. This limits the predictive content of the DMI route.","section":"Sec. III, Figs. 2 and 3"}],"minor_comments":[{"comment":"Figure 1 uses 'τ = 1, 2, ...' but τ is defined as T/Tc, and the curves are not explained in the caption; clarify which reduced temperatures are shown and specify the units of the axes.","section":"Fig. 1 and Sec. II"},{"comment":"The sentence following Eq. (5) mislabels the free-energy terms: it refers to the second term as 'the energy of the electric polarization with susceptibility χE', but in Eq. (5) the second term is the DMI term DfL and the polarization energy appears only through F_E. This should be reworded.","section":"Sec. III, after Eq. (5)"},{"comment":"The units listing for γ ('in units of inverse polarization') and the parenthetical '(in ,' before 'the second is the energy' appear incomplete; add the missing units for A and χE.","section":"Sec. III, units paragraph"},{"comment":"In the finite-difference description, the identification of nodes (i,0) with (i,N_z) and (i,N_z+1) with (i,2) is stated without a figure reference; label the computational grid nodes as referenced in Fig. 4 so that the periodic boundary conditions are unambiguous.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The central symmetry argument appears sound, and the paper contains a clean analytic solution in Sec. II plus a reproducible numerical scheme in Appendix B. The main risk is overclaiming detectability on the basis of an unconstrained coupling constant. If the authors revise the quantitative claims to be explicitly conditional on γχE, or provide a material-specific estimate, the paper could be publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a decent phenomenological paper with a genuinely new qualitative claim — collinear bulk magnets can develop a surface-bound electric polarization via the magnetization gradient — but the quantitative detectability claim outruns the evidence. The GL derivation in Sec. II is clean and the Mills profile is used properly. The explicit p(z) profile peaked a few correlation lengths from the surface is a real addition, and the authors are honest in footnote 34 that the coupling constant is the weak point. The DMI sections assemble Bogdanov–Rössler machinery into a second mechanism; the numerics are plausible but under-documented.\n\nSoft spots. Eq. (4) appears to have a factor-of-two error: taking the derivative of m = m∞ tanh((z+a0)/(√2 ξ)) gives ∇(m²) coefficient √2 m∞²/ξ, not 1/√2. The pmax formula in the text uses the correct √2 coefficient, so Eq. (4) and the text disagree internally. Easy fix, but it should be flagged.\n\nThe bigger issue is the detectability claim. The paper itself gives γχE ranging from 10^-23 to 10^-13 sm/A and says γχE > 10^-19 needed for measurable polarization. That makes the Discussion assertion \"detectable within current accuracy\" a statement about the upper end of a ten-order-of-magnitude window, not a prediction. No specific material is evaluated, and γ is never computed or constrained. The symmetry argument survives, but the amplitude claim does not. This is a conditional result, and the paper should say so in the abstract.\n\nThe DMI sections use arbitrary dimensionless parameters (A=D=k=γ=1, χE=1.5) and no convergence study. The mathematics in Appendix B is standard finite differences with Newton's method, so the method is sound, but the link to real materials is thin.\n\nOn the positive side: no fitting to targets, citations are used appropriately, and self-citation (Ref. 16) appears in the context of distinguishing from phase-dislocation mechanisms, not self-promotion. The central qualitative mechanism is credible.\n\nWho this is for: experimentalists working on interfacial DMI and multiferroic heterostructures, and theorists wanting a simple GL framework for surface magnetoelectric effects. It deserves a serious referee — referee time is worthwhile, with the main requests being correction of Eq. (4), softening the detectability claim, and better documentation of the numerical parameters.","headline":"Worth a referee's time, but the detectability claim is not backed by the paper's own coupling constant ranges.","tokens_in":12616,"tokens_out":2479,"would_cite":false,"duration_ms":22987,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["75.85.+t","75.70.-i","77.80.-e"],"model":"deepseek-v4-flash","headline":"A Ginzburg-Landau analysis shows that surfaces make collinear magnets ferroelectric even when their bulk is not.","keywords":["surface multiferroicity","Ginzburg-Landau free energy","magnetoelectric coupling","collinear magnetism","Dzyaloshinskii-Moriya interaction","thin films","multilayers","ferroelectric polarization"],"falsifier":"A first-principles calculation of $\\gamma\\chi_E$ for a clean (001) surface of a simple collinear magnet such as bcc iron would settle the quantitative claim: if the computed value is below about $10^{-19}\\,\\mathrm{sm/A}$, the polarization predicted by the paper's formula $|p_{\\max}|\\sim\\chi_E|\\gamma|m_\\infty^2/\\xi$ is too small to detect, and the central measurable prediction fails.","tokens_in":11611,"feed_emoji":"🧲","tokens_out":12672,"duration_ms":123564,"temperature":0.7,"pith_summary":"This paper tries to establish that a surface alone can make a bulk collinear ferromagnet or antiferromagnet multiferroic. In the Ginzburg-Landau description, magnetization must fall to a boundary value at the surface, creating a gradient of $m^2$; the magnetoelectric term $p\\cdot\\gamma\\nabla(m^2)$ then forces a local electric polarization $p(z)$ peaked a few correlation lengths from the surface. A second route appears when surface Dzyaloshinskii-Moriya interactions wind the magnetization direction, giving a polarization tied to $d\\varphi/dz$. Using reported coupling values, the predicted polarizations range from $10^{-10}$ to $10^4\\,\\mu\\mathrm{C/cm}^2$, and the paper argues they are detectable with current experimental techniques if the coupling products exceed about $10^{-19}$ to $10^{-17}\\,\\mathrm{sm/A}$. This matters because surfaces are everywhere in magnetic devices, so nominally non-multiferroic magnets could host local ferroelectric order at their boundaries.","feed_headline":"Surface alone can make collinear magnets ferroelectric","feed_subtitle":"Magnetization gradients at a boundary couple to electric polarization, creating a measurable surface multiferroic response","key_machinery":"The key machinery is the Ginzburg-Landau free energy with boundary conditions, augmented by a linear magnetoelectric term and a quadratic electric term. The object that carries the first mechanism is the dimensionless magnetization profile $f(z)=\\tanh((z+a_0)/(\\sqrt{2}\\xi))$, whose gradient $\\nabla(m^2)$ plays the role of the inversion-symmetry-breaking field; substituting it into $p=-\\chi_E\\nabla(\\gamma m^2)$ yields the polarization. For the second mechanism, the central object is the Lifshitz invariant $m_x\\,dm_y/dz-m_y\\,dm_x/dz$ (and its three-dimensional generalizations), multiplied by a surface profile $f(z)$ representing the spatially inhomogeneous DMI strength; minimizing the coupled equations gives $p_z=-D\\gamma\\chi_E f(z)\\,d\\varphi/dz$ and the skyrmion-like expression $p=\\chi_E D f\\gamma\\cos^2\\theta|\\partial\\theta/\\partial z|$. The DMI profile and the spin-angle gradients carry the boundary effect into the polarization.","core_discovery":"At the level of a Ginzburg-Landau free energy, the paper shows that a collinear ferromagnet or antiferromagnet with a (001) surface has an equilibrium magnetization $m(z)=m_\\infty\\tanh((z+a_0)/(\\sqrt{2}\\xi))$, whose square has a nonzero gradient in a boundary layer. Adding the magnetoelectric coupling $F_{ME}=p\\cdot(\\gamma\\nabla(m^2)+\\gamma'[m(\\nabla\\cdot m)-(m\\cdot\\nabla)m]+\\dots)$ plus the electric term $p^2/(2\\chi_E)$ and minimizing gives $p(z)=-\\chi_E\\gamma\\nabla(m^2)$, namely $p(z)=-\\frac{m_\\infty^2}{\\sqrt{2}\\xi}\\chi_E\\gamma\\,\\tanh\\!\\big(\\frac{z+a_0}{\\sqrt{2}\\xi}\\big)\\big[1-\\tanh^2\\!\\big(\\frac{z+a_0}{\\sqrt{2}\\xi}\\big)\\big]$, peaked at $z=\\sqrt{2}\\xi\\tanh^{-1}(1/\\sqrt{3})$. The second route holds when a surface-induced DMI exists: the Lifshitz-invariant coupling yields $p_z=-D\\gamma\\chi_E f(z)\\,d\\varphi/dz$ for in-plane spin rotation, and $p=\\chi_E D f\\gamma\\cos^2\\theta|\\partial\\theta/\\partial z|$ for localized three-dimensional textures, giving the same qualitative surface-bound polarization. The paper claims both mechanisms produce polarization detectable with current techniques for realistic parameter ranges.","pith_inferences":["The symmetry argument implies surface-bound polarization should be generic to any magnetoelectric material, so published surface-sensitive data on nominally non-multiferroic magnets may already contain an overlooked electric signal.","Temperature sweeps provide a direct test: as $T\\to T_c$ from below, $\\xi$ diverges, so the predicted peak moves farther from the surface while $|p_{\\max}|\\sim 1/\\xi$ diminishes; a measured profile with this scaling would confirm the mechanism.","Engineering the interfacial DMI, for instance through heavy-metal/ferromagnet bilayers, should turn the chirality of surface spin textures into a voltage-addressable polarization, linking skyrmion chirality to ferroelectric switching.","Because the coupling $\\gamma$ is material-specific and not derived here, a first-principles calculation of $\\gamma$ for a clean Fe or Ni surface would settle whether the predicted polarization is observable or purely formal."],"forward_implications":["A clean surface of any collinear ferromagnet or antiferromagnet with nonzero $\\gamma\\chi_E$ should show a polarization layer of thickness set by $\\xi$, even when the bulk is centrosymmetric and nonpolar.","The predicted maximum polarization scales as $|p_{\\max}|\\sim \\chi_E|\\gamma| m_\\infty^2/\\xi$; for the reported coupling range it spans $10^{-10}$ to $10^2\\,\\mu\\mathrm{C/cm}^2$, and becomes measurable once $\\gamma\\chi_E\\gtrsim10^{-19}\\,\\mathrm{sm/A}$.","Temperature controls the profile: $\\xi$ grows near $T_c$, so the polarization peak moves away from the surface and weakens as the transition is approached from below.","In thin films and multilayers, an engineered interfacial DMI creates a second route: polarization is proportional to the surface spin rotation rate $d\\varphi/dz$, so chiral textures (including skyrmion-like states) carry a local electric polarization.","The existence of a surface magnetoelectric response means thin-film devices based on nominally non-multiferroic magnets may display parasitic or exploitable ferroelectric responses at interfaces."],"supporting_citations":[{"why":"Supplies the molecular-field Ginzburg-Landau equation and boundary condition that produce the tanh magnetization profile used throughout the first mechanism.","marker":"28"},{"why":"Supplies the Lifshitz-invariant DMI formulation, the spatially varying DMI profile $f(z)$, and the magnetic solutions used for the K>0 case.","marker":"38"},{"why":"Establishes the spin-current mechanism by which spin rotation induces polarization, the bulk analogue that the DMI surface route extends.","marker":"11"},{"why":"Gives the Ginzburg-Landau Lifshitz-invariant treatment of spiral multiferroics that the surface DMI term generalizes.","marker":"12"},{"why":"Provides the experimentally observed inverse effect, surface-induced magnetization in a ferroelectric, which this paper's surface polarization is the reciprocal of.","marker":"46"},{"why":"Provide the reported range of magnetoelectric coupling constants used to estimate observable polarization amplitudes.","marker":"29–32"},{"why":"Supplies the typical magnetic correlation length $\\xi\\approx10\\,\\mathrm{nm}$ used for the numerical estimates.","marker":"33"}],"fun_headline_variants":["Surfaces turn collinear magnets ferroelectric","Magnetic surfaces spark multiferroic polarization","Boundary layers induce ferroelectricity in magnets","Surface gradients make magnets ferroelectric","Collinear magnets get ferroelectric at surfaces"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the magnetic materials in question possess a nonzero coupling between magnetization gradients and electric polarization, with strength $\\gamma\\chi_E$ at least around $10^{-19}\\,\\mathrm{sm/A}$; if that coupling is absent or smaller, the surface polarization may be theoretically present but experimentally invisible.","fun_headline_variants_meta":{"raw":{"variants":["Surfaces turn collinear magnets ferroelectric","Magnetic surfaces spark multiferroic polarization","Boundary layers induce ferroelectricity in magnets","Surface gradients make magnets ferroelectric","Collinear magnets get ferroelectric at surfaces"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000269,"raw_usage":{"total_tokens":1633,"prompt_tokens":965,"completion_tokens":668,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":581,"completion_tokens_details":{"reasoning_tokens":598}},"tokens_in":581,"tokens_out":668,"duration_ms":6567,"temperature":1.0,"reasoning_tokens":598,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:44:42.390767+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A first-principles calculation of $\\gamma\\chi_E$ for a clean (001) surface of a simple collinear magnet such as bcc iron would settle the quantitative claim: if the computed value is below about $10^{-19}\\,\\mathrm{sm/A}$, the polarization predicted by the paper's formula $|p_{\\max}|\\sim\\chi_E|\\gamma|m_\\infty^2/\\xi$ is too small to detect, and the central measurable prediction fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the molecular-field Ginzburg-Landau equation and boundary condition that produce the tanh magnetization profile used throughout the first mechanism."},{"cited_title":"Emori et al","cited_arxiv_id":null,"evidence_quote":"Supplies the Lifshitz-invariant DMI formulation, the spatially varying DMI profile $f(z)$, and the magnetic solutions used for the K>0 case."},{"cited_title":"Katsura, N","cited_arxiv_id":null,"evidence_quote":"Establishes the spin-current mechanism by which spin rotation induces polarization, the bulk analogue that the DMI surface route extends."},{"cited_title":"Mostovoy, Phys","cited_arxiv_id":null,"evidence_quote":"Gives the Ginzburg-Landau Lifshitz-invariant treatment of spiral multiferroics that the surface DMI term generalizes."},{"cited_title":"Hu, L.-Q","cited_arxiv_id":null,"evidence_quote":"Supplies the typical magnetic correlation length $\\xi\\approx10\\,\\mathrm{nm}$ used for the numerical estimates."}],"review_version":1}