Negative Effective Mass in Plasmonic Systems.

We report the negative effective mass (density) metamaterials based on the electro-mechanical coupling exploiting plasma oscillations of a free electron gas. The negative mass appears as a result of the vibration of a metallic particle with a frequency of ω, which is close the frequency of the plasma oscillations of the electron gas m 2 relative to the ionic lattice m 1 . The plasma oscillations are represented with the elastic spring k 2 = ω p 2 m 2 , where ω p is the plasma frequency. Thus, the metallic particle vibrated with the external frequency ω is described by the effective mass m e f f = m 1 + m 2 ω p 2 ω p 2 - ω 2 , which is negative when the frequency ω approaches ω p from above. The idea is exemplified with two conducting metals, namely Au and Li.


Introduction
Metamaterials are recently developed artificial materials demonstrating properties that are not found in naturally occurring materials [1,2]. The domain of metamaterials covers a broad diversity of fields in physics and engineering: electromagnetics, acoustics, mechanics and thermodynamics. In metamaterials, both electric permittivity and magnetic permeability may be negative at certain frequencies [1][2][3]. Moreover, they may be tuned in a broad range of values [4]. The electromagnetic metamaterials are usually synthesized by embedding various constituents/inclusions with novel geometrical shapes and forms in some host media [2]. Various types of electromagnetic composite media, such as double-negative materials, chiral materials and omega media have been studied by various research groups worldwide [1][2][3][4].
A relatively new kind of metamaterials are acoustic metamaterials [5][6][7][8][9][10]. Acoustic metamaterial, in which both the effective density and bulk modulus are simultaneously negative, in the true and strict sense of an effective medium have been reported [5]. Acoustic metamaterials demonstrating the negative Poisson's ratio have been discussed [10]. Acoustic metamaterials demonstrate a potential to be perfect absorbers of mechanical vibrations [11] and also as materials enabling the focusing of ultrasound [7]. The present paper introduces the negative effective mass metamaterials based on mechano-electromagnetic coupling. The idea of the negative effective mass (density) acoustic metamaterials was demonstrated and discussed in [12,13]. We propose to exploit the plasma oscillations of the electron gas [14] in the development of metamaterials with the negative effective mass (density) [13,15]. The applications of the negative mass (density) materials include: acoustic tunneling through narrow channels, control of the radiation field, perfect transmission through sharp corners and power splitting as discussed in [16]. Elastic wave control and seismic wave protection with acoustic metamaterials possessing the negative mass (density) is considered in [17].

Negative Effective Mass and Plasma Oscillations in Metals
The mechanical model, giving rise to the negative effective mass effect, is depicted in Figure 1A. The core with mass m 2 is connected internally through the spring with constant k 2 to a shell with mass Materials 2020, 13, 1890; doi:10.3390/ma13081890 www.mdpi.com/journal/materials m 1 . The system is subjected to the external sinusoidal force (t) =Fsinωt. If we solve the equations of motion for the masses m 1 and m 2 and replace the entire system with a single effective mass m e f f , we obtain [12,13,15]: where ω 0 = k 2 m 2 . Clearly, when the frequency ω approaches ω 0 from above the effective mass m e f f will be negative [12,13,15]. Now consider the electro-mechanical analogy of the aforementioned model, giving rise to the negative effective mass. Consider a cubic metal particle, seen as ionic lattice m 1 . containing the Drude-Lorenz free electrons gas possessing a total mass of m 2 = m e nV, where m e = 9.1 × 10 −31 kg is the mass of electron, n is the concentration (number density) of the electron gas and V is the volume of the particle [14,18,19]. Electron gas is free to oscillate with the plasma frequency [14,15]. The mechanical model, giving rise to the negative effective mass effect, is depicted in Figure 1A. The core with mass is connected internally through the spring with constant to a shell with mass . The system is subjected to the external sinusoidal force ( ) = . If we solve the equations of motion for the masses and and replace the entire system with a single effective mass , we obtain [12,13,15]: Clearly, when the frequency approaches from above the effective mass will be negative [12,13,15]. Now consider the electro-mechanical analogy of the aforementioned model, giving rise to the negative effective mass. Consider a cubic metal particle, seen as ionic lattice containing the Drude-Lorenz free electrons gas possessing a total mass of = , where = 9.1 × 10 kg is the mass of electron, n is the concentration (number density) of the electron gas and V is the volume of the particle [14,18,19]. Electron gas is free to oscillate with the plasma frequency = [14,15]. is the plasma frequency (the left sketch). The equivalent mechanical scheme of the system (right sketch).
We exposed the entire metal particle to the external sinusoidal force ( ) = . The effective mechanical scheme of the metallic particle is shown in Figure 1B (the right sketch) and it coincides exactly, giving rise to the negative effective mass, supplied in this case by: where is the mass of the ionic lattice, is the total mass of the electronic gas and = ; it is seen that it may be negative when the frequency approaches from above. The negativity of the effective mass appears as a result of the attempt to use a single mass to represent a two mass system comprising masses , , as noted in [13]. Considering ≪ 1 yields: Core with mass m 2 is connected internally through the spring with k 2 to a shell with mass m 1 . The system is subjected to the sinusoidal force F(t) =Fsinωt. (B) Free electrons gas m 2 is embedded into the ionic lattice m 1 ; ω p is the plasma frequency (the left sketch). The equivalent mechanical scheme of the system (right sketch).
We exposed the entire metal particle to the external sinusoidal force F(t) =Fsinωt. The effective mechanical scheme of the metallic particle is shown in Figure 1B (the right sketch) and it coincides exactly, giving rise to the negative effective mass, supplied in this case by: where m 1 is the mass of the ionic lattice, m 2 is the total mass of the electronic gas and k 2 = ω 2 p m 2 ; it is seen that it may be negative when the frequency ω approaches ω p from above. The negativity of the effective mass appears as a result of the attempt to use a single mass m e f f to represent a two mass system comprising masses m 1 , m 2 , as noted in [13]. Considering It is clear from Equation (3) that the effective dimensionless mass m e f f m 1 +m 2 m e f f m 1 depends only on the ratio m 2 m 1 ; thus, it is independent on the metallic particle size. Thus, for the purposes of calculation, m 2 is taken as the mass of electron m e , and m 1 is the mass of the atom of metal (see Table 1). The dependence of the dimensionless effective mass m e f f /(m 1 + m 2 ) on the dimensionless frequency ω/ω p for two model metals Li and Au is plotted in Figure 2 (the data relevant to these metals is supplied in Table 1). The macro-scale values of the "plasma spring" constant k 2 10 2 N m are noteworthy. The dependencies of the dimensionless effective mass m e f f /(m 1 + m 2 ) on the dimensionless difference ω−ω p ω p = ∆ω ω p calculated for Li and Au are presented in Figures 3 and 4.
It is clear from Equation (3) that the effective dimensionless mass ≅ depends only on the ratio ; thus, it is independent on the metallic particle size. Thus, for the purposes of calculation, is taken as the mass of electron , and is the mass of the atom of metal (see Table 1). The dependence of the dimensionless effective mass /( + ) on the dimensionless frequency / for two model metals Li and Au is plotted in Figure 2 (the data relevant to these metals is supplied in Table 1). The macro-scale values of the "plasma spring" constant

Negative Mass and Low Frequency Plasmons in 1D Metallic Meso-Structures
The plasma oscillations shown in Figure 1 demonstrate the negative mass in the vicinity of the plasma frequency, which is on the order of magnitude of ≅ 10 Hz , which is very high. However, this frequency may be decreased very strongly for meso-structures built of thin metallic wires, as demonstrated in [20]. Depression of the plasma frequency into the far infrared and even GHZ band becomes possible due to the mutual inductance that appear in the periodic arrays built of thin metallic wires [20]. We consider the 1D lattice built of the metallic wires with diameter 2r connected with springs , as depicted in Figure 5. The effective (pseudo) density of electrons in the metamaterial lattice shown in Figure 5 is given by [20]: where n is the concentration of the free electron gas supplied in Table 1 for Li and Au.

Negative Mass and Low Frequency Plasmons in 1D Metallic Meso-Structures
The plasma oscillations shown in Figure 1 demonstrate the negative mass in the vicinity of the plasma frequency, which is on the order of magnitude of ω p 10 16 Hz , which is very high. However, this frequency may be decreased very strongly for meso-structures built of thin metallic wires, as demonstrated in [20]. Depression of the plasma frequency into the far infrared and even GHZ band becomes possible due to the mutual inductance that appear in the periodic arrays built of thin metallic wires [20]. We consider the 1D lattice built of the metallic wires with diameter 2r connected with springs k 1 , as depicted in Figure 5. The effective (pseudo) density of electrons in the metamaterial lattice shown in Figure 5 is given by [20]: where n is the concentration of the free electron gas supplied in Table 1 for Li and Au.
Materials 2020, 13, x FOR PEER REVIEW 5 of 7 Figure 5. One-dimensional lattice built of metallic wires 2r connected with springs . The separation between wires is a.
The pseudo-mass of electrons in such matrices is given by [18]: where n is the concentration of the free electron gas supplied in Table 1. The value expressed by Equation (5) is called in [20] as the "effective mass"; however, in our paper the notion of the "effective mass" is already ascribed to the mass of the vibrated element, given by Equation (1). Thus, we call the value expressed by Equation (5) the "pseudo-mass", and the effective density of electrons expressed by Equation (4) Substituting the aforementioned numerical parameters yields effective plasma frequencies of the lattices built from Au and Li wires * = 4.6 × 10 Hz; * = 5.2 × 10 Hz, which are much smaller that the aforementioned values of the "true" plasma frequencies.
The spring constants corresponding to aforementioned plasma frequencies are already small and equal k2(Li) = 2.4 × 10 −9 N/m, k2(Au) = 1.9 × 10 −9 N/m. The optical and acoustical branches of the The pseudo-mass of electrons in such matrices is given by [18]: where n is the concentration of the free electron gas supplied in Table 1. The value expressed by Equation (5) is called in [20] as the "effective mass"; however, in our paper the notion of the "effective mass" is already ascribed to the mass of the vibrated element, given by Equation (1). Thus, we call the value expressed by Equation (5) the "pseudo-mass", and the effective density of electrons expressed by Equation (4) we label as the "pseudo-density". Assuming r = 1.0 × 10 −6 m; a = 5.0 × 10 −3 m we estimate m Li 6.4 × 10 −27 kg; m Au 8.1 × 10 −27 kg. Equations (4) and (5) enable calculation of the effective pseudo-plasma frequencies ω * p for Au and Li according to Equation (6): Substituting the aforementioned numerical parameters yields effective plasma frequencies of the lattices built from Au and Li wires ω * Au p = 4.6 × 10 10 Hz; ω * Li p = 5.2 × 10 10 Hz, which are much smaller that the aforementioned values of the "true" plasma frequencies.
The spring constants k 2 corresponding to aforementioned plasma frequencies are already small and equal k 2 (Li) = 2.4 × 10 −9 N/m, k 2 (Au) = 1.9 × 10 −9 N/m. The optical and acoustical branches of the longitudinal modes propagation in the 1D lattice, depicted in Figure 5, should be elucidated. It should be emphasized that the ensembles of metallic wires, shown schematically in Figure 5, will not demonstrate simultaneously the negative mass (density) and the negative refraction effects [20,21]. This is due to the fact that the negative refraction becomes possible below the plasma frequency ω p [20,21]; contrastingly, the effect of the negative mass in our model emerges when the frequency ω approaches ω p from above; thus, the creation of material demonstrating the negative density and dielectric constant simultaneously remains challenging. A more comprehensive approach should consider inevitable losses resulting in the decay of plasmons [22], consequently influencing the effect of the negative mass considerably, as discussed in [23].