Abstract
Wavelength determines the length scale of the cross section when electromagnetic waves are scattered by an electrically small object. The cross section diverges for resonant scattering, and diminishes for nonresonant scattering, when wavelength approaches infinity. This scattering law explains the colour of the sky as well as the strength of a mobile phone signal. We show that such wavelength scaling comes from the conical dispersion of free space at zero frequency. Emerging Weyl systems, offering similar dispersion at nonzero frequencies, lead to new laws of electromagnetic scattering that allow cross sections to be decoupled from the wavelength limit. Diverging and diminishing cross sections can be realized at any target wavelength in a Weyl system, providing the ability to tailor the strength of wave–matter interactions for radiofrequency and optical applications.
Introduction
Electromagnetic scattering is a fundamental process that occurs when waves in a continuum interact with an electrically small scatterer. Scattering is weak under nonresonant conditions; an example is Rayleigh scattering, which is responsible for the colours of the sky. Conversely, scattering becomes much stronger with resonant scatterers, which have an internal structure supporting localized standing waves, such as antennas, optical nanoresonators, and quantum dots. Resonant scatterers have wide application because the resonance allows physically small scatterers to capture wave energy from a large area. As such, large electromagnetic cross sections, σ, are always desirable: a larger σ value means, for example, stronger mobile phone signals^{1} and higher absorption rates for solar cells^{2}.
The maximum cross section of resonant scattering is bounded by the fundamental limit of electrodynamics. One might be tempted to enlarge the scatterer to increase the cross section, but this strategy only works for nonresonant scattering, or electrically large scatterers. In resonant scattering, physical size only affects spectral bandwidth, while the limit of cross section is determined by the resonant wavelength λ as^{3}:
The directivity D describes the anisotropy of the scattering; D = 1 for isotropic scatterers. Equation 1 shows that an atom^{4} can have a \(\sigma _{{\mathrm{max}}}\) similar to that of an optical antenna^{5}, despite the subnanometre size of the atom. This also means that optical scatterers cannot attain resonant cross sections as large as those of radiofrequency (RF) antennas, due to the smaller wavelengths involved.
Overcoming the limit just described has farreaching implications for RF and optoelectronic applications. Many efforts have been devoted to realizing this goal, including the use of enhanced directivity D ^{6,7,8,9,10}, degenerate resonances^{11}, decreased dielectric constants^{12} \({\it{\epsilon }}\) and materials with negative refractive index^{13}. While these approaches exploit certain tradeoffs to slightly increase the prefactor in Eq. (1), the fundamental limit of λ ^{2} remains, which can be proven directly from Maxwell’s equations without requiring specific scatterer details^{5}. Until now, extremely large cross sections have only been obtained at long wavelengths near the DC frequency.
Here we show that the scattering laws in Weyl systems^{14} allows the cross section to be decoupled from the wavelength limit. This opens a new path to realizing strong wave–matter interaction, providing potential benefits to RF and optoelectronic devices that rely on resonant scattering. Moreover, the scattering effect discovered here is equally applicable to acoustic or electronic waves^{15,16,17}.
Results
Length scale of the cross section
In free space, the dispersion directly leads to the wavelength limit of the cross section, which is shown by Eq. (1). Specifically, the DC point, which is located at the apex of the conical dispersion relation as shown in Fig. 1a, gives rise to the diverging cross sections at low frequencies. The apex of the conic dispersion can be also realized at other spectral regimes, such as that shown in Fig. 1b. All the special scattering properties associated with the DC point can be reproduced, resulting in diverging cross sections and exceptionally strong light–matter interactions at high frequencies. Weyl points^{14, 18,19,20,21,22,23,24,25,26,27,28}, the threedimensional (3D) analogy of Dirac points, have recently been shown to exhibit such conical dispersion relations. The scattering laws in Weyl systems allows the cross section to be decoupled from the wavelength limit. Extraordinarily large cross section can be realized even for very small scatterers inside a Weyl photonic crystal, which is schematically illustrated in Fig. 1c. Unlike the lensing effect that can concentrate incident waves to a fixed focus, the scatterer concentrates incident waves to itself no matter where it is placed inside the Weyl photonic crystal.
Conservation law of resonant scattering
To illustrate the underlying physics of the relation between the dispersion and cross section, we will first show a conservation law of resonant scattering. We start by considering the resonant cross section of a dipole antenna (Fig. 2a). Without losing generality, we only discuss the scattering cross section, assuming zero absorption. Similar conclusions can be drawn for the absorptive case, with the maximum absorption cross section \(\frac{1}{4}\) that of the scattering cross section^{29}.
The dipole antenna is anisotropic due to its elongated shape, so the cross section \(\sigma (\theta ,\varphi )\) depends on the incident direction of the wave. At a normal direction, when \(\theta = \pi /2\), it reaches its maximum value^{5} of \(\sigma _{{\mathrm{max}}} = 3\lambda ^2/2\pi\). Along the axial direction, when \(\theta = 0\), the cross section vanishes. While it is straightforward to calculate the cross section of a dipole antenna, it is not immediately apparent why the cross section follows the λ ^{2} rule and diverges around the DC frequency.
Figure 2b shows the realspace representation of \(\sigma (\theta ,\varphi )\) at the resonant frequency ω _{0}. We can also represent \(\sigma (\theta ,\phi )\) in momentum space as \(\sigma ({\mathbf{k}})\), with k located on the isosurface defined by \(\omega = \omega _0\). As shown in Fig. 2c, the isosurface is a sphere with \({\mathbf{k}} = \omega _0/c\). To visualize the momentumspace representation, \(\sigma ({\mathbf{k}})\) is indicated by the colour intensity on the isosurface in Fig. 2d. As we show in Supplementary Note 1, the resonant cross section satisfies the following conservation law:
The integration is performed on the isosurface \(\omega \left( {\mathbf{k}} \right) = \omega _0\). Our proof is based on quantum electrodynamics^{30}, so it applies to classical scatterers such as antennas, as well as to quantum scatterers such as electronic transitions that absorb and emit light. More importantly, the continuum, in which the scatterer is embedded, does not need to be free space; it can be anisotropic materials, or even photonic crystals^{31}, as long as a welldefined dispersion relation \(\omega = \omega ({\mathbf{k}})\) exists.
Equation 2 dictates the scaling of σ with respect to the resonant frequency ω _{0} of the scatterer. As ω _{0} decreases, the area of the isosurface shrinks. To conserve the value of the integration over a smaller isosurface, the cross section \(\sigma \left( {\mathbf{k}} \right)\) must increase accordingly. For example, Fig. 2e and f show the isosurfaces of dipole antennas with resonant frequencies at \(\omega _0/2\) and \(\omega _0/3\), respectively. The cross sections, indicated by colour intensity, must increase proportionally to maintain a constant integration over these smaller isosurfaces. This can also be seen by the average cross section in the momentum space:
Here, \(S \equiv {\iint} {ds}\) is the area of the isosurface. When approaching the apex of the conical dispersion relation, the isosurface diminishes, i.e., \(S \to 0\). As a result, the cross section diverges around the DC point (Fig. 1a).
Resonant scattering in Weyl photonic crystal
Recent demonstrations of Weyl points in photonic crystals^{18} show that the dispersion of a 3D continuum can exhibit conical dispersion at any designed frequency (Fig. 1b). The Hamiltonian for the continuum around the Weyl point is given by \({\cal H}\left( {\mathbf{q}} \right) = v_xq_x\sigma _x + v_yq_y\sigma _y + v_zq_z\sigma _z\), where \(\sigma _{x,y,z}\) are Pauli matrices. The momentum \({\mathbf{q}}\left( {q_x,q_y,q_z} \right) = {\mathbf{k}}  {\mathbf{k}}_{{\mathrm{Weyl}}}\) defines the distance to the Weyl point in the momentum space, and q = 0 at the Weyl point. The linear dispersion described by this Hamiltonian produces an ellipsoidal isosurface that encloses the Weyl point. The isosurface shrinks to a point at the Weyl frequency \(\omega _{{\mathrm{Weyl}}}\). Scattering properties associated with the DC point are carried to high frequencies within Weyl medium, resulting in exceptionally strong resonant scattering. As illustrated in Fig. 1b, the average cross section scales as \(\bar \sigma \sim \frac{1}{{\left( {\omega _0  \omega _{{\mathrm{Weyl}}}} \right)^2}}\) around the Weyl point. This allows high frequency resonant scatterers, such as atoms and quantum dots, to attain large cross sections, potentially at a macroscopic scale.
We now demonstrate a specific example of resonant scattering in a Weyl photonic crystal^{22}. We consider a localized resonant scatterer in an infinitely large photonic crystal. The simulations are performed in two steps. First, we numerically calculate the eigenmodes of the resonant frequency in the Brillouin zone by using MIT photonic bands^{32} (MPB). Next, we use each eigenmode as excitation and numerically calculate the scattering cross section by using the quantum scattering theory we developed recently^{33}, which is described in detail in Supplementary Note 2.
The structure of the Weyl photonic crystal consists of two gyroids, as shown in Fig. 3a. The magenta gyroid is defined by the equation \(f\left( {\mathbf{r}} \right) >1.1\), where \(f\left( {\mathbf{r}} \right) = \sin (2\pi x/a)\cos (2\pi y/a) + \sin (2\pi y/a)\cos (2\pi z/a) + \sin (2\pi z/a)\cos (2\pi x/a)\) and a is the lattice constant. The yellow gyroid is the spatial inversion of the magenta one. The two gyroids are filled with a material with a dielectric constant of \({\it{\epsilon }} = 13\). To obtain Weyl points, we add four air spheres to the gyroids to break the inversion symmetry (Fig. 3a). These spheres are related by an \(S_4(z)\) transformation. The resulting band structure has four isolated Weyl points at the same frequency of \(\omega _{{\mathrm{Weyl}}} = 0.5645\,(2\pi c/a)\). Figure 3b illustrates conical dispersion in the \(2k_z =  k_x  k_y\) plane.
The resonant scatterer is a quantum twolevel system (TLS) embedded in the above photonic crystal. We numerically solve the scattering problem using quantum electrodynamics^{33, 34}. The Hamiltonian is \(= H_{{\mathrm{pc}}} + H_{{\mathrm{TLS}}} + H_{\mathrm{I}}\). The first two terms, \(H_{{\mathrm{pc}}} = \mathop {\sum}\nolimits_{\mathbf{q}} {\hbar \omega _{\mathbf{q}}c_{\mathbf{q}}^\dagger c_{\mathbf{q}}}\) and \(H_{{\mathrm{TLS}}} = \hbar \omega _0b^\dagger b\), are the Hamiltonian of the photons and the TLS, respectively^{35}. Here, \(\hbar\) is the reduced Planck constant, b ^{†} and b are the raising and lowering operators for the quantum dot, respectively, and \(c_{\mathbf{q}}^\dagger\)and c _{ q } are the bosonic creation and annihilation operators of the photons, respectively. The Lamb shift^{35} is incorporated into the resonant frequency ω _{0}. The third term, \(H_{\mathrm{I}} = i\hbar \mathop {\sum}\nolimits_{\mathbf{q}} {g_{\mathbf{q}}(c_{\mathbf{q}}^\dagger b  c_{\mathbf{q}}b^\dagger )}\), is the interaction between the TLS and the radiation. The coupling coefficient is \(g_{\mathbf{q}} = {\mathbf{d}} \cdot {\hat{\mathbf E}}_{\mathbf{q}}\sqrt {\omega _0/2\hbar {\it{\epsilon }}_0L^3}\), where d is the dipole moment, \({\it{\epsilon }}_0\) is vacuum permittivity, \({\hat{\mathbf E}}_{\mathbf{q}}\) is the unit polarization vector of the photons and L ^{3} is the quantization volume.
Scattering inside photonic crystals is much more complex than that in free space. The continuum is highly dispersive, anisotropic and nonuniform, and thus the scattering cross section depends strongly on the location of the scatterer and orientation within the photonic crystal.
As an example, we consider a TLS with a transition frequency ω _{0} slightly below the Weyl frequency: \(\omega _0  \omega _{{\mathrm{Weyl}}} =  0.0005\,(2\pi c/a)\). The calculated cross section \(\sigma \left( {\mathbf{q}} \right)\) (see Supplementary Note 2 for derivation) is plotted on the isosurface as shown by Fig. 3e, f. As expected, it strongly depends on the incident wavevector q. In addition, \(\sigma \left( {\mathbf{q}} \right)\) varies greatly at different locations, as shown by comparing Fig. 3e and Fig. 3f. Despite all these differences, when integrated over the isosurface, \({\iint} {\sigma \left( {\mathbf{q}} \right)ds}\) always results in the same constant: \(16\pi ^2\). We perform the integration for a TLS at 20 different locations, all with the same constant, as shown in Fig. 3d.
As the transition frequency of the TLS approaches the Weyl point, i.e., \(\omega _0 \to \omega _{{\mathrm{Weyl}}}\), the isosurface shrinks in size, as illustrated by the insets of Fig. 4a. The conservation law leads to an increasing \(\sigma ({\mathbf{q}})\), as shown by stronger colours. Near the Weyl frequency (black dashed line), the average cross sections \(\bar \sigma\) is enhanced by three orders of magnitude compared to that in free space, eventually diverging at the Weyl point (Fig. 4a). The analytical prediction from Eq. 2 and the area of the isosurface agree very well with predictions from numerical simulation (circles in Fig. 4a).
Frequency dependence of resonant scattering
A Weyl point greatly enhances the cross section at the resonant frequency. However, it comes at the price of suppressed cross section away from the resonant frequency. Next, we discuss the spectral feature of the average cross section \(\bar \sigma (\omega )\) for a given TLS. The spectral dependence is shown in Supplementary Note 3 as
Here, \(\wp\) is a constant that depends on the local electric field at the position of the TLS, but does not vary significantly with frequency. In order to derive Eq. 4, we use the fact that the spontaneous decay rate is proportional to \(\left( {\omega _0  \omega _{{\mathrm{Weyl}}}} \right)^2\) (Supplementary Eq. 32). At the resonance when \(\omega _0  \omega = 0,\) the average cross section scales as \(1/\left( {\omega _0  \omega _{{\mathrm{Weyl}}}} \right)^2\), which increases as the resonant frequency moves closer to the Weyl point. However, away from the resonance when \(\omega _0  \omega  \gg \omega _0  \omega _{{\mathrm{Weyl}}}\), Eq. 4 reduces to \(\bar \sigma \left( {\omega _0} \right)\sim \left( {\omega _0  \omega _{{\mathrm{Weyl}}}} \right)^2/\left( {\omega _0  \omega } \right)^2\), which shows that being close to the Weyl point suppresses the cross section. In Fig. 4b, we calculate the spectra for three different TLSs with their transition frequencies approaching the Weyl point (black dashed line). While the peak value of the cross section grows, the full width at half maximum of the spectrum decreases. The spectral integration of the cross section remains around a constant (Supplementary Note 3).
Rayleigh or nonresonant scattering
The nonresonant scattering of electrically small objects in free space follows the Rayleigh scattering law. In great contrast to resonant scattering, the cross section of Rayleigh scattering scales as \(\sigma \sim \omega ^4\), and diminishes at the DC frequency, as shown in Fig. 5a. This property is also closely related to the isosurface and can be carried to high frequency at the Weyl point (see proof in Supplementary Note 5). While the resonant cross section diverges, the nonresonant cross section diminishes at Weyl points, as illustrated in Fig. 5b.
Using perturbation theory and the firstorder Born approximation^{36}, the Rayleigh scattering cross section can be shown as (details in Supplementary Note 5):
Here, \({\mathbf{u}}_{\mathbf{k}}({\mathbf{r}})\) is the eigenmode of the Weyl photonic crystal associated with wavevector k. \({\mathbf{k}}_{\mathrm{s}}\) and \({\mathbf{k}}_{{\mathrm{inc}}}\) are the wavevectors of the scattered and incident eigenmodes, respectively. V is a tensor for the scattering potential of the Rayleigh scatterer. The integral is proportional to the area of the isosurface S. At the DC point, the isosurface shrinks to a point with \(S = 0\), and the cross section of Rayleigh scattering is zero. Similarly, around the Weyl point, the area of the isosurface \(S\sim {\mathrm{\Delta }}\omega ^2 = \left( {\omega  \omega _{{\mathrm{Weyl}}}} \right)^2\). The cross section scales as \(\sigma \sim \omega ^2\Delta \omega ^2\), and is zero at the Weyl point (Fig. 5b).
To validate our theoretical prediction above, we numerically calculate the Rayleigh scattering cross section of a small dielectric sphere embedded in the same Weyl photonic crystal. We use MPB^{32} to obtain the eigenmode of the Weyl photonic crystal, then numerically calculate the Rayleigh scattering cross section using the normalmode expansion^{37,38,39} (see more details in Supplementary Note 6). The calculated Rayleigh scattering cross sections around the Weyl point are obtained by averaging 100 scatterers at random locations within one unit cell of the photonic crystal, and are plotted in Fig. 5c as red circles. They are normalized by \(\sigma _{\mathrm{R}}\), the scattering cross section of the same scatterer in free space. The blue dashed line indicates the scaling law \(\sigma \sim \omega ^2\Delta \omega ^2\). The calculated Rayleigh scattering cross section agrees with the theoretical prediction well, with a vanishing cross section observed at the Weyl point. For highfrequency devices, such as integrated waveguides and laser cavities, Rayleigh scattering caused by interface roughness degrades performance and increases noise. The combination of suppressed Rayleigh scattering and enhanced resonant scattering could make Weyl media attractive for these optoelectronic applications.
Discussion
Moreover, the conservation law of resonant scattering can also be extended to lowerdimensional space. In twodimensional photonic crystals, diverging cross sections can be realized at Dirac points; an example is provided in Supplementary Note 4. To further generalize the findings in this paper, we may not necessarily need conical dispersion. Quadratic dispersion found around the band edges of photonic crystals also provides shrinking isosurfaces. However, this is less useful in practice because the zero group velocity at the band edge makes it difficult to obtain propagating waves^{40} in the presence of disorders. In addition, coupling into such media is difficult due to the large impedance mismatch.
As a final remark, the transport properties of electrons around Dirac and Weyl points has also been studied in the past few years^{15,16,17}. Some of the observations are consistent with the physics of photon scattering shown in this paper. Here, we explicitly show the general conservation law of cross section and its connection to the dispersion relation. We expect that similar conclusions can be drawn for both electrons and phonons. It provides useful insight to understand general scattering physics beyond Dirac and Weyl systems.
In conclusion, large resonant cross sections are of great practical importance. They are only achievable with long resonant wavelengths when the frequency approaches the DC point. As shown in this work, the dispersion, rather than the wavelength, is responsible for the cross section. As a result, Weyl points, which can achieve similar conic dispersion as that around the DC point, lead to the diverging resonant cross section at any desired frequency. The exceptionally strong resonant scattering is also accompanied by diminishing nonresonant scattering, which is also similar to that around the DC point. Since Weyl points can be realized at any frequency, we can effectively decouple the cross section and the wavelength. It opens up possibilities for tailoring wave–matter interaction with extraordinary flexibility, which also can be extended to acoustic and electronic wave scattering.
Data availability
The data that support the finding of this study are available from the corresponding author upon reasonable request.
References
 1.
