Plasmons in Finite Spherical Electrolyte Systems: RPA Effective Jellium Model for Ionic Plasma Excitations

Plasmons are fundamental collective excitations in many particle charged systems like in free electron liquid in metals, high energy nuclear plasma in solar core or in fusion devices, in ion gas in ionosphere or in intra- and inter-galactic gas clouds. Plasmons play a central role also in small systems, in particular in metallic nanoparticles and in their arrays allowing for subdiffraction light manipulation. In analogy to metallic nanoparticles, we have developed description of the soft plasmonics in finite electrolyte systems confined in micrometer scale by insulating membranes. Plasmon-type excitations in such finite ionic systems are determined via originally formulated theoretical model allowing to describe surface and volume plasmons in confined geometry of the ion liquid. Size-effect for attenuation of surface plasmons in the finite electrolyte system is described and its various regimes are identified. The cross-over in the plasmon damping system-size-dependence is demonstrated including scattering of ions and their energy losses via irradiation. The plasmon resonances in ion systems replicate the metal cluster plasmon phenomena, though in distinct energy and size scale related to larger ion mass and lower ion concentration (in low energy plasma) in comparison to electrons in metals. The possibility for tuning plasmon resonances in finite ionic systems in a wide range by changing system size, ion, and electrolyte parameters is demonstrated.


