{"id":"b6fc266c-86c4-49ba-8ca0-3d16e10ca74b","arxiv_id":"2501.05025","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Monolayer NiBr2 has a cycloidal spin order with polarization from the gKNB mechanism, while monolayer NiI2 needs an additional p-d hybridization term proportional to sin(4πq) to explain its polarization.","lead":"This paper uses computer simulations to explain why two atomically thin magnetic insulators, NiBr2 and NiI2, develop an electric polarization when their magnetic spins form spiral patterns. It identifies two competing microscopic mechanisms and shows how to tell them apart by how the polarization changes with spin-orbit coupling and spiral wavelength.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The NiI2 decomposition in Eq. 13 is a two-parameter fit: the fitted A does not match the independently calculated gKNB amplitude, so the quantitative separation of PKM and Ppd is not microscopically anchored.","rationale":"The reader's conditional verdict identifies the same core problem: the quantitative decomposition for NiI2 is fit-driven rather than derived from independent microscopic calculations. My stress-test sharpens this by showing that the fitted A is not merely an amplitude absorbed from the residual, but actually conflicts with the gKNB value that the paper itself computes from the M matrix. This makes Eq. 13 a phenomenological fit with two free parameters, not a physically constrained decomposition into PKM and Ppd. The NiBr2 gKNB reproduction is strong and provides independent support for that part of the paper, and the qualitative prediction of a cycloidal ground state for NiBr2 is credible. However, the abstract's claim that the p-d hybridization and gKNB contributions in NiI2 are 'quantitatively evaluated' is not supported by the evidence presented. The concern is load-bearing because it directly targets the central quantitative result, but it does not invalidate the qualitative mechanism picture, so the reader's CONDITIONAL verdict remains appropriate. My proposed test would settle the issue by forcing the fit to use the independently computed gKNB amplitude and checking whether the p-d term is still needed and what its magnitude becomes. This is a concrete, feasible check that can be done with the data already in the paper.","tokens_in":14275,"tokens_out":7646,"duration_ms":75036,"concrete_test":"Fix A in Eq. 13 to the gKNB value obtained from the reported M matrix for NiI2 (A_gKNB = |P(M1)|/√3 at q = 0.2) and refit only B to the q-dependent DFT data. If P(q) = A_gKNB sin(2πq) + B sin(4πq) reproduces the DFT points within numerical error, the original two-parameter fit was over-parameterized and the extracted B is not unique; if it fails, the residual must be explained by additional q-dependence beyond the gKNB term. Report the residual and the implied B in either case.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim for NiI2 rests on Eq. 13, P_total(q) = A sin(2πq) + B sin(4πq), with A and B obtained by fitting to DFT. The problem is that A is not the value predicted by the gKNB model. From the M1 matrix in Eq. 10 and Eq. 5 for Q = (q,0,0), the gKNB amplitude is |P|/√3 ≈ 1.32×10^-13/1.73 ≈ 0.76×10^-13 C/m (or, if the quoted value is the [010] component rather than the magnitude, A ≈ 0.88×10^-13 C/m). The fitted A = 1.14×10^-13 C/m is 30–50% larger. Thus the 'PKM' term in Eq. 13 is not the gKNB contribution computed in the paper; it is a free parameter that absorbs part of what is later attributed to p-d hybridization. Likewise, B is not computed from the extended Arima model; only the sin(4πq) functional form is borrowed. The λSOC test in Eq. 14 uses the same fitted A and B, so it is a self-consistency check rather than an independent determination. The qualitative statement that a second harmonic in q is present is plausible, but the claimed quantitative separation of PKM and Ppd is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates the microscopic origin of magnetoelectric polarization in monolayer NiBr2 and NiI2 using first-principles DFT calculations. For NiBr2, the authors predict a cycloidal magnetic ground state, show that the total electric polarization is linear in the spin-orbit coupling strength λSOC, and demonstrate that the generalized Katsura-Nagaosa-Balatsky (gKNB) model, including first through third nearest-neighbor spin dimers, quantitatively reproduces the DFT polarization (Eq. 8). For NiI2, they adopt a proper-screw helical state with Q=(q,0,0), find that the gKNB model alone cannot describe the q- and λSOC-dependence of the polarization, and propose an extension of Arima's higher-order p-d hybridization mechanism. The total polarization is written as P_total(q) = A sin(2πq) + B sin(4πq) (Eq. 13), with A and B determined by fitting the DFT data, and the λSOC dependence is described by P'_total(λSOC) = (A sin 2πq)λSOC + (B sin 4πq)λSOC² (Eq. 14).","tokens_in":14605,"tokens_out":5647,"duration_ms":50695,"significance":"If the claims were fully established, the paper would provide a valuable validation of the gKNB model in a two-dimensional material and a quantitative decomposition of two spin-driven ferroelectric mechanisms. The NiBr2 result is a genuine quantitative success: the gKNB model with first through third nearest-neighbor dimers reproduces the DFT polarization within a few percent and captures the observed linear λSOC scaling. The NiI2 analysis, however, is not on the same footing. The central quantitative separation of the 'PKM' and 'p-d' contributions rests on a two-parameter fit whose amplitudes are not tied to the independently calculated gKNB matrix elements or to a microscopic derivation of the p-d term. The paper is therefore significant mainly for the NiBr2 case and as a qualitative suggestion that an additional q-dependent mechanism operates in NiI2.","major_comments":[{"comment":"The fitted amplitude A in Eq. (13) is inconsistent with the gKNB amplitude that the same paper calculates from the first-nearest-neighbor M matrix. From Eq. (10) and the expression P_i^tot = (√3/2 A, −3/2 A, 0) with A = (M11 − M22) sin(2πq), the reported value |P(M1)| ≈ 1.32×10^−13 C/m at q = 0.2 implies A_KNB ≈ 1.32×10^−13 / (√3 sin 72°) ≈ 0.80×10^−13 C/m. The fitted A = 1.14×10^−13 C/m is roughly 40% larger. Since the paper does not compute second- or third-neighbor M matrices for NiI2, the sin(2πq) term cannot be identified with the microscopic gKNB contribution; the fitted A absorbs part of whatever additional physics is present. Consequently, the paper's quantitative separation of PKM and Ppd is not microscopically anchored.","section":"§IIID, Eq. (13)"},{"comment":"The λSOC-dependence test in Eq. (14) uses the same fitted A and B obtained from the q-dependence fit in Eq. (13). The agreement with the DFT data in Fig. 4 is therefore a self-consistency check of the two-parameter representation, not an independent verification that the linear term is the gKNB/Kaplan-Mahanti mechanism and the quadratic term is p-d hybridization. In particular, any mechanism contributing a sin(2πq) shape would be absorbed into A, and any mechanism with a sin(4πq) shape would be absorbed into B, so the λSOC scaling test does not discriminate between microscopic origins.","section":"§IIID, Eq. (14)"},{"comment":"The extension of Arima's p-d hybridization model from a proper-screw state with Q=(q,q,0) to the present case Q=(q,0,0) is asserted rather than derived. Equation (12), P_pd = C sin(θ1 − θ2) sin(4πq), is introduced without a microscopic Hamiltonian or a symmetry analysis specific to the NiI2 structure, and the constant C is never computed or constrained. Thus the 'extended p-d hybridization mechanism' is not independently tested; only its functional form sin(4πq) is borrowed, and its amplitude B is fixed by fitting the residual DFT polarization. The existence and magnitude of the p-d contribution therefore remain unverified.","section":"§IIID, Eq. (12) and surrounding text"}],"minor_comments":[{"comment":"The caption refers to a 'red dashed green line' which appears to be a typo; presumably one line is red dashed and another is green.","section":"§IIID, Fig. 3 caption"},{"comment":"The term 'no-substitution method' used for the gKNB matrix element calculation is not defined or referenced; please provide a citation or a brief explanation.","section":"§II, 'no-substitution method'"},{"comment":"The word 'repectively' is misspelled; it should be 'respectively'.","section":"Table II caption"},{"comment":"In Eq. (9), the Taylor expansion expression is written with 'O(q)^4' but it should be 'O(q^4)' or 'O(q)^4' is unconventional; please clarify the order symbol.","section":"§IIIC, Eq. (9)"}],"recommendation":"major_revision","confidential_remarks":"The NiBr2 part of the paper is solid and could support a good publication on its own. The NiI2 analysis as written oversells the quantitative decomposition: the key amplitudes are fit parameters, and the fitted gKNB amplitude does not match the computed first-neighbor value. I would like to see either (i) a full gKNB calculation for NiI2 including the relevant higher-neighbor dimers, with a direct comparison to the fitted A, or (ii) a reframing of the NiI2 result as a phenomenological fit with an explicit statement that the p-d attribution is speculative. The extension of Arima's model also needs a proper symmetry-based justification for Q=(q,0,0) before the quantitative claim can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about arXiv:2501.05025. First, the NiBr2 half is a credible quantitative success: the gKNB M-matrices through third-neighbor dimers sum to (0.00, -2.76, 0.07) x 10^-13 C/m, matching the DFT value of about 2.67 x 10^-13 C/m, and the predicted cycloidal ground state is a concrete, testable claim. Second, the NiI2 half, which the abstract sells as a quantitative separation of mechanisms, is a two-parameter fit. Eq. 13, P(q)=A sin(2pi q)+B sin(4pi q), gets A and B from DFT data; the fitted A is 1.14 x 10^-13 C/m, while the gKNB M1 matrix implies A about 0.76-0.88 x 10^-13 C/m depending on how you interpret the [010] component. That is a 30-50% gap, so the 'PKM' curve in Fig. 3(a) is not the computed gKNB contribution. B is also fitted; the Arima-style derivation gives only the sin(4pi q) form, not the amplitude. The lambda_SOC curve from Eq. 14 uses those same fitted parameters, so it is a consistency check rather than an independent prediction.\n\nWhat is genuinely new: the extension of Arima's p-d hybridization mechanism from Q=(q,q,0) to the experimentally relevant Q=(q,0,0) proper-screw, and the finding that a quadratic-in-lambda_SOC component can serve as evidence for p-d hybridization. Those are useful ideas.\n\nThe soft spot is the gap between the abstract's claim of quantitative evaluation and what the fit actually establishes. The authors do not compute the p-d amplitude from a microscopic Hamiltonian, and they omit second- and third-neighbor dimers from the NiI2 gKNB baseline, which could shift the expected A. The qualitative story, that there is a sin(4pi q) component and a lambda_SOC^2 piece, is plausible and well supported by the two-harmonic fit. But the numbers are not microscopically determined.\n\nWho should read it: anyone working on NiI2 or NiBr2 monolayers, or on fit-based analyses of spin-driven polarization. It deserves a serious referee. My recommendation: send to review, but require a major revision that either computes the p-d amplitude independently, includes the higher-neighbor gKNB terms for NiI2, or explicitly reframes the NiI2 contributions as fit-based estimates rather than quantitatively evaluated mechanisms.","headline":"NiBr2 half is solid; the NiI2 'quantitative separation' is a two-parameter fit without microscopic anchors.","tokens_in":15226,"tokens_out":5228,"would_cite":true,"duration_ms":47115,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In a single layer of the magnet NiI2, the electric polarization induced by spiral spin order combines two distinct mechanisms that scale differently with spin-orbit coupling, while single-layer NiBr2 is described by one mechanism.","keywords":["type-II multiferroics","monolayer NiI2","monolayer NiBr2","spin-spiral ferroelectricity","gKNB model","p-d hybridization mechanism","spin-orbit coupling","first-principles calculations"],"falsifier":"Compute the polarization of monolayer NiI2 at fixed $\\boldsymbol{Q}=(0.2,0,0)$ for a series of spin-orbit strengths from $\\lambda_{\\mathrm{SOC}}=0$ to 2 and plot $P/\\lambda_{\\mathrm{SOC}}$ against $\\lambda_{\\mathrm{SOC}}$; if the ratio is flat within numerical error, the claimed quadratic p-d term is absent, and fitting $P(q)$ over many $q$ values should likewise return $B$ consistent with zero.","tokens_in":13976,"feed_emoji":"🧲","tokens_out":12731,"duration_ms":111609,"temperature":0.7,"pith_summary":"This paper tries to establish where the electric polarization in the monolayer magnets NiBr2 and NiI2 comes from. Using density-functional calculations, it predicts that NiBr2 has a cycloidal spin ground state whose polarization is linear in the spin-orbit coupling strength and is reproduced by the generalized Katsura-Nagaosa-Balatsky (gKNB) model once first-, second-, and third-neighbor spin dimers are included. For the proper-screw spiral of NiI2 (spins rotating in the plane perpendicular to the propagation vector), the paper finds that the gKNB term is not enough: both the wave-vector dependence and the spin-orbit dependence require a second contribution, which it attributes to an extended p-d hybridization mechanism. The central quantitative result is $P(q) = A\\sin(2\\pi q) + B\\sin(4\\pi q)$ for NiI2, where the second term is quadratic in spin-orbit coupling and is about 17% of the first at the experimental wave vector. If correct, this gives a way to separate spin-current and hybridization mechanisms in other triangular-lattice helimagnets.","feed_headline":"NiI2's polarization splits into two spin mechanisms","feed_subtitle":"Calculations split NiI2's polarization into a linear and a quadratic spin-orbit part; NiBr2 stays purely linear","key_machinery":"The central machinery is the gKNB magnetoelectric tensor, a $3\\times3$ matrix $M_{ij}$ assigned to each Ni-Ni dimer so that the pair contributes $M_{ij}(\\mathbf{S}_i\\times\\mathbf{S}_j)$ to the polarization; the paper computes these matrices from density-functional calculations for first-, second-, and third-neighbor dimers and applies the sum to the spiral spin configurations. The complementary piece is the p-d hybridization mechanism, in which the covalency of each metal-ligand bond is modulated by the two neighboring bonds sharing the same ligand; extended from the proper-screw case $\\boldsymbol{Q}=(q,q,0)$ to NiI2's $\\boldsymbol{Q}=(q,0,0)$, it gives $P_{pd}=C\\sin(\\theta_1-\\theta_2)\\sin(4\\pi q)$ on a Ni2I4 cluster. The two mechanisms are separated by their signatures: the gKNB term is proportional to $\\sin(2\\pi q)$ and linear in the spin-orbit coupling $\\lambda_{\\mathrm{SOC}}$, while the p-d term is proportional to $\\sin(4\\pi q)$ and quadratic in $\\lambda_{\\mathrm{SOC}}$.","core_discovery":"The paper's central claim is that the electric polarization of the NiI2 monolayer is a sum of two spin-orbit-mediated mechanisms with different wave-vector and spin-orbit dependence. In NiBr2, the cycloidal state's polarization is fully explained by the gKNB spin-current model once first-, second-, and third-neighbor Ni dimers are included; the sum reproduces the density-functional result as $(0.00, -2.76, 0.07)\\times10^{-13}\\ \\mathrm{C/m}$, and the polarization is strictly linear in $\\lambda_{\\mathrm{SOC}}$. In NiI2, a proper-screw state with $\\boldsymbol{Q}=(q,0,0)$ gives a total polarization that deviates from the gKNB-only prediction, so the authors extend the p-d hybridization mechanism to this propagation direction and obtain $P_{pd}=C\\sin(\\theta_1-\\theta_2)\\sin(4\\pi q)$. Combining the mechanisms yields $P_{\\mathrm{total}}(q)=A\\sin(2\\pi q)+B\\sin(4\\pi q)$ with $A=1.14\\times10^{-13}\\ \\mathrm{C/m}$ and $B=3.15\\times10^{-14}\\ \\mathrm{C/m}$, and a spin-orbit dependence $P'_{\\mathrm{total}}(\\lambda_{\\mathrm{SOC}})=A\\sin(2\\pi q)\\,\\lambda_{\\mathrm{SOC}}+B\\sin(4\\pi q)\\,\\lambda_{\\mathrm{SOC}}^2$ that the authors report reproduces the direct density-functional data.","pith_inferences":["Inference: because the amplitude $B$ of the p-d term is set by fitting the residual density-functional polarization rather than derived, the model's robust prediction is the functional form; computing $C$ independently from a cluster or Wannier model would test whether the fitted $B$ is physically reasonable.","Inference: the $\\sin(4\\pi q)$ form follows from three metal-ligand bonds sharing one ligand, so similar higher-harmonic polarization terms should appear in other triangular-lattice helimagnets whose proper-screw order runs along a nearest-neighbor direction.","Inference: a quadratic-in-spin-orbit polarization could in principle also be produced by higher-order spin interactions in an effective spin model; distinguishing that description from the p-d hybridization picture would require independent calculation of the four-spin exchange constants.","Inference: the paper's symmetry argument that anisotropic symmetric exchange cannot contribute to the polarization implies a sharp selection rule that could be tested by reversing the spiral helicity and checking that the polarization flips sign through the gKNB channel."],"forward_implications":["Measurements of electric polarization versus spiral wave vector in monolayer NiI2 should show a two-harmonic profile, $\\sin(2\\pi q)$ plus $\\sin(4\\pi q)$, rather than a single harmonic.","The quadratic-in-spin-orbit part of NiI2's polarization grows relative to the linear part as $q$ moves toward 0.25, so the balance between the two mechanisms is tunable through the magnetic period.","Because the halogen spin-orbit coupling dominates in both compounds, substituting Br for I changes not only the magnitude but also the ratio of the two NiI2-style contributions.","The same decomposition is claimed to carry over to bulk NiI2 and NiBr2, so a $\\sin(4\\pi q)$ component should appear in the polarization of bulk NiI2 as its spiral wave vector is varied."],"supporting_citations":[{"why":"Supplies the gKNB magnetoelectric-tensor method used to compute spin-dimer polarization.","marker":"[43]"},{"why":"Supplies the p-d hybridization mechanism whose proper-screw form is extended to the Q=(q,0,0) case.","marker":"[62]"},{"why":"Supplies the Kaplan-Mahanti contribution that produces the diagonal matrix elements and the linear spin-orbit term.","marker":"[60]"},{"why":"Supplies the original spin-current expression for the off-diagonal gKNB term.","marker":"[59]"},{"why":"Supplies the intrasite versus intersite classification showing the intrasite p-d term vanishes in NiX2.","marker":"[17]"},{"why":"Provides the experimental proper-screw wave vector Q=(0.2203,0,0) used to set up the NiI2 monolayer.","marker":"[24]"},{"why":"Provides experimental evidence of a multiferroic state in monolayer NiI2 that motivates the modeling.","marker":"[21]"},{"why":"Provides the bulk NiBr2 cycloidal multiferroic phase used as a comparison for the monolayer ground state.","marker":"[18]"},{"why":"Provides the frustration condition and the cycloidal-versus-proper-screw analysis used for NiBr2.","marker":"[19]"}],"fun_headline_variants":["NiI2 polarization splits into linear and quadratic spin-orbit parts","NiI2's polarization: two mechanisms, NiBr2 stays linear","Microscopic origin: NiI2 has two spin-orbit polarization paths","NiI2 polarization: linear plus quadratic spin-orbit, NiBr2 linear only"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the p-d hybridization formula extended to the $\\boldsymbol{Q}=(q,0,0)$ propagation direction is correct, since its amplitude is fixed by fitting the residual density-functional polarization rather than derived from a microscopic Hamiltonian.","fun_headline_variants_meta":{"raw":{"variants":["NiI2 polarization splits into linear and quadratic spin-orbit parts","NiI2's polarization: two mechanisms, NiBr2 stays linear","Microscopic origin: NiI2 has two spin-orbit polarization paths","NiI2 polarization: linear plus quadratic spin-orbit, NiBr2 linear only"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000734,"raw_usage":{"total_tokens":3379,"prompt_tokens":1141,"completion_tokens":2238,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":757,"completion_tokens_details":{"reasoning_tokens":2156}},"tokens_in":757,"tokens_out":2238,"duration_ms":14213,"temperature":1.0,"reasoning_tokens":2156,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:20:09.619185+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the polarization of monolayer NiI2 at fixed $\\boldsymbol{Q}=(0.2,0,0)$ for a series of spin-orbit strengths from $\\lambda_{\\mathrm{SOC}}=0$ to 2 and plot $P/\\lambda_{\\mathrm{SOC}}$ against $\\lambda_{\\mathrm{SOC}}$; if the ratio is flat within numerical error, the claimed quadratic p-d term is absent, and fitting $P(q)$ over many $q$ values should likewise return $B$ consistent with zero.","supporting_citations":[{"cited_title":"Xiang, E","cited_arxiv_id":null,"evidence_quote":"Supplies the gKNB magnetoelectric-tensor method used to compute spin-dimer polarization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the p-d hybridization mechanism whose proper-screw form is extended to the Q=(q,0,0) case."},{"cited_title":"Bolens, Theory of electronic magnetoelectric coupling in d5 mott insulators, Physical Review B 98, 125135 (2018)","cited_arxiv_id":null,"evidence_quote":"Supplies the Kaplan-Mahanti contribution that produces the diagonal matrix elements and the linear spin-orbit term."},{"cited_title":"Katsura, N","cited_arxiv_id":null,"evidence_quote":"Supplies the original spin-current expression for the off-diagonal gKNB term."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the intrasite versus intersite classification showing the intrasite p-d term vanishes in NiX2."},{"cited_title":"Layer thickness crossover of type-II multiferroic magnetism in NiI2","cited_arxiv_id":"2307.10686","evidence_quote":"Provides the experimental proper-screw wave vector Q=(0.2203,0,0) used to set up the NiI2 monolayer."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides experimental evidence of a multiferroic state in monolayer NiI2 that motivates the modeling."},{"cited_title":"Tokunaga, D","cited_arxiv_id":null,"evidence_quote":"Provides the bulk NiBr2 cycloidal multiferroic phase used as a comparison for the monolayer ground state."},{"cited_title":"In-plane strain tuning multiferroicity in monolayer van der Waals NiI$_{2}$","cited_arxiv_id":"2209.12392","evidence_quote":"Provides the frustration condition and the cycloidal-versus-proper-screw analysis used for NiBr2."}],"review_version":1}