Fujimoto, K. & James, J. R. Mobile Antenna Systems Handbook (Artech House, Norwood, 2001).
 2.
Nozik, A. J. Quantum dot solar cells. Phys. E LowDimens. Syst. Nanostructures 14, 115–120 (2002).
 3.
Balanis, C. A. Antenna Theory: Analysis and Design, 3rd edn (WileyInterscience, Hoboken, 2005).
 4.
CohenTannoudji, C., DupontRoc, J., Grynberg, G. & Thickstun, P. AtomPhoton Interactions: Basic Processes and Applications (Wiley Online Library, New York, 1992).
 5.
Jackson, J. D. Classical Electrodynamics, 3rd edn (Wiley, New York, 1998).
 6.
Yagi, H. Beam transmission of ultra short waves. Proc. Inst. Radio Eng. 16, 715–740 (1928).
 7.
Kosako, T., Kadoya, Y. & Hofmann, H. F. Directional control of light by a nanooptical Yagi–Uda antenna. Nat. Photonics 4, 312–315 (2010).
 8.
Curto, A. G. et al. Unidirectional emission of a quantum dot coupled to a nanoantenna. Science 329, 930–933 (2010).
 9.
Coenen, T., Bernal Arango, F., Femius Koenderink, A. & Polman, A. Directional emission from a single plasmonic scatterer. Nat. Commun. 5, 3250 (2014).
 10.