Introduction
Recent experimental and theoretical investigations of plasmon oscillations in metallic nanoparticles have focused attention on their fundamental character and numerous applications. In particular, the so-called plasmon effect in solar cells modified on the nanoscale with surface-deposited metallic particles has led to improvements of their efficiency [1][2][3][4][5][6]. The surface plasmon oscillations in these particles play a mediating role in harvesting energy from sunlight because of the particle's radiative properties. The radiated energy from plasmon oscillations is preferable for transport applications. As was observed experimentally and predicted theoretically, irradiation losses of plasmon energy are strongly sensitive to the size of metallic nanoparticles [7,8].
The strong irradiation of plasmon oscillations in metallic nanoparticles also plays a major role in the construction of plasmonic waveguides with high transference efficiency. Several experimental studies [9,10] have indicated that periodic linear structures of metallic nanoparticles serve as efficient plasmon waveguides with low damping [11][12][13]. The wavelengths of plasmon polaritons propagating in such structures are typically shorter than wavelengths of light with the same frequency by one or two orders of magnitude, enabling avoidance of diffraction limits in light circuits [14][15][16]. This avoidance enables construction of plasmonic opto-electronic nanodevices not available when using only light waveguides limited by diffraction constraints. The efficient energy transfer in plasmonic waveguides is also supported by the radiative properties of surface plasmons in metallic nanocomponents.
Accelerating charges irradiate electromagnetic waves, and the related energy loss can be accounted for by an effective electric field that hampers charge movement called as the Lorentz friction [17,18]. In the case of an oscillating dipole such as the dipole-type surface plasmons in a metallic cluster, the Lorentz friction force is proportional to the third-order time-derivative of the dipole [17]. Let us emphasize here that the strong irradiation of surface plasmons in metallic clusters is exclusively present in sufficiently large metallic particles (larger than ca. 15 nm in diameter, for Au or Ag). Ultra-small metallic clusters with diameters of 2-10 nm do not exhibit irradiation efficiency as high as that of nanospheres with radii a > 15 nm, mostly because of the much lower number of electrons in so small clusters in comparison to the number of electrons in larger nanospheres. In particular, extensive attention has been focused on large nanoparticles of noble metals (gold, silver, and copper) because their plasmon resonances are located within the visible region of the electromagnetic spectrum.
Plasmons in metallic nanostructures focused interest on similar local charge density oscillations in amorphous materials or in other unconventional systems [19][20][21] including also ionic systems and interaction of ions with metallic plasmons [22] as well as soft flexible photonic crystals with wide applications [23]. Combining metallic nanoparticles with host materials whose dielectric properties can be tuned by means of an external control is one route to create active plasmonics. By exploiting the outstanding properties of self-organizing materials, so-called smart matter, including plasmons interacting with light, a bridge between two branches of physics: 'hard matter' and 'soft matter' can be built. The soft plasmonics and plasmonic crystals have found already many applications in photonics and in sensing [24]. Some examples of natural and man-made plasma besides of the metal-electron-plasma can be listed as follows: laboratory gas discharge, e.g., in vacuum tubes, spark gaps, welding arcs, and neon or fluorescent lights, controlled thermonuclear fusion experiments, e.g., in tokamak, earth ionosphere that is partially photo-ionized by solar ultraviolet radiation, sun's core, where fusion of hydrogen to form helium generates the sun's heat, solar wind, i.e., the wind of plasma that blows off the sun, interstellar, and intergalactic ionic gas medium. The plasmon oscillationslocal fluctuations of the charge density-have the frequency proportional to the square root of the charge carrier density and this frequency varies between 10 18 1/s (for solar core), across 10 8 1/s (ionosphere) to 10 2−4 1/s (interstellar and intergalactic ionized gas). The plasmon excitations in the ionosphere have, in particular, a great importance in radio-communication and in over-horizon radar techniques. Much attention, also in experiments with ionized fullerene gas, has been recently paid to electrostatic waves in so-called pair-ion-plasma consisting of only positive-and negative-charged equal mass particles with a time-space parity kept because the mobility of the equal mass particles in electromagnetic fields is the same [25,26]. The magnetohydrodynamics of such a system may have importance in understanding of cosmic ion systems in electromagnetic field and in fusion plasma, previously analyzed upon the Vlasov kinetic equation [27]. Despite the ionic plasmons are thus well recognized in various large and open systems [28,29], their counterparts in small confined electrolytes are not penetrated as of yet, but probably offer a rich physics as they might replicate the plasmonics in small metallic particles. Micrometer scale of electrolyte confinement is frequent in bio-cell organization, where a cytoplasm containing various ions is separated from surroundings by lipid cellular membrane. The local charge fluctuations in such ionic finite systems can be important in biophysical phenomena of communication and signaling as well as in local energy transport. An example is the newly developed plasmon-polariton model of the so-called saltatory-conduction in neuron long axon [30].
An interesting question arises as to whether similar to complicated metallic plasmonic effects are possible with ionic carriers instead of electrons. Many finite ionic systems in the form of electrolyte enclosed by membranes are found in biological structures. The question then arises as to what role plasmonic phenomena would play in such structures and whether the radiative properties of plasmon fluctuations would be as significant in ionic systems as in metals. One can reasonably expect that ionic plasmon effects would be located in different energy and wavelength scale regions compared to those of metallic systems due to the larger mass of ions and smaller concentration than for metal electrons. The ionic soft plasmonics could be linked to the functionality of biological systems in which electricity is of an ionic rather than an electronic character, such as cell signaling, membrane transfer, and nerve-cell conductivity.
Ionic systems are much more complicated in comparison to a metal crystal structure with free electrons. Therefore, identification of appropriate model simplifications is of primary significance to properly describe collective charge excitations in electrolytes along with keeping analogy with metal plasmonics.
In the present paper, we will consider a finite spherical ionic system (e.g., liquid electrolyte confined within a spherical membrane) to identify plasmonic excitation. We will determine their energies with respect to various parameters of the ionic system, with special attention paid to the irradiation properties of ionic plasmons.
The paper is organized as follows. In the first paragraph, the effective model for plasmon fluctuations in finite spherical binary electrolyte system is formulated via introducing of specially defined two ion-jellium components (for electrolytes, the jellium is fictitious and an auxiliary model construction unlike to metals where the jellium is the real rigid crystal positive core). The model is utilized in the next paragraph to identify surface and volume ionic plasmon excitations in finite spherical electrolyte system ranged by neutral dielectric membrane. Both self-frequencies and the attenuation rates are next estimated for ionic surface and volume plasmons in the considered system, for various ion and electrolyte concentration parameters and system sizes. The radiative properties of the ion surface dipole plasmons are next examined in details which are of importance for the energy and information transfer in ionic systems in possible application to electrolyte finite components in biological cellular structures and their electrical functioning.

