{"id":"042e0c80-3f77-441a-b713-e4f714ce66ff","arxiv_id":"2411.09876","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For 3C-SiC, diagonalizing the atomistic dielectric response reveals hidden nonplanar optical eigenwaves and an effective screening picture yielding the local-field-corrected macroscopic dielectric constant.","lead":"Using a quantum-mechanical atomistic model, the authors compute deeply subwavelength optical eigenmodes inside a silicon carbide crystal and visualize them in real space. They find hidden, nonplanar optical waves that classical refractive index models miss, and use these modes to re-derive the crystal's local-field-corrected dielectric response.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central quantitative claims—eigenvalue spectrum, mode weights, and the local-field-induced lowering of epsilon_eff—rest on an unvalidated empirical-pseudopotential RPA dielectric matrix that disagrees with the measured dielectric constant; no convergence study is supplied, so the core numbers…","rationale":"The reader identified the unvalidated dielectric-matrix calculation as the weakest assumption, and I agree. The central claim is not just that eigenmodes exist—it is that specific eigenvalues, weights, and the resulting effective dielectric constant (ε_eff = 4.908) reveal a field-dependent refractive index. Those numbers are direct outputs of the RPA matrix, and the large discrepancy with the measured dielectric constant of 3C-SiC is a concrete, falsifiable warning that the input is not reliable. Without a convergence study or an independent benchmark, the quantitative conclusions in Fig. 4 and the non-uniqueness claim are conditional. The longitudinal-only scope is a real limitation but is explicitly acknowledged by the authors; the missing numerical validation is the more load-bearing issue. The reader's CONDITIONAL verdict is appropriate, so I recommend no change.","tokens_in":9897,"tokens_out":9833,"duration_ms":110962,"concrete_test":"Recompute the eigen-decomposition using a converged first-principles plane-wave RPA dielectric matrix with the same crystal structure and q = 2π/a (0.1,0,0): increase the k-mesh to at least 25×25×25, increase the plane-wave cutoff until ε00 converges to within 1%, include at least 40 conduction bands, and compare the resulting ε_eff and the weights w0, w4 with the paper's values. Also compare ε00(q→0) with the measured high-frequency dielectric constant of 3C-SiC (≈6.5). If the converged ε_eff differs from 4.908 by more than ~10%, or the mode weights shift materially, the quantitative central claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All of the paper's quantitative results—the eigenmode eigenvalues πλ, the weights wλ, and the central value ε_eff = 4.908 (vs. ε00 = 5.435)—are derived from the RPA dielectric matrix of Eq. S1. That matrix is computed with a 51×51 reciprocal-space cutoff, a 15×15×15 k-mesh, 4 valence and 20 conduction bands, and empirical pseudopotentials, and no convergence study is reported. The most direct available benchmark, the measured high-frequency electronic dielectric constant of 3C-SiC, is about 6.5; the paper's ε00 = 5.435 is ~17% lower and ε_eff = 4.908 is ~25% lower. This is a serious red flag: if the dielectric matrix systematically underestimates screening, then the eigenvalues, the mode weights, and the claimed local-field-induced lowering of the dielectric constant could all be artifacts of the unconverged or inaccurate calculation. The qualitative existence of higher-order nonplanar eigenmodes may survive, but the specific quantitative claims in Fig. 4 and the field-dependence conclusion require accurate eigenvalues and weights. The longitudinal-only treatment is a further scope limitation: the claim that this constitutes a 'quantum generalization of the refractive index' needs justification for transverse optical fields, since the calculation treats only longitudinal response.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper introduces a 'deep microscopic optical band structure' for crystalline solids by eigendecomposing the RPA dielectric matrix in reciprocal space, and applies it to 3C-SiC. The dominant eigenmode is shown to be approximately planewave-like while higher-order eigenmodes exhibit rapidly varying polarization textures. The paper derives an effective screening strength as a weighted average of the eigenvalues, and reports that for a macroscopic planewave the effective dielectric constant is ε_eff = 4.908, lower than the head element ε_00 = 5.435. The authors interpret this as evidence that the refractive index is not a unique material property but depends on the spatial structure of the optical field.","tokens_in":1305,"tokens_out":1539,"duration_ms":62191,"significance":"The conceptual framework of decomposing the dielectric response into microscopic eigenwaves is attractive and may offer new insight into nonlocal and local-field effects in crystals. The visualizations of the eigenmodes are evocative. However, the central quantitative claims are not supported by a convergence study or an independent benchmark, and the connection to the transverse refractive index is not established. The qualitative existence of higher-order eigenmodes is robust because any non-diagonal dielectric matrix has multiple eigenmodes, but the specific eigenvalues, weights, and the reported lowering of the dielectric constant require validation. The paper is a promising application of an existing formalism (ref. 26) rather than a fully validated new prediction.","major_comments":[{"comment":"The quantitative results—the eigenvalues πλ, the weights wλ, and the central values ε_eff = 4.908 and ε_00 = 5.435—are computed from a single RPA dielectric matrix with a 51×51 reciprocal-space cutoff, a 15×15×15 k-mesh, 4 valence and 20 conduction bands, and empirical pseudopotentials. No convergence study is reported, and there is no comparison to the measured high-frequency electronic dielectric constant of 3C-SiC (~6.5) or to converged first-principles results. Since the claimed local-field-induced lowering of the dielectric constant depends on the difference between ε_eff and ε_00, the accuracy of the underlying matrix is load-bearing. The authors should provide convergence tests in both the reciprocal-space cutoff and the k-mesh, and benchmark ε_00 against experiment or established DFT calculations.","section":"S1 (Eq. S1) and Fig. 4"},{"comment":"The eigen-decomposition in Eq. (2) is performed on the longitudinal dielectric matrix, and the paper explicitly limits the treatment to longitudinal fields. Yet the abstract and conclusions claim a 'quantum generalization of refractive index', which is a property of transverse optical fields. The paper does not justify the connection between the longitudinal eigenmodes and the transverse refractive index; the standard equivalence in the long-wavelength limit for cubic crystals is not discussed. Unless this connection is established, the claim should be reframed as a generalization of longitudinal dielectric response, or the transverse treatment should be included.","section":"Sec. II and Sec. V"},{"comment":"The effective screening formula ε_eff^(-1) = Σ_λ w_λ π_λ^(-1) for a macroscopic planewave is exactly (ε^(-1))_00, the head of the inverse dielectric matrix, which is the standard definition of the macroscopic dielectric function including local-field effects. The paper should state this relation explicitly, because otherwise the derivation appears circular: the eigenvalues and weights are constructed from the same matrix whose macroscopic response they are supposed to explain. With that clarification, the mathematical content is sound, but the predictive claim that 'the classical refractive index is not a unique quantity' (Sec. IV) is not demonstrated: no example of a structured field (e.g., a point source or a tightly focused beam) is given in which the effective screening is actually lowered, nor is the coupling of such fields to the high-order eigenmodes computed. Without such an example, this conclusion is speculative.","section":"Sec. IV, Eq. (4)"}],"minor_comments":[{"comment":"The notation n_hat(q+g) is used in Sec. III without definition; it is defined later in Eq. (S3) of the SI, but should be defined at first use in the main text.","section":"Sec. III"},{"comment":"The figures lack scale bars and clear units; the 'a.u.' intensity scale is not defined.","section":"Fig. 3 and SI figures"},{"comment":"The SI text contains placeholder citation markers '[? ? ?]' in the reference to Adler-Wiser theory and to the symmetrized Coulomb kernel; these need to be replaced with proper references.","section":"SI S1"},{"comment":"The paper states that the eigenwaves with λ = 2 and 3 are degenerate and related by a four-fold improper rotation, but this is not shown explicitly; a symmetry analysis or a statement about the irreducible representations would strengthen the claim.","section":"Sec. III"}],"recommendation":"major_revision","confidential_remarks":"The paper overlaps with the authors' own recent preprint (ref. 26), which introduces the theoretical framework. The current manuscript is an application and extension to SiC, but the lack of validation of the quantitative numbers is the main concern. The absence of a reproducibility statement and the placeholder references in the SI suggest the manuscript is not yet in a polished state. I would not recommend acceptance until the convergence study and the benchmark are provided."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful thing here is the material-specific application: the first computed deep microscopic optical band structure and eigenwave visualizations for a real crystal, 3C-SiC. The pictures are genuinely helpful for seeing what 'nonlocal response at the lattice level' means, and the eigenwave decomposition of the macroscopic dielectric constant into mode weights is a clean way to present local-field effects. The math is straightforward Hermitian linear algebra and it is executed correctly.\n\nThe soft spots are real, and they sit exactly where the reader and the stress-test note put them. The central numbers—the eigenvalues, the weights, and especially epsilon_eff = 4.908 versus epsilon_00 = 5.435—all come from an empirical-pseudopotential RPA dielectric matrix computed with a 51x51 cutoff and a 15x15x15 k-mesh, and no convergence study is reported. The measured electronic dielectric constant of 3C-SiC is around 6.5; getting 5.435 at the head of the matrix is a ~17% miss, and the headline 4.908 is ~25% low. That does not automatically sink the paper, but it means the quantitative results are not currently validated against any independent benchmark, and the claimed local-field lowering could be partly an artifact of an unconverged or inaccurate calculation. The authors should be asked to show convergence with respect to matrix size and k-mesh, and to compare their epsilon_00 to experiment or to a well-tested DFT calculation.\n\nThe other issue is scope: the calculation treats only the longitudinal response, and the title 'quantum generalization of refractive index' is too strong for that. The paper does acknowledge this in the outlook, but the framing in the abstract and introduction overreaches. The effective screening is also, as the reader notes, a weighted average of eigenvalues of the very same matrix that defines the macroscopic element, so the specific value 4.908 is a basis re-expression rather than an independent prediction. That is fine as an interpretation, but it should not be sold as a new physical effect without more validation.\n\nWho gets value from this? People working on local-field effects in semiconductors, EELS interpretation, and nonlocal nanophotonics will appreciate the visualization and the modal picture. But I would not cite the quantitative numbers until the convergence and benchmark questions are answered.\n\nRecommendation: send it to peer review, but insist on a major revision: a convergence study, a direct comparison of the computed macroscopic dielectric constant with measured or high-level theoretical values, and a more careful statement about the longitudinal-only regime. The qualitative framework is sound enough to deserve referee time.","headline":"A well-illustrated application of the authors' own eigen-decomposition framework to 3C-SiC, but the quantitative claims are not yet backed by converged or benchmarked calculations.","tokens_in":10704,"tokens_out":1391,"would_cite":false,"duration_ms":16822,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["78.20.Ci","71.45.Gm"],"model":"deepseek-v4-flash","headline":"Inside SiC, the refractive index stops being a single number.","keywords":["refractive index","microscopic optical band structure","dielectric matrix","nonlocality","silicon carbide","optical eigenmodes","local-field effects","atomistic electrodynamics"],"falsifier":"Repeat the dielectric-matrix calculation with a larger reciprocal-space cutoff and denser k-mesh (for instance a 101×101 matrix and a 25×25×25 k-mesh) and compare the resulting mode-weighted ε_eff to the experimental electronic dielectric constant of 3C-SiC; if the converged value moves substantially away from 4.908 toward the measured value, the quantitative claim that local fields lower the dielectric constant to this degree is called into question, although the qualitative existence of nonplanar higher-order eigenwaves could still survive.","tokens_in":9671,"feed_emoji":"💎","tokens_out":5048,"duration_ms":48473,"temperature":0.7,"pith_summary":"The paper argues that the classical refractive index is only an approximation that holds for uniform planewaves in a homogeneous medium. At the atomistic level, a crystal's optical response is a nonlocal, multi-modal object: diagonalizing the inverse dielectric matrix yields a deep microscopic optical band structure with many eigenwaves, each carrying its own screening strength. For silicon carbide, the authors compute these eigenwaves and show that the dominant one resembles a planewave, while higher-order ones are nonplanar and inhomogeneous even near the optical limit. They reconstruct the macroscopic dielectric constant as a mode-weighted average of screening strengths, obtaining ε_eff = 4.908, and conclude that the refractive index of a crystal depends on the spatial structure of the field inside it.","feed_headline":"Inside SiC, the refractive index stops being a single number","feed_subtitle":"A quantum decomposition of the crystal's dielectric response reveals hidden optical modes that a planewave never excites.","key_machinery":"The central object is the deep microscopic optical band structure, obtained by the eigen-decomposition of the dielectric matrix: ε̄(q,ω) = Σ_λ |E_λ⟩ π_λ ⟨E_λ|. Here π_λ(q,ω) is the microscopic optical index (whose screening strength is 1 − $π_λ^{{-1}}$) and |E_λ⟩ is the microscopic optical eigenwave. Below the bandgap the dielectric matrix is Hermitian, so the eigenwaves form an orthonormal basis; this lets an arbitrary external field be decomposed into these eigenwaves, and the effective screening strength follows as $ε_eff^{{-1}}$ = Σ_λ w_λ $π_λ^{{-1}}$ with weights w_λ = Σ_g |E_{λg}† E_g|^2. This decomposition is what turns the refractive index from a scalar into a mode-dependent quantity.","core_discovery":"The central discovery is that the dielectric response of a crystal decomposes into a discrete set of microscopic optical eigenwaves, and that the effective screening felt by any field is a weighted average over these eigenwaves, with weights determined by the field's overlap with each mode. In SiC, the eigenvalue decomposition of the Hermitian dielectric matrix at q = 2π/a (0.1,0,0) and ω = 10 THz yields a dominant eigenwave (λ = 0) that is nearly a planewave, plus higher-order eigenwaves (λ = 1,2,3,4) with rapidly varying polarization textures at the lattice level. The macroscopic dielectric constant ε_eff = 4.908 is obtained from the eigenmode overlaps and screening strengths 1 − $π_λ^{{-1}}$, and it is lower than the value ε_00 = 5.435 that ignores local-field effects, consistent with the known tendency of local fields to lower the dielectric constant. The paper therefore claims that the refractive index is not a uniquely determined material constant but depends on the structure of the light inside the medium.","pith_inferences":["If the mode-weighted picture is correct, far-field measurements (which approximate planewaves) and near-field probes (electron beams, optical antennas) will systematically report different effective refractive indices for the same crystal; this could be tested by comparing electron energy loss spectra with ellipsometry.","The nonplanar higher-order eigenwaves imply that a point emitter inside a crystal radiates into a structured near field that no local dielectric constant can describe, which would affect predicted Purcell factors and radiation patterns of quantum emitters in SiC.","One testable extension is to compute the same band structure with a larger reciprocal-space cutoff and denser k-mesh and compare the resulting ε_eff to the measured electronic dielectric constant of 3C-SiC; the quantitative claim about the lowering of the dielectric constant stands only if the converged value remains close to the mode-weighted result."],"forward_implications":["Near the Γ-point, higher-order eigenwaves remain nonplanar and inhomogeneous, so atomistic nonlocal optical response persists even in the small-momentum limit.","Structured or localized light that couples strongly to high-order eigenwaves will experience a lower effective screening strength than a macroscopic planewave.","The paper predicts hidden optical eigenwaves that do not couple to macroscopic planewaves, and suggests they may carry topological optical invariants.","The macroscopic dielectric constant is lowered by local-field effects: for SiC at 10 THz, the mode-weighted value is 4.908 versus 5.435 without local fields.","The framework is intended to extend beyond SiC to other materials, including two-dimensional materials and transverse excitations."],"supporting_citations":[{"why":"Introduces the deep microscopic optical band structure and the eigen-decomposition of the dielectric matrix that the present paper applies to SiC.","marker":"[26]"},{"why":"Provides the quantum theory of the dielectric constant in real solids that underlies the random-phase approximation used for the dielectric matrix.","marker":"[21]"},{"why":"Establishes the dielectric constant with local-field effects, the framework the paper extends to a full eigenmode decomposition.","marker":"[22]"},{"why":"Supplies the symmetry-based indicators used to argue that the dielectric matrix and its eigenwaves carry crystal-symmetry information.","marker":"[30]"},{"why":"Gives the first-principles dielectric-matrix eigen-equation that the paper recasts as its microscopic optical band structure.","marker":"[31]"},{"why":"Motivates the multi-modal, multiple-scattering picture of the optical response of an atomic crystal, which the paper realizes through eigen-decomposition.","marker":"[25]"}],"fun_headline_variants":["Refractive index is a set of eigenwaves, not a single number","Hidden lattice-scale optical waves govern crystal light bending","Crystal optics demystified: nonlocal eigenmodes decide refraction","SiC reveals picoscale waves that redefine macroscopic dielectric"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative results rest on the accuracy of the empirical-pseudopotential RPA dielectric matrix computed on a 51×51 reciprocal-space grid with a 15×15×15 k-mesh, and the paper reports no convergence study and does not benchmark its ε_eff = 4.908 (or ε_00 = 5.435) against the measured electronic dielectric constant of 3C-SiC (about 6.5).","fun_headline_variants_meta":{"raw":{"variants":["Refractive index is a set of eigenwaves, not a single number","Hidden lattice-scale optical waves govern crystal light bending","Crystal optics demystified: nonlocal eigenmodes decide refraction","SiC reveals picoscale waves that redefine macroscopic dielectric"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1377,"prompt_tokens":933,"completion_tokens":444,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":549,"completion_tokens_details":{"reasoning_tokens":374}},"tokens_in":549,"tokens_out":444,"duration_ms":5214,"temperature":1.0,"reasoning_tokens":374,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:12:20.194140+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the dielectric-matrix calculation with a larger reciprocal-space cutoff and denser k-mesh (for instance a 101×101 matrix and a 25×25×25 k-mesh) and compare the resulting mode-weighted ε_eff to the experimental electronic dielectric constant of 3C-SiC; if the converged value moves substantially away from 4.908 toward the measured value, the quantitative claim that local fields lower the dielectric constant to this degree is called into question, although the qualitative existence of nonplanar higher-order eigenwaves could still survive.","supporting_citations":[{"cited_title":"Unraveling Optical Polarization at Deep Microscopic Scales in Crystalline Materials","cited_arxiv_id":"2407.15189","evidence_quote":"Introduces the deep microscopic optical band structure and the eigen-decomposition of the dielectric matrix that the present paper applies to SiC."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the quantum theory of the dielectric constant in real solids that underlies the random-phase approximation used for the dielectric matrix."},{"cited_title":"Wiser, Dielectric constant with local field effects included, Phys","cited_arxiv_id":null,"evidence_quote":"Establishes the dielectric constant with local-field effects, the framework the paper extends to a full eigenmode decomposition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the symmetry-based indicators used to argue that the dielectric matrix and its eigenwaves carry crystal-symmetry information."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the first-principles dielectric-matrix eigen-equation that the paper recasts as its microscopic optical band structure."}],"review_version":1}