pith. sign in

arxiv: 1907.02376 · v1 · pith:C7XNWZLJnew · submitted 2019-07-04 · ⚛️ physics.optics · physics.app-ph

Toroidal eigenmodes in all-dielectric metamolecules

Anna C. Tasolamprou , Odysseas Tsilipakos , Maria Kafesaki , Costas M. Soukoulis , Eleftherios N. Economou This is my paper

Pith reviewed 2026-05-25 09:11 UTC · model grok-4.3

classification ⚛️ physics.optics physics.app-ph
keywords toroidal dipolemetamoleculesall-dielectricpolaritonic rodseigenmodescoupled mode theorycollective resonancesmagnetic dipole
0
0 comments X

The pith

All-dielectric metamolecules with an odd number of rods support toroidal dipole modes that are more spectrally isolated from neighboring resonances than those with even numbers.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines electromagnetic resonant modes supported by polaritonic rods placed at the vertices of regular polygons using finite-element eigenvalue simulations, coupled mode theory, and a lumped wire model. These cyclic metamolecules support a toroidal dipole mode formed from the magnetic dipole mode of each rod, along with other collective modes whose frequencies oscillate around the single-rod resonance because the constituent modes are leaky. Odd-numbered rod ensembles produce larger frequency separation between the toroidal mode and its nearest neighbor than even-numbered ensembles. This separation, together with the toroidal mode's low quality factor, allows the usually silent toroidal dipole to couple more readily to external fields and to interact with the electric dipole. A reader would care because toroidal excitations have been hard to observe and control in conventional structures.

Core claim

Systems of polaritonic rods at the vertices of canonical polygons support an unconventional toroidal dipole mode consisting of the magnetic dipole mode in each rod. The resonant frequencies of all collective modes oscillate about the single-rod magnetic dipole resonance owing to the leaky character of the constituent modes. Ensembles with an odd number of rods exhibit larger frequency separation between the toroidal mode and its spectral neighbor than ensembles with an even number of rods; this increased isolation, combined with the low quality factor of the toroidal mode, favors coupling of the toroidal dipole to the outside world.

What carries the argument

The toroidal dipole mode formed by collective reorganization of magnetic dipole currents in the rods, captured by finite-element eigenvalue simulations and coupled mode theory.

If this is right

  • The toroidal dipole becomes more accessible for external excitation and interaction with the electric dipole in odd-rod arrangements.
  • All collective mode frequencies oscillate about the isolated-rod magnetic dipole frequency as rod separation varies.
  • The low quality factor of the toroidal mode assists its coupling to incident fields.
  • Coupled mode theory accounts for the current reorganizations that produce the observed collective modes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • These polygon-based structures could serve as building blocks for devices that selectively respond to toroidal excitations.
  • The odd-versus-even separation pattern may appear in other multipolar collective modes when similar rod arrangements are used.
  • Exploring non-canonical polygons or mixed rod sizes could test whether the isolation advantage generalizes.

Load-bearing premise

The finite-element eigenvalue simulations and coupled mode theory accurately capture the reorganization of currents and the resulting collective modes without significant numerical artifacts or invalid approximations for the leaky constituent modes.

What would settle it

Fabricate and measure the scattering spectra of odd-numbered and even-numbered rod clusters at the same separation; the toroidal resonance should show a measurably larger gap to its nearest neighbor in the odd case.

Figures

Figures reproduced from arXiv: 1907.02376 by Anna C. Tasolamprou, Costas M. Soukoulis, Eleftherios N. Economou, Maria Kafesaki, Odysseas Tsilipakos.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematics of the cyclic metamolecules under study [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. First four eigenmodes and corresponding eigenfre [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Results for [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Results for [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (a) TE [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. (a,b) Mode frequencies and quality factors as a funct [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. (a) TE [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Resonant frequency evolution with rod separation [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Resonant frequency and half-power bandwidth of the t [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Polarization current distribution for selected TE [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Toroidal mode: magnetic field lines and magnetic [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. TE [PITH_FULL_IMAGE:figures/full_fig_p014_12.png] view at source ↗
read the original abstract

We present a thorough investigation of the electromagnetic resonant modes supported by systems of polaritonic rods placed at the vertices of canonical polygons. The study is conducted with rigorous finite-element eigenvalue simulations. To provide physical insight, the simulations are complemented with coupled mode theory (the analog of LCAO in molecular and solid state physics) and a lumped wire model capturing the coupling-caused reorganizations of the currents in each rod. The systems of rods, which form all-dielectric cyclic metamolecules, are found to support the unconventional toroidal dipole mode, consisting of the magnetic dipole mode in each rod. Besides the toroidal modes, the spectrally adjacent collective modes are identified. The evolution of all resonant frequencies with rod separation is examined. They are found to oscillate about the single-rod magnetic dipole resonance, a feature attributed to the leaky nature of the constituent modes. Importantly, we observe that ensembles of an odd number of rods produce larger frequency separation between the toroidal mode and its neighbor than the ones with even number of rods. This increased spectral isolation, along with the low quality factor exhibited by the toroidal mode, favors the coupling of the commonly silent toroidal dipole to the outside world, rendering the proposed structure a prime candidate for controlling the observation of toroidal excitations and their interaction with the usually present electric dipole

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper investigates electromagnetic resonant modes in all-dielectric metamolecules consisting of polaritonic rods arranged at the vertices of canonical polygons. Rigorous finite-element eigenvalue simulations are complemented by coupled-mode theory (analogous to LCAO) and a lumped wire model to identify collective modes, including the toroidal dipole mode formed from the magnetic dipole resonances of individual rods. The evolution of resonant frequencies with rod separation is tracked, revealing oscillations about the single-rod magnetic dipole resonance attributed to the leaky nature of the modes. A key observation is that odd-numbered rod ensembles exhibit larger frequency separation between the toroidal mode and its nearest neighbor compared to even-numbered ensembles; combined with the toroidal mode's low quality factor, this is argued to enhance coupling to external fields and interaction with the electric dipole.

Significance. If the reported odd/even distinction in spectral isolation holds under scrutiny, the work supplies a concrete design rule for rendering toroidal dipole excitations observable in all-dielectric structures, where they are otherwise dark. The integration of eigenvalue simulations with analytic models for current reorganization provides useful physical insight, and the emphasis on leaky-mode effects leading to oscillatory frequency shifts is a distinctive feature that may generalize to other metamolecule systems.

major comments (1)
  1. [Numerical methods and results sections on frequency evolution] The central claim of larger toroidal-neighbor frequency separation for odd versus even rod counts (abstract and results on frequency evolution with separation) rests entirely on finite-element eigenvalue simulations of leaky constituent modes. The manuscript must demonstrate that the reported separations are robust against numerical artifacts by providing explicit convergence data with respect to mesh density, PML implementation, and computational domain truncation; without this, small shifts could erase or reverse the odd/even distinction, especially given that separations oscillate around the single-rod resonance.
minor comments (2)
  1. Clarify in the main text which specific polygons (e.g., triangle, square, pentagon) are simulated, as the abstract refers only to 'canonical polygons' without enumeration.
  2. [Coupled mode theory] The coupled-mode theory section would benefit from an explicit statement of the coupling coefficients used and how they are extracted from the simulations.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the detailed review and constructive feedback on our manuscript. We agree that explicit convergence data are necessary to support the numerical claims regarding frequency separations in leaky modes. We address the major comment below and will incorporate the requested material in the revised manuscript.

read point-by-point responses

  1. Referee: [Numerical methods and results sections on frequency evolution] The central claim of larger toroidal-neighbor frequency separation for odd versus even rod counts (abstract and results on frequency evolution with separation) rests entirely on finite-element eigenvalue simulations of leaky constituent modes. The manuscript must demonstrate that the reported separations are robust against numerical artifacts by providing explicit convergence data with respect to mesh density, PML implementation, and computational domain truncation; without this, small shifts could erase or reverse the odd/even distinction, especially given that separations oscillate around the single-rod resonance.

    Authors: We agree that convergence studies are essential for validating the reported odd/even distinction in frequency separation, given the oscillatory behavior and leaky character of the modes. In the revised manuscript we will add an appendix (or subsection in the numerical methods) that presents explicit convergence tests. These will include: (i) mesh density refinement (showing frequency shifts stabilize below a chosen tolerance), (ii) PML thickness and conductivity sweeps, and (iii) computational domain truncation (varying the outer boundary distance). Our existing internal checks indicate that the larger separation for odd rod counts persists across these refinements and is not an artifact; the new data will be included to make this transparent to readers. revision: yes

Circularity Check

0 steps flagged

No circularity: claims rest on independent finite-element simulations and coupled-mode analysis

full rationale

The paper reports resonant frequencies and odd/even separation trends directly from eigenvalue simulations of the metamolecule geometries, supplemented by standard coupled-mode theory. No fitted parameters are renamed as predictions, no self-citations supply load-bearing uniqueness theorems, and no ansatz or definition reduces the reported separation to an input by construction. The derivation chain is self-contained and externally falsifiable via the simulation method itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard electromagnetic simulation methods and the validity of coupled-mode approximations for leaky modes; no free parameters or invented entities are evident from the abstract.

axioms (1)
  • domain assumption Finite-element eigenvalue simulations and coupled mode theory (analog of LCAO) can reliably identify and explain collective toroidal and adjacent modes in rod polygons.
    Abstract states that simulations are complemented with these methods to provide physical insight.

pith-pipeline@v0.9.0 · 5783 in / 1217 out tokens · 35146 ms · 2026-05-25T09:11:46.519455+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages

  1. [1]

    Y. B. Zel’dovich, Sov. Phys. JETP 6, 1184 (1958)

    work page 1958

  2. [2]

    C. S. Wood, S. C. Bennett, D. Cho, B. P. Masterson, J. L. Roberts, C. E. Tanner, and C. E. Wieman, Science 275, 1759 (1997)

    work page 1997

  3. [3]

    W. C. Haxton and C. E. Wieman, Annu. Rev. Nucl. Sci. 51, 261 (2001)

    work page 2001

  4. [4]

    Ceulemans, L

    A. Ceulemans, L. F. Chibotaru, and P. W. Fowler, Phys. Rev. Lett. 80, 1861 (1998)

    work page 1998

  5. [5]

    Y. F. Popov, A. M. Kadomtseva, D. V. Belov, G. P. Vorob’ev, and A. K. Zvezdin, J. Exp. Theor. Phys. Lett. 69, 330 (1999)

    work page 1999

  6. [6]

    Kl¨ aui, C

    M. Kl¨ aui, C. A. F. Vaz, L. Lopez-Diaz, and J. A. C. Bland, J. Phys. Condens. Matter 15, R985 (2003)

    work page 2003

  7. [7]

    I. I. Naumov, L. Bellaiche, and H. Fu, Nature 432, 737 (2004)

    work page 2004

  8. [8]

    A. K. Zvezdin, V. V. Kostyuchenko, A. I. Popov, A. F. Popkov, and A. Ceulemans, Phys. Rev. B 80, 172404 (2009)

    work page 2009

  9. [9]

    Ungur, S

    L. Ungur, S. K. Langley, T. N. Hooper, B. Moubaraki, E. K. Brechin, K. S. Murray, and L. F. Chibotaru, J. Am. Chem. Soc. 134, 18554 (2012)

    work page 2012

  10. [10]

    G. N. Afanasiev and Y. P. Stepanovsky, J. Phys. A 28, 4565 (1995)

    work page 1995

  11. [11]

    V. M. Dubovik and V. V. Tugushev, Phys. Rep. 187, 145 (1990)

    work page 1990

  12. [12]

    E. E. Radescu and D. H. Vlad, Phys. Rev. E 57, 6030 (1998)

    work page 1998

  13. [13]

    E. E. Radescu and G. Vaman, Phys. Rev. E 65, 046609 (2002)

    work page 2002

  14. [14]

    A. T. G´ ongora and E. Ley-Koo, Rev. Mex. Fis. E 52, 188 (2006)

    work page 2006

  15. [15]

    G. N. Afanasiev, J. Phys. D: Appl. Phys 34, 539 (2001)

    work page 2001

  16. [16]

    Sawada and N

    K. Sawada and N. Nagaosa, Phys. Rev. Lett. 95, 237402 (2005)

    work page 2005

  17. [17]

    M. M. Hingorani and M. O’Donnell, Nat. Rev. Mol. Cell Biol. 1, 22 (2000)

    work page 2000

  18. [18]

    A. A. Simpson, Y. Tao, P. G. Leiman, M. O. Badasso, Y. He, P. J. Jardine, N. H. Olson, M. C. Morais, S. Grimes, D. L. Anderson, T. S. Baker, and M. G. Rossmann, Nature 408, 745 (2000)

    work page 2000

  19. [19]

    C. M. Ho and R. J. Scherrer, Phys. Lett. 722, 341 (2013)

    work page 2013

  20. [20]

    H. Wang, E. Yan, E. Borguet, and K. B. Eisenthal, Chem. Phys. Lett. 259, 15 (1996)

    work page 1996

  21. [21]

    F. Wise, L. Qian, and X. Liu, J. Nonlinear Opt. Phys. Mater. 11, 317 (2002)

    work page 2002

  22. [22]

    S. A. Maier and H. A. Atwater, J. Appl. Phys. 98, 011101 (2005)

    work page 2005

  23. [23]

    J. Shan, J. I. Dadap, I. Stiopkin, G. A. Reider, and T. F. Heinz, Phys. Rev. A 73, 023819 (2006)

    work page 2006

  24. [24]

    S. Lal, S. Link, and N. J. Halas, Nat. Photon. 1, 641 (2007)

    work page 2007

  25. [25]

    Kujala, B

    S. Kujala, B. K. Canfield, M. Kauranen, Y. Svirko, and J. Turunen, Opt. Express 16, 17196 (2008)

    work page 2008

  26. [26]

    Kaelberer, V

    T. Kaelberer, V. A. Fedotov, N. Papasimakis, D. P. Tsai, and N. I. Zheludev, Science 330, 1510 (2010)

    work page 2010

  27. [27]

    ¨Og¨ ut, N

    B. ¨Og¨ ut, N. Talebi, R. Vogelgesang, W. Sigle, and P. A. Van Aken, Nano Lett. 12, 5239 (2012)

    work page 2012

  28. [28]

    Z.-G. Dong, J. Zhu, J. Rho, J.-Q. Li, C. Lu, X. Yin, and X. Zhang, Appl. Phys. Lett. 101, 144105 (2012)

    work page 2012

  29. [29]

    Z.-G. Dong, P. Ni, J. Zhu, X. Yin, and X. Zhang, Opt. Express 20, 13065 (2012)

    work page 2012

  30. [30]

    Y. Fan, Z. Wei, H. Li, H. Chen, and C. M. Soukoulis, Phys. Rev. B 87, 115417 (2013)

    work page 2013

  31. [31]

    Savinov, V

    V. Savinov, V. A. Fedotov, and N. I. Zheludev, Phys. Rev. B 89, 205112 (2014)

    work page 2014

  32. [32]

    S.-H. Kim, S. Oh, K.-J. Kim, J.-E. Kim, H. Y. Park, O. Hess, and C.-S. Kee, Phys. Rev. B 91, 035116 (2015)

    work page 2015

  33. [33]

    D. W. Watson, S. D. Jenkins, J. Ruostekoski, V. A. Fe- dotov, and N. I. Zheludev, Phys. Rev. B 93, 125420 (2016)

    work page 2016

  34. [34]

    Huang, W

    Y.-W. Huang, W. T. Chen, P. C. Wu, V. A. Fedotov, N. I. Zheludev, and D. P. Tsai, Sci. Rep. 3, 1237 (2013)

    work page 2013

  35. [35]

    A. M. Zagoskin, A. Chipouline, E. Ilichev, J. R. Johans- son, and F. Nori, Sci. Rep. 5, 16934 (2015). 16

    work page 2015

  36. [36]

    A. A. Basharin, M. Kafesaki, E. N. Economou, C. M. Soukoulis, V. A. Fedotov, V. Savinov, and N. I. Zhe- ludev, Phys. Rev. X 5, 011036 (2015)

    work page 2015

  37. [37]

    A. E. Miroshnichenko, A. B. Evlyukhin, Y. F. Yu, R. M. Bakker, A. Chipouline, A. I. Kuznetsov, B. Luk’yanchuk, B. N. Chichkov, and Y. S. Kivshar, Nat. Commun. 6, 8069 (2015)

    work page 2015

  38. [38]

    W. Liu, J. Zhang, B. Lei, H. Hu, and A. E. Mirosh- nichenko, Opt. Lett. 40, 2293 (2015)

    work page 2015

  39. [39]

    K. C. Huang, M. L. Povinelli, and J. D. Joannopoulos, Appl. Phys. Lett. 85, 543 (2004)

    work page 2004

  40. [40]

    Yannopapas and E

    V. Yannopapas and E. Paspalakis, J. Opt. 12, 104017 (2010)

    work page 2010

  41. [41]

    A. S. Barker, A. A. Ballman, and J. A. Ditzenberger, Phys. Rev. B 2, 4233 (1970)

    work page 1970

  42. [42]

    J. A. Stratton, Electromagnetic theory (John Wiley & Sons, 1941)

    work page 1941

  43. [43]

    S. Fan, P. R. Villeneuve, J. Joannopoulos, M. Khan, C. Manolatou, and H. Haus, Phys. Rev. B 59, 15882 (1999)

    work page 1999

  44. [44]

    C. M. Gentry and M. A. Popovi´ c, Opt. Lett. 39, 4136 (2014)

    work page 2014

  45. [45]

    H. A. Haus and W. Huang, Proc. IEEE 79, 1505 (1991)

    work page 1991

  46. [46]

    Popovi´ c, C

    M. Popovi´ c, C. Manolatou, and M. Watts, Opt. Express 14, 1208 (2006)

    work page 2006

  47. [47]

    Sauvan, J

    C. Sauvan, J. P. Hugonin, I. S. Maksymov, and P. Lalanne, Phys. Rev. Lett. 110, 237401 (2013)

    work page 2013

  48. [48]

    S. G. Johnson, M. Ibanescu, M. Skorobogatiy, O. Weis- berg, J. Joannopoulos, and Y. Fink, Phys. Rev. E 65, 066611 (2002)

    work page 2002