Liu, J., Zhou, M., Ying, L., Chen, X. & Yu, Z. Enhancing the optical cross section of quantum antenna. Phys. Rev. A 95, 013814 (2017).
 11.
Ruan, Z. & Fan, S. Superscattering of light from subwavelength nanostructures. Phys. Rev. Lett. 105, 013901 (2010).
 12.
Zhou, M., Shi, L., Zi, J. & Yu, Z. Extraordinarily large optical cross section for localized single nanoresonator. Phys. Rev. Lett. 115, 023903 (2015).
 13.
Miroshnichenko, A. E. NonRayleigh limit of the LorenzMie solution and suppression of scattering by spheres of negative refractive index. Phys. Rev. A 80, 013808 (2009).
 14.
Weyl, H. Elektron und gravitation. I. Z. Für Phys. 56, 330–352 (1929).
 15.
Nandkishore, R., Huse, D. A. & Sondhi, S. L. Rare region effects dominate weakly disordered threedimensional Dirac points. Phys. Rev. B 89, 245110 (2014).
 16.
Tabert, C. J., Carbotte, J. P. & Nicol, E. J. Optical and transport properties in threedimensional Dirac and Weyl semimetals. Phys. Rev. B 93, 085426 (2016).
 17.
Holder, T., Huang, C.W. & Ostrovsky, P. Electronic properties of disordered Weyl semimetals at charge neutrality. Preprint at https://arxiv.org/abs/1704.05481 (2017).
 18.
Lu, L., Fu, L., Joannopoulos, J. D. & Soljačić, M. Weyl points and line nodes in gyroid photonic crystals. Nat. Photonics 7, 294–299 (2013).
 19.
Dubček, T. et al. Weyl points in threedimensional optical lattices: synthetic magnetic monopoles in momentum space. Phys. Rev. Lett. 114, 225301 (2015).
 20.
BravoAbad, J., Lu, L., Fu, L., Buljan, H. & Soljačić, M. Weyl points in photoniccrystal superlattices. 2D Mater. 2, 034013 (2015).
 21.
Chen, W.J., Xiao, M. & Chan, C. T. Photonic crystals possessing multiple Weyl points and the experimental observation of robust surface states. Nat. Commun. 7, 13038 (2016).
 22.
Wang, L., Jian, S.K. & Yao, H. Topological photonic crystal with equifrequency Weyl points. Phys. Rev. A 93, 061801 (2016).
 23.
Lin, Q., Xiao, M., Yuan, L. & Fan, S. Photonic Weyl point in a twodimensional resonator lattice with a synthetic frequency dimension. Nat. Commun. 7, 13731 (2016).
 24.
Gao, W. et al. Photonic Weyl degeneracies in magnetized plasma. Nat. Commun. 7, 12435 (2016).
 25.
Xiao, M., Lin, Q. & Fan, S. Hyperbolic Weyl point in reciprocal chiral metamaterials. Phys. Rev. Lett. 117, 057401 (2016).
 26.
Xu, Y. & Duan, L.M. TypeII Weyl points in threedimensional coldatom optical lattices. Phys. Rev. A 94, 053619 (2016).
 27.
Oono, S., Kariyado, T. & Hatsugai, Y. Section Chern number for a threedimensional photonic crystal and the bulkedge correspondence. Phys. Rev. B 94, 125125 (2016).
 28.
Noh, J. et al. Experimental observation of optical Weyl points and Fermi arclike surface states. Nat. Phys. 13, 611–617 (2017).
 29.
Ruan, Z. & Fan, S. Temporal coupledmode theory for light scattering by an arbitrarily shaped object supporting a single resonance. Phys. Rev. A 85, 043828 (2012).
 30.
Berestetskii, V. B., Pitaevskii, L. P. & Lifshitz, E. M. Quantum Electrodynamics, 2nd edn, Vol. 4 (ButterworthHeinemann, London, 1982).
 31.
Yablonovitch, E. Inhibited spontaneous emission in solidstate physics and electronics. Phys. Rev. Lett. 58, 2059–2062 (1987).
 32.
Johnson, S. G. & Joannopoulos, J. D. Blockiterative frequencydomain methods for Maxwell’s equations in a planewave basis. Opt. Express 8, 173–190 (2001).
 33.
Liu, J., Zhou, M. & Yu, Z. Quantum scattering theory of a singlephoton Fock state in threedimensional spaces. Opt. Lett. 41, 4166–4169 (2016).
 34.
Shen, J.T. & Fan, S. Coherent single photon transport in a onedimensional waveguide coupled with superconducting quantum bits. Phys. Rev. Lett. 95, 213001 (2005).
 35.
Zhu, S.Y., Yang, Y., Chen, H., Zheng, H. & Zubairy, M. S. Spontaneous radiation and Lamb shift in threedimensional photonic crystals. Phys. Rev. Lett. 84, 2136–2139 (2000).
 36.
Sakurai, J. J. & Napolitano, J. Modern Quantum Mechanics (AddisonWesley, Boston, 2011).
 37.
Dowling, J. P. & Bowden, C. M. Atomic emission rates in inhomogeneous media with applications to photonic band structures. Phys. Rev. A 46, 612–622 (1992).
 38.
Sakoda, K. & Ohtaka, K. Optical response of threedimensional photonic lattices: solutions of inhomogeneous Maxwell’s equations and their applications. Phys. Rev. B 54, 5732–5741 (1996).
 39.
Chigrin, D. N. Radiation pattern of a classical dipole in a photonic crystal: photon focusing. Phys. Rev. E 70, 056611 (2004).
 40.
Mookherjea, S., Park, J. S., Yang, S.H. & Bandaru, P. R. Localization in silicon nanophotonic slowlight waveguides. Nat. Photonics 2, 90–93 (2008).
Acknowledgements
M.Z., L.Y. and Z.Y. acknowledge the finanical support of DARPA (YFA17 N660011714049 and DETECT program). L.L. was supported by the National key R\&D Progam of China under Grant No. 2017YFA0303800, 2016YFA0302400 and supported by the National Natural Science fundation of China (NSFC) under Project No. 11721404. L.S. and J.Z. were supported by 973 Program (2015CB659400), China National Key Basic Research Program (2016YFA0301100) and NSFC (11404064).
Author information
Affiliations
Contributions
M.Z. developed the theoretical formalism, performed the analytic calculations and performed the numerical simulations. All authors contributed to the final versions of the manuscript. Z.Y. supervised the project.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhou, M., Ying, L., Lu, L. et al. Electromagnetic scattering laws in Weyl systems. Nat Commun 8, 1388 (2017). https://doi.org/10.1038/s41467017015330
Received:
Accepted:
Published:
Further reading

Experimental observation of photonic nodal line degeneracies in metacrystals
Nature Communications (2018)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.