Signatures of Complex Optical Response in Casimir Interactions of Type I and II Weyl Semimetals

The Casimir interaction is induced by electromagnetic fluctuations between objects and it is strongly dependent upon the electronic and optical properties of the materials making up the objects. Here we investigate this ubiquitous interaction between Weyl semimetals, a class of 3D systems with low energy linear dispersion and nontrivial topology due to symmetry conditions and stemming from separated energy cones. A comprehensive examination of all components of the bulk conductivity tensor as well as the surface conductivity due to the Fermi arc states in real and imaginary frequency domains is presented using the Kubo formalism for Weyl semimetals with different degree of tilting of their linear energy cones. The Casimir energy is calculated using a generalized Lifhsitz approach, for which electromagnetic boundary conditions for anisotropic materials were derived and used. We find that the Casimir interaction between Weyl semimetals is metallic-like and its magnitude and characteristic distance dependence can be modified by the degree of tilting and chemical potential. The nontrivial topology plays a secondary role in the Casimir interaction of these 3D materials and thermal fluctuations are expected to have similar effects as in metallic systems.

L ight-matter interactions are at the forefront of current fundamental and applied research. The exchange of zero-point energy modes of fluctuating electromagnetic fields between objects with finite boundaries gives rise to the ubiquitous Casimir interaction 1 . This type of force is strongly dependent on the geometry of the objects as well as their electronic and optical response 2 . The Casimir force is especially pronounced at micro and submicrometer separations and it can limit the operation of nano and micro electronic devices due to unwanted sticktion and adhesion 3 . Recent developments have shown that it is beneficial to explore the materials aspect of Casimir interactions in order to probe novel physics 4 . This in turn can be of great importance to find ways to control the magnitude and sign of this ubiquitous force, which may be used to improve the performance of tiny devices.
Materials with Dirac spectrum have been a fruitful platform for Casimir force discoveries. Perhaps the most studied system in this regard has been graphene whose reduced dimensionality and linear energy band structure led to many unusual functionalities in this interaction, including much enhanced thermal effects [5][6][7][8][9][10][11][12] . Additionally, Casimir force phase transitions driven by external fields have been predicted in related 2D materials, such as silicene, germanene, and stanene 13 . Moreover, repulsive and force quantization effects have been described in other materials, such as 3D topological insulators and Chern insulators [14][15][16][17] . Weyl semimetals (WSMs) are 3D materials and they have analogous to graphene linear energy dispersion band structure near the Fermi level. WSMs are also characterized by linear energy cones, whose nodes can be viewed as magnetic monopoles in reciprocal space, associated with Berry curvature [18][19][20] . Tilting of the Weyl cones is also possible as one distinguishes between type I and type II WSMs depending on the degree of tilting [21][22][23] . Exotic surface states, which are associated with the nontrivial topology in these materials have also been demonstrated 19,20 .
All of these unique electronic structure properties have a direct consequence in the optical response of Weyl semimetals. The interplay between the 3D linear energy dispersion, separations of the Weyl cones, and their tilting results in an anisotropic optical conductivity tensor with distinct features in the longitudinal and Hall components, as well as the existence of surface conductivity due to the special surface Fermi arc states. In this work, using an effective model for the linear band structure and the Kubo formalism a comprehensive examination of the different optical conductivity components is presented for type I and II WSMs. The Casimir interaction is calculated using a generalized Lifshitz approach for which the boundary conditions for anisotropic media are resolved via the Fresnel reflection matrix. We focus on the purely quantum mechanical and purely thermal regimes, which helps understand how different features in the optical conductivity affect the interaction. Although the optical response is complex, the Casimir interaction is found to be dominated by the bulk diagonal conductivity components, while the surface states and Hall conductivity play a secondary role. Unlike 2D Dirac materials for which several unusual functionalities were found due to their non-trivial topology, the Casimir interaction in WSMs is similar to the one for metallic systems. We find that the magnitude of the interaction is strongly affected by the cutoff energy for the linear band dispersion, the chemical potential, and the cone tilting.

Results and discussion
Electronic properties of Weyl semimetals. The essential electronic properties of Weyl materials with linearly dispersing bands crossing at the Fermi level can be captured by the following Hamiltonian per Weyl node 19,24 with eigenenergies This is a minimal model with a pair of Weyl nodes specified by their chirality η = ±1. Here v F is the Fermi velocity, σ = (σ x , σ y , σ z ) are the Pauli matrices, k is the 3D wave vector, and 1 is the 2 × 2 identity matrix. The nontrivial topology corresponding to the above Hamiltonian is determined by the separation between the two crossing points in the band structure. For systems with broken Inversion symmetry this separation, denoted as b 0 , marks the distance between the Weyl points in energy space, while for systems with broken time reversal symmetry b gives this separation in momentum space. These two cases are generalized in Eq. (1) by using the four-vector b μ = (b 0 , b).
In addition to their separation, the Weyl cones can also be tilted, which is determined by the tilt velocity v t . One distinguishes between v t /v F < 1 and v t /v F > 1 cases corresponding to type I and type II Weyl semimetals, respectively. The v t = 0 situation corresponds to the simplest type I material with cones in an upright position (no tilting). In type II WSM the tilt is large enough so the cones are tipped over transforming the point-like Fermi surface, typical for type I WSM, to a different shape. Due to the v t /v F > 1 condition, there is also a finite density of states at the Fermi level characterized by electron and hole pockets. Let us note that the nontrivial topology is irrelevant of the cone tilting and it is present in both types of semimetals. Nevertheless, many properties are expected to be affected by v t and its relative value with respect to v F , giving rise to distinct features in the optical response, anisotropic chiral anomaly, and anomalous Hall effects among others [24][25][26][27][28] .
In this work we consider the separation between the cones as described by b μ ¼ ðb 0 ; bk z Þ and we choose the cone tilting with respect to the z-axis in momentum space with v t ¼ v tkz . As a result, the above Hamiltonian and its eigenergies and eigenspinors can be simplified to y q , and jM η j ¼ jK η j þ K η z . Also, Ψ = v t /v F denotes the tilting ratio. In Fig. 1a, we show representative band structures for type I WSM whose cones are separated in momentum space, which corresponds to materials with broken Time Reversal symmetry. Fig. 1b shows type II WSM with cones separated along the energy direction, which reflects systems with broken Inversion symmetry. Due to the tilted cones in type II WSM free carriers at the Fermi level are also available as schematically shown in Fig. 1c.
Optical response of Weyl semimetal. Using the eigenstates and energies from Eqs. (4, 5) the conductivity tensor for both types of WSMs is obtained using the Kubo formalism and it is found of the following form Each bulk component is decomposed into interband and intraband terms, such that σ To obtain the total conductivity corresponding to the Hamiltonian in Eq. (3) summation of the contributions from the Weyl cones must be done, such that σ = σ (+) + σ (−) . Detailed calculations and explicit expressions for zero temperature, T = 0, and no spatial dispersion, q = 0, for the Real and Imaginary parts of the interband and intraband conductivity components are given in the Supplementary Note 1.
Our results indicate that both types of materials have anisotropic optical properties with distinct diagonal components. Eqs. (3)(4)(5) correspond to the 3D bulk band structure, which is responsible for σ zz,b ≠ σ xx,b = σ yy,b . In addition, the nontrivial topology of the separated Weyl cones results in a surface conductivity (discussed later), which gives rise to an anisotropic surface conductivity tensor. The bulk components have an explicit dependence upon the cone tilting and the cutoff energy ϵ c , as found by several other authors [29][30][31][32][33] . The cutoff energy serves as a limit beyond, which the Dirac bands are not of linear dispersion anymore 34,35 . This is different than the situation in graphene materials, where the optical response is found to be independent of such a cutoff 36 . The unique role of ϵ c is thus related to the interplay between the linear band dispersion and dimensionality of the material.
The general expressions for the different conductivity components are quite complicated (see Supplementary Note 1). Thus in Fig. 2 we show the nonzero components for a single cone for both types of WSMs and their dependence on frequency. For ω → 0, Re σ ½ xx;b and Re σ ½ zz;b of type I WSM tend to zero when μ = 0 (Fig. 2a), while these components exhibit Drude-like behavior for type II WSM due to the finite carrier concentration originating from the severe tilting (Fig. 2c). The Drude-like behavior is also seen in Fig. 2e, g, in which intraband contributions from the non-zero chemical potential are included. We find that Re σ ½ xx;b and Re σ ½ zz;b exhibit a linear dependence with frequency until 2ϵ c , however there is some nonlinearity in the ðΩ L ¼ 2μ jΨÀ1j ; Ω U ¼ 2μ jΨþ1j Þ region in Fig. 2e, g 30 . The real part of the Hall conductivity shows a discontinuous behavior as ℏω approaches ϵ c . The imaginary parts of the different components are also given as a function of frequency. The Im σ ½ xx;b and Im σ ½ zz;b are discontinuous as ℏω approaches ϵ c (Fig. 2b, d, f, h). At small frequency the imaginary parts are mainly determined by the Drude-like behavior 29,31,32 .
In addition to the graphical representation of the WSMs optical response in Fig. 2, analytical expressions in different limits for the plasma frequencies and Hall response can be found. Specifically, in both cases the intraband conductivities are of Drude type with plasma frequencies ω I;II P explicitly dependent on the WSM type, as shown in Supplementary Note 1. In addition, the plasma frequencies for WSM type II are found to diverge as a function of the cutoff, consistent with other works 29,32 . For example, in the limit of large ϵ c the xx-component is obtained as whereμ ¼ jμ À hb 0 j. Similar second power law and logarithimic divergencies upon ϵ c are present for the zz-component of the bulk intraband conductivity. It is also interesting to examine the behavior of the plasma frequencies as a function of chemical potential. For type I WSM limμ !0 ω I P;xxðzzÞ ¼ 0, which shows vanishing intraband conductivity as the chemical potential passes through the Fermi level. For type II WSM, one finds that Thus even ifμ ¼ 0, there is a nonzero intraband conductivity. This correlates with the severe tilting in the WSM from Ψ > 1 which results in a band crossing at the Fermi level, meaning that there are always free carriers available in the system (Fig. 1c). We further study the Hall conductivity for each type of WSM. The zero frequency limit, which is a measure of the anomalous Hall effect, is obtained as The first term in the above expressions is common to both types of WSMs. It is proportional to the distance b between the cones in momentum space and it reflects the nontrivial topology of the materials. The second term in each σ IðIIÞ; η xy comes from the tilting of the cones and the chemical potential scaled by the separation in energy spaceμ ¼ jμ À hb 0 j and it does not have topological origin. The Hall conductivity of type II WSM also shows a logarithmic divergence with the energy cutoff, while σ I;η xy is finite at the large ϵ c limit. Our results are consistent with calculations reported by other authors 31,33 . Let us note, however that the total Hall conductivity for both types of Weyl semimetals is σ IðIIÞ xy ðω ¼ 0Þ ¼ αc_b 4π 2 as the second terms in Eqs. (8,9) cancel out upon addition due to the opposite in sign cone chiralities.
Surface conductivity due to Fermi arcs. One of the hallmarks of WSMs is their Fermi arc surface states. Other materials support Fermi surface states that have circular or deformed circular shapes 37,38 . In a Weyl semimetal occupying semi-infinite space with an abrupt boundary, however, only half of this usual shape resides on the surface giving rise to a Fermi arc, which essentially starts and ends at the surface projections of the bulk Weyl points, as schematically shown in Fig. 1a. The Fermi arcs are also a signature of the nontrivial bulk topology of WSMs regardless of tilting and their existence is determined by the bulk-surface correspondence [39][40][41][42][43] . In the case of Weyl cones separated in momentum space the Fermi arc can be represented as a single Dirac fermion (similar to graphene), whose conductivity in the low frequency is of Drude-like type and it is determined by the distance between the cones 44 . Specifically, for the case of b μ ¼ ð0; bk z Þ at the x = 0 surface one finds there is only one nonzero conductivity component At the y = 0 boundary, the nonzero surface conductivity is σ xx,s with the same value as the above equation. There is no surface conductivity at the z = 0 boundary since the projection of the vector b μ ¼ ð0; bk z Þ is zero in this case 44 . Thus this optical response is associated with the WSM boundary due to the distance between the cones b and it is highly anisotropic.
Optical response at imaginary frequency. For the Casimir interaction the conductivity components in imaginary frequency domain ω = iξ are needed. Thus it is convenient to use the Kramers-Krönig relations for which, Realizing that Im½σ ij ðiξÞ ¼ 0, the analytical continuation of Re½σ ij ðiξÞ to all positive imaginary frequencies in the upper half of the complex plane gives the final expressions for the conductivity components. By taking ξ → ξ + ℏ/τ finite dissipation can be included in the optical response (Γ = ℏ/τ).
In Fig. 3 results for some functional dependences of the conductivity components (shown explicitly in Supplementary Note 2), commensurate with the ones in Fig. 2, for a single cone at imaginary frequency ω = iξ are given. In all cases σ inter xx;b and σ inter zz;b increase linearly in the small ξ region (Fig. 3a-d). The Hall conductivity has a very weak dependence upon ξ and in both WSMs the ξ = 0 limit is governed by the analytical expressions in Eqs. (6,7). There is a discontinuous jump at Ψ = 1, which marks the transition between type I and type II WSMs as shown in Fig. 3e. As Ψ increases for type II WSM the conductivity is determined primarily by the intraband contributions due to increasing the number of free carriers from the severe tilting even if μ = 0. Also, σ ij,b as a function of μ is depicted for type I WSM in Fig. 3e at ξ = 0. This behavior is consistent with our analytical considerations and it shows that the response is mainly determined by the intraband contributions as the chemical potential grows.
Many of the topologically nontrivial features in Weyl semimetals are related to the separation of their energy cones in momentum space. In order to better understand this issue in the context of the Casimir interaction, we consider some limiting cases with respect to b μ ¼ ð0; bk z Þ. The so obtained results are then compared with the numerical calculations, in which the full optical response of the materials is taken into account.
Casimir interaction for coinciding Weyl cones. At this point, we consider non-tilted Weyl cones with b = 0. Such a situation with The legend given above the panels denotes the total σ xx,b (full black) and σ zz,b (full red), interband contributions σ inter xx;b (dashed black) and σ inter xx;b (dashed red), the bulk Hall conductivity σ xy,b (full green) and the topological Hall conductivity σ T (dashed gray). The following parameters were used in the calculations: ϵ c = 10 eV, coinciding cones corresponds to a 3D Dirac semimetal 21 . In this case, one finds that the interband conductivity components in the imaginary frequency domain are of the following form This isotropic response simplifies the reflection matrices entering in Eq. (20) and some asymptotic behavior of the Casimir energy can be found analytically. For example, the quantum mechanical (T = 0) Casimir energy in the limit of large cutoff frequency is obtained as the energy for perfect metals, such that E M A ¼ À π 2 _c 720d 3 . Since ϵ c is finite, the limiting cases of large and small separations with respect to the characteristic distance d 0 ¼ _c , and Ei(x) is the exponential integral evaluated at x. Thus the small d limit results in a d −2 scaling law for the energy. Due to the weak logarithmic dependence on d in X, however, the large d limit is essentially similar to the one for perfect metals with a d −3 scaling law in the energy.
In Fig. 4a we show the numerical calculations for the quantum mechanical interacting energy of materials with optical response from Eq. (12). These results indicate that the energy cutoff practically determines the characteristic behavior of the Casimir energy through the distance d 0 , which is in agreement with Eqs. (13,14). Specifically, for ϵ c = 1 eV one finds that d 0 = 100 nm, while for ϵ c = 5 eV, d 0 = 20 nm. Since the quantum mechanical Casimir interaction is domated by Matsubara frequencies with ℏξ < ϵ c , the onset of a metallic-like interaction is strongly related to the range of validity of the linear band structure approximation controlled by ϵ c . For separations smaller than d 0 one needs to include bands beyond the linear energy approximation studied here. Such an explicit dependence upon the cutoff energy is not present in 2D materials with Dirac energy spectra, where the optical response and Casimir energy are independent of the bandwidth of the linear bands, except at very small distances 13,15,36 .
The role of cone separation in the interaction. The parameter b enters into the reflection matrices in a complicated way, which makes it difficult to obtain analytical expressions for the Casimir energy for the considered materials. To get some idea about the role of the cone separation, we consider the situation when there is no surface conductivity and the bulk conductivity tensor contains only off-diagonal components with σ η xy;b ¼ σ T ¼ e 2 b 4π 2 , which capture the anomalous Hall effect Eqs. (8,9). We find that in the limiting case of small distance, the Casimir energy is where b 1,2 are the cone separation parameters for the two interacting materials. In the large distance limit, the energy for perfect metals is obtained. Equation 15 shows a d −1 scaling law in the short distance separation and the α 2 dependence in the numerator indicates that the energy is of rather small magnitude. It is also recognized that the product b 1 b 2 can be positive or negative given that the cone separation can be in positive or negative domains in momentum space 21,22 . This behavior is reminiscent of the Casimir interaction in other topological materials, whose response is dominated by the Hall conductivity [13][14][15] . The positive sign of the product of the cone separation parameters, which are essentially proportional to the constant anomalous Hall conductivity, is associated with repulsion as found in 3D topological insulators and 2D Chern insulators. However, since b 1,2 are continuous, the Casimir energy is not quantized as found in other topological materials due to the discrete nature of their Hall conductivities. Additionally, the Casimir energy in topological and Chern insulators has d −3 scaling law, while in Weyl semimetals there is a crossover in distance dependence from d −1 at small separations to d −3 at large separations.
The numerical results for the interaction energy when only the bulk Hall conductivity σ η 4π 2 is taken into account are shown in Fig. 4b. The repulsive interaction at small separations is dictated by the sign of the b 1 b 2 product according to Eq. (15), while at large separations the energy approaches the limit of perfect metals, which is in agreement with the above discussed analytical expressions. The strongest repulsion is found for b ⋅ d (×10 3 )~1/2, thus characteristic values b = 5-20 nm −1 correspond to d = 0.1-0.025 microns. Nevertheless, this repulsive interaction is at least two orders of magnitude smaller than E M in the sub-micron and micron separation region (depending on b values), as suggested by Fig. 4b. In Fig. 4c, d the numerical results for the normalized Casimir energy as a function of separation calculated by taking into account the full optical conductivity tensor of dissipative materials are shown for two values of the cutoff energy. The corresponding energies obtained when only the diagonal components of the conductivity with and without dissipation are taken into account, are also shown. These results reveal the surprising role of dissipation: while the low and high d limits in E are best captured for calculations with the diagonal conductivity (no Hall components), the low d limit is better described when Γ = 0, while the higher d range is better described when dissipation is finite. Figure 4c, d further shows that the role of cone separation is to bring E closer to E M in the intermediate distance range and and all repulsive effects are washed out by the dominant diagonal components of the optical response.
At this point, one might ask how the surface conductivity due to the Fermi arcs affects the Weyl Casimir interaction. For the chosen b μ ¼ ð0; bk z Þσ zz,s = 0 as discussed earlier, thus to probe this property WSMs separated along the y-axis or x-axis must be considered with the appropriate Fresnel matrices given in the Supplementary Material. Although the numerical calculations are technically difficult we find that the surface conductivity plays even a lesser effect in the Casimir energy as compared to the one of the Hall conductivity as discussed above. The surface conductivity increases the metallic-like nature of the interaction, but given the smallness of the effect further results are not given.

