Next Article in Journal
Casimir Puzzle and Casimir Conundrum: Discovery and Search for Resolution
Next Article in Special Issue
Jet-Induced Feedback in the [O III] Lines of Early Evolution Stage Active Galactic Nuclei
Previous Article in Journal
Blocks Size Frequency Distribution in the Enceladus Tiger Stripes Area: Implications on Their Formative Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

On the Heating of AGN Magnetospheres

1
School of Physics, Free University of Tbilisi, 0183 Tbilisi, Georgia
2
E. Kharadze Georgian National Astrophysical Observatory, 0301 Abastumani, Georgia
3
Institute for Fusion Studies, The University of Texas at Austin, Austin, TX 78712, USA
*
Author to whom correspondence should be addressed.
Universe 2021, 7(4), 83; https://doi.org/10.3390/universe7040083
Submission received: 26 February 2021 / Revised: 23 March 2021 / Accepted: 29 March 2021 / Published: 31 March 2021
(This article belongs to the Special Issue Panchromatic View of the Life-Cycle of AGN)

Abstract

:
The Langmuir–Landau-Centrifugal Drive (LLCD), which can effectively “convert” gravitational energy into particles, is explored as a driving mechanism responsible for the extreme thermal luminosity acquired by some active galactic nuclei (AGN). For this purpose, we consider equations governing the process of heating of AGN magnetospheres. In particular, we examine the Fourier components of the momentum equation, the continuity equation and the Poisson equation in the linear approximation and estimate the growth rate of the centrifugally excited electrostatic waves and the increment of the Langmuir collapse. It is shown that the process of energy pumping is composed of three stages: in the first stage the energy is efficiently transferred from rotation to the electrostatic modes. In due course of time, the second regime-the Langmuir collapse-occurs, when energy pumping is even more efficient. This process is terminated by the Landau damping, when enormous energy is released in the form of heat. We show that the magnetospheres of the supermassive black holes with luminosities of the order of 10 45 46 erg/s can be heated up to 10 6 10 K.

1. Introduction

Rotating magnetospheres, with relatively strong magnetic fields, surrounding active galactic nuclei (AGN), are believed to be the regions where cosmic rays might originate [1,2]. These particles are highly relativistic and one of the open problems in modern astrophysics is how do they attain such high energies.
Several mechanisms have been suggested to explain acceleration of leptons or hadrons to ultra-high energies. In particular, in the so-called Fermi class (and its modifications) of processes [3,4,5], the acceleration occurs in relatively strong magnetic fields in the magnetospheres of astrophysical objects. The Fermi processes, it turns out, are efficient only if the particles are already relativistic. For overall efficiency, therefore, some kind of a pre-acceleration is desirable [1].
In this paper, however, we explore the Langmuir-Landau-Centrifugal Drive (LLCD) mechanism that makes use of the plasma collective phenomena manifested though the Langmuir wave that will serve as a conduit for transferring gravitational energy in the rotating magnetosphere of an AGN to particle kinetic /thermal energy. Although LLCD has been discussed and developed in several papers [2,6,7,8], we will summarize its salient qualitative features in the introduction.
Since in the AGN magnetospheres magnetic fields are of the order of 10 1 4 G, the charged particle flows will follow co-rotating magnetic field lines (frozen-in condition), experiencing very strong relativistic centrifugal force in the light cylinder (LC) zone (area where the linear velocity of rotation exactly equals the speed of light) [9]. If the plasma number density is high enough to provide the screening of longitudinal (along the magnetic field lines) electric field up to the LC area, the centrifugal mechanism becomes very efficient [10]. The strong relativistic differential centrifugal force leads to several consequences, for instance, the direct acceleration of electrons to very high energies (Lorentz factors of the order of 10 7 ) [1,11]. Although the LLCD mechanism has been shown to accelerate particles to even higher energies, we in this paper emphasize its possible role in heating of the magnetosphere.
The first step in the LLCD process is the parametric excitation of Langmuir waves through charge separation caused by differential rotation of different charge species (electrons and protons ) endowed with different Lorentz factors. This process, examined for millisecond pulsars [12] and AGN [13], proves to be very effective in pumping enormous energy (from rotation, which in turn is gravitationally driven) into what could be called centrifugally driven electrostatic waves.
In the second stage of LLCD, the electrostatic waves Landau damp and transfer energy from centrifugally amplified modes to particle kinetic energy. The efficacy of this wave-particle energy transfer has been demonstrated for millisecond and newly born pulsars [8,14]; it was shown that LLCD could energize particles to the tune of 10 18 eV.
In the study of the role of LLCD in the AGN magnetospheres [2], and SgrA (located in the center of the Milky Way) [7], we found that the parametrically amplified Langmuir waves, before undergoing standard Landau damp, go through another intermediate boosting via Langmuir collapse. One learns from the classical work by [15] that the Langmuir turbulence in the non-linear regime creates unstably excited caverns or the low-density regions. As a result, the high frequency pressure pulls the particles from the cavern, which in turn, provokes highly unstable collapse leading to fast amplification of the electrostatic field. The corresponding problem has been studied for three-dimensional geometry by [16] considering spectra of Langmuir turbulence. The termination by means of the Landau damping of the collapse process has been numerically studied by [17]. This last stage is the most important one in energy pumping process, because the whole energy of Langmuir waves transfers to particle kinetic energies. Applying this process to AGN [2] and the central black hole of our galaxy [7], it has been found that protons might achieve extremely high energies of the order of Z e V (in AGNs) and P e V (in SgrA ).
A principal consequence of Langmuir collapse might be not only the particle acceleration in a particular direction but also the efficient heating of the sustaining plasma. It is this heating that is posited here to be the cause of the observed intrinsic brightness (high temperatures 10 6 9 K) of some AGNs [18].
The paper is organized as follows: in Section 2, the essential formulation of LLCD is recalled, in Section 3, the theoretical model is applied to AGN heating problem, and in Section 4, our principal results are summarized.

2. The Theoretical Model of LLCD

The plasma under consideration lies in the magnetosphere of a supermassive black hole of typical mass M = 10 8 9 × M ( M 2 × 10 33 g is the solar mass), and corresponding angular velocity of rotation
Ω a c 3 G M 10 3 a M 8 rad / s 2 ,
where c is the speed of light, G 6.67 × 10 8 dyne-cm 2 /g 2 is the gravitational constant and M 8 M / ( 10 8 M ) and 0 < a 1 are respectively dimensionless mass and a dimensionless parameter measuring the rate of rotation.
If the plasma in the near vicinity of a supermassive black hole is anchored by a magnetic field strong enough to maintain frozen-in condition, i.e., the charged particles will follow the field lines. In our model, we assume that the AGN radiation energy density and magnetic field energy density are of the same orders of magnitude (approximate equipartition). Then, on the light cylinder surface, where the relativistic effects of rotation are the most efficient, the strength of magnetic field is given by [11],
B 2 L R l c 2 c 870 × a × L 45 1 / 2 × M 8 G ,
where R l c = c / Ω is the light cylinder radius and L 45 = L / 10 45 erg / s is the dimensionless luminosity. As has been explained in detail by [1], the charged particles co-rotate with the field lines until the Lorentz factor reaches its maximum value, γ m a x e B / ( 2 m p Ω c ) 2 / 3 1.2 × 10 7 L 45 1 / 3 (here m p is the proton’s mass). Due accretion, the nearby area of AGN contains soft thermal photons which, by means of the IC scattering, might affect the process of acceleration. IC operates primarily on electrons while it is strongly suppressed for protons (by a factor 10 13 ) and does not limit the maximum energy protons may acquire. In the framework of our approach the magnetic field lines almost up to the LC are supposed to be approximately straight. In the rotating magnetosphere this means that the particles in the rotating frame of reference move along straight trajectories.
However, it has been shown that on the LC itself the field lines are swept back, lagging behind the rotation [19] but up to the mentioned zone the field lines are quasi straight.
The original content of the LLCD mechanism, the parametric generation of centrifugally driven Langmuir waves, is contained in the set of linearized fluid equations (in Fourier space) coupled to the Poisson equation [2,8,14]:
p β t + i k υ β 0 p β = υ β 0 Ω 2 r β p β + e β m β E ,
n β t + i k υ β 0 n β , + i k n β 0 υ β = 0
i k E = 4 π β n β 0 e β ,
where Equation (3) is the equation of motion in the presence of a centrifugal force, Equation (4) represents the continuity equation, and Equation (5) is the Poisson equation for the electric field E. In the preceding equations, β is the species index (either electrons or protons), p β is the first order dimensionless momentum ( p β p β / m β ), k represents the wave number of the excited mode, υ β 0 ( t ) c cos Ω t + ϕ β is the zeroth order velocity and r β ( t ) c Ω sin Ω t + ϕ β is the radial coordinate [2], e β is the particle’s charge and n β and n β 0 are, respectively the perturbed and unperturbed Fourier components of the number density.
The first term on the righthand side of Equation (3) represents the relativistic analogue of the centrifugal force. Since the centrifugal force differentiates between different species (electrons and protons), the resulting charge separation excites the Langmuir instability.
Notice that due to the time dependence of the centrifugal force, the preceding system cannot be “mode analyzed” by the standard techniques. However, we can extract almost as much information from the system as in a typical system, by the method developed in [14]. The ansatz
n β = N β e i V β k Ω sin Ω t + ϕ β ,
converts the aforementioned set of governing equations to a pair of coupled equations (time dependence is left only in χ )
d 2 N p d t 2 + ω p 2 N p = ω p 2 N e e i χ ,
d 2 N e d t 2 + ω e 2 N e = ω e 2 N p e i χ ,
where ω e , p 4 π e 2 n e , p / m e , p γ e , p 3 is the relativistic plasma frequency of the corresponding specie, γ e , p is the Lorentz factor, χ = b cos Ω t + ϕ + , b = 2 c k Ω sin ϕ and 2 ϕ ± = ϕ p ± ϕ e .
A quasi “dispersion relation” for the electrostatic wave is derived after a formal Fourier transform (For the detailed derivation please see the paper by [14] of Equations (7) and (8)
ω 2 ω e 2 ω p 2 J 0 2 ( b ) = ω p 2 μ J μ 2 ( b ) ω 2 ( ω μ Ω ) 2 ,
where J μ ( x ) is the Bessel function. For resonant modes, ω = μ Ω + Δ = ω r + Δ , Δ ω r , we find that the effective dispersion relation is contained in
Δ 3 = ω r ω p 2 J μ r ( b ) 2 2 ,
that is readily solved to yield the instability growth rate
Γ = 3 2 ω e ω p 2 2 1 3 J μ ( b ) 2 3 .
where μ = ω e / Ω .

3. Discussion

Let us now examine the relevance of this instability as an agent for energy transfer in the AGN context. Since this instability is driven by differential rotation between the electron and ion fluids, it falls in the general class of two stream instabilities. We now show that for the magnetospheric plasma, consisting of relativistic electrons and protons, the instability growth rates are “large” implying that Langmuir wave generation is very efficient.
From Equation (10), one can straightforwardly show that for electrons with Lorentz factors γ 1 = 7 × 10 3 , and protons with the same value γ 2 = 10 4 , the timescale of energy pumping from rotation to electrostatic waves, τ 1 / Γ , varies in the range 600–2000 s. On the other hand, the kinematic timescale (escape timescale) of co-rotating particles (electrons and protons) equals P / 4 = π / ( 2 Ω ) 1.5 × 10 4 s. As is evident, the instability timescale is much less than the kinematic timescale, indicating extremely high efficiency of centrifugally driven electrostatic waves.
These relatively high amplitude Langmuir waves will, now, induce high frequency pressure pushing the particles out from the perturbed zone [15]. In the resulting flow density areas (caverns), the penetrating waves will amplify pressure, which in turn, augments the process of pulling out of the particles resulting in the Langmuir collapse.
We assume that the kinetic and potential energies of plasmons inside the caverns are of the same orders of magnitude [20]
k 2 λ D 2 δ n n 0 ,
where λ D k B T 0 / ( 4 π n 0 e 2 ) is the Debye length scale, k B 1.38 × 10 16 erg K 1 is the Boltzmann constant and T 0 is the temperature provided by the accretion process. Since the perturbation density is much less than the unperturbed density, δ n n 0 , energy of plasmons is almost constant
d r E 2 = c o n s t .
It is useful to recognize that k defines the system length scale (k l 1 ). Equation (13), then, implies that the electrostatic energy density behaves as E 2 l q where q denotes the dimensionality of the process. Coupled with Equation (12), we find that the high frequency pressure, P h f E 2 δ n / ( 24 π k 2 λ D 2 n 0 ) E 2 l q [20], suppresses the thermal pressure, P t h = k B T 0 δ n δ n k 2 l 2 only for three-dimensional geometry. Inside the magnetosphere, however, the effective geometry is one-dimensional ( q = 1 ) since the particles are in the frozen-in condition and follow the field lines. In this region, thus, the collapse is impossible ( P h f ∼l 1 and P t h ∼l 2 ); it can be realized only in the outer regions of the magnetospheres, for distances exceeding the LC radius.
Zakharov showed that for three-dimensional geometry, q = 3 , the driven electrostatic field and the corresponding length scale of the cavern behave as [15]
E E 0 t 0 t 0 t
l l 0 t 0 t 0 t 2 / 3 ,
where t 0 is the time when the electric field fully collapses, and E 0 4 π n e Δ r exp Γ P / 4 is the electrostatic field before the collapse starts (the initial field amplified by means of the electrostatic instability) and Δ r R l c / ( 2 γ m a x ) represents a length scale in the LC zone where the process of energy pumping occurs [2] and we have taking into account that the escape time of particles equals P / 4 [21]. It is clear from the preceding equations that in due course of time, the length scale of the cavern goes to zero, whereas the electrostatic field asymptotically increases.
This process is terminated by means of Landau damping, when l reaches the dissipation length scale, l d 2 π λ D ,20]. From Equations (14) and (15) one can show that the electrostatic field will be boosted by the factor Δ r / l d 3 / 2 . Correspondingly, a certain fraction, κ , of the electrical energy converts to heat, increasing the temperature of the ambient plasma,
2 κ π R l c l d H E 2 8 π 4 σ c T 4 π R 0 2 H ,
where (in the framework of the equipartition approach) it has been assumed that the energy pumped by means of the collapse nearby the LC zone in a thin layer, l d (being of the order of ( 3 5 ) × 10 3 cm), is uniformly distributed in a cylindrical area with typical radius, R 0 ; the energy is emitted away as a black body radiation. Here we assume that the typical heated radius is of the order of 0.01 pc (an approximate outer radius of the accretion disc for the given mass of the BH [22], and H represents the height of the corresponding cylinder. From Equation (16) one can, straightforwardly, obtain the increased temperature of the AGN ambient
T c κ l d R l c E 2 16 π σ R 0 2 1 / 4 .
The proposed heating mechanism is efficient, if the corresponding temperature exceeds that of the accretion disk temperature, which is estimated as [22]
T 0 3 G M M ˙ 8 π σ R 3 1 / 4 R R l c 3 / 4 1 R R l c 1 / 4 ,
where R 2 G M / c 2 is the Schwarzschild radius of the supermassive black hole, M ˙ = L / ( ϵ c 2 ) is the accretion rate and ϵ 1 describes efficiency of the accretion process.
In Figure 1 we show the dependance of T 0 and T on the bolometric luminosity of AGN. The set of parameters is: a = 0.1 , ϵ = 0.1 , κ = 0.5 , M 8 = 1 , γ 1 = 7 × 10 3 , γ 2 = 10 4 , and R 0 0.01 pc. As is evident from the plots, the centrifugally driven heating mechanism becomes efficient for highly luminous AGN. In particular, for AGN with L > > 3 × 10 44 erg/s the temperature becomes much higher than the initial value ( T 0 ∼10 5 K) acquired in the accretion processes. On the other hand, one can straightforwardly check that the heating process is not sensitive with the rotation rate, which is a direct result of the power 1 / 4 (see Equation (17)). Here we assumed the equipartition distribution of energy: half of the pumped energy goes to acceleration and half-to the heating process.
The heating mechanism considered in the present work is so efficient that it can heat the system to blackbody radiation temperatures 10 9 K, which have been observed in the X-ray corona of AGNs [18].
It is worth noting that during the acceleration the particles might lose energy by means of the several cooling mechanisms. This problem is considered in detail by [2]. The synchrotron process almost from the very beginning of motion becomes insignificant because the corresponding radiation is so efficient that the particles very soon transit to the ground Landau level, follow the magnetic field lines and do not emit any more in the synchrotron regime.
The curvature radiation does not impose any significant constraints as well. In particular, on the last stage of LLCD the particles achieve such high energies that the magnetic field becomes dependent on plasma energy density, the particles move along straight trajectories and as a result the curvature radiation is terminated.
The inverse Compton (IC) scattering for the considered energies occurs in the Klein–Nishina regime, when the corresponding timescale is a continuously increasing function of the proton energy [23], therefore this mechanism does not limit maximum achievable energies.
Another mechanism, which potentially might limit the process of proton acceleration is the photo-pion energy losses. The corresponding timescale is a continuously decreasing function of energy and for the extremely high values 10 19 20 eV is of the order of 10 7 s. On the other hand, the timescale of energy pumping into Langmuir waves is by several orders of magnitude less than the aforementioned value. The collapse in turn, since characterized by blowing up (see Equations (14) and (15)) is even more efficient than any other mechanisms. Therefore, the photo-pion cooling does not impose any constraints on proton acceleration and consequently on the heating process.

4. Summary

By examining the LLCD mechanism we have considered the system of equations composed by the Euler equation, continuity equation and the Poisson equation respectively. It has been shown that by means of the centrifugal force, the Langmuir waves parametrically amplify. As a result, the excited electrostatic waves efficiently pump energy from rotation.
Amplification of electrostatic field on the next stage is provided by means of the Langmuir collapse, which terminates on the Debye length-scales, resulting in the heating of the surrounding medium. It has been shown that for luminosities greater than 2 × 10 45 erg/s the mentioned process can provide temperatures in the following interval 10 6 9 K.

Author Contributions

Z.O. and S.M. developed the analytical model, obtained the results and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are available in the article and can be accessed via a DOI link.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
LLCDThe Langmuir-Landau-Centrifugal Drive
AGNactive galactic nuclei
LCLight cylinder
ICinverse Compton

References

  1. Rieger, F.M.; Mannheim, K. Particle acceleration by rotating magnetospheres in active galactic nuclei. Astron. Astrophys. 2000, 353, 473. [Google Scholar]
  2. Osmanov, Z.; Mahajan, S.; Machabeli, S.; Chkheidze, N. Extremely efficient Zevatron in rotating AGN magnetospheres. Mon. Not. R. Astron. Soc. 2014, 445, 4155–4160. [Google Scholar] [CrossRef] [Green Version]
  3. Fermi, E. On the Origin of the Cosmic Radiation. Phys. Rev. 1949, 75, 1169. [Google Scholar] [CrossRef]
  4. Bell, A.R. The acceleration of cosmic rays in shock fronts—I. Mon. Not. R. Astron. Soc. 1978, 182, 147–156. [Google Scholar] [CrossRef] [Green Version]
  5. Bell, A.R. The acceleration of cosmic rays in shock fronts—II. Mon. Not. R. Astron. Soc. 1978, 182, 443–455. [Google Scholar] [CrossRef]
  6. Osmanov, Z. Relativistic Effects of Rotation in γ-ray Pulsars—Invited Review. arXiv 2021, arXiv:2101.05584. [Google Scholar]
  7. Osmanov, Z.; Mahajan, S.; Machabeli, S. On the Rotationally Driven Pevatron in the Center of the Milky Way. Astrophys. J. 2017, 835, 164. [Google Scholar] [CrossRef] [Green Version]
  8. Osmanov, Z.; Mahajan, S.; Machabeli, S.; Chkheidze, N. Millisecond newly born pulsars as efficient accelerators of electrons. Nat. Sci. Rep. 2015, 5, 14443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Gold, T. Rotating Neutron Stars and the Nature of Pulsars. Nature 1969, 221, 25. [Google Scholar] [CrossRef]
  10. Rieger, F.M. Nonthermal Processes in Black Hole-Jet Magnetospheres. Int. J. Mod. Phys. D 2011, 20, 1547–1596. [Google Scholar] [CrossRef] [Green Version]
  11. Osmanov, Z.; Rogava, A.S.; Bodo, G. On the efficiency of particle acceleration by rotating magnetospheres in AGN. Astron. Astrophys. 2007, 470, 395–400. [Google Scholar] [CrossRef] [Green Version]
  12. Machabeli, G.; Osmanov, Z.; Mahajan, S. Parametric mechanism of the rotation energy pumping by a relativistic plasma. Phys. Plasmas 2005, 12, 062901. [Google Scholar] [CrossRef] [Green Version]
  13. Osmanov, Z. Centrifugally driven electrostatic instability in extragalactic jets. Phys. Plasmas 2008, 15, 032901. [Google Scholar] [CrossRef] [Green Version]
  14. Mahajan, S.; Machabeli, G.; Osmanov, Z.; Chkheidze, N. Ultra High Energy Electrons Powered by Pulsar Rotation. Nat. Sci. Rep. 2013, 3, 1262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Zakharov, V.E. Collapse of Langmuir Waves. Sov. J. Exp. Theory Phys. 1972, 35, 908–914. [Google Scholar]
  16. Galeev, A.A.; Sagdeev, R.Z.; Shapiro, V.D.; Shevchenko, V.I. Langmuir turbulence and dissipation of high-frequency energy. Sov. J. Exp. Theory Phys. 1977, 46, 711. [Google Scholar]
  17. Degtiarev, L.M.; Zakharov, V.E.; Rudakov, L.I. Dynamics of Langmuir collapse. Sov. J. Plasma Phys. 1976, 2, 240–246. [Google Scholar]
  18. Beckmann, V.; Shrader, C. Active Galactic Nuclei; Wiley-Vch: New York, NY, USA, 2012. [Google Scholar]
  19. Osmanov, Z. Efficiency of the centrifugally induced curvature drift instability in AGN winds. Astron. Astrophys. 2008, 490, 487–492. [Google Scholar] [CrossRef] [Green Version]
  20. Artsimovich, L.A.; Sagdeev, R.Z. Plasma Physics for Physicists; Atomizdat: Moscow, Russia, 1979. [Google Scholar]
  21. Machabeli, G.Z.; Rogava, A.D. Centrifugal force: A gedanken experiment. Phys. Rev. A 1994, 50, 98. [Google Scholar] [CrossRef] [PubMed]
  22. Carroll, B.W.; Ostlie, D.A. An Introduction to Modern Astrophysics and Cosmology; Pearson: San Francisco, CA, USA, 2010. [Google Scholar]
  23. Blumenthal, G.R.; Gould, R.J. Bremsstrahlung, Synchrotron Radiation, and Compton Scattering of High-Energy Electrons Traversing Dilute Gases. Rev. Mod. Phys. 1970, 42, 237. [Google Scholar] [CrossRef]
Figure 1. Behavior of T (solid line) and T 0 (dashed line) with respect to the bolometric luminosity of AGN. The set of parameters is: a = 1 , ϵ = 0.1 , κ = 0.5 , M 8 = 1 , γ 1 = 7 × 10 3 , γ 2 = 10 4 and R 0 0.01 pc.
Figure 1. Behavior of T (solid line) and T 0 (dashed line) with respect to the bolometric luminosity of AGN. The set of parameters is: a = 1 , ϵ = 0.1 , κ = 0.5 , M 8 = 1 , γ 1 = 7 × 10 3 , γ 2 = 10 4 and R 0 0.01 pc.
Universe 07 00083 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Osmanov, Z.; Mahajan, S. On the Heating of AGN Magnetospheres. Universe 2021, 7, 83. https://doi.org/10.3390/universe7040083

AMA Style

Osmanov Z, Mahajan S. On the Heating of AGN Magnetospheres. Universe. 2021; 7(4):83. https://doi.org/10.3390/universe7040083

Chicago/Turabian Style

Osmanov, Zaza, and Swadesh Mahajan. 2021. "On the Heating of AGN Magnetospheres" Universe 7, no. 4: 83. https://doi.org/10.3390/universe7040083

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop