{"id":"168d3a73-55a1-4ea8-95c7-a1dc5f1172c6","arxiv_id":"2607.13053","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A complete, orthonormal multipole basis including inter-shell operators is constructed for the local density matrix and applied to show that the electric toroidal monopole moment of tellurium flips sign between the two enantiomers.","lead":"This theory paper builds a complete set of local electron \"pattern rulers\" — multipole operators — for atoms in crystals, including operators that mix different electron shells that earlier methods left out. The authors use the new tool to compute the electric toroidal monopole, a chiral scalar that changes sign between the left- and right-handed forms of tellurium, proposing it as a local readout of chirality.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"ETM moment is a projection-dependent quantity; a single untested DFT value does not establish it as a robust local chirality measure.","rationale":"The reader's weakest_assumption is exactly the projection truncation dependence of the ETM value. This is the load-bearing concern because the mathematical completeness proof is internal and correct; the physical significance of the application depends on the stability of a single computed number. The paper itself states that only one radial projector per l-shell is used, and the radial dependence would enter the formalism as a 'scalar factor' — that factor is precisely what makes the ETM sensitive to the projector choice. No convergence tests are reported, and the data are not yet released. Therefore the conditional verdict is appropriate. The sign inconsistency in eqs. (20)–(23) is a real but minor normalization issue that does not affect the completeness argument or the symmetry-enforced sign change. My concrete test would settle whether the ETM magnitude is meaningful or an artifact.","tokens_in":23855,"tokens_out":18916,"duration_ms":162372,"concrete_test":"Recompute w^{ν=0,110}_{s,p,0} for both enantiomers of tellurium using the same VASP settings but (i) a PAW projector radius increased by 20% and (ii) l_max=2 (including d states in the local subspace). If in either case the sign flips or the magnitude changes by more than a factor of 2, the ETM moment is not a robust chirality descriptor and the application claim is unsupported. If the sign remains opposite and the magnitude varies by <10%, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim has two parts: (1) a complete orthogonal real Hermitian multipole basis spanning the truncated local operator space, and (2) the electric toroidal monopole (ETM) moment computed for tellurium as a local chirality descriptor. The mathematical construction in Sec. II and App. A is solid: the dimension counting and 3j orthogonality relations convincingly demonstrate completeness within the finite subspace defined by one radial projector per l-shell. The weakest load-bearing point is the physical interpretation of the ETM moment. In the DFT application, the density matrix is projected onto PAW projectors with one radial function per l-shell and l_max=1 (s and p only). The ETM operator W^{ν=0,110}_{s,p,0} is an inter-shell operator whose expectation value involves the s–p block of the projected density matrix. That block is proportional to the radial overlap <R_s|r|R_p>, which depends on the choice of projector radius and the PAW method. The paper reports w = ±2.6×10^{-4} without any test against projector radius, k-mesh, cutoff, functional, or inclusion of higher l states. If the magnitude or even the sign changes under such variations, the ETM is an artifact of the truncation, not a robust measure of chirality. The Data Availability statement says files will be public only upon publication, so the result is not currently independently reproducible.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a complete orthogonal real Hermitian multipole basis for the local single-site density matrix, extending the fixed-shell multipole formulations of van der Laan/Thole and Bultmark et al. to include inter-shell operators. The construction proceeds by decomposing the orbital and spin operator spaces under SO(3), coupling the resulting spherical tensors, and then forming real Hermitian combinations. The basis is classified by spatial parity and time-reversal symmetry. Orthogonality and completeness are addressed in Appendix A4, including a dimension count of 4(l_max+1)^4 for the real Hermitian operator space. As an application, the authors compute the electric toroidal monopole (ETM) moment w^{ν=0,110}_{s,p,0} for the two enantiomers of trigonal tellurium from DFT density matrices and report equal magnitudes and opposite signs on the two enantiomers.","tokens_in":24033,"tokens_out":24829,"duration_ms":244833,"significance":"The mathematical construction is a useful and mostly self-contained extension of the standard multipole basis: it supplies an explicit, symmetry-adapted, complete basis for the full local Hermitian operator space, including inter-shell blocks, and reduces to the known fixed-shell formalism when l'=l. The proofs in Appendix A are explicit, and the dimension-counting argument is convincing. The application to tellurium is conceptually attractive, but the reported ETM value is a single DFT number with no robustness checks. Because the sign change between enantiomers is enforced by spatial parity, the physically informative content of Table II is the magnitude and its stability; the manuscript currently does not establish that stability.","major_comments":[{"comment":"The main application is not yet supported by sufficient evidence. The ETM moment is an expectation value of an inter-shell operator in a local subspace defined by one radial PAW projector per angular-momentum shell and l_max=1 (Sec. II). Its magnitude depends on the s–p block of the projected density matrix and on the radial matrix element between the s and p projector functions. No tests are reported against the PAW sphere/projector radius, augmentation settings, energy cutoff, k-mesh, exchange-correlation functional, or inclusion of higher-l states. Since the sign change between enantiomers is fixed by parity, a single value ±2.6×10^-4 does not by itself demonstrate that the ETM is a robust chirality descriptor. Please add convergence/robustness tests or explicitly reframe the Te calculation as an illustration of the formalism rather than a validation of the ETM as a measure.","section":"Sec. IV.B / Table II"},{"comment":"The concluding claim that the ETM 'acts as an atomic-scale measure for structural chirality' overreaches the presented evidence. It is a symmetry statement that a parity-odd time-reversal-even scalar changes sign under exchange of enantiomorphic structures; the nontrivial content of the calculation is the nonzero amplitude and its stability under the local-truncation and PAW choices, neither of which is tested. The conclusion should be tempered pending the robustness analysis described above.","section":"Sec. V / Abstract"}],"minor_comments":[{"comment":"Please specify the PAW projector/augmentation-sphere parameters and the precise way the one-center density matrix is extracted from VASP. This is needed both for reproducibility and for assessing the sensitivity of the inter-shell moments.","section":"Sec. IV.A"},{"comment":"The notation l1,l2 in W^{νkpr,real}_{l1,l2,t} after defining l'=min and l=max is a little confusing; state explicitly that the operator is associated with the (l',l) block with l'≤l.","section":"Sec. II / Eq. (60)"},{"comment":"The table entries '≠0' and '0' should be read as 'for the given ν' and 'for the complementary ν'; a sentence clarifying that the two rows are complementary rather than contradictory would help avoid misinterpretation.","section":"Sec. III / Table I"},{"comment":"The caption mentions blue inward-pointing and red outward-pointing hedgehogs, but the figure shows arrows around one Te atom; clarify the relation between arrow direction and the sign of w^{ν=0,110}_{s,p,0}.","section":"Fig. 1"},{"comment":"The statement that input files and data 'will be publicly available upon publication' means the current result cannot be independently reproduced. Consider providing the relevant computational parameters and data in the supplementary material or a repository at the time of submission.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The mathematical construction is solid and likely publishable, but the physical application needs more work before the ETM claim can be accepted. I would recommend major revision: the central derivation is defensible, but the application section requires either robustness/convergence tests or a clearly reduced claim. The paper may be a good fit for the journal once the application is strengthened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: the mathematics is the real contribution; the tellurium calculation is a plausible demo but not yet a validated quantitative result.\n\nWhat's new: a complete orthogonal real Hermitian multipole basis for the local single-particle operator space, extending the fixed-shell formalism to inter-shell (l≠l') operators. The construction is careful — they separate orbital and spin spaces, use Wigner–Eckart with the van der Laan normalization, couple spin and orbital parts, then build real Hermitian combinations. I checked the orthogonality and the dimension count N=4(l_max+1)^4, which matches dim End(H). The parity and time-reversal selection table is useful and follows cleanly. This fills a real gap, since standard tools like multipyles only do fixed-shell.\n\nThe paper does well: the derivation is explicit and self-contained, no fitting or external parameters, and the reduction to the old formalism for l'=l is checked. The appendices give enough detail to re-derive from scratch. The Te application demonstrates that the ETM picks out chirality: the sign flips between enantiomers, which is symmetry-enforced, so it's a valid consistency check.\n\nWhere it's soft: the application is under-supported. One DFT calculation with a single k-mesh, cutoff, functional, and one projector per l-shell (l_max=1) gives w=±2.6×10^-4, with no test of how this depends on PAW projector radius or inclusion of d states. The operator involves the s–p overlap <R_s|r|R_p>, which is radial-dependent, so the magnitude is not obviously robust. The sign flip is guaranteed by symmetry; the quantitative claim as a 'local measure of chirality' needs convergence checks or a second chirality observable. This is the weakest load-bearing point. Also, the data availability statement (upon publication) means the numbers are not independently checkable now.\n\nMinor but real: eqs. (20)–(22) identify n_lkl' with a stretched 3j symbol, but that symbol can be negative (e.g. (1 1 1; -1 0 1) = -1/√6) while the closed-form square root is positive. The normalization in (23) might still work if n is allowed to be signed, but the equations as printed are inconsistent. This should be fixed in revision; it's the kind of thing a careful reader will trip on.\n\nWho this is for: people using multipole decompositions in DFT, especially toroidal/chiral observables and inter-shell physics. It deserves a serious referee. I'd send it to review with a request to strengthen the application (convergence tests or softened claims) and fix the sign.","headline":"Solid mathematical extension of the multipole formalism; the tellurium demo is thin and needs convergence work, but the core deserves refereeing.","tokens_in":24694,"tokens_out":4759,"would_cite":true,"duration_ms":51286,"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":"This paper constructs a complete orthogonal basis of real Hermitian multipoles for the local single-particle density matrix, extends it to inter-shell operators, and uses the resulting electric toroidal monopole to distinguish the two enant","keywords":["multipole decomposition","local density matrix","inter-shell operators","time-reversal symmetry","electric toroidal monopole","chirality","tellurium","spherical tensor operators"],"falsifier":"Compute the electric toroidal monopole moment in trigonal tellurium with different projector radii (e.g., 1.0 to 1.5 bohr), denser k-meshes, higher energy cutoffs, and with d states included; if the sign difference between the two enantiomers disappears or flips under any of these changes, the claim that the ETM provides a robust local measure of chirality fails.","tokens_in":23585,"feed_emoji":"🌀","tokens_out":5230,"duration_ms":45259,"temperature":0.7,"pith_summary":"Fixed-shell multipole decompositions of the local density matrix only resolve operators within a single orbital l; the authors show that adding inter-shell operators — those mapping between different l — completes the basis of all Hermitian single-particle operators on the site. They build the complete orthogonal real Hermitian basis from coupled orbital and spin tensors, classify every operator by spatial parity and time-reversal, and prove the dimension count matches the full operator space. In the new classification the lowest odd-parity, time-reversal-even pseudoscalar is the electric toroidal monopole, a local chiral scalar. Applying the decomposition to trigonal tellurium, they compute this monopole moment at each Te site and find it is equal in magnitude and opposite in sign in the two enantiomers, making it a candidate atomic-scale measure of chirality. The reader should care because the completed basis gives a symmetry-resolved language in which any local observable can be matched to exactly the multipole moments that can contribute to it.","feed_headline":"Toroidal monopole flips sign between tellurium enantiomers","feed_subtitle":"New multipole basis gives atoms a local handedness score that flips between mirror-image crystals.","key_machinery":"The machinery is the coupled tensor operator W^{kpr}_t(l',l) = n^{-1}_{kpr} (-1)^{k+p} Σ_{q,y} (-1)^{-q-y} [k r p; -q t-y] C^k_q(l',l) ⊗ σ^p_y, built from the orbital tensor C^k_q(l',l) (normalized so the stretched matrix element is 1) and the spin tensor σ^p_y (identity and Pauli matrices in spherical form). Hermitian adjoints relate opposite (l',l) blocks, so the authors form real components and then '+/–' Hermitian combinations; time-reversal projection selects definite parity ν and the notation W^{νkpr,real}_{l1,l2,t} with l1≤l2. Completeness is checked by counting: the number of operators equals 4(l_max+1)^4, the dimension of End(H), and orthogonality follows from the orthogonality of W","core_discovery":"The central claim is that the time-reversal even/odd real Hermitian multipole operators W^{νkpr,real}_{l1,l2,t} form a complete orthogonal real basis of the Hermitian local single-particle operator space, extending the van der Laan–Thole/Bultmark fixed-shell basis to inter-shell blocks (l1 ≠ l2). The construction starts from the tensor product structure End(H) ≅ (⊕ V_l ⊗ V*_{l'}) ⊗ V_s ⊗ V*_s, decomposes orbital and spin operator spaces separately under SO(3), couples them with Clebsch–Gordan coefficients, and then forms real Hermitian combinations. The authors also prove that with one radial projector per angular-momentum shell and truncation at l_max, the counting of operators matches dim","pith_inferences":["If the ETM sign is robust to the projection details, it could serve as a first-principles descriptor for structural chirality in other helical or screw-axis crystals, and possibly for chirality-dependent transport such as the current-induced magnetoelectric response in tellurium.","The basis completeness argument suggests a natural test: recomputing the same Te ETM with d and f shells included should preserve sign and change magnitude; if the s–p-only result vanished at higher l_max, the reported descriptor would be an artifact of truncation rather than a stable chiral invariant.","Because the ETM is the lowest-rank parity-odd time-reversal-even pseudoscalar, one might expect it to correlate with the sign of natural optical activity or circular dichroism in chiral crystals; computing the rotational-strength-weighted spectra for the two enantiomers would give a direct observable test."],"forward_implications":["Any symmetry-resolved observable O with definite parity and time-reversal character can only couple to multipole moments of the matching class, so the full basis lets one read off which local degrees of freedom a given measurement or response probes.","Fixed-shell multipole analyses (van der Laan–Thole, Bultmark et al.) are recovered exactly for l'=l; inter-shell operators such as the electric toroidal monopole are the new, previously missing channels.","The electric toroidal monopole moment computed from the s–p block of the density matrix distinguishes the two enantiomers of trigonal tellurium by sign while preserving equal magnitude at all three Te sites.","The same decomposition can be applied to any site-projected DFT density matrix, turning the local density matrix into a complete set of symmetry-labeled amplitudes comparable across materials."],"fun_headline_variants":["Multipole basis reveals handedness in tellurium crystals","Complete multipole basis exposes toroidal monopole in chiral Te","Inter-shell multipoles enable toroidal monopole measurement","Time-reversal symmetry classifies multipoles, finds Te handedness","Electric toroidal monopole distinguishes tellurium enantiomers"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the projected local Hilbert space used in the application — one radial projector per angular-momentum shell and only s and p orbitals — captures a physically meaningful electric toroidal monopole moment; if the projected density matrix were strongly basis-dependent, the reported ±2.6×10^{-4} values would be artifacts of the projection rather than a robust chirality descriptor.","fun_headline_variants_meta":{"raw":{"variants":["Multipole basis reveals handedness in tellurium crystals","Complete multipole basis exposes toroidal monopole in chiral Te","Inter-shell multipoles enable toroidal monopole measurement","Time-reversal symmetry classifies multipoles, finds Te handedness","Electric toroidal monopole distinguishes tellurium enantiomers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000209,"raw_usage":{"total_tokens":1249,"prompt_tokens":751,"completion_tokens":498,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":495,"completion_tokens_details":{"reasoning_tokens":414}},"tokens_in":495,"tokens_out":498,"duration_ms":4613,"temperature":1.0,"reasoning_tokens":414,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T08:52:59.829345+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the electric toroidal monopole moment in trigonal tellurium with different projector radii (e.g., 1.0 to 1.5 bohr), denser k-meshes, higher energy cutoffs, and with d states included; if the sign difference between the two enantiomers disappears or flips under any of these changes, the claim that the ETM provides a robust local measure of chirality fails.","supporting_citations":[],"review_version":1}