Next Article in Journal
Rabi Frequency Management of Collapsing Quasi-Two-Dimensional Bose-Einstein Condensates with Pseudospin-1/2
Previous Article in Journal
Electrostatic Simulations for the DUNE ND-GAr Field Cage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Protoneutron Star Magnetized Envelops in Neutrino Energy Spectra

by
Vladimir N. Kondratyev
1,2,*,
Tamara D. Lobanovskaya
2 and
Dimash B. Torekhan
2
1
Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research (JINR), 141980 Dubna, Russia
2
Nuclear Physics Department, Dubna State University, University Str., 19, 141982 Dubna, Russia
*
Author to whom correspondence should be addressed.
Submission received: 21 November 2021 / Revised: 3 April 2022 / Accepted: 6 April 2022 / Published: 12 April 2022

Abstract

:
The neutrino dynamics in hot and dense magnetized matter, which corresponds with protoneutron star envelopes in the core collapse supernova explosions, is considered. The kinetic equation for a neutrino phase space distribution function is obtained, taking into account inelastic scattering by nuclear particles. The transfer component in a momentum space using transport properties is studied. The energy transfer coefficient is shown to change from positive to negative values when the neutrino energy exceeds four times the matter temperature. In the vicinity of a neutrino sphere, such effects are illustrated to lead to the energy strengthening in the neutrino spectra. As this paper demonstrates, such a property is favorable for the possibility of observing supernova neutrino fluxes using Large Volume Neutrino Telescopes.

1. Introduction

Type II supernovae (SNe) represent one of the most energetic events, as they can, plausibly, give high energy cosmic rays, sites for synthesis of heavy nuclides (e.g., e-, s-, p-, r-process atomic nuclei), renew other nuclear components, and so on. Moreover, the SN explosion mechanism constitutes an important issue that is still under debate. In particular, the mechanism of energy transfer to all-star matter (initially bound) represents the main SN problem. Neutrino flux and/or magnetic pressures plausibly make a key contribution to explosive shock wave formation, in addition to core bounce and thermal pressure, respectively.
Being in the vicinity of neutrino sphere strong convection [1], and/or magnetorotational instability [2], can trigger a turbulent dynamo leading to a magnetic field amplification of up to tens of teratesla (TT). Prompt SN explosions can be associated with magnetic pressure contributing significant energy into the stellar matter, and bringing the predominant mechanism of shock wave formation, respectively. Another possible explosion mechanism is related to neutrino heating and reviving the stopped shock wave using the neutrinos and antineutrinos that are emitted by a cooling protoneutron star [3,4,5]. Neutrino absorption of nucleons gives rise to the increasing temperature and pressure of the matter behind the shock that starts to expand, pushing the shock forward. The efficiency of such a delayed CCSN mechanism depends on neutrino luminosity and the hardness of their spectra [4,5].
Since neutrinos or magnetic pressure are capable of making a significant contribution to the mechanism of supernova explosions, an analysis of neutrino transport in supernova matter, taking into account magnetic effects, represents an important issue. In addition, the possible magnetic influence on neutrino spectra are crucial for the interpretation of the r- and neutrino processes that can be also affected by magnetic fields [6,7,8,9]. As has been recently shown [10,11], the energy exchange in neutrino–nuclear scattering can be noticeably enhanced due to the magnetic field. In the next section we consider neutrino kinetics, paying particular attention to the transfer in momentum space, on the basis of an energy transfer cross section [11]. As is illustrated in Section 3, such an approach gives a clear picture of the influence of inelastic scattering on neutrino spectra. Possibilities to observe the possible effects of this using Large Volume Neutrino Telescopes are analyzed in Section 4. Conclusions are given in Section 5.

2. Neutrino Kinetics

To describe neutrino kinetics, we use quite a general kinetic equation for the phase–space distribution function f(r, p, l)
df d l = f l + z f r + p l f p = Λ + St [ f ] ,
where a distance passed by neutrino l = c t with speed of light c and time t, r and p representing the partial derivatives with respect to the spatial, r, and momentum coordinates, p = z E/c, with the unit vector z defining the direction of neutrino momentum at energy E. Here, we take into account that since supernova neutrinos possess typical energies in the MeV range, much larger than the experimental rest-mass limit for active flavors <1 eV, they essentially propagate with the speed of light c. The momentum derivative in Equation (1), p / l accounts for an energy and momentum exchange during the neutrino collisions with environment particles. On the right hand side of Equation (1), Λ stands for all the rates of neutrino production, absorption, and annihilation, whereas the term St[f], accounts for fluctuations in scattering processes.

