Fig. 1. EFS of magnetized plasma at f=0.4  THz, for (a) B=4  T and (b) B=0.4  T. The magnetic field is along the z direction and the parameters used in the calculation are the plasma frequency ωp/2π=0.5  THz, ε∞=16, and γ=5×1010  rad/s.

Fig. 2. (a) Top: 2D dispersion relation curve for f=0.4  THz on the kx=0 plane, and the selected wavevector value as marked on the curve: (ky,kz)∈{(3,5.68),(4,6.69),(5,7.80),(6,8.97),(7,10.18)}. Bottom: schematic of the 2D cavity array. (b) In-plane electric field Ey distribution of the cavities at the resonance frequency with Ez incidence, the cavity sizes are from left to right (ly,lz)∈{(125,66),(94,56),(75,48),(63,42),(54,37)}  μm, corresponding to the wavevector coordinate in (a); the resonant mode orders used to determine the cavity sizes are (m,n)=(1,1). (c) The transmission spectra of different cavities around the resonance frequency corresponding to mode orders (m,n)=(1,1).

Fig. 3. (a) Resonance frequency shift with increasing mode number m when n is set to 1 for a cavity size of P1(200  μm,100   μm). (b) The resonance frequency shift with n increasing and m is set to 1 for a cavity size of P2(150   μm,200   μm). (c) The transmission spectrum for the case of (a). (d) The transmission spectrum for the case in (b). (e), (f) The configuration and the electric field distribution of different resonance frequencies in (c) and (d), with separate polarizations of Ey and Ez.

Fig. 4. (a) Transmission spectra for a fixed cavity size of (ly,lz)=(125 μm,66 μm) under different magnetic field strengths, with the same configuration as in Fig. 2. (b) The dependence of the cavity resonance frequency over the magnetic field for mode (1,1), where the solid line is calculated by using Eqs. (1) and (3), and the squares are obtained from the full wave simulation.