Casimir interaction and cone tilting.
In what follows, we further investigate how other characteristics affect the interaction between Weyl materials. More specifically, by utilizing the full optical response dissipative type I and type II WSMs separated along the z-axis are considered and the Casimir energy upon cutoff energy, degree of tilting, and chemical potential is calculated. For this purpose, the generalized Lifshitz formula in Eq. (20), while the optical response and Fresnel matrices can be found in the Supplementary Material. These results help us find ways to modulate the Weyl Casimir interaction as well as recognize signatures of the optical response and compare with the above discussed analytical expressions.
The scaling dependence shown in Fig. 5a indicates that Casimir energy increases as d is increased, although for the case of Ψ = 0, μ = 0 this increase is at a much smaller scale as compared to the other displayed cases. This behavior is reminiscent of the interaction between Drude-like metals controlled by the magnitude of their plasma frequencies. The larger the plasma frequencies, the closer the interaction is to the one for perfect metals. In the case of larger Ψ and μ, the intraband conductivities dominate the Casimir interaction due to their large ω P , while for Ψ = 0 and μ = 0 the optical response is dominated by the interband components and E/E M is a small fraction.
The behavior of E∕E M vs. d for Ψ = 0, μ = 0 in Fig. 5a is actually very similar to the analytical results of the Casimir energy for coinciding Weyl cones whose d 0 = 20 nm when ϵ c /E b = 5 (Fig. 4a). Indeed, E∕E M at large distances (d ≫ d 0 ) is practically determined by the prefactor in Eq. (13), which indicates that type I WSM with μ = 0 behaves as a Dirac material with coinciding cones in the displayed region of distance separation. For the other cases in Fig. 5a, however, the nonzero plasma frequency elevates the role of the intraband conductivity and the interaction is essentially similar to the one of Drude-like metals 2 . Given that the cutoff energy plays a prominent role in ω II P Eqs. (6, 7), Fig. 5b shows how the Casimir interaction changes as a function of ϵ c for type I and II WSMs. The much more pronounced enhancement of E for type II WSMs as ϵ c is increased is due to the large intraband optical conductivity and specifically the strong dependence of its plasma frequency on the energy cutoff (Eqs. (6, 7)). Figure 5c displays how the Casimir energy changes as a function of degree of cone tilting. The change of an almost flat E∕E M to a linear-like behavior at Ψ~1 signals a transition in the Casimir energy from type I to type II WSM. Again, the much stronger interaction for type II WSMs as compared to type I WSMs in Fig. 5c gives another perspective of the much enhanced role of the intraband optical conductivity for Ψ > 1. Figure 5d further shows that the chemical potential is also an effective parameter in controlling the magnitude of the Casimir energy for both types of WSMs, although the degree of tilting and distance separation can affect the degree of change. Additionally, increasing the Fermi velocity of the linear bands in Eq. (4) can increase the magnitude of Casimir energy although the characteristic behavior in terms of distance and other parameters are preserved (Fig. 5e).
Our numerical calculations indicate that the b parameter plays a rather small role and the Casimir interaction is dominated by the contributions from the diagonal components in the optical tensor. This can be easily understood in the context of the above discussions for the limiting cases of coinciding Weyl cones and materials with only Hall conductivity components. In fact, taking the case of non-tilted Weyl cones and expanding the Fresnel reflection matrices for small σ xy due to σ xy ≪ σ xx = σ zz , we find that the first order correction (which contains b) in the small distance limit between Weyl semimetals with μ = 0 is the same as Eq. (15). Given the much longer range due to the d −1 dependence and the α 2 in the numerator, it becomes clear that first order correction is much smaller than the dominant Casimir energy due to the diagonal components of the bulk conductivity response. Thermal Casimir energy. Thermal fluctuations in the Casimir energy are captured in the n = 0 Matsubara term of Eq. (20) and they are expected to dominate at room temperature at separations on the sub-micron and micron scales for many materials 2,4 . To study the thermal effects in the Casimir interaction of WSMs, we calculate the conductivity at finite T in the limit ξ → 0 + using the results in the Supplementary Note 2 and the Maldague's formula 13,45,46 . The interband conductivity is linear in ξ in the limit of small frequency and it is well approximated as where γ is the Euler's constant. The intraband conductivity at finite T and ξ → 0 has the typical Drude-like expression where ω I;II P are the plasma frequencies for type I and type II WSMs. Since the optical response at ξ → 0 + is dominated by the intraband contribution, the thermal Casimir energy is found as where ζ(s) is the Riemann zeta function. Thus the thermal effects in the WSM Casimir interaction are of the usual form of a Drude metal. Experimental measurements of the Casimir interaction in Weyl semimetals can be realized using a setup based on an atomic force microscope, which has been used for several other materials such as noble metals, ferromagnets, oxides, and graphene covered substrates 12,47,48 . Such a setup allows measuring a Casimir force whose magnitude is a fraction of the one of perfect metals and probe separations in the submicron-micron range. By taking type I and II Weyl semimetals one can potentially measure the interaction as a function of distance and get information on the role of cone tilting. By doping the materials or applying external gate voltage, one can also probe the role of the chemical potential. For distance separations in the submicron-micron range the linear band approximation dominates the Casimir interaction in most Weyl semimetals which ensures valid comparison with our theory. Due to the similarity of Weyl materials with metallic systems it is expected that thermal effects will not be important at such distances and the atomic force microscope technology at room temperature 12,47,48 will capture the quantum mechanical limit of the interaction. In this study, we present a comprehensive investigation of the bulk and surface conductivity components by distinguishing between type I and type II WSMs, which are then applied to the quantum mecahnical and thermal limits of the Casimir interaction. Our results indicate that although there is a complex interplay between the electronic and optical properties of these materials, the interaction exhibits a behaviors similar to the one of metallic-like systems due to the dominant role of the diagonal bulk conductivity components. The quantum mechanical energy can be modulated as a function of chemical potential and cone tilting as captured by the μ and Ψ dependences in the plasma frequency for type I and type II materials. The explicit ϵ c dependence in the bulk conductivity renders strong dependence of the energy upon the cutoff corresponding to the validity of the linear energy dispersion. For separations larger than the characteristic distance d 0 set by ϵ c the linear energy dispersion will dominate the interaction, while for larger than d 0 separations other energy bands are expected to be important and have to be considered in the electronic and optical properties of the materials. Unlike topological and Chern insulator, where the nontrivial topology can result in significant repulsive and even quantization effects, the anomalous Hall conductivity plays a rather small role in the quantum interaction in Weyl materials. We also find that because of the dominant diagonal bulk conductivity, thermal fluctuations are expected to affect the Casimir interaction in a similar way as the case of metallic systems. Our investigation is a testament that a thorough understanding of the fundamental electronic and optical properties of materials is necessary in order to make progress towards other research areas, such as light-matter interactions and the Casimir force as a particular example.

Methods
Kubo formalism. The optical response of a given material is directly related to its electronic structure. The dynamical optical conductivity tensor is a fundamental property, which is a necessary ingredient in understanding light-matter interactions, such as the Casimir force. The components of this tensor can be calculated using the Kubo formalism 49 . For a given cone specified by its chirality η, these are σ η ij ðω; qÞ ¼ Àie 2 Here ju λ;η k > and E λ;η k (λ; λ 0 ¼ ± correspond to the electron and hole states) are the eigenspinors and eigenenergies for the Hamiltonian in Eqs. (3)(4)(5) and v η ¼ ∇ k H η s =_ ¼ ηv F ðσ x ; σ y ; σ z Þ is the velocity operator. Also, the Fermi-Dirac distribution function is n F ðE λ;η k Þ ¼ 1=ðe ðE λ k ÀμÞ=k B T þ 1Þ and τ is the relaxation time.
Lifshitz formalism. The Casimir interaction can be calculated using the Lifshitz approach 2,4 , however for lossy materials with reflection matrices with complex frequency dependence this formalism must be applied carefully. Specifically, one must account for contributions from propagating and evanescent waves, which can arise from frequencies above and below the light cone, respectively 50 . Thus, for semi-infinite planar objects taken to be separated by a distance d along the z-axis, the interaction energy per unit area A can be written as Here R j = R j (k n , k x , k y ) for j = 1, 2 are the 2 × 2 Fresnel reflection matrices for the interacting objects, where k n ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi k 2 x þ k 2 y þ ξ 2 n =c 2 q with ξ n = 2πnk B T∕ℏ being the Matsubara frequencies (n = 0, 1, 2, … and k ∥ = (k x , k y )). Also, R 0 j ¼ R j ðÀk n ; k x ; k y Þ denotes the reflection matrix for scattering with opposite direction 17 . The generalized Lifshitz expression from Eqs. (20) is necessary for the interaction involving Weyl semimetals. From our subsequent calculations we find that the reflection matrices may not be real in some imaginary frequency ranges. This is attributed to the interplay between propagating and evanescent wave stemming from the anisotropic optical response, which motivates the generalization of the Lifshitz formalism 50 .
The reflection matrices depend explicitly on the optical response of the involved materials and they are obtained by applying electromagnetic boundary conditions for semi-infinite 3D objects with planar surfaces separated at a finite distance. Resolving the reflection matrices is technically difficult as each type of Weyl semimetal has anisotropic bulk and surface optical response as well as nonzero bulk Hall response. Nevertheless, the electric and magnetic fields can be expanded in terms of the suitable for planar geometry orthogonal functions M k ðrÞ ¼

Data availability
The data that support the plots in the paper and other findings of this study are available from the corresponding author upon reasonable request.