3. Neutrino Dynamics in Magnetized Neutrino Sphere

In this work, we concentrate on neutrino kinetics in a region of neutrino decoupling from matter for dynamo active SNe. In such neutrino sphere regions, the neutrino dynamics change from diffusive to free streaming. Spectra of neutrinos emerging from a protoneutron star can be parameterized by the following equation:
W ( E , T ) = E 2   d   f ( r ,   p ,   l ) ~   E a exp { ( 2 + α ) E / E av }
Here, Ω denotes the solid angle of vector z, E a v is an average energy, and α is a numerical parameter describing the amount of spectral pinching; the value α = 2 corresponds to a Maxwell–Boltzmann spectrum, and α = 2.3 to a Fermi–Dirac distribution with zero chemical potential.
Beyond the protoneutron star surface in a neutrino sphere region, it is impossible to maintain both the chemical equilibrium between neutrinos and stellar matter and diffusion; however, a noticeable energy exchange between neutrinos and strongly magnetized stellar material can affect neutrino spectra. Since heavy-leptonic neutrinos only interact with star matter through a neutral current, they are energetically less coupled to stellar plasma than the electron flavor when neutral and charged currents are involved; therefore, the heavy-leptonic neutrinos break out of thermal equilibrium in the energy sphere, which is significantly deeper inside the nascent protoneutron star than the transport sphere (near to neutrino sphere), where the transition from diffusion to free flow breaks in. Within the scattering atmosphere, respective heavy-leptonic neutrinos still frequently collide with neutrons and protons. As is demonstrated in Section 3.2 and Section 3.3, in this case, magnetic effects noticeably enhance the energy exchange in neutrino–nucleon scattering due to the neutral current.

3.1. Neutrino Sphere Properties

Matter in the neutrino sphere region corresponds to a moderate density n ~ 0.1–10 Tg cm−3 (1 Tg = 1012 g) and temperature T ~ 5–10 MeV. We assume strong fluctuations of temperature T and density n in this region since it meets a strong convection of matter and corresponds with the vicinity of the bifurcation point for stellar material between a collapse of the central compact object and supernova ejecta. Figure 1a shows the Fermi energy of nucleons   E F N and electrons E F e versus a beta equilibrium parameter Ye at density n = 1 Tg/cm3. One sees that at the realistic numbers of the beta equilibrium parameter Ye ~ 0.2–0.3, these values for nucleons (protons-   E F p ~ 0.6 MeV, neutrons-   E F n ~ 1.1 MeV) and electrons ( E F e ~ 35 MeV) are small and large when compared with temperature, respectively. Therefore, nucleon components, with   E F N « T, represent non-degenerate gas, whereas an electron gas, with E F e T, is strongly degenerated. As a consequence, a neutrino–electron scattering cross section is strongly suppressed because of the Pauli principle. Such a blocking effect also leads to the actual termination of a charged current component in neutrino–nucleon scattering. Magnetization gives rise to the effective increase in Fermi energy and further diminution of respective scattering. Moreover, the corresponding mean free path (mfp) rises up to 10ths km at considered densities; therefore, we hereafter neglect the right hand side of Equation (1). On the contrary, neutrino–nucleon scattering due to a neutral current component can be considered an independent process with corresponding mfp lf = (NN σGT0)−1 ~ 100 m. Here, Ni = ni/mi represents the number density of i-th nuclear particle (N denotes nucleon), with mass mi and contribution ni to the total mass density n, and σGT0 denotes the respective cross section σGT0  10−40 cm2 (E/37 MeV)2; see [4,5].
The energy transfer cross section in a magnetized neutrino sphere for ν + Nν′ + N′ scattering has been considered by Kondratyev et al. [11]. We briefly recall that this value is defined as S1i = −∫dε ε (dσiνν/dε) with an energy transfer ε differential cross section (dσiνν/dε). The nucleon energy levels with spin magnetic moments directed along (spin up) and opposite (spin down) the magnetic field are split on the value Δ = |gα| μNH ≡ |gα|ωL because of an interaction with field H. Here, μN denotes the nuclear magneton, ωL = μNH gives the Larmour frequency, and gα is the nucleon g-factor. At the temperature T for neutral GT0 neutrino nucleon scattering, the energy transfer cross section reads [11]
S1 ≈ σGT0 Δ (2δE − (1 + δ2E) th(δT/2))|Δ<E,T
σ GT 0   Δ   ( 2 δ E   δ T / 2 )  
where δE = Δ/E, δT = Δ/2 T and th(x) is the hyperbolic tangent. For a magnetized nucleon gas, the average transferred energy in inelastic neutrino scattering is evaluated as being <δE> = −<ε> ≈ S1/σGT0. This value depends on temperature T, splitting Δ in the nucleon gas, and the energy of the incoming neutrino E. Figure 1b shows the average transferred energy <δE> for a magnetized nucleon gas in units of Δ as a function of neutrino energy E at various temperatures T. As is seen in Figure 1b, the average transferred energy varies from a positive value for a hot nucleon gas to a negative value for a cold system. This change corresponds to a transition from exo-energetic to endo-energetic neutrino scattering, which occurs under the conditions corresponding to the temperature TE/4, see [11]. The physical reason of such a transition is obviously because of the decreasing thermal population of the upper split energy level of the nucleons, which results in a suppression of the contribution of GT0 transitions from this level to the underlying level. The condition of such a transition from one mode to another is well described by the relation of E ≈ 4T, which is independent on splitting value Δ.

3.2. Energy Transfer in Neutrino Spectra

Making use of Equation (4), we determine an energy transfer coefficient as
E l = i N i S i E ( 1 E 4 T ) / l t ,
where the sum goes over nuclear particles and the energy transfer length reads lt = E2/(∑i Ni σiGT02 ≈ 100 m (3 MeV/Δ)2 (10 Tg cm−3/n). Neglecting the right hand side of Equation (1) for the uniform flux z  f / r = 0, the solution of Equation (1) is given by replacing E with the solution of Equation (5); see Appendix A, i.e.,
EνelE (el + (1 − el) E/4T)−1
with el = exp{l/lt}.
Figure 2 shows the energy transfer effect in neutrino energy spectra during evolution in the vicinity of the neutrinospheric region. The Maxwell–Boltzmann distribution corresponds to α = 2 and Eav = 10 MeV in Equation (2) and is taken as the initial one. One sees that the energy transfer effect leads to an increase in neutrino energy at the maximum distribution. When neutrino path l approaches a mean energy transfer length lt, we obtain a spreading in distribution W(E) with the maximum point increasing in a near linear manner with a growing el. Such an acceleration is particularly effective at larger gas temperatures.
The case of a single effective collision in Figure 2 clarifies the obtained results. In this case, the relationship between corresponding exo- and endo-energetic regimes is determined by a ratio of the occupation of respective nucleon levels and neutrino phase space volume in the exit channel (i.e., Exp{δT} (1 − δE)2θ(1 − δE)/(1 + δE)2, with step function θ(x)). When this ratio is less than 1, the number of endo-energetic collisions is larger than the number of exo-energetic collisions, and vice versa; therefore, for neutrino dynamics in magnetized nucleon gas, it is preferable that a change of acceleration will occur, and that the stopping regimes correspond to the condition δT = −2 ln{(θ(1 − δE)(1 − δE)/(1 + δE)}|δET<1 ≈ 4δE. Such a switch in neutrino dynamics is also displayed in Equations (3) and (4). For large energies E > 4T, such effects proceed faster, compared with small energies E < 4T, because of the strong energy dependence of the energy transfer coefficient.

3.3. Fluctuation Effects in Energy Spectra

The strong convection in the vicinity of the neutrino sphere and bifurcation point gives rise to large fluctuations in the properties of respective stellar materials. We average the results of the energy spectra modification over the fluctuations. For temperature T, we assume a uniform distribution in a range from 5 to 10 MeV, which is independent from density fluctuations. As is seen in Figure 3a, the maximum distribution W(E) is shifted towards larger energies, approaching the region of 10–20 MeV. The properties of such averaged energy distribution resemble the results of 10 MeV, which support an effective acceleration mechanism at higher temperatures.

4. Strengthening of Spectrum Hardness and Large Neutrino Detector Sensitivity

The spectrum hardness can be characterized quantitatively using the average energy of particles:
<E> = ∫dE W(E) Eel Eav (1 + (el − 1) 3Eav/16T)−1
In Equation (7), we used the Saddle Point method justified at el ~ 1. As can be seen in Figure 3b, the value <E> increases almost linearly with the increasing path length l. One sees an effective energy absorption by neutrinos from magnetized stellar material. Such a feature can result in double (or multiple) peaked SN light curves. The growing hardness of neutrino energy spectra makes it a favorable feature of the considered effect for detector sensitivity.
Strongly variable transient particle fluxes can be detected using large-volume neutrino telescopes: KM3NeT [12], Baikal-GVD [13], and so on. Sensitivity to neutrinos on a scale of 10 MeV can be achieved by observing a collective increase, using multiple detectors to measure the rate of counting coincidences. A sharp increase in the spatially uniform neutrino flux Φ(t) is associated with a CCSN infall phase that occurs during the half a second [14] that defines the observation time. Assuming a spherically uniform neutrino emission, one obtains Φ(t) = L(t)/4πd2 with the neutrino luminosity L, and the distance to the source d. The CCSN neutrino detection rate, rSN(t), can be evaluated as
rSN(t) = φ Φ(t) Σi ni ∫dε W(ε) σi(ε),
where the sum index i ∈ {p, e, 16O} represents the most important target components producing energetic charged particles (i.e., e+/−, ni is the number of targets, σi(ε) the total interaction cross section for the given target i, and W(ε) gives the energy spectrum from Equation (3)). The detector efficiency φ corresponds to the ratio between the number of detected events and the number of interacting neutrinos per unit of water [15].
At time t of time interval δt, the probability of a triggered detector p = r(t) δt/Π, where Π is the number of detectors and total the detection rate r = rSN + rB includes the background event rate rB. The multiple coincidences of k detectors meet with the probability given by Poissonian law pk/k! ep. In this case, the ratio signal/background is given by (1 + rSN/rB) k  (1 + k rSN/rB). Evidently, the k-fold coincidence enhances by a factor of k, which is the detection sensitivity for a weak SN neutrino signal. When the condition (k rSN/rB) ≫ 1 corresponds to value k approaching ten, then an excess of hundreds or thousands in terms of the total number of detectors Π is required for d~10 kpc, respectively.

5. Conclusions

We considered energy transfer for neutrino nuclear scattering in strong magnetic fields, which may plausibly arise in supernovae, and its respective effect in neutrino energy spectra. Nuclear magnetization is shown to bring new, neutral-current induced reaction channels, giving additional noticeable mechanisms to the dynamic of neutrinos being weakly coupled with matter. The energy transfer coefficient in kinetic equations is demonstrated to change from positive to negative values with increasing neutrino energy. For magnetized nondegenerate nucleon gas, such cross over between acceleration and stopping regimes occurs when neutrino energy is about factor four of the gas temperature, while the nucleon Larmour frequency is sufficiently small. This switching in dynamical properties originates from the detailed balance principle, and a difference of phase space volume for neutrinos in the final channel when scattering spin-up and spin-down nucleons, and independent on splitting value ∆. Consequently, such a property is insensitive to the magnetization geometry. The respective acceleration and/or stopping rates are determined by the product of splitting ∆ and scattering cross section σGT0 in the nucleon gas. For realistic properties of stellar material (see Section 3.1), such neutrino–nuclear scattering effects result in an increase in the hardness of neutrino energy spectra. Since electronic neutrinos decouple from matter at the neutrino-sphere, and thereafter, they experience a few (single in average) effective collisions, the corresponding acceleration effect is relatively small. Beyond the energy sphere, the dynamics of heavy-leptonic neutrinos is mainly governed by collisions with nucleons. Within the scattering atmosphere (up to the neutrino-sphere) these collisions are frequent enough to maintain spatial diffusion for heavy-leptonic neutrinos. The significant completed path l, within a magnetized region of a star, leads to a considerable acceleration effect in the case of a heavy-leptonic component. The strengthening of neutrino hard energies is favorable for supernova neutrino observations that use Large Volume Neutrino Telescopes. In this case, the CCSN neutrino flux is revealed from an increase in the counting rate of multiple detectors that are coincidentally triggered. As illustrated, the k-fold coincidence enhances the detection sensitivity for a weak SN neutrino signal by factor k. Finally, we notice that such strong magnetization also arises in neutron star mergers, magnetar crusts, and heavy-ion collisions.

Author Contributions

Conceptualization and methodology, V.N.K.; formal analysis, investigation, V.N.K., D.B.T. and T.D.L.; writing—original draft preparation, V.N.K.; writing—review and editing, V.N.K., D.B.T. and T.D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

To find neutrino energy El after passing a distance l in neutrinospheric region we rewrite Equation (5) as
dEl/[El (1 − El/4T)] = dl/lt, El|l=0 = E
The solution is given by
El/(El − 4T) = [E/(E − 4T)]exp{l/lt},
and leads to Equation (6).

References

  1. Mösta, P.; Ott, C.D.; Radice, D.; Roberts, L.; Schnetter, E.; Haas, R. A large-scale dynamo and magnetoturbulence in rapidly rotating core-collapse supernovae. Nature 2015, 528, 376–379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Guilet, J.; Müller, E. Numerical simulations of the magnetorotational instability in protoneutron stars—I. Influence of buoyancy. Mon. Not. R. Astron. Soc. 2015, 450, 2153–2171. [Google Scholar] [CrossRef]
  3. Bethe, H.A.; Wilson, H.A. Revival of a stalled supernova shock by neutrino heating. Astrophys. J. 1985, 295, 14. [Google Scholar] [CrossRef]
  4. Bethe, H.A. Supernova mechanisms. Rev. Mod. Phys. 1990, 62, 801. [Google Scholar] [CrossRef]
  5. Janka, H.-T.; Melson, T.; Summa, T. Physics of core-collapse supernovae in three dimensions: A sneak preview. Annu. Rev. Nucl. Part. Sci. 2016, 66, 341–375. [Google Scholar] [CrossRef] [Green Version]
  6. Kondratyev, V.N.; Korovina, Y.V. Synthesis of chemical elements in dynamo active supernovae. JETP Lett. 2015, 102, 131. [Google Scholar] [CrossRef]
  7. Kondratyev, V.N. Properties and composition of magnetized nuclei. Particles 2020, 3, 272–277. [Google Scholar] [CrossRef] [Green Version]
  8. Kondratyev, V.N. Neutron capture reactions in strong magnetic fields of magnetars. Phys. Rev. C 2004, 69, 038801. [Google Scholar] [CrossRef]
  9. Kondratyev, V.N. R-Process with magnetized nuclei at dynamo-explosive supernovae and neutron star mergers. Universe 2021, 7, 487. [Google Scholar] [CrossRef]
  10. Horowitz, C.J. Weak magnetism for antineutrinos in supernovae. Phys. Rev. D 2002, 65, 043001. [Google Scholar] [CrossRef] [Green Version]
  11. Kondratyev, V.N.; Dzhioev, A.A.; Vdovin, A.I.; Cherubini, S.; Baldo, M. Energy exchange in neutrino nuclear scattering. Phys. Rev. C 2019, 100, 045802. [Google Scholar] [CrossRef]
  12. Adrián-Martínez, S.; Ageron, M.; Aharonian, F.; Aiello, S.; Albert, A.; Ameli, F.; Anassontzis, E.; Andre, M.; Androulakis, G.; Anghinolfi, M.; et al. Letter of intent for KM3NeT 2.0. J. Phys. G Nucl. Part. Phys. 2016, 43, 084001. [Google Scholar] [CrossRef]
  13. Avrorin, A.D.; Avrorin, A.V.; Aynutdinov, V.M.; Bannash, R.; Belolaptikov, I.A.; Brudanin, V.B.; Budnev, N.M.; Doroshenko, A.A.; Domogatsky, G.V.; Dvornický, R.; et al. Baikal-GVD. EPJ WoC 2017, 136, 04007. [Google Scholar] [CrossRef] [Green Version]
  14. The Garching Core-Collapse Supernova Research. Available online: https://wwwmpa.mpa-garching.mpg.de/ccsnarchive/ (accessed on 7 April 2022).
  15. Aiello, S.; Albert, A.; Garre, S.A.; Aly, Z.; Ambrosone, A.; Ameli, F.; Andre, M.; Androulakis, G.; Anghinolfi, M.; Anguita, M.; et al. The KM3NeT potential for the next core-collapse supernova observation with neutrinos. Eur. Phys. J. C 2021, 81, 445. [Google Scholar] [CrossRef]
Figure 1. Neutrino sphere properties: (a) Fermi Energy at density n = 1 Tg cm−3 versus beta equilibrium parameter Ye. (b) The neutrino energy dependence of an average transferred energy <δE> during inelastic scattering at temperatures T = 5 and 10 MeV, corresponding to the curves 1 and 2. The approximate relationship in Equations (3) and (4) coincide within the line widths.
Figure 1. Neutrino sphere properties: (a) Fermi Energy at density n = 1 Tg cm−3 versus beta equilibrium parameter Ye. (b) The neutrino energy dependence of an average transferred energy <δE> during inelastic scattering at temperatures T = 5 and 10 MeV, corresponding to the curves 1 and 2. The approximate relationship in Equations (3) and (4) coincide within the line widths.
Particles 05 00011 g001
Figure 2. Neutrino energy spectra in arbitrary units for α = 2, Eav = 10 MeV, and el = 1 (i.e., l = 0)—curves 1, 2—curves 2, and 3—curves 3 at T = (a) 5 MeV (b) 10 MeV. The curves 4 correspond to a respective single effective collision event at ∆ = 2 MeV.
Figure 2. Neutrino energy spectra in arbitrary units for α = 2, Eav = 10 MeV, and el = 1 (i.e., l = 0)—curves 1, 2—curves 2, and 3—curves 3 at T = (a) 5 MeV (b) 10 MeV. The curves 4 correspond to a respective single effective collision event at ∆ = 2 MeV.
Particles 05 00011 g002
Figure 3. (a) Neutrino energy spectra in arbitrary units for α = 2, Eav = 10 MeV, and el = 1 (i.e., l = 0)—curve 1, 2—curve 2, and 3—curve 3 average over T from 5 MeV to 10 MeV. (b) The average neutrino energy depending on the path length l. For α = 2, Eav = 10 MeV, T = 5 MeV—curve 1 and 10 MeV—curve 2. Curve 3 corresponds to the average over temperature T in the range of 5 to 10 MeV.
Figure 3. (a) Neutrino energy spectra in arbitrary units for α = 2, Eav = 10 MeV, and el = 1 (i.e., l = 0)—curve 1, 2—curve 2, and 3—curve 3 average over T from 5 MeV to 10 MeV. (b) The average neutrino energy depending on the path length l. For α = 2, Eav = 10 MeV, T = 5 MeV—curve 1 and 10 MeV—curve 2. Curve 3 corresponds to the average over temperature T in the range of 5 to 10 MeV.
Particles 05 00011 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kondratyev, V.N.; Lobanovskaya, T.D.; Torekhan, D.B. Effect of Protoneutron Star Magnetized Envelops in Neutrino Energy Spectra. Particles 2022, 5, 128-134. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5020011

AMA Style

Kondratyev VN, Lobanovskaya TD, Torekhan DB. Effect of Protoneutron Star Magnetized Envelops in Neutrino Energy Spectra. Particles. 2022; 5(2):128-134. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5020011

Chicago/Turabian Style

Kondratyev, Vladimir N., Tamara D. Lobanovskaya, and Dimash B. Torekhan. 2022. "Effect of Protoneutron Star Magnetized Envelops in Neutrino Energy Spectra" Particles 5, no. 2: 128-134. https://0-doi-org.brum.beds.ac.uk/10.3390/particles5020011

Article Metrics

Back to TopTop