Optical Tellegen metamaterial with spontaneous magnetization

The nonreciprocal magnetoelectric effect, also known as the Tellegen effect, promises a number of groundbreaking phenomena connected to fundamental (e.g., electrodynamics of axion and relativistic matter) and applied physics (e.g., magnetless isolators). We propose a three-dimensional metamaterial with an isotropic and resonant Tellegen response in the visible frequency range. The metamaterial is formed by randomly oriented bi-material nanocylinders in a host medium. Each nanocylinder consists of a ferromagnet in a single-domain magnetic state and a high-permittivity dielectric operating near the magnetic Mie-type resonance. The proposed metamaterial requires no external magnetic bias and operates on the spontaneous magnetization of the nanocylinders. By leveraging the emerging magnetic Weyl semimetals, we further show how a giant bulk effective magnetoelectric effect can be achieved in a proposed metamaterial, exceeding that of natural materials by almost four orders of magnitude.

Most materials interact with light in the optical spectral range through electric polarization and are characterized by permittivity ε.The effect of magnetization in such materials, described by permeability µ, at optical frequencies is negligible, however, it can be enhanced using the metamaterial paradigm [1].In addition to the electric polarization and magnetization, materials operating in the linear regime may host bianisotropic phenomena, i.e., when magnetization can be induced by the electric component of light and polarization can be generated by the magnetic component.There are two types of bianisotropic phenomena [2,3]: reciprocal, commonly attributed to chirality, and nonreciprocal, usually referred to as the magnetoelectric, Tellegen, or magnetochiral effect.Optical chirality, occurring in various materials exhibiting reciprocity but broken parity symmetry, has been extensively studied in the literature [4].In contrast, the nonreciprocal magnetoelectric (ME) phenomena (not to be confused with the magneto-optical effect) are extremely rare and less explored as they require materials with both broken parity symmetry and reciprocity [5,6].
The existence of the ME effect was first conceived conceptually by B. Tellegen in a composite formed by particles bearing parallel or anti-parallel static electric and magnetic dipoles [7].Independently, in the physics community, the magnetoelectric (equivalent to Tellegen) effect was conceptualized by L. Landau [8] and predicted in Cr 2 O 3 by I. Dzyaloshinskii [9], followed by experimental verifications [10,11].These developments led to exten- * These authors contributed equally.
sive research on ME materials and multiferroics (ferroelectric ferromagnets) [12].As magnetism is strong at microwave frequencies, ME effects are particularly pronounced in this spectral range [13].Low-frequency ME materials are thus considered for applications in spintronics, next-generation data storage, and ultralow power logic memory devices [12].
Since the late 1990-s, several attempts to move from natural ME materials towards artificial ME composites were made using metamaterials approach.However, meta-atoms designs were still limited to the microwave frequency range [14][15][16][17].Furthermore, inside such threedimensional metamaterials, the bulk isotropic Tellegen effect necessarily vanishes unless a complex local magnetic biasing is applied.Recently, two alternative ideas were conceptualized suggesting the bulk ME effect using temporally modulated permittivity and permeability at optical frequencies [18] and a multilayer metamaterial in the out-of-plane antiferromagnetic configuration [19].However, the practical realization of both material types is highly complex, let alone in the optical regime most probably unfeasible with the current state of nanotechnology.
The observed ME phenomena at optical frequencies in known materials are negligible and thus remain elusive for realistic applications [20,21].The exception is the topological ME effect in certain topological insulators [22][23][24].This effect is linked to surface states and exhibits opposite signs on either side of the surface.Consequently, its experimental investigation using light is challenging due to the canceling contributions from both interfaces of a thin film sample.Classical (i.e., non-topological) optical ME effect is vanishingly weak because both cyclotron and Larmor frequencies of electrons are typically in the GHz range [6].Even at cryogenic temperatures already at 1 THz, the largest reported bulk ME parameter drops to χ = 0.02 [25].Despite the absence of practical materials with a measurable bulk optical ME effect, the ME phenomena have been extensively theoretically studied over the last two decades.It was predicted that ME phenomena could lead to novel nonreciprocal light reflection and transmission [26], photonic topological edge states [27], PT-symmetric ME energy density [28], synthetic movement [29], directional dichroism [30,31], spin-dependent thermal radiation [32], persistent planar heat current [33], to name a few.Several phenomena, such as Fresnel drag for light [34] and photonic gauge fields [35], in principle, cannot occur in natural materials and require some structural engineering.
The possibility of an isotropic ME effect has particular importance for fundamental physics.It was pointed out in seminal works [36,37] that electrodynamic equations of isotropic ME materials have the same form as those of an axion medium.Axion is an elementary particle in the quantum field theory that was theoretically predicted to resolve two grand problems in modern physics, namely, the strong charge conjugation and parity (CP) problem [38], and the problem of the existence of the dark matter [39].Therefore, in materials science and condensed matter physics, there has been intense interest in exploring exotic concepts of axion electrodynamics with practical ME materials [40,41].Currently, the two suggested material classes with the ME effect for exploring axion electrodynamics are conventional crystalline [42] and topological insulators [22,43].However, in addition to the aforementioned limitations, their ME effect is non-tunable (due to the fixed crystal lattice), making the probing of the effective axion field technologically complicated.
Here we theoretically propose a three-dimensional metamaterial with an isotropic ME effect in the visible spectral range.Its meta-atom is a nanocylinder with a ferromagnetic nanodisk in a single-domain magnetic state and a high-permittivity dielectric nanodisk supporting the magnetic Mie-type resonance.The 3D metamaterial is formed by randomly distributing (positionand orientation-wise) the nanocylinders in a host medium such as water or polymer.Due to the optimized magnetic shape anisotropy the ferromagnetic nanodisks exhibit spontaneous magnetization and the ME effect (without the need for an external magnetic bias).Random orientation of the nanocylinders in metamaterial ensures the disappearance of other types of bulk bianisotropic effects while preserving the bulk Tellegen effect, enhanced by the Lorenz-Mie-type resonances [44] in the dielectric parts of the nanocylinders [45,46].Such a bias-free optical Tellegen metamaterial represents both qualitative (allowing tunability in space, time, and frequency) and quantitative (resonant enhancement) advancement towards real-izing ME phenomena in realistic materials.We demonstrate that by using conventional materials such as cobalt and silicon, we achieve two orders of magnitude improvement of the ME effect compared to other known natural materials at room temperatures.Finally, we show that by using the emerging magnetic Weyl semimetals as the component of the meta-atoms in a metamaterial, we can enhance the ME effect almost four orders of magnitude compared to natural materials, while making bulk ME parameter the same order as the effective permittivity and permeability of the metamaterial.

Comparison of gyrotropic, chiral, and Tellegen effects
In the approximation of linear response and weak spatial dispersion, the most general isotropic materials are characterized by the following constitutive relations that relate the displacement fields D and B to the electric and magnetic fields E and H [20]: Here, c stands for the speed of light, ε, µ, κ, and χ are the scalar permittivity, permeability, chirality parameter, and Tellegen parameter of a continuous medium, respectively.In the case of a metamaterial, ε, µ, κ, and χ in eqs.(1) are the effective material parameters.In what follows, we assume time-harmonic oscillations in the form e +iωt .Among the four scalar material parameters, only one, namely, the Tellegen parameter, corresponds to a nonreciprocal effect, that is an effect where the Lorentz reciprocity is broken [6].
We devise the Tellegen metamaterial by designing a single meta-atom, analyzing its response inside a bulk composite of equivalent meta-atoms and finding the effective material parameters of the composite.We assume a random orientation and arrangement of the meta-atoms in the metamaterial.While this eases up the future large-scale nanofabrication, it also prompts exploiting anisotropic meta-atoms providing more design degrees of freedom.Indeed, a collection of randomly oriented anisotropic meta-atoms produces an isotropic metamaterial since any anisotropic electromagnetic effects in the individual meta-atoms will be compensated in such metamaterial [2, p. 148].For simplicity, we assume that individual meta-atoms have uniaxial (transversely isotropic) symmetry.This enables a variety of electromagnetic effects (see below) and is more accessible experimentally.Without losing the generality, we select the symmetry axis to be along the z-axis.For light incident on the meta-atom along its symmetry axis (E z = H z = 0), the induced electric p and magnetic m dipole moments can p (2)   p (1)   E f (1) p (1)  p (1)   p (2)  p (2)   E b (1) m (2)   m (2)   m (2) m (1)   m (1)   m (1)    be expressed as [3] p where I is the two-dimensional unit matrix and J is the antisymmetric two-dimensional (skew-symmetric) matrix with only nonzero components J yx = −J xy = 1.We limit our discussion to the dipolar moments since the higher-order multipoles in individual meta-atoms do not contribute to the bulk magnetoelectric effect when they are randomized in a metamaterial.[47, p. 142, 227].In eqs.( 2 Here, the orientation of the external static magnetic field (or the material own magnetization) is assumed to be along the +z-direction.Such double gyrotropic response is referred to as bi-gyrotropic [50, Sec.2.3].Polarizabilities α χs = α xx χ = α yy χ and α χa denote the symmetric (isotropic) and antisymmetric (anisotropic) Tellegen response of the meta-atom, whereas α κs = α xx κ = α yy κ and α κa describe its isotropic and anisotropic chirality [51].
Since the isotropic Tellegen and the other three mentioned effects, described by α ge , α gm , and α κs , show up only in the cross-polarized scattered fields, we analyze separately their characteristic signatures.Using eqs.(2), we qualitatively illustrate in Fig. 1 how the four selected effects influence light scattering by a uniaxial meta-atom.Light scattering with the same (co-) polarization as that of the incident light is not shown for clarity.In each subfigure, we depict dipole moments induced by the incident fields E i and H i and their radiation patterns for two scenarios: light incident along the −z-direction marked with a superscript ' (1) ' (red arrows) and the +z-direction marked with a superscript ' (2) ' (blue arrows).Subscripts 'f' and 'b' refer to the forward and backward scattering, respectively.The black arrows denote the wave vectors k of the incident and scattered waves.As is seen in Fig. 1(a), gyroelectric response of a meta-atom results in a cross-polarized (with respect to the incident electric field) electric dipole that radiates symmetrically crosspolarized light in the backward (magneto-optical Kerr effect) and forward (magneto-optical Faraday effect) directions.Notably, the light scattering is independent of the incidence direction.According to Fig. 1b, a gyromagnetic meta-atom, likewise, scatters cross-polarized light with the same intensity into backward and forward directions, but the phase of light now depends on the illumination.
Chirality in a meta-atom is responsible for reciprocal cross-polarization conversion in transmission (Fig. 1c).Here, the radiation patterns for the two illuminations are oppositely oriented but both correspond to a Huygens'type pattern with no backward and maximum forward radiation [52].The sense of polarization rotation (clockwise or anti-clockwise with respect to the propagation direction) is independent of the illumination direction confirming reciprocity of a chiral meta-atom [3].
Finally, a uniaxial meta-atom with the Tellegen response behaves somewhat oppositely to a chiral metaatom: regardless of the illumination direction, the induced dipole moments do not scatter light in the forward direction (Fig. 1d).Instead, they scatter asymmetrically backwards.Note that this does not violate the optical theorem [53] since the zero forward scattering occurs for the cross-polarized portion of scattered light.Thus, the nonreciprocity in a Tellegen meta-atom manifests itself in the cross-polarized reflection but does not appear for transmitted light.Compared to gyroelectric and gyromagnetic effects where both Kerr and Faraday rotations occur, the Tellegen effect hosts only the former.Due to this feature, even bulk Tellegen material slabs do not produce cross-polarized light in transmission and can be probed only in reflection [20].
By comparing eqs.( 1) and ( 2), we observe that in a bulk isotropic metamaterial with randomly oriented meta-atoms with statistical averaging of their scattered fields only the effects characterized by α e , α m , α χs , and α κs survive [2, p. 148].Therefore, to design a purely Tellegen metamaterial, we must ensure that all the metaatoms have zero α κs polarizability.

Design of a Tellegen meta-atom
The challenge to reach an isotropic Tellegen metamaterial is to produce a space-inhomogeneous external magnetization (i.e., each meta-atom should be magnetized externally) with a sufficiently strong field amplitude.This is a bottleneck of all the known microwave designs of Tellegen meta-atoms [14][15][16][17].If the external magnetization is uniform in space, the bulk Tellegen effect would be fully compensated because magnetization in each metaatom will be random with respect to its geometry.In other words, isotropy of the metamaterial forbids uniform external bias fields.Another challenge is that at optical frequencies the magneto-optical effects exploited for the generation of the effective Tellegen response are inherently very weak [6] since their resonances occur in the microwave region for a typical magnetic flux density in the order of a few Tesla.
In order to address both challenges, we propose a metaatom with the geometry shown in Fig. 2a.The proposed meta-atom is a nanocylinder with one part ferromagnet or ferrimagnet (ferrite) and one part highpermittivity semiconductor or dielectric.As an example, we select cobalt as a ferromagnet having one of the largest magneto-optical parameters Q = |ε xy /ε xx | ∼ 0.02 in the visible frequency range [54].The ferromagnet must be in a single-domain state [55,Sec. 7], which puts a constraint on its size and aspect ratio δ = h Co /D.Although nanoscale ferromagnets with high aspect ratios are superior for attaining the single-domain configuration, they are usually challenging to fabricate.We thus keep δ = 1.There is an upper limit on the diameter of the ferromagnetic nanocylinder in a singledomain state since the magnetostatic energy scales up faster with its size (proportional to its volume) than the energy required to create a domain wall (proportional to the nanocylinder cross-section) [56, Sec.17.5].For a chosen cobalt nanocylinder the maximum diameter is D max = 96 nm (see Sec. 1 in Supplementary Information [57]).There is also a lower size limit as the thermal energy k B T might become sufficiently large to flip the magnetization produced by the magnetocrystalline and shape anisotropies in the smaller nanocylinder [58] (here k B is the Boltzmann constant and T is the temperature of the nanocylinder).We estimate the minimum diameter of a cobalt cylinder to be D min = 7.4 nm at the room temperature of T = 300 K (see Sec. 2 in Supplementary Information [57]).Importantly, these theoretical limits are rather close to the experiment [59].Here we choose the cobalt nanocylinder having D = D max ) to reach the optical resonance in the visible spectrum for the whole meta-atom.
The ferromagnetic nanocylinder in a single-domain state possesses permanent static magnetization close to the saturation magnetization of cobalt M s = 1.44 × 10 6 A/m.The magnetization direction (up or down) determines the sign of the effective Tellegen polarizability α χs .It can be defined upon the fabrication.Fixed magnetization along the nanocylinder axis can be obtained through the magnetocrystalline, shape, and/or the perpendicular magnetic anisotropies of the nanocylinder [55, Sec.7], [60].While the shape anisotropy is generated when the nanocylinder has a high aspect ratio (δ > 1), the magnetocrystalline anisotropy is an intrinsic property due to the cobalt lattice.We assume the magnetocrystalline anisotropy to be parallel to the nanocylinder's axis, thus, the coercive field can reach up to H c = 6.8 kOe (see Sec. 3 in Supplementary Information [57]).Such a coercive field ensures that the magnetization in each meta-atom is securely fixed to its geometry, which in turn means that a randomly oriented mixture of the meta-atoms will have an uncompensated bulk Tellegen response.Addressing the generally weak magneto-optical response at optical frequencies, we note that in the singledomain state, the magnetization is the highest M ≈ M s , leading to the highest magneto-optical effect for a given ferromagnet [6].Additionally, the bottom nanocylinder made of a high-permittivity semiconductor (crystalline Si, Fig. 2a) allows us to substantially enhance the Tellegen polarizability α χs of the meta-atom since it operates at the magnetic Mie-type resonance [52].
To provide a qualitative picture of the magnetoelectric effect in the meta-atom, we excite it by a standing wave with the electric field polarized along the y-axis and with the node of the magnetic field located at the center of the meta-atom (Fig. 2a).Note that such illumination has a different wavevector direction k i compared to the one in Fig. 1(d).We choose it here for simplifying the analysis as it minimizes the effects of spatial dispersion in the meta-atom.While in the dipolar approximation, there are several moments excited in the meta-atom by the electric field, in our discussion and in the figure, for clarity, we focus only on the moments that lead to the ME effect.Due to the antisymmetric (gyroelectric) part of the nanocylinder permittivity tensor, the incident electric field induces an oscillating dipole moment p orthogonal to both incident field and magnetization M s [6].This dipole moment generates the curl of the high-frequency magnetic field inside the semiconductor part of the nanocylinder according to Ampere's law.This magnetic field, in turn, induces a non-zero optical magnetic dipole moment m inside the semiconductor nanocylinder.The chosen size of the semiconductor nanocylinder supporting magnetic Mie-type optical mode ensures this induced magnetic dipole moment is brought to the resonance.To support this qualitative description, we plot the magnetic field distribution (complex amplitude) at the cross-section of meta-atom at resonance (λ = 750 nm) (Fig. 2a) obtained via full-wave simulations.A magneticdipole-like field distribution indeed emerges at the interface between the two nanocylinders.The Mie-type resonant mode in Si is hybridized due to the close proximity of the cobalt nanocylinder.Thus, the incident ypolarized optical electric field induces the y-directed optical magnetic moment in the meta-atom, which is the essence of the optical ME (Tellegen) effect.Conversely, an incident optical magnetic field would generate an electric dipole moment.
The semiconductor nanocylinder of the meta-atom has the additional role of breaking the parity symmetry that a pure cobalt nanocylinder would have.The magnetic point group of the meta-atom C ∞v (C ∞ ) [61] includes only the symmetry operations m x , m y , and 4 z .Here, m x denotes a mirror symmetry with respect to the x axis followed by a time reversal and 4 z denotes a fourfold rotational symmetry [62].Thus, the meta-atom has a parity-time (PT) symmetry required for Tellegen metamaterial constituent.Due to its magnetic point group, the meta-atom has only 9 non-zero independent polarizability components: α e , α zz e , α m , α zz m , α ge , α gm , α χs , α zz χ , and α κa (see Sec. 4 in Supplementary Information [57]).The meta-atom cannot possess chirality due to the symmetry, that is α κs = 0.Moreover, as discussed above, the gyrotropy effects will disappear upon statistical averaging of the scattered fields in a metamaterial with randomly oriented meta-atoms.