Fluctuations of Charge Density in a Finite Spherical Ionic System
For a simple two-component ionic system, we address the water solution of ions of both signs, creating an electrolyte with balanced total charge, enclosed in finite size spherical volume ranged by electrically neutral insulating membrane. At equilibrium, also local charge cancellation holds. Both types of ionic carriers can form, however, density fluctuations, resulting in disruption of the local electric equilibrium. The total charge conservation and neutrality condition require that any density fluctuation of negative charges must be accompanied by equivalent fluctuation in the positive ions, possibly in another even distant place of the system, and vice versa. Therefore, we effectively deal with density fluctuations of ions, positive and negative (but always mutually compensated), with respect to a uniform fictitious background charge distributions of the opposite sign in some analogy to the jellium model in metals. In the case of two component electrolyte, these auxiliary uniform background charge distributions cancel mutually themselves and do not modify the system and its energy. Each of these auxiliary ionic jelliums has the total charge equaled to the sum of charges of all ions of the opposite sign. Thus, the opposite oscillations of both types of ions are represented here as the equivalent sum of two ion-jellium systems also simultaneously oscillating.

Definition of the Model
To develop a model according to these guidelines, let us consider a spherical system with a radius a and a balanced total charge of both sign ions with uniform equilibrium density distributions n +(−) (r) = n (a − r) (where (r) is the Heaviside step function). The equilibrium density of the charged liquid, denoted by n, will be treated as a parameter and n = ηN 0 will be assumed, where η is the molarity of the electrolyte within the sphere and N 0 is the one-molar electrolyte concentration of ions.
The Hamiltonian for the simplest two-component ion system has the form, where q −(+) , m −(+) , and N −(+) are the charge, mass, and total number of the −(+) ions, respectively. Indices i and j are introduced to distinguish two sorts of ions. To analyze this complicated system, we propose the following approximation. Assume, for simplicity, q − = −q + = q, N − = N + = N, and m − = m + = m (the generalization to distinct charges and masses of both sign ions is straightforward), and let us add and subtract the same terms (the last four terms in the following expression for the Hamiltonian), In this way, we have formally introduced a jellium of spherical shape for both types of ions, with the density n ideally compensating opposite charges of uniformly distributed ions, n(r) = n (a − r), a is the sphere radius. Assuming now, upon a rough approximation that, one can can separate the Hamiltonian (1) into the sum H ions =Ĥ − +Ĥ + , wherê The latter term in the right-hand side of Eq. 4 corresponds to the interaction between ions of the same sign, whereas the second term in the first sum describes the interaction of these ions with the jelliums of opposite sign (ε is the dielectric constant of the electrolyte medium). Because of the separation of the Hamiltonian (1), one can consider a single Hamiltonian (4).
The question is the applicability of the condition (3). This condition means that the mutual interaction between two sign ion fluctuations is equal to the sum of interactions of these fluctuations with the fictitious jelliums of opposite signs ideally neutralized themselves. For small fluctuations, one can argue that the total energy of the ion interaction is not changed by the approximation (3) and it may be used to the assessment of energy scale of ion fluctuations. The advantage of such an approach is the close analogy to the description of plasmons in metals, including the direct definition of the shape of the system with the explicit rigid jellium form. Let us emphasize that the apparent decoupling of both sign ion fluctuations via Eq. 3 is not complete in fact. Each fluctuation of ion density beyond the uniform equilibrium distribution, let say of negative ions, produces noncompensated positive charging of the fictitious jellium being in fact the fluctuation of positive ions in the real system. The decoupling corresponds thus to duplication of the ion-jellium dynamics description for both ion types, without, however, the change of interaction energy, for small fluctuations at least. Both the ion-jellium dynamics describe the same fluctuations of ion densities in terms of opposite sign ions, which are actually coupled in the binary electrolyte. Such a picture is of particular usefulness for the case of the same charges and masses of ions, whereas it worsens with rising differences between ion parameters, when asymmetry between separate ion-jellium oscillations grows. The alternative way to introduce the jellium model for ions in the electrolyte might be the definition of an effective ion (with effective charge and mass in a far analogy to the two body problem) comprising all the two component dynamics described as fluctuations with respect to the single opposite sign jellium. Such a model would diminish, however, the level of degrees of freedom in comparison to the real system and in opposition to the two component ion-jellium model. Therefore, we will develop and test the two component ion-jellium model resulted from the approximate Eq. 3 aiming on even rough assessment of the energy scale of ion plasmon fluctuations. Upon the proposed model, the equilibrium ion density determines the bulk plasmon frequency for the ion system (for each type of ions) according to the formula analogous to bulk metal [31], ω 2 p = 4πnq 2 m , where n and m are the equilibrium uniform concentration and the mass of ions with charge q, respectively. Because m is larger than the electron mass, m e , and the ion concentration is usually smaller than that of electrons in metals, ω p can be considerably reduced, even by several orders of magnitude. Noticeably, for electrons in metals, ω p 10 eV and typically falls in the ultraviolet region. In an ionic system, the plasmon frequency ω p can be much lower: in the infrared or even lower-energy regions.
The form of the Hamiltonian (4) allows for its utilization in the scheme applied to electrons in metals [31][32][33]. A local density of ions can be written, analogous to the semiclassical Pines-Bhom random-phase approximation (RPA) of electrons in metals [31,32], in the following form: where r j denotes the coordinate of the j − th ion and the Dirac delta semiclassically fixes the j − th ion position; ion (t) denotes the ion wave-function corresponding to the Hamiltonian (4). The Fourier picture of the local density of ions has the form: where the operatorρ(k) = j e −ik·r j .
Using the aforementioned notation, one can rewriteĤ ion in the following form, analogous to the case for metallic plasmons [31][32][33]: whereñ(k) = d 3 rn(r)e −ik·r is the Fourier picture of the jellium distribution (in the derivation of Eq. 7 we have taken into account that 4π Utilizing this form of the effective ion Hamiltonian, one can write out the dynamic equation in Heisenberg representation for the ion density fluctuations ([..] denotes the commutator), which attains the following form: where δρ(k) =ρ(k) −ñ(k) describes the operator of local ion density fluctuations with respect to the equilibrium uniform density. Averaging over the quantum states | ion >, we obtain the following equation for the ion density fluctuations: For small k, as with the semiclassical approximation for electrons [31,33], the contributions of the second and third components of the first term on the right-hand side of Eq. 10 can be neglected as small in comparison to the first component (with the lowest power of k). The third term in the right-hand side of Eq. 10 is also small (and thus negligible), involving a product of two δρ (which we assumed to be small, δρ/n << 1). This approach corresponds to the RPA formulated for bulk metal [31,32]. Within the RPA, Eq. 10 takes the following shape: and, because of spherical symmetry, Equation 11 can be rewritten in the position representation: t) . (12) In the case of metals, the Thomas-Fermi formula is used to assess the averaged kinetic energy [32]: This formula, however, refers to fermionic and degenerate quantum systems, such as electrons in metals. For ionic systems, such an estimation of kinetic energy is inappropriate because the ion concentration is usually much lower than the concentration of electrons in metals and because the system is not degenerate even if the ions are fermions. The Maxwell-Boltzmann distribution should be applied instead of the Fermi-Dirac or Bose-Einstein distribution. Independent of fermionic or bosonic ion statistics, the Maxwell-Boltzmann distribution allows for an estimation of the averaged kinetic energy of ions located inside a sphere of radius a in the following form: where k is the Boltzmann constant and T is the temperature. For ionic molecules with 3D or linear shapes, the inclusion of rotational degrees of freedom results in the factor 6kT 2 or 5kT 2 , respectively, rather than 3kT 2 for the point-like ion model.
Using the formula (14) and taking into account that ∇ (a − r) = − r r δ(a − r), Eq. 12 can be rewritten in the following manner: In this formula, ω p is the bulk ion-plasmon frequency, The solution of Eq. 15 can be decomposed into two components related to the distinct domains inside the sphere and on the sphere surface: The domains correspond to volume and surface excitations, respectively. These two parts of local ion density fluctuations satisfy the equations (according to Eq. 15), and (here = 0+), The Dirac delta in Eq. 18 results from the derivative of the Heaviside step function, the ideal jellium charge distribution. In Eq. 18, an infinitesimal shift, = 0+, is introduced to fulfill the requirements of the Dirac delta definition (its singular point must be an inner point of an open subset of the domain). This shift is only of a formal character and does not reflect any asymmetry.
The electric field due to surface charges is zero inside the sphere and therefore cannot influence the volume excitations. Conversely, the volume charge fluctuation-induced electric field can excite surface fluctuations. Therefore, the equation for volume plasmons is independent of surface plasmons, whereas the volume plasmons contribute to the equation for the surface plasmons.
The problem of separation between surface and volume plasmons has been thoroughly analyzed for metal clusters and has been identified as particularly significant for very small clusters. In the size scale of 1-3 nm for metallic clusters, the effect of so-called electron spillout beyond the jellium edge is important and causes a fuzzy surface resulting in the coupling of volume and surface plasmon oscillations. Direct numerical simulations of time-dependent local density approximation (TDLDA) [34,35] have verified that the volume-surface excitation mishmash gradually disappears in larger metallic clusters [34,35], which supports the accuracy of the semiclassical RPA description, within which volume plasmons are separated from surface ones. The role of spill-out effect diminishes gradually with growing sphere size as the ratio of surface to volume falls down and spill-out becomes negligible in the range of several nanometers for metals and similarly for large ionic spheres. Moreover, for the electrolyte system confined by the insulating membrane, the spill-out of ions is irrelevant. Therefore, the RPA description of ionic density fluctuations is a proper model. Moreover, the analytical RPA semiclassical picture in the form of an oscillator equation allows for convenient inclusion of damping effects, which is especially important since the plasmon damping caused by irradiation losses turns out to be an overwhelming physical property of plasmons in the case of large metallic nanospheres [7,33] (with a > 15 nm for Au or Ag) as well as of large ionic systems, as demonstrated in the paragraph "Damping of Plasmon Oscillations in Ionic Systems."

