REVIEW 3 major objections 5 minor 38 references
Primordial Media: the shrouded realm of composite materials
T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Composites with nonlocal components do not converge to the local effective medium as their parts shrink; instead they pass through a 'primordial metamaterial' regime and only deep sub-nonlocality reproduces homogenization, with different…
desk verdict The noncommutation claim is real but only shown inside a truncated quadratic model; worth refereeing, but the universal phrasing outruns the derivation. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the nonlocal wave equation (Eq. (3)) for the TM-polarized magnetic field $B_y$ in a uniaxial medium whose constitutive law (Eq. (1)) expresses nonlocality as a single, position-dependent coefficient $\alpha(z)$ multiplying the second $z$-derivative of $E_z$. This equation is fourth order in $z$ and therefore supports two TM waves per direction; their dispersion is the quartic (Eq. (11)) $\alpha \tilde{k}_z^4 + (\epsilon_{zz} - \alpha\epsilon_{\perp})\tilde{k}_z^2 + \epsilon_{\perp}(\tilde{k}_x^2 - \epsilon_{zz}) = 0$, whose second root is the additional wave that defines the primordial scale $l_p$. Smoothly varying fields require continuity of $D_z$, $E_x$, $E_z$, and $\alpha\,\partial E_z/\partial z$, which supply the additional boundary conditions and guarantee Poynting-flux conservation. Homogenizing Eq. (3) under slow field variation yields the nonlocal effective parameters of Eq. (7), and a 4x4 transfer-matrix method built on the same dispersion produces the transmission and Bloch-mode results that map the regimes.
What would settle it
Measure TM transmission at fixed incidence through a plasmonic-dielectric stack of fixed total thickness and fill fraction while the individual layer thickness $d$ is tuned across the nonlocality scale $l_p$ (e.g., from $\lambda_p/20$ down to $\lambda_p/10^4$). The paper predicts the single local effective-medium transmission dip splits into multiple minima as $d$ passes through $l_p$, and that at the smallest $d$ the spectrum re-converges to a single dip positioned according to Eq. (7), not Eq. (6). Observing instead a single dip that always follows the local effective medium would falsify the central claim.
Extended reading notes
Core claim
The central claim is that the response of a composite made of nonlocal components is not the local effective-medium response plus a small correction; it is a different object entirely. For a layered uniaxial medium with spatial dispersion along the stacking direction, the paper derives (Eq. (7)) the nonlocal effective-medium parameters $\epsilon^{\rm nl}_{\perp} = \langle \epsilon_{\perp} \rangle$, $\epsilon^{\rm nl}_{zz} = \langle \epsilon_{zz} \rangle$, and $\tilde{\alpha} = \langle 1/\alpha \rangle^{-1}$, in contrast to the local results $\epsilon^l_{zz} = \langle 1/\epsilon_{zz} \rangle^{-1}$. Because these two sets of parameters do not agree, a transition realm must exist between them; the paper identifies this realm as the primordial metamaterial, where the component size $d$ is comparable to the nonlocality scale $l_p$ and the optical response is set by the interference of additional (nonlocal) waves. The paper further shows the two limits $\alpha \to 0$ and $d \to 0$ do not commute, and derives the required additional boundary conditions from the macroscopic continuity of $E_z$ and $\alpha \partial E_z/\partial z$ rather than from a microscopic model.
Load-bearing premise
The derivation rests on the constitutive model in Eq. (1), where all nonlocality is a single scalar $\alpha(z)$ multiplying the second $z$-derivative of $E_z$; if real materials have higher-order or off-axis nonlocal terms, the two new regimes and the noncommutation of limits need not survive.
Editorial extensions
If this is right
- Deeply subwavelength multilayers made of nonlocal components must be described by the nonlocal effective medium of Eq. (7), not the local one; routine use of local Maxwell-Garnett or hyperbolic effective-medium theory in this regime would be qualitatively wrong.
- The region $d \sim l_p$ offers a new design space of 'primordial metamaterials' with highly oscillatory internal fields and multiple resonances, already reported experimentally (Ref. [24]).
- In materials with very long nonlocality (phononic, excitonic, $l_p\sim\lambda$), the local effective medium regime may be unreachable, leaving only photonic-crystal and primordial behavior.
- The additional boundary conditions derived from macroscopic continuity remove the need for quantum-mechanical input in layered nonlocal structures, making the formalism directly usable in design codes.
Reading between the lines
- Because Eq. (7) makes $\tilde{\alpha} = \langle 1/\alpha \rangle^{-1}$, the composite's nonlocality is governed by its most local component, an arithmetic-harmonic duality that may generalize to random composites; we suspect an effective-medium theory for disordered nonlocal stacks would take the same form in the homogenized limit.
- The noncommutation of limits is a hallmark of a singular perturbation: the nonlocal term raises the differential order, so a perturbative expansion around $\alpha=0$ cannot be uniform in $d$. An asymptotic matching across the primordial regime could yield quantitative criteria for when local effective-medium theory applies.
- A controlled experiment sweeping $d$ through $l_p$ with in-situ transmission and phase measurements would map the predicted regime boundaries; a natural next step is to test whether the primordial resonance splitting survives in disordered or 3D composite geometries.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies electromagnetic composites whose constituents have weak spatial dispersion, modeled in Eq. (1) by a quadratic nonlocal term proportional to α(z) ∂_z² E_z. From Maxwell equations the authors derive a fourth-order master equation for the magnetic field, Eq. (3), and use it to identify boundary conditions (continuity of E_x, D_z, E_z, and α ∂_z E_z) for nonlocal multilayers. For very fine layers they propose a nonlocal effective medium with parameters ε_⊥^nl = ⟨ε_⊥⟩, ε_zz^nl = ⟨ε_zz⟩, and ᾱ = ⟨1/α⟩^{-1} (Eq. (7)), which differs from the local effective medium (Eq. (6)). The paper claims that the limits of vanishing nonlocality (α→0) and vanishing layer thickness (d→0) do not commute, and that this gives rise to two new regimes: 'primordial metamaterials' (d∼l_p) and 'nonlocal effective medium' (d≪l_p). The claims are supported by transfer-matrix calculations for a plasmonic/dielectric multilayer.
Significance. If the central claim is correct, the paper would constitute a significant conceptual correction to effective-medium theory: nanostructured composites with nonlocal components need not converge to the local effective medium as the components shrink. The specific prediction that ε_zz^nl = ⟨ε_zz⟩ rather than the local harmonic average ⟨1/ε_zz⟩^{-1}, together with the identification of a primordial regime, is falsifiable and design-relevant for plasmonic, phononic, and excitonic composites. The paper also provides a useful explicit transfer-matrix framework for layered nonlocal media and derives additional boundary conditions from macroscopic electrodynamics, and it makes sample code available (Ref. [38]). However, the results presently rest on a quadratic truncation of the nonlocal response and on the authors' own numerical implementation, so the strength of the claim is not yet matched by independent verification.
major comments (3)
- [§II, Eq. (1) and Appendix Eq. (11)] The constitutive model (1) truncates the nonlocal response at the second derivative, justified by inversion symmetry and weak spatial dispersion. However, the primordial regime is defined by d∼l_p, where the additional wave has k_z∼l_p^{-1} (from Eq. (11) with α∼(l_p/λ)^2). At this wavevector the dimensionless expansion parameter is k_z l_p∼O(1), so the next allowed term β(z)∂_z^4 E_z, with β∼l_p^4, contributes to the dispersion at the same order as the retained α term. Consequently the noncommutation claim is established only within the quadratic model, not for the full nonlocal response of the physical systems cited (Fermi gas, phonon-polaritons, excitons). The authors should either provide quantitative bounds on the omitted terms in the primordial regime or redo the central derivation with a nonlocal kernel that is not quadratically truncated.
- [§III, Eq. (7)] Equation (7) is introduced with 'a similar procedure' applied to nonlocal composites, but the homogenization step from the master equation (3) to the nonlocal effective-medium equation (5) is not shown. The local result (6) follows from continuity of E_x and D_z across the unit cell; the nonlocal case additionally requires handling the α∂_z E_z continuity condition and a two-scale asymptotic expansion. Without this derivation the reader cannot verify the striking claim ε_zz^nl=⟨ε_zz⟩ rather than the local harmonic average, nor the harmonic average of 1/α. Please provide the explicit reduction, including the conditions under which field variations are slow compared with d in the nonlocal effective-medium regime.
- [§IV, Figs. 3–5 and Appendix B] The evidence for the two new regimes is generated entirely with the authors' own transfer-matrix code based on the same constitutive model (1), and the primordial regime is linked to the authors' prior work [24] rather than to an independent measurement reported here. The central claim would be materially strengthened by at least one cross-check: e.g., comparison with a finite-difference or finite-element solution of the full Maxwell-nonlocal equations, or with the experimental data of Ref. [24]. Without such a check it is difficult to exclude the possibility that the predicted regime structure is an artifact of the quadratic truncation combined with the particular mode-matching scheme.
minor comments (5)
- [Throughout] The manuscript contains several typographical errors: 'vaccuum' in Fig. 1, 'qunatum' in §II, 'metamateirals' in §VII, 'Permagon' in Refs. [25] and [28], and the GitHub URL in Ref. [38] contains a space ('primordial metamaterials').
- [§II, Eq. (1) and §IV, Eq. (8)] The relation between the abstract coefficient α(z) in Eq. (1) and the material parameters used in the numerical example (Fermi velocity, plasma frequency, scattering rates) is not given explicitly; please provide the mapping α = 12π² v_f²/(5c²) for the degenerate Fermi gas or state the formula used to generate Fig. 2.
- [§II, paragraph after Eq. (3)] The sentence 'It can be shown that the expression in square brackets in Eq. (3) is proportional to E_z' is a key step for the additional boundary conditions, but the proportionality is not demonstrated. A short derivation or an explicit expression would help the reader verify the ABCs.
- [Appendix B] The material matrix \(\hat N_l\) in Eq. (12) has rows labeled by E_x, B_y, α ∂E_z/∂z, and E_z, but the boundary conditions are stated to require continuity of E_x, D_z, E_z, and α ∂E_z/∂z. Since D_z is proportional to 1/ε_⊥ ∂B_y/∂z, please clarify how the rows in Eq. (12) are ordered relative to these four quantities.
- [§VI] The discussion of the noncommutation of limits would benefit from a more precise statement of what 'α→0' means when the nonlocal scale l_p is finite; as written, the limit is taken at fixed d and λ, but the dependence of α on the material parameters and on frequency is not made explicit in the limiting statement.
Circularity Check
No significant circularity: the central derivation is self-contained from the stated constitutive model, and same-author citations are not load-bearing.
full rationale
The paper's derivation chain is self-contained. Equation (3) follows from Maxwell's equations plus the explicitly stated constitutive relation in Eq. (1), which is an input assumption rather than a fitted or imported result. The effective-medium parameters in Eq. (7) are obtained by homogenizing Eq. (3) under the stated slow-variation condition; they are not fitted to the numerical data. The claimed noncommutation of the limits alpha -> 0 and d -> 0 is a mathematical consequence of the different homogenized expressions in Eq. (6) versus Eq. (7) and is exhibited directly in the paper. Same-author citations, including [24] for the experimental realization of primordial metamaterials and [38] for sample code, are used as supporting evidence or implementation details, not as load-bearing steps in the derivation; no claim in the derivation chain reduces to those citations. The quadratic truncation of the nonlocal constitutive model is an explicit modeling assumption whose validity at k*l_p ~ 1 is a legitimate correctness concern, but it is not circularity. Accordingly, no circular step is identified and the overall circularity score is 0.
Assumptions & free parameters
free parameters (4)
- alpha(z) nonlocality coefficient =
not independently measured; implied l_p from v_f/c ~ 3e-3 in numerics
- dielectric nonlocality strength =
1e-6 (5+1i) k^2 c^2 / omega^2
- tau, local scattering rate =
0.1 omega_p
- tau_nl, nonlocal scattering rate =
0.2 omega_p
assumptions (5)
- standard math Maxwell equations and harmonic time and space dependence
- domain assumption Nonlocal response has the self-adjoint quadratic form Eq.(1)
- ad hoc to paper Weak spatial dispersion: higher-order derivatives in D(E) can be neglected even when d is comparable to l_p
- domain assumption Homogenization by averaging continuous field components across the unit cell remains valid for nonlocal media
- domain assumption Uniaxial geometry and TM-only analysis are sufficient for the general regime classification
Cite this review
Pith. "Pith review of Primordial Media: the shrouded realm of composite materials." pith.science (2026). https://pith.science/paper/RIBS5FU7
@misc{pith2026250711373,
author = {Pith},
title = {Pith review of: Primordial Media: the shrouded realm of composite materials},
year = {2026},
howpublished = {\url{https://pith.science/paper/RIBS5FU7}},
note = {Machine review of arXiv:2507.11373}
}
read the original abstract
Electromagnetic composites (metamaterials) recently underwent explosive growth fueled, in part, by advances in nanofabrication. It is commonly believed that as the size of the components decreases, the behavior of a composite converges to the response of a homogeneous material. Here we show that this intuitive understanding of the electromagnetic response of composite media is fundamentally flawed, even at the qualitative level. In contrast to the well-understood local effective medium response, the properties of nanostructured composites can be dominated by electromagnetic nonlocality. We demonstrate that the interplay between the nonlocality and the structural inhomogeneity introduces two fundamentally new electromagnetic regimes, primordial metamaterials and nonlocal effective medium. We develop an analytical description of these regimes and show that the behavior of metamaterials in the limits of vanishing nonlocality and of vanishing component size do not commute. Our work opens a new dimension in the design space of nanostructured electromagnetic composites.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[24]
A. Ware, J. LaMountain, R. White, S. Bank, D. Wasserman, E. Narimanov, and V. Podolskiy, Primordial metamaterials, ArXiv , 2505.21359 (2025)
arXiv 2025
-
[38]
Primordial metmaterials github, https://github.com/viktor- podolskiy/primordial metamaterials. 20
- [1]
-
[2]
J. B. Pendry, Negative refraction makes a perfect lens, Phys. Rev. Lett.85, 3966 (2000)
2000
-
[3]
N. Engheta, Circuits with light at nanoscales: Optical nanocircuits inspired by metamaterials, Science317, 1698 (2007), https://www.science.org/doi/pdf/10.1126/science.1133268
-
[4]
A. Poddubny, I. Iorsh, P. Belov, and Y. Kivshar, Hyperbolic metamaterials, Nat. Photon7, 948 (2013)
work page 2013
-
[5]
Shalaev, Optical negative index metamaterials, Nat
V. Shalaev, Optical negative index metamaterials, Nat. Photon1, 41 (2007)
work page 2007
-
[6]
V. Veselago and E. Narimanov, The left hand of brightness: past, present and future of negative index materials, Nat.Mater5, 759 (2006)
work page 2006
Show all 38 references
-
[7]
Milton,The theory of composites(Cambridge U
G. Milton,The theory of composites(Cambridge U. Press, 2002)
2002
-
[8]
Elser, V
J. Elser, V. A. Podolskiy, I. Salakhutdinov, and I. Avrutsky, Nonlo- cal effects in effective-medium response of nanolayered metamaterials, Ap- 17 plied Physics Letters90, 191109 (2007), https://pubs.aip.org/aip/apl/article- pdf/doi/10.1063/1.2737935/14374786/191109 1 online.pdf
2007 doi
-
[9]
J. Kern, S. Großmann, N. V. Tarakina, T. H¨ ackel, M. Emmerling, M. Kamp, J.-S. Huang, P. Biagioni, J. C. Prangsma, and B. Hecht, Atomic-scale confinement of resonant optical fields, Nano Letters12, 5504 (2012), pMID: 22984927, https://doi.org/10.1021/nl302315g
2012 doi
-
[10]
D. O. Sigle, J. Mertens, L. O. Herrmann, R. W. Bowman, S. Ithurria, B. Dubertret, Y. Shi, H. Y. Yang, C. Tserkezis, J. Aizpurua, and J. J. Baumberg, Monitoring morphological changes in 2d monolayer semiconductors using atom-thick plasmonic nanocavities, ACS Nano9, 825 (2015), ...
2015 doi
-
[11]
C. S. Levin, C. Hofmann, T. A. Ali, A. T. Kelly, E. Morosan, P. Nordlander, K. H. Whitmire, and N. J. Halas, Magnetic-plasmonic core-shell nanoparticles, ACS Nano3, 1379 (2009), pMID: 19441794, https://doi.org/10.1021/nn900118a
2009 doi
-
[12]
C. P. Byers, H. Zhang, D. F. Swearer, M. Yorulmaz, B. S. Hoener, D. Huang, A. Hog- gard, W.-S. Chang, P. Mulvaney, E. Ringe, N. J. Halas, P. Nordlander, S. Link, and C. F. Landes, From tunable core-shell nanoparticles to plasmonic drawbridges: Ac- tive control of nanoparticle ...
2015 doi
-
[13]
Fischer and M
J. Fischer and M. Wegener, Three-dimensional optical laser lithography beyond the diffraction limit, Laser & Photonics Reviews7, 22 (2013), https://onlinelibrary.wiley.com/doi/pdf/10.1002/lpor.201100046
2013 doi
-
[14]
J. J. Baumberg, J. Aizpurua, M. H. Mikkelsen, and D. R. Smith, Extreme nanophotonics from ultrathin metallic gaps, Nature materials18, 668 (2019)
2019
-
[15]
A. J. Muhowski, E. Simmons, K. Li, E. E. Narimanov, V. A. Podolskiy, and D. Wasserman, Extending plasmonic response to the mid-wave infrared with all-epitaxial composites, Opt. Lett.47, 973 (2022)
2022
-
[16]
I. A. Larkin and M. I. Stockman, Imperfect perfect lens, Nano Letters5, 339 (2005), pMID: 15794622, https://doi.org/10.1021/nl047957a
2005 doi
-
[17]
S. Peng, J. M. McMahon, G. C. Schatz, S. K. Gray, and Y. Sun, Reversing the size-dependence of surface plasmon resonances, Proceedings of the National Academy of Sciences107, 14530 (2010), https://www.pnas.org/doi/pdf/10.1073/pnas.1007524107
2010 doi
-
[18]
F. J. Garc ´ ıa de Abajo, Nonlocal effects in the plasmons of strongly interacting nanopar- 18 ticles, dimers, and waveguides, The Journal of Physical Chemistry C112, 17983 (2008), https://doi.org/10.1021/jp807345h
2008 doi
-
[19]
Nordlander and J
P. Nordlander and J. C. Tully, Energy shifts and broadening of atomic levels near metal surfaces, Phys. Rev. B42, 5564 (1990)
1990
-
[20]
Sp¨ agele, M
C. Sp¨ agele, M. Tamagnone, D. Kazakov, M. Ossiander, M. Piccardo, and F. Capasso, Multi- functional wide-angle optics and lasing based on supercell metasurfaces, Nature Communica- tions12, 3787 (2021)
2021
-
[21]
Malek, A
S. Malek, A. Overvig, A. Alu, and Y. N., Multifunctional resoannt wavefront-shaping meta- optics based on multilayer and multi-perturbation nonlocal metasurfaces, Light: Sci. and Appl.11, 246 (2022)
2022
-
[22]
Y. Zhou, S. Guo, A. C. Overvig, and A. Al` u, Multiresonant nonlocal metasurfaces, Nano Letters23, 6768 (2023), pMID: 37307588, https://doi.org/10.1021/acs.nanolett.3c00772
2023 doi
-
[23]
R. J. Pollard, A. Murphy, W. R. Hendren, P. R. Evans, R. Atkinson, G. A. Wurtz, A. V. Zayats, and V. A. Podolskiy, Optical nonlocalities and additional waves in epsilon-near-zero metamaterials, Phys. Rev. Lett.102, 127405 (2009)
2009
-
[25]
Landau, E
L. Landau, E. Lifshitz, and L. Pitaevskii,Electrodynamics of Continuous Media(Permagon Press, Oxford, 1984)
1984
-
[26]
Agranovich and V
V. Agranovich and V. Ginzburg,Crystal Optics with Spatial Dispersion and Excitons(Nauka, Moscow, 1979)
1979
-
[27]
Pekar, The theory of electromagnetic waves in a crystal in which excitons are produced, Sov
S. Pekar, The theory of electromagnetic waves in a crystal in which excitons are produced, Sov. Phys. JETP6, 785 (1958)
1958
-
[28]
Landau, E
L. Landau, E. Lifshitz, and L. Pitaevskii,Physical Kinetics(Permagon Press, Oxford, 1984)
1984
-
[29]
Kiselev, B
V. Kiselev, B. Razbirin, and I. Uraltsev, Interference states of excitons, observation of addi- tional waves, Sov. Phys. JETP18, 504 (1973)
1973
-
[30]
J. J. Hopfield and D. G. Thomas, Theoretical and experimental effects of spatial dispersion on the optical properties of crystals, Phys. Rev.132, 563 (1963)
1963
-
[31]
G. S. Agarwal, D. N. Pattanayak, and E. Wolf, Electromagnetic fields in spatially dispersive media, Phys. Rev. B10, 1447 (1974)
1974
-
[32]
Wangberg, J
R. Wangberg, J. Elser, E. E. Narimanov, and V. A. Podolskiy, Nonmagnetic nanocomposites 19 for optical and infrared negative-refractive-index media, Journal of the Optical Society of America B23, 498 (2006)
2006
-
[33]
A. J. Hoffman, L. Alekseyev, S. S. Howard, K. J. Franz, D. Wasserman, V. A. Podolskiy, E. E. Narimanov, D. L. Sivco, and C. Gmachl, Negative refraction in semiconductor metamaterials, Nature Materials6, 946 (2007)
2007
-
[34]
Avrutsky, I
I. Avrutsky, I. Salakhutdinov, J. Elser, and V. Podolskiy, Highly confined optical modes in nanoscale metal-dielectric multilayers, Physical Review B75, 241402 (2007)
2007
-
[35]
V. A. Podolskiy and E. E. Narimanov, Strongly anisotropic waveguide as a nonmagnetic left- handed system, Physical Review B—Condensed Matter and Materials Physics71, 201101 (2005)
2005
-
[36]
P. Yeh, A. Yariv, and C.-S. Hong, Electromagnetic propagation in periodic stratified media. I. General theory, J. Opt. Soc. Am.67, 423 (1977)
1977
-
[37]
Langevin, P
D. Langevin, P. Bennet, A. Khaireh-Walieh, P. Wiecha, O. Teytaud, and A. Moreau, Pymoosh: a comprehensive numerical toolkit for computing the optical properties of multilayered struc- tures, J. Opt. Soc. Am. B41, A67 (2024)
2024
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.