Bulk optical Tellegen metamaterial and giant Tellegen response
By using the full-wave simulations, we calculate the traces of the electric, magnetic, and magnetoelectric polarizability tensors.Next, we use the Maxwell Garnett mixing rule to determine the macroscopic effectivemedium parameters of a bulk Tellegen metamaterial (Fig. 2b) consisting of a random arrangement of the designed meta-atoms (see Methods and Sec. 5 of Supplementary Information [57]).For simplicity, we assume the host medium to be air (a more practical scenario with a dielectric host medium would have qualitatively similar results with the only change of increased effective permittivity).We choose the volumetric concentration of the meta-atoms N = 1.45 × 10 20 m −3 , which is equivalent to the volume fraction N V = 0.15 (ratio between the volume occupied by the meta-atoms and the total volume of the metamaterial).Figures 2(c)-(d dominant electric polarizability α e (see Sec. 5 of Supplementary Information [57]) and relatively high density of the meta-atoms.In contrast, µ eff ≈ 1 due to the small magnetic polarizability of the meta-atoms.Although it is resonant, the peak appears flattened when compared to the relative permeability of air.Likewise, the effective Tellegen parameter is resonant and reaches χ eff ∼ 10 −3 .Such a high value is two orders of magnitude greater than the observed bulk ME effect in natural materials at room temperatures [63,64].We stress that the proposed ME metamaterial has important fundamental advantages over natural materials.By adjusting its meta-atom geometry and constituent materials we are potentially able to tune its resonant wavelength within the optical range.Moreover, the possibility of structuring the metamaterial (making it nonuniform in space and/or time) has important implications for generating and analyzing the effective nonuniform axion field [40,41] and related effects [37,43,[65][66][67].Nonuniformity in the axion field is essential since this field enters into the Maxwell equations only through derivatives with respect to time or space [68].
Although the cobalt-based Tellegen metamaterial (Fig. 2b) readily provides a record-high bulk ME response, we further extend our design approach to other magneto-optical materials.Specifically, the anomalously strong magneto-optical response was recently predicted and experimentally confirmed in the so-called magnetic Weyl semimetals [69,70].Weyl semimetals exhibit con-duction and valence bands that touch at discrete points (Weyl nodes) in momentum space, i.e., the Brillouin zone.Each Weyl node acts as a source or sink of Berry curvature, a quantity in quantum mechanics that describes the phase evolution of quantum states.Weyl nodes come in pairs of opposite chirality (analogous to the concept of 'handedness') and are topologically protected [71], meaning that they cannot be annihilated without a significant change in the system's properties.In the case of magnetic Weyl semimetals, the Weyl nodes of opposite chiralities are at the same energy but separated from one another in the momentum space.Weyl semimetals have been recently put forward as a means to generate compact optical isolators, orbital angular momentum detectors, and nonreciprocal thermal emitters among many others [68].
We consider the topological semimetal EuCd 2 As 2 experimentally characterized in [72].For the calculation of its antisymmetric dielectric tensor, we use the same material parameters as in [73].The frequency dispersion of the diagonal and off-diagonal permittivity components of the Weyl semimetal is plotted in Sec.6 of Supplementary Information [57] at room temperature.For simplicity, in our qualitative analysis, we do not consider the Fermi arcs existing in Weyl semimetals (unique surface states that connect the projections of different Weyl points onto the surface Brillouin zone).The giant magneto-optical response of the Weyl semimetal occurs at the mid-and far-infrared frequency ranges.The ge-ometry of the designed Weyl-based Tellegen meta-atom is shown in Fig. 3a.Similarly to the cobalt-based metaatom, it is built out of the upper magneto-optic (Weyl semimetal) and lower high-permittivity (amorphous silicon [74]) nanocylinders.The diameter D Weyl = 530 nm of the Weyl nanocylinder with unit aspect ratio is chosen under the upper limit (around 1 µm) for the singledomain magnetization in similar semimetals [75].Figure 3a also features the optically-resonant magnetic field distribution (complex amplitude) when the meta-atom is illuminated by an incident standing wave.Similarly to the case in Fig. 2a, one observes the magnetic-dipole-like field distribution.Next, we determine the six polarizability components of the Weyl meta-atom (see Methods and Sec.7 in Supplementary Information [57]) and, using the Maxwell Garnett mixing rule, find the effective material parameters of Weyl-based mid-infrared Tellegen metamaterial.The meta-atoms are randomly oriented and distributed within the host medium (air) (see Fig. 3b).We choose the volumetric concentration of the meta-atoms N = 1 × 10 18 m −3 , equivalent to the volume fraction of N V = 0.12.The effective material parameters of the metamaterial are depicted in Fig. 3c-e.We note the presence of multiple Lorentzian resonances in all the parameters as the considered frequency range is below the plasma frequency of the Weyl semimetal.Due to the anomalous magneto-optical response of the semimetal, the Tellegen bulk parameter χ eff at 5.2 µm reaches the value of 0.36, four orders of magnitude higher than any known material.Remarkably, at the same wavelength, the effective permittivity and permeability are weakly resonant at ε eff = 1.28 − i1.04 and µ eff = 0.96 − i0.05, comparable in amplitude to the Tellegen parameter.This giant ME parameter potentially elevates the utility of Tellegen media, ushering it into practical applications.
To demonstrate one of the implications of the giant Tellegen response, we consider the magneto-optical Kerr effect for plane waves incident on the planar interface between air and designed Weyl-based metamaterial [6].As was discussed above, the Tellegen effect leads to the Kerr rotation and ellipticity that are independent of the illumination side of the isotropic metamaterial.Figure 3f depicts both quantities in the studied wavelength range.The Kerr rotation and ellipticity follow the Kramers-Kronig relation and exhibit resonance at the wavelength of 5.2 µm, reaching the record values in the considered optical range (for typical magneto-optical materials, the Kerr rotation does not exceed 1 deg in the optical frequency range [50,Sec. 1.4]).

Interplay between Tellegen and gyrotropic responses
Next, we analyze the contribution of the Tellegen and gyrotropic effects in the optical Tellegen metamaterial and visualize how different structural blocks in the designed meta-atom affect the Tellegen response.We do this on the example of Weyl-based Tellegen meta-atoms (Fig. 3a).As mentioned above, the gyrotropic effect present in each meta-atom cannot contribute to the macroscopic gyrotropy of the metamaterial due to meta-atoms' random orientations.Next, we confirm this statement.Since the full-wave simulations of a threedimensional metamaterial are very challenging (the simulation domain is large, and one needs a large set of such simulations for obtaining meaningful results using statistical averaging), we instead analyze a single layer of metaatoms (metasurface).We further simplify the problem by assuming a periodic arrangement of the meta-atoms along the x-and y-directions with the period of 2 µm (see Methods).Each supercell includes four meta-atoms (Fig. 4a).By changing the mutual orientation of the four meta-atoms within the supercell, we emulate both oriented and random arrangements of the same meta-atoms as in a bulk metamaterial.Following the results shown in Fig. 1, we can define the two auxiliary parameters R χs and R ge that allows us to separate the Tellegen and gyroelectric (same applies to gyromagnetic) effects from one another.Although both effects lead to the crosspolarized reflection, the phases of such reflections for two opposite illuminations (1) and ( 2) are equal in the case of gyrotropic and opposite in the case of the Tellegen effect.Therefore, we define the gyroelectric reflection parameter as R ge = R xy and R xy stand for complex-valued cross-polarized reflection coefficients from the array of meta-atoms (the reference plane is located at the geometric center of the array, i.e. at z = 0) when the incident wave propagates along the −z-and +z-directions, respectively.
First, the four meta-atoms are all oriented along the z axis (Fig. 4a).This is equivalent to an oriented (anisotropic) three-dimensional Tellegen metamaterial.We then calculate the reflection parameters R ge and R χs for the forward and backward incidences (Fig. 4a).In this case, in addition to the Tellegen, the array exhibits an uncompensated and dominant gyroelectric response.However, by making the metamaterial become isotropic through randomizing its meta-atoms (having an equal number of meta-atoms looking up and down inside the supercell), the gyroelectric effect R ge vanishes, while the Tellegen response R χs is nearly unchanged due to the local magnetization in the meta-atoms (Fig. 4b).Thus, such an isotropic configuration represents a purely Tellegen metamaterial.
Next, we check how the metamaterial behaves if the magnetization in each meta-atom is independent of its orientation, i.e. equal for all Tellegen meta-atoms.This occurs when the meta-atom is not in the single-domain magnetic state and the magnetization can be viewed as global (Fig. 4c).Although each meta-atom has the Tellegen response, it vanishes for the metamaterial as a whole (R χs = 0).This explains why despite that the design of individual Tellegen meta-atoms is known for a long time [14][15][16][17], bulk Tellegen metamaterials based on them could not be proposed.However, the single-domain magnetic state in a meta-atom alone is not enough for the bulk Tellegen response.Additionally, we need to break the inversion (parity) symmetry of the meta-atom.Here we have a Si nanocylinder on one side of the meta-atom that does this.Once we remove the silicon nanocylinders, the Tellegen response disappears both in oriented (see Fig. 4d) and random arrangements (see Fig. 4e).

DISCUSSION
We proposed a three-dimensional (bulk) Tellegen (ME) metamaterial with an isotropic ME effect in the optical regime.The individual meta-atoms exhibit spontaneous magnetization and, therefore, do not require an external bias field.We further demonstrate that using an emerging class of magnetic Weyl semimetals, the effective ME response can reach values on par with those of permittivity and permeability.Realized experimentally, such a giant ME response would enable developments and applications in the field that so far has remained purely theoretical [26][27][28][29][31][32][33].The proposed concept can be further extended to create anisotropic Tellegen metamaterials, in particular, to the "synthetically moving" media [3,30,34,35], and enable their implementation in the optical regime.An important fundamental implication of the proposed Tellegen metamaterial is in the electrodynamics of hypothetical axion matter.The axion electrodynamics is characterized by an additional term L θ = θ(r, t)e 2 E • B/(2πh) appearing in the Maxwell Lagrangian [68] (θ is the axion field, e is the electron charge, and h is the Planck constant).Applying the duality transformations to these equations, one can obtain the conventional Maxwell equations, but the constitutive relations get the new form, precisely as in Eq. ( 1) [68,76].This circumstance provides us with the opportunity to explore light propagation in arbitrary axion fields θ(r, t) by exploring light propagation in the Tellegen metamaterial with the Tellegen parameter χ(r, t) that can be straightforwardly tuned experimentally by design.This is in sharp contrast with recently proposed condensed-matter platforms emulating axion electrodynamics.These materials can model only a constant axion field due to their fixed crystalline structures (topological and conventional insulators [22,42,43]) or the axion field with a specific linear dependence on time and coordinate (bulk Weyl semimetals [68]).This obviously limits their ability to fully reproduce the behavior of a true axion field.Therefore, the proposed class of designer Tellegen metamaterials can offer novel avenues towards the first experimental implementations of various concepts related to axion matter, including dyon quasiparticles [65], anyon statistics [65], the Witten effect [37], and axionic polaritons [43,66].

Polarizabilities calculation
We calculate the polarizabilities of the Tellegen metaatoms depicted in Figs.2a and 3a using the full-wave simulations.The frequency dispersions for cobalt (diagonal and off-diagonal components) and crystalline silicon are taken from the experimental data in [54] and [77], respectively.We employ the rigorous approach based on the calculation of the high-frequency polarizations and magnetizations induced in the meta-atom under given excitation [78].To independently determine different polar-izability components, the meta-atom is illuminated by a set of specific standing waves [79].The calculated polarizabilities can be found in Secs. 5 and 7 of Supplementary Information [57].By calculating the traces of the three polarizability tensors 2α e +α zz e , 2α m +α zz m , and 2α χs +α zz χ and using the Maxwell Garnett mixing rule [2, Sec.7.2.2],we determine the macroscopic effective-medium parameters of a bulk Tellegen metamaterial consisting of the designed meta-atoms.

Simulation
The responses of the two-dimensional arrays of Tellegen meta-atoms in Fig. 4 were simulated using commercial COMSOL Multiphysics software based on the Frequency Domain Source Sweep (FDSS) solver.The peri-odic boundary conditions were used to create the effect of an infinite simulation domain for the structure.The polarization of incident light was defined in the incident port, and the reflection and transmission spectra were calculated using the scattering parameters of COMSOL.

FIG. 1 .
FIG. 1. Characteristic signatures in the cross-polarized light scattering of different electromagnetic effects in a meta-atom with a uniaxial symmetry along the z-axis.(a) Gyroelectric magneto-optical effect.The effect can be obtained at microwaves using meta-atoms with magnetized electron plasmas and at optical frequencies using ferromagnetic meta-atoms.(b) Gyromagnetic magneto-optical effect that can be induced in the meta-atom at microwaves by exploiting ferrites.(c) Chirality effect that can be obtained in reciprocal meta-atoms with broken parity symmetry.(d) Tellegen effect that can be generated only by nonreciprocal meta-atoms with broken parity symmetry and reciprocity.Meta-atoms in (a), (b), and (d) have static magnetizations along the +z-direction.
), α e = α xx e = α yy e and α m = α xx m = α yy m are the electric and magnetic polarizabilities describing reciprocal responses, while α ge and α gm are gyroelectric (due to the cyclotron orbiting of free electrons [48, p. 571]) and gyromagnetic (due to the precession of electron spins [49, p. 454]) polarizabilities responsible for magneto-optical effects.

FIG. 2 .
FIG. 2. (a) Uniaxial Tellegen meta-atom operating in the visible.The upper nanocylinder is a ferromagnet (Co) in the singledomain state.Saturation magnetization direction Ms defines the sign of the Tellegen polarizability.The lower nanocylinder is a high-index dielectric (Si) operating at the magnetic Mie-type resonance (hCo = 96 nm, hSi = 43 nm, DCo = 96 nm).Incident y-polarized electric field Ei in the standing-wave configuration induces, among a magnetic dipole moment m along the y-direction (for clarity, not all induced moments are shown).The field map shows the magnetic field distribution inside the meta-atom at resonance (λ = 750 nm).(b) Schematics of the bulk optical Tellegen metamaterial with cobalt-silicon meta-atoms.The yellow arrows depict the orientations of the local magnetization in single-domain ferromagnetic nanocylinders.Effective parameters of the metamaterial in (b): (c) relative permittivity, (d) relative permeability, and (e) bulk Tellegen parameter.

FIG. 3 .
FIG. 3. (a) Tellegen meta-atom for mid-infrared wavelengths.The upper nanocylinder consists of a ferromagnetic (Weyl semimetal EuCd2As2) in the single-domain magnetization state.The lower nanocylinder is amorphous silicon.The dimensions of the meta-atom are h Weyl = 530 nm, hSi = 26 nm, and D Weyl = 530 nm.Incident y-polarized electric field induces a magnetic dipole moment along the y-direction.The field map plots the magnetic field distribution inside the meta-atom at resonance (λ = 5 µm).(b) Schematics of the bulk mid-infrared Tellegen metamaterial with Weyl-silicon meta-atoms.The yellow arrows depict the orientations of the local magnetization in single-domain ferromagnetic nanocylinders.(c)-(e) Effective parameters of the metamaterial in (b).(f) Kerr rotation and Kerr ellipticity calculated for plane waves incident on the planar interface between air and metamaterial in (b).
xy , while the Tellegen reflection parameter is given by R χs = R

FIG. 4 .
FIG. 4. Periodic two-dimensional arrays of meta-atoms equivalent in their optical reflection properties to bulk metamaterials.In each illustration, only one supercell is shown.(a) Oriented and (b) effectively random arrays of Tellegen meta-atoms, equivalent to anisotropic and isotropic Tellegen metamaterials, respectively.The macroscopic Tellegen response vanishes if the magnetization in the meta-atoms is global which is the case in (c) or if the meta-atoms have inversion symmetry which is the case in (d) or (e).