Solution of RPA Equations: Volume and Surface Ionic Plasmon Frequencies
Equations 17 and 18 are solved for metallic nanospheres [33], and these solutions can be directly applied to ionic systems. To briefly summarize this analysis, we represent both parts of the plasma fluctuation as follows: with initial conditions F (r, t)| t=0 = 0, σ( , t)| t=0 = 0, ( is the spherical angle), F (r, t)| r=a = 0, ρ(r, t)d 3 r = N (neutrality condition). With the above initial and boundary conditions and taking advantage of the spherical symmetry, we write the time-dependent parts of the ion concentration fluctuations in the form [33] (cf. Appendix): and where j l (ξ ) = π 2ξ I l+1/2 (ξ ) is the spherical Bessel function, Y lm ( ) is the spherical function (some examples are presented in Fig. 1), ω li = ω p 1 +  The function F (r, t) describes volume plasmon oscillations, whereas σ ( , t) describes surface plasmon oscillations. We emphasize that the first term in Eq. 21 corresponds to the surface self-oscillations, whereas the second describes the surface oscillations induced by the volume plasmons. The frequencies of the surface self-oscillations are equal to which, for l = 1, is a dipole-type surface oscillation frequency, described for metallic nanospheres by Mie [36],

Ionic Surface Plasmon Frequencies for Nanospheres Embedded in a Dielectric Medium
The influence of dielectric surroundings (generally distinct from the inner dielectric of the ionic system) on plasmons in this system can now be included. Let us assume that ions on the surface (r = a+, i.e., r ≥ a, r → a) interact with Coulomb forces renormalized by the relative dielectric constant ε 1 > 1 (distinct from the ε of the inner medium). Thus, a small modification of Eq. 18 is in order: (note that Eq. 17 is not affected by the outer medium). The solution of Eq. 23 is of the same form as that of Eq. 18, but with renormalized surface plasmon frequencies:

Damping of Plasmon Oscillations in Ionic Systems
The presented above semiclassical RPA treatment of plasmon excitations in finite ion systems does not account for plasmon damping. The damping of plasmon oscillations can be included in a phenomenological manner by the addition of an attenuation term to the plasmon dynamic equations,  Fig. 2 The spherical Bessel functions J l (r) for l = 0, 1, 2 displaying possible charge density fluctuations in the sphere along the radius r for volume plasmon modes. The angular distribution of these modes is governed by the real spherical functions Y lm ( ), similar to the surface plasmon modes (cf. Fig. 1 right). Right: the exemplary temperature dependence of the self-frequencies of the volume plasmon modes ω li , li = 11, 12, 21 for a dilute electrolyte n 10 14 1/m 3 and ion mass ∼ 10 4 m e , a ∼ 50μm i.e., the term − 2 τ 0 ∂δρ(r,t) ∂t added to the right-hand sides of both Eqs. 17 and 18, taking into account their oscillatory form. The introduced damping ratio 1 τ 0 accounts for ion scattering losses and can be approximated, in analogy to metallic systems, by the inclusion of energy dissipation caused by its irreversible transformation into heat via various microscopic channels, similar to Ohmic resistivity [37]: where a is the sphere radius, v = 3kT m is the mean velocity of ions, and λ b is the ion mean free path in the bulk electrolyte material (comprising effects of scattering of ions on other ions, on solvent particles and admixtures). The second term in Eq. 25 accounts for the scattering of ions on the boundary of the finite ionic sphere of radius a, where the constant C is of the order of unity and reflects the type of scattering of ions by the boundary of the sphere dependent of microscopic particularities of the membrane [37].
To explicitly express a forcing field that moves ions in the system, the inhomogeneous time-dependent term should be added to the homogeneous equations (17) and (18). The forcing field may be a time-dependent electric field (e.g., the electric component of the incident e-m wave which may excite plasmons). Similarly as for metallic nanospheres, the surface plasmon resonant wavelength highly exceeds the system dimension in the case of finite ionic systems and the e-m forcing field is essentially space-homogeneous along the whole sphere. Such a perturbation fulfills the so-called dipole approximation requirements and excites only surface dipole plasmons, i.e., the mode with l = 1, which can be described by the function Q 1m (t) (l = 1 and m are angular momentum numbers related to the assumed spherical symmetry). The corresponding dynamical equation for the surface plasmons reduced to only mode Q 1m (t) has the following form: where ω 1 = ω p √ 3ε 1 (a dipole surface plasmon frequency, ε 1 is the dielectric susceptibility of the system surroundings). Because only Q 1m modes contribute to the plasmon response to the homogeneous electric field, the effective ion density fluctuation has the form [33] δρ(r, t) = ⎧ ⎨ ⎩ 0, r < a, where Y lm ( ) is the spherical function with l = 1. One can also explicitly calculate the dipole D(t) corresponding to surface plasmon oscillations given by Eq. 27: The dipole D(t) satisfies the equation (from Eq. 26) Noticeably, the dipole (28) scales as the system volume, ∼ a 3 , indicating that all ions actually contribute to the surface plasmon oscillations. This observation is connected with the fact that the surface modes correspond to uniform translation-type oscillations of ions in the system when the charge of ions inside the sphere is exactly compensated by oppositely signed ions, whereas unbalanced charge density occurs only on the surface, despite all the ions oscillating. For volume plasmons, non-compensated charge density fluctuations are present inside the sphere because volume plasmon modes have compressional character with unbalanced charge fluctuations along the system radius.
The scattering effects accounted for by the approximate formula (25) cause damping of plasmons and are especially strong for small systems because of the nanosphere-edge scattering contribution, which is proportional to 1 a . The significance of this term, however, decreases with increasing radius. We will show that radiation losses (due to Lorentz friction) scales initially as a 3 and that, for increasing a, these irradiative energy losses quickly dominate plasmon attenuation. Because of the opposite size dependencies of the scattering and irradiation contributions to plasmon damping, we observe a cross-over of damping with respect to the size, as depicted in Fig. 3. In addition, the radius a * for which the total attenuation rate of surface plasmons is mini- , can also be determined. The system sizes a * for two distinct ionic systems are listed in Table 1. The radiative energy loss of the oscillating dipole is expressed by the Lorentz friction [17], i.e., the effective electric field slowing the motion of charges: Hence, we can rewrite Eq. 29 to include the Lorentz friction term:  Fig. 3 The cross-over in the system size-dependence of the surface plasmon damping rate for T = 300 K, m = 10 4 m e , q = 3e, n = 10 −2 N 0 (where N 0 is the concentration of one molar electrolyte) (left) and for n = 10 −3 N 0 (right); in the region close to the cross-over, the perturbative treatment for Lorentz friction well coincides with the exact approach or for E = 0, The perturbation method can be applied for a solution of Eq. 32 when the right-hand side of this equation is treated as a small perturbation. In the zeroth step of the perturbation, we have ∂ 2 Hence, for the first step of the perturbation, we substitute the latter formula into the right-hand side of Eq. 32, i.e., Within the first step of perturbation, the Lorentz friction can be included in the total attenuation rate 1 τ . Nevertheless, this approximation is justified only for sufficiently small perturbations, i.e., when the second term in Eq. 34, which is proportional to a 3 , is sufficiently small to fulfill the perturbation restrictions. The related limiting value,ã, of the ionic system size depends on the ion concentration, charge, mass, and dielectric susceptibility, as is exemplified in the following paragraph.
The solution of Eq. 33 is of the form D(t) = Ae −t/τ cos(ω 1 t + φ), where ω 1 = ω 1 1 − 1 (ω 1 τ ) 2  Radius for maximal damping a * * from maximum of I m 2 given by Eq. 36 8 × 10 −6 m 2 5 × 10 −6 m a red shift to the plasmon resonance because of a strong (∼ a 3 ) increase of attenuation caused by the irradiation. The Lorentz friction term in Eq. 34 dominates the plasmon damping for a * < a <ã because of this a 3 dependence (cf. Fig. 3). Plasmon damping grows rapidly with a, which results in the pronounced red shift of the resonance frequency.

Exact Inclusion of Lorentz Damping to the Attenuation of Ionic Dipole Surface Plasmons
Let us now consider the dynamic equation for surface plasmons in an ionic spherical system, Eq. 32, with the Lorentz friction term, but without applying the perturbation method for a solution. To compare various contributions to Eq. 32, we change to a dimensionless variable t → t = ω 1 t. Equation 32 becomes the form, When solving Eq. 35 by perturbation, we obtain a renormalized attenuation rate for an effective damping term, . This term quickly reaches unity, for which the oscillator falls into the over-damped regime. For the system parameters assumed for Fig. 3, the attenuation rate reaches unity at 25.5 and 8 μm for n = 10 −3 N 0 and n = 10 −2 N 0 , respectively. At these values of a, the frequency ω 1 = ω 1 1 − 1 (ω 1 τ ) 2 goes to zero, which indicates an apparent artifact of the perturbation method. To verify the exact damped frequency behavior in this system, one must solve the dynamical equation (35) without any approximations. As this equation is a third-order linear differential equation, its solution takes the form ∼ e i t , with analytical expressions for three possible values of the exponent: In Fig. 4, we have plotted the damping rate (I m 2 ) and the self-frequency (Re 2 , in right panel the corresponding wavelength is visualized) with respect to the system radius a. For comparison, the approximate perturbative solutions are also plotted (in blue line whereas the exact solution of Eq. 35 in red line). The blue line ends at a limit when the attenuation rate within the perturbation approach reaches the critical value 1 (then λ → ∞). For an accurate solution of Eq. 35, this singular behavior disappears and the oscillating solution exists for larger a as well.
We note that the red shift of the plasmon resonance is strongly overestimated in the framework of the perturbative approach to Lorentz friction unless a <ã, whereã is sensitive to ionic system parameters and especially to the ion concentration (demonstrated in Fig. 4).
We emphasize that the Eq. 35 has, in general, two types of particular solutions: e i t , with complex self-frequencies . The solutions given by 2 and 3 are of the damped oscillator type (i 2 and i 3 are mutually conjugated, thus ble and exponentially rising solution (negative imaginary solution). This unstable solution is a well-known artifact in Maxwell electrodynamics (cf. e.g., $ 75 in [17]) and corresponds to infinite self-acceleration of the free charge due to the Lorentz friction force (i.e., the singular solution of the equation mv = const.×v, which is associated with a formal renormalization of the field-mass of the charge: infinite for a point-like charge and canceled in an artificial manner by an arbitrarily assumed negative infinite non-field mass, resulting in the ordinary mass of, e.g., an electron, although not defined in a properly mathematical way). This unphysical singular particular solution ( 1 ) should thus be discarded. The other oscillatory solution resembles that of the ordinary damped harmonic oscillator, though with a distinct attenuation rate and frequency. This attenuation rate and frequency are expressed by analytical formulae for 2 (or equivalently, 3 ) by Eq. 36 and are calculated for various values of a and compared with the corresponding quantities found within the perturbation approach. This comparison is presented in Fig. 4. This comparison reveals that the application of the perturbation approach leads to a high overestimation of the damping rate for a >ã. Therefore, we conclude that the use of the approximate formula for Lorentz friction damping in Eq. 34 is justified up to a ã, whereas for a >ã, this approximate expression causes a large error, rapidly growing with a, in comparison to the exact solution. The valuẽ a < a limit sharply depends on ionic system parameters and can be approximated byã a limit 2 , where a limit denotes here the sphere radius at which the oscillatory perturbation solution terminates (cf. Fig. 4).
It must be emphasized, however, that in the electrolyte systems there occur additional excitations of the system components which may mix with plasmon oscillations and even could wash out plasmonic effects in similar time scale of oscillations. The most pronounced such an effect is the strong absorption of microwave e-m radiation by rotational degrees of freedom of water dipole molecules. In wide region of GHz, water molecules in liquids and in other materials heavily absorb energy, rotate, and dissipate it next as the heat to the surroundings (for utilization of this effect in microwave ovens typically is used 2.5 GHz e-m wave). Besides of rotational excitations, a separate water molecule exhibits the vibration resonant absorption at ca. 100 THz in vapor phase, whereas in a liquid phase at ca. 50 THz. Therefore, the plasmon frequencies of GHz order and of several THz order frequencies inconveniently are placed in the region of water absorption: the micro-wave rotational one at wide range between 1 and 300 GHz and the vibrational one at ca. 50 THz. To avoid the overlap between plasmon frequencies and absorption frequencies of water, one can address the soft plasmonics of water electrolytes toward lower (of MHz range frequencies) via reduction of the ion concentration and taking advantage of significant increase of water dielectric permittivity in MHz range, lowering plasmon frequency (the relative permittivity of water is ε 80 for frequencies in the MHz range [38], although for higher frequencies, beginning at approximately 10 GHz, this value decreases to approximately 1.7, corresponding to the optical refractive index of water, η √ ε 1 = 1.33).
Moreover, in non-spherical geometry of an elongate cylinder or strongly prolate ellipsoid, the lowering of longitudinal surface plasmon frequency might be very effective [39,40], and this frequency can be pushed below the inconvenient GHz window for even larger ion concentrations. Note that such a situation can occur in biological electrolyte systems, e.g., in thin and long electrolyte neuron cords [30]. The developed above model of plasmons allowed for description of plasmon-polaritons in ionic micro-chains of periodically confined electrolytes which has been next utilized to elucidate of previously not known and longsearched mechanism of so-called saltatory conduction in myelinated axons in peripheral nervous system and in the white substance of central nervous system. The transfer of action potential in long axons covered with white lipid myelin (and therefore the name of white substance in brain and in spinal cord) thick sheath exhibits much larger velocity (100 times larger) than the diffusive velocity of ions in neuron cytoplasm. The mechanism of this acceleration is not known and according to observed jumps of the signal between neighboring small intervals in myelin layer called as Ranvier nodes, the related behavior has been called as the saltatory conduction. This jumping of the signal across the myelinated sectors of length typically of 100 μm between consecutive Ranvier nodes accelerates the transfer of action potential and causes that so-called firing of the axon is quick, which in turn is essential for communication and functioning of the body. Upon the assumption of the plasmon-polariton mechanism of neuron firing in the case of periodically myelinated axons, a very good coincidence of signal velocity for realistic parameters of neuron cord size, internal cytoplasm electrolyte molarity, and thickness of the myelin sheath can be achieved [30].

Conclusions
In summary, we can state that in finite ionic systems, one can observe plasmons similar to those in metallic nanoparticles. The structure of ionic surface and volume plasmons is analogous to similar properties of electronic plasmons in metallic spherical systems, albeit with a significant red shift of resonance energy corresponding to the far larger mass of ions compared to that of the electrons and the lower concentration of ions in electrolytes compared to the concentration of electrons in metals. Thus, corresponding to the resonance energy, wavelength is shifted to deep-infrared or even longer wavelengths depending on ion concentration. The typical cross-over in the plasmon damping size dependence for metal clusters between the scattering, Ohmic-type energy dissipation, and the radiative losses is observable in spherical ionic systems. This dependence is similar to the size-dependence in metals, though shifted toward the micrometer scale for ions instead of toward the nanometer scale for metals. Of particular interest is the high irradiation regime for dipole plasmons in ionic systems, with potential applications for signaling and energy transfer. The initial strong enhancement of the efficiency of the Lorentz friction with increasing radius of the electrolyte sphere is observed on the micrometer scale with typical a 3 radius dependence above some threshold that depends on the electrolyte parameters. At a certain value of the radius (which also strongly depends on the ion system parameters), this enhancement saturates and the radiative losses slowly diminish, which allows for the definition of the most convenient size of the finite electrolyte system for optimizing radiation-mediated energy transport efficiency, preferring the highest radiation losses. The last two terms of the right-hand side of Eq. 44 can be transformed using the generating function for Legendre polynomials [41], where a = a r r and cosγ = a·r 1 ar 1 . Employing Eq. 48, the last two terms in Eq. 44 can be transformed as follows: