pith. sign in

arxiv: 2604.25013 · v1 · submitted 2026-04-27 · 🌌 astro-ph.SR · math-ph· math.MP

The Singular Behaviour of Ambipolar Diffusion Revealed by 1D Cartesian Solutions

F. Moreno-Insertis , E. R. Priest , D. N\'obrega-Siverio This is my paper

Pith reviewed 2026-05-08 01:24 UTC · model grok-4.3

classification 🌌 astro-ph.SR math-phmath.MP
keywords ambipolar diffusionmagnetic null pointsstagnation point floweigenmodesflux annihilationMHD code testingsolar chromosphere
0
0 comments X

The pith

Ambipolar diffusion forms three distinct layers that transfer magnetic flux uniformly around a null point.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives exact solutions for ambipolar diffusion in a simplified one-dimensional setup near magnetic null points. It identifies a steady stagnation-point flow in which magnetic flux moves inward at a constant rate: first by advection from far away, then through an ambipolar-diffusion layer where the field strength rises as the cube root of distance from the null, and finally through an inner Ohmic layer where the field is linear and flux is destroyed. Separate eigenmode solutions of the diffusion equation are found, both symmetric and antisymmetric, that contain sharp current sheets; time-dependent runs show that these modes relax toward the simplest allowable state. The same solutions are used to verify that an MHD code reproduces the expected behavior to high accuracy.

Core claim

A stagnation-point flow solution was found with a uniform flux transfer rate across three regions: an external advection region; an internal ambipolar diffusion region with magnetic profile B ∝ x^(1/3); and an innermost Ohmic region with B ∝ x; in the latter, flux annihilation occurs at a rate imposed by the advection. Both symmetric and antisymmetric eigenmode solutions to the ambipolar diffusion problem are found with sharp current sheets at the internal nulls. The time evolution of the eigenmodes shows how higher-order or perturbed modes evolve toward the lowest-order allowable eigenmodes.

What carries the argument

The self-similar ansatz applied to the one-dimensional ambipolar diffusion equation, solved with phase-plane methods to obtain steady stagnation-point profiles and nonlinear eigenmodes.

Load-bearing premise

The one-dimensional Cartesian self-similar form and the restriction that the flow stays steady or follows only the identified eigenmodes continue to describe the physics when the setup is made three-dimensional and the ionization fraction is allowed to vary.

What would settle it

A simulation or observation that measures the magnetic-field profile in a partially ionized plasma near a null and finds no interval where field strength scales as distance to the power one-third would falsify the three-layer ambipolar solution.

Figures

Figures reproduced from arXiv: 2604.25013 by D. N\'obrega-Siverio, E. R. Priest, F. Moreno-Insertis.

Figure 1
Figure 1. Figure 1: A sketch of the solution for steady-state ambipolar diffusion join￾ing B = B1 at x = −L0 to B = B0 at x = L0 for three different values of B1, namely, 2B0/3 (top curve), −2B0/3 (middle curve), and −B0 (bot￾tom curve) view at source ↗
Figure 2
Figure 2. Figure 2: The solution to the mixed Ohmic and ambipolar diffusion problem of Eqs. (11)-(12), in which the inset shows the whole range −L0 ≤ x ≤ L0) and the main frame shows just the central domain −0.018L0 ≤ x ≤ 0.018L0)]. Green: the exact solution (16); blue: the purely ambipolar solution (17); dashed red: the linear, purely Ohmic solution with slope at the origin equal to that of the full solution (19); pink: the … view at source ↗
Figure 3
Figure 3. Figure 3: Upper panel: a schematic of a stagnation-point flow configu￾ration, showing magnetic field lines (light-headed arrows) for a one￾dimensional current sheet together with a stagnation-point flow (solid￾headed arrows). The ideal, ambipolar and Ohmic regions are indicated. Middle and lower panels: solution of the dimensionless stagnation flow equation (Eq. 29) with logarithmic (middle panel) and linear axes (l… view at source ↗
Figure 4
Figure 4. Figure 4: (Top) The fundamental (i.e., ZKBP) and first four symmetric har￾monics as functions of the stretched coordinate ξˆ = H 1/2 ξ. The eigen￾values Kˆ = K/H are indicated. (Bottom) Spatial profile of the dimen￾sionless ambipolar diffusive flux for the eigenfunctions shown in the top panel. The locations of the null crossings of the solution are marked with vertical dashed lines. The boundary conditions (41) ser… view at source ↗
Figure 5
Figure 5. Figure 5: (Top) The fundamental (i.e., dipole) and first three antisymmet￾ric harmonics as functions of the stretched coordinate ξˆ = H 1/2 ξ. The eigenvalues Kˆ = K/H are indicated. (Bottom) Spatial profile of the di￾mensionless diffusive flux D for the eigenfunctions shown in the top panel. The locations of the null crossings of the solution are marked with vertical dashed lines. 4. Mode stability: a few time-depe… view at source ↗
Figure 6
Figure 6. Figure 6: Magnetic field evolution starting from a superposition of the symmetric first harmonic and a zero net-flux finite perturbation which does not add any extra internal null. The panels show four selected in￾stants (top: initial state; bottom: asymptotic evolution). Note the dif￾ferent ranges of the ordinate axes in the different panels. The evolving self-similar first symmetric harmonic is shown as a dashed l… view at source ↗
Figure 7
Figure 7. Figure 7: Selected stages in the evolution of the numerical solution of the ambipolar diffusion equation with initial condition having the shape of the 7th symmetric self-similar harmonic and with the numerical mesh containing 2048 grid points. Red: the exact solution for the 7th self￾similar harmonic. Black: actual numerical solution. Cyan (only in the three lowermost panels): exact solution for the 1st self-simila… view at source ↗
Figure 8
Figure 8. Figure 8: Selected stages in the evolution of the numerical solution de￾scribed in Sect. 4.2. Panel (a): initial condition (black solid curve), re￾sulting from the addition of the exact first symmetric self-similar har￾monic (red dashed) and a finite antisymmetric perturbation of zero net flux (green dashed). In panels (d) - (f), the exact self-similar dipole har￾monic evolving towards the final shape of the numeric… view at source ↗
Figure 9
Figure 9. Figure 9: A comparison of the Bifrost runs (black thick curve) with the theoretical solutions (red curves), showing the initial condition (left panel) and the final configuration for nx = 512 (central panel). The relative mutual deviation of the maxima between the theoretical and Bifrost solutions for four resolution levels is given in the right panel view at source ↗
Figure 10
Figure 10. Figure 10: Time-dependence of the maximum of the magnetic field (left panel), the location of the singular points [middle panel; upper curve for L(t), lower curve for Lcs(t)] and the value of the integral R B(x, t) dx in the central lobe of the solution (right panel) for the solution of the first harmonic, illustrated in view at source ↗
Figure 11
Figure 11. Figure 11: A comparison of the Bifrost solution (black thick curve) for the third harmonic with the theoretical solution (red curve), showing the initial (left) and final (middle) profiles for a resolution of nx = 512. Right panel: the relative mutual deviation of the maxima between the theoretical and Bifrost solutions for the four resolution levels used in the test. state at time t = 100 (τs = 1206) for the lowest… view at source ↗
read the original abstract

Aims. We seek to (a) study 1D Cartesian ambipolar diffusion near null points; (b) characterise the nonlinear eigenmodes for ambipolar diffusion; (c) propose tests for ambipolar diffusion solvers in MHD codes. Methods. (a) Direct analysis is used to find analytical solutions for ambipolar diffusion. (b) To study the eigenmodes, we solve the ODE for self-similar solutions of the 1D ambipolar diffusion equation using phase-plane techniques. We also solve the general time-dependent 1D problem for initial conditions of interest. (c) We test the Bifrost code by trying to reproduce the behaviour of the eigenmodes. Results. (a) A stagnation-point flow solution was found with a uniform flux transfer rate across three regions: an external advection region; an internal ambipolar diffusion region with magnetic profile B propto x**(1/3); and an innermost Ohmic region with B propto x; in the latter, flux annihilation occurs at a rate imposed by the advection. (b) Both symmetric and antisymmetric eigenmode solutions to the ambipolar diffusion problem are found with sharp current sheets at the internal nulls. The time evolution of the eigenmodes (pure or perturbed) is probed, showing how higher-order eigenmodes, or perturbed ones, evolve in time towards the lowest-order allowable eigenmodes. (c) The Bifrost code reproduces the behaviour of the eigenmodes with excellent accuracy. Conclusions. Stagnation-point configurations exist with ambipolar diffusion carrying magnetic flux in an inner layer and serving as an intermediary between the external advection and an Ohmic-diffusion core around the null. Our tests are compatible with the hypothesis that zero-flux higher harmonics of the self-similar equation evolve toward either the first symmetric or antisymmetric harmonic. The self-similar solutions can serve as strong tests for ambipolar diffusion solvers in general MHD codes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript derives analytical 1D Cartesian solutions for ambipolar diffusion near magnetic null points under a stagnation-point flow, identifying a matched three-region structure with uniform flux transfer: an external advection-dominated region, an ambipolar-diffusion interior with B ∝ x^{1/3}, and an innermost Ohmic core with B ∝ x in which annihilation is imposed by the outer advection. It further obtains symmetric and antisymmetric nonlinear eigenmodes of the self-similar ambipolar ODE via phase-plane methods, demonstrates their time-dependent relaxation toward the lowest-order allowable modes, and reports that the Bifrost code reproduces the eigenmode evolution with high accuracy.

Significance. If the derivations hold, the work supplies exact, parameter-free analytical profiles and eigenmodes that constitute strong, reproducible benchmarks for ambipolar-diffusion solvers in MHD codes. The explicit three-region flux-transfer construction and the documented relaxation of higher harmonics toward the fundamental modes provide concrete, falsifiable predictions for flux transport near nulls in partially ionized plasmas.

minor comments (3)
  1. [§2.2] §2.2, Eq. (7): the self-similar ansatz for the velocity field is introduced without an explicit statement of the assumed functional form for v(x,t); adding this would make the reduction to the ODE fully transparent.
  2. [Figure 4] Figure 4: the time-evolution panels for perturbed eigenmodes would be clearer if the initial perturbation amplitude and the final relaxed state were indicated by overlaid reference curves.
  3. [§5] §5: the claim that the solutions 'serve as strong tests' would be strengthened by a short table comparing the analytically predicted decay rates of the first few eigenmodes with the numerically measured rates from Bifrost.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive and accurate summary of our manuscript, the assessment of its significance as providing benchmarks for ambipolar-diffusion solvers, and the recommendation for minor revision. No specific major comments were listed in the report for us to address point by point.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper obtains its central results (stagnation-point three-region solution with B ∝ x^{1/3} in the ambipolar layer and B ∝ x in the Ohmic core, plus symmetric/antisymmetric eigenmodes) by direct analytical integration of the governing ODEs and phase-plane analysis of the self-similar ambipolar equation, without any fitted parameters, self-referential definitions, or load-bearing self-citations. The time-dependent relaxation toward lowest-order eigenmodes is shown by explicit numerical integration of the 1D PDE. The Bifrost reproduction constitutes an external numerical check independent of the analytic construction. No step in the derivation chain reduces to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on the standard MHD equations with an ambipolar diffusion term and the assumption of self-similar 1D Cartesian solutions; no free parameters are fitted to data, no new entities are postulated, and no ad-hoc axioms beyond standard mathematics are introduced.

axioms (1)
  • domain assumption The governing equations admit self-similar solutions that can be reduced to an ODE solvable by phase-plane analysis.
    Invoked to obtain the eigenmodes and stagnation-point flow.

pith-pipeline@v0.9.0 · 5672 in / 1383 out tokens · 31920 ms · 2026-05-08T01:24:25.926973+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

71 extracted references · 71 canonical work pages

  1. [1]

    D., Botha, G

    Arber, T. D., Botha, G. J. J., & Brady, C. S. 2009, ApJ, 705, 1183

    work page 2009

  2. [2]

    D., Haynes, M., & Leake, J

    Arber, T. D., Haynes, M., & Leake, J. E. 2007, ApJ, 666, 541

    work page 2007

  3. [3]

    L., Alexeev, I., Collados, M., et al

    Ballester, J. L., Alexeev, I., Collados, M., et al. 2018, Space Sci. Rev., 214, 58

    work page 2018

  4. [4]

    1952, Prikl

    Barenblatt, G. 1952, Prikl. Mat. i Mekh, 16, 67

    work page 1952

  5. [5]

    Bauer, S. J. 1973, Physics of planetary ionospheres (Springer, Berlin)

    work page 1973

  6. [6]

    Bernis, F., Hulshof, J., & Vázquez, J. 1993, J. reine angew Math., 435, 1

    work page 1993

  7. [7]

    Braginskii, S. I. 1965, Reviews of Plasma Physics, 1, 205

    work page 1965

  8. [8]

    & Zweibel, E

    Brandenburg, A. & Zweibel, E. G. 1994, Astrophys J Lett, 427, L91

    work page 1994

  9. [9]

    Cally, P. S. & Khomenko, E. 2018, Astrophys. J., 856, 20

    work page 2018

  10. [10]

    2024, A&A, 685, A2

    Cheng, G., Ni, L., Chen, Y ., & Lin, J. 2024, A&A, 685, A2

    work page 2024

  11. [11]

    Cheung, M. C. M. & Cameron, R. H. 2012, ApJ, 750, 6

    work page 2012

  12. [12]

    2023, ApJ, 952, 21

    Chouliaras, G., Syntelis, P., & Archontis, V . 2023, ApJ, 952, 21

    work page 2023

  13. [13]

    1957, Magnetohydrodynamics, Interscience tracts on physics and astronomy (Interscience Publishers)

    Cowling, T. 1957, Magnetohydrodynamics, Interscience tracts on physics and astronomy (Interscience Publishers)

    work page 1957

  14. [14]

    Crutcher, R. M. 1999, ApJ, 520, 706 De Pontieu, B., Martínez-Sykora, J., & Chintzoglou, G. 2017, Astrophys. J. Letts., 849, L7

    work page 1999

  15. [15]

    Draine, B. T. 1986, Mon. Not. Roy. Astron. Soc., 220, 133

    work page 1986

  16. [16]

    Draine, B. T. & McKee, C. F. 1993, Ann. Rev. Astron. Astrophys., 31, 373 González-Morales, P. A., Khomenko, E., Downes, T. P., & de Vicente, A. 2018, A&A, 615, A67

    work page 1993

  17. [17]

    Goodman, M. L. 2000, Astrophys. J., 533, 501

    work page 2000

  18. [18]

    Goodman, M. L. 2004, Astron. Astrophys., 416, 1159

    work page 2004

  19. [19]

    1979, Quarterly of Applied Mathematics, 37, 259

    Grundy, R. 1979, Quarterly of Applied Mathematics, 37, 259

    work page 1979

  20. [20]

    V ., Carlsson, M., Hansteen, V

    Gudiksen, B. V ., Carlsson, M., Hansteen, V . H., et al. 2011, A&A, 531, A154

    work page 2011

  21. [21]

    & Zweibel, E

    Heitsch, F. & Zweibel, E. G. 2003, ApJ, 583, 229

    work page 2003

  22. [22]

    & André, P

    Hennebelle, P. & André, P. 2013, Astron. Astrophys., 560, A68

    work page 2013

  23. [23]

    2019, Physics of Plasmas, 26, 082902

    Hillier, A. 2019, Physics of Plasmas, 26, 082902

    work page 2019

  24. [24]

    Hillier, A. S. 2024, Philosophical Transactions of the Royal Society of London Series A, 382, 20230229

    work page 2024

  25. [25]

    S., Snow, B., & Luna, M

    Hillier, A. S., Snow, B., & Luna, M. 2024, Philosophical Transactions of the Royal Society of London Series A, 382, 20230230

    work page 2024

  26. [26]

    Hulshof, J. 1991, J. Math. Anal. Appl., 157, 75

    work page 1991

  27. [27]

    & Vázquez, J

    Kamin, S. & Vázquez, J. 1991, SIAM J. Math. Anal., 22, 34

    work page 1991

  28. [28]

    2023, A&A, 673, A66

    Keppens, R., Popescu Braileanu, B., Zhou, Y ., et al. 2023, A&A, 673, A66

    work page 2023

  29. [29]

    & Collados, M

    Khomenko, E. & Collados, M. 2012, Astrophys. J., 747, 87

    work page 2012

  30. [30]

    Khomenko, E., Collados, M., & Díaz, A. J. 2016, Astrophys. J., 823, 132

    work page 2016

  31. [31]

    2014, Astron

    Khomenko, E., Díaz, A., de Vicente, A., Collados, M., & Luna, M. 2014, Astron. Astrophys., 565, A45

    work page 2014

  32. [32]

    2018, A&A, 618, A87

    Khomenko, E., Vitas, N., Collados, M., & de Vicente, A. 2018, A&A, 618, A87

    work page 2018

  33. [33]

    & Roberts, B

    Kumar, N. & Roberts, B. 2003, Solar Phys., 214, 241

    work page 2003

  34. [34]

    Kunz, M. W. & Mouschovias, T. C. 2010, MNRAS, 408, 322

    work page 2010

  35. [35]

    Leake, J. E. & Arber, T. D. 2006, Astron. Astrophys., 450, 805

    work page 2006

  36. [36]

    Leake, J. E. & Linton, M. G. 2013, Astrophys. J., 764, 54

    work page 2013

  37. [37]

    E., Lukin, V

    Leake, J. E., Lukin, V . S., & Linton, M. G. 2013, Physics of Plasmas, 20, 061202

    work page 2013

  38. [38]

    E., Lukin, V

    Leake, J. E., Lukin, V . S., Linton, M. G., & Meier, E. T. 2012, Astrophys. J., 760, 109

    work page 2012

  39. [39]

    E., et al

    Li, Z.-Y ., Banerjee, R., Pudritz, R. E., et al. 2014, Protostars and Planets VI, 173 Mac Low, M.-M. & Klessen, R. S. 2004, Reviews of Modern Physics, 76, 125 Martínez-Sykora, J., De Pontieu, B., Carlsson, M., et al. 2017a, Astrophys. J., 847, 36 Martínez-Sykora, J., De Pontieu, B., De Moortel, I., Hansteen, V . H., & Carlsson, M. 2018, ApJ, 860, 116 Mart...

    work page 2014

  40. [40]

    & Spitzer, L

    Mestel, L. & Spitzer, L. 1956, Mon. Not. Roy. Astron. Soc., 116, 503

    work page 1956

  41. [41]

    & Kruger, C

    Mitchner, M. & Kruger, C. H. 1973, Partially Ionized Gases, Wiley series in plasma physics (Wiley)

    work page 1973

  42. [42]

    R., & Hood, A

    Moreno-Insertis, F., Nóbrega-Siverio, D., Priest, E. R., & Hood, A. W. 2022, Astron. Astrophys., 662, A42

    work page 2022

  43. [43]

    Mouschovias, T. C. & Paleologou, E. V . 1979, ApJ, 230, 204

    work page 1979

  44. [44]

    & Li, Z.-Y

    Nakamura, F. & Li, Z.-Y . 2008, ApJ, 687, 354

    work page 2008

  45. [45]

    2022, A&A, 665, A116

    Ni, L., Cheng, G., & Lin, J. 2022, A&A, 665, A116

    work page 2022

  46. [46]

    2015, ApJ, 799, 79

    Ni, L., Kliem, B., Lin, J., & Wu, N. 2015, ApJ, 799, 79

    work page 2015

  47. [47]

    I., & Schmieder, B

    Ni, L., Lin, J., Roussev, I. I., & Schmieder, B. 2016, ApJ, 832, 195 Nóbrega-Siverio, D., Martínez-Sykora, J., Moreno-Insertis, F., & Carlsson, M. 2020a, A&A, 638, A79 Nóbrega-Siverio, D., Moreno-Insertis, F., Martínez-Sykora, J., Carlsson, M., &

    work page 2016

  48. [48]

    Parker, E. 1973, J. Plasma Phys., 9, 49

    work page 1973

  49. [49]

    Parker, E. N. 1963, Astrophys. J. Suppl, 8, 177

    work page 1963

  50. [50]

    Pattle, R. E. 1959, The Quarterly Journal of Mechanics and Applied Mathemat- ics, 12, 407

    work page 1959

  51. [51]

    Pfaff, R. F. 2012, Space Sci. Rev., 168, 23 Popescu Braileanu, B. & Keppens, R. 2021, A&A, 653, A131

    work page 2012

  52. [52]

    1982, Solar Magnetohydrodynamics (Dordrecht: D

    Priest, E. 1982, Solar Magnetohydrodynamics (Dordrecht: D. Reidel)

    work page 1982

  53. [53]

    Priest, E. R. 2014, Magnetohydrodynamics of the Sun (Cambridge, UK: Cam- bridge University Press)

    work page 2014

  54. [54]

    & Przybylski, D

    Rempel, M. & Przybylski, D. 2021, ApJ, 923, 79

    work page 2021

  55. [55]

    S., Ballai, I., Khomenko, E., & Collados, M

    Ruderman, M. S., Ballai, I., Khomenko, E., & Collados, M. 2018, Astron. As- trophys., 609, A23

    work page 2018

  56. [56]

    Sakai, J. I. & Smith, P. D. 2009, Astrophys. J. Letts., 691, L45

    work page 2009

  57. [57]

    Shu, F. H. 1983, ApJ, 273, 202

    work page 1983

  58. [58]

    Snow, B., Botha, G. J. J., McLaughlin, J. A., & Hillier, A. 2018, A&A, 609, A100

    work page 2018

  59. [59]

    J., Ballester, J

    Soler, R., Díaz, A. J., Ballester, J. L., & Goossens, M. 2012, ApJ, 749, 163

    work page 2012

  60. [60]

    Soler, R., Oliver, R., & Ballester, J. L. 2010, Astron. Astrophys., 512, A28

    work page 2010

  61. [61]

    & Priest, E

    Sonnerup, B. & Priest, E. 1975, J. Plasma Phys., 14, 283

    work page 1975

  62. [62]

    1962, Physics of Fully Ionized Gases (New York: Interscience)

    Spitzer, L. 1962, Physics of Fully Ionized Gases (New York: Interscience)

    work page 1962

  63. [63]

    2025, A&A, 700, A152 Vázquez, J

    Tritsis, A. 2025, A&A, 700, A152 Vázquez, J. 2007, The Porous Medium Equation. Mathematical Theory, Oxford Mathematical Monographs (Oxford: Clarendon Press)

    work page 2025

  64. [64]

    P., & De Pontieu, B

    Vranjes, J., Poedts, S., Pandey, B. P., & De Pontieu, B. 2008, Astron. Astrophys., 478, 553

    work page 2008

  65. [65]

    M., Martínez-Sykora, J., Hansteen, V

    Wargnier, Q. M., Martínez-Sykora, J., Hansteen, V . H., & Pontieu, B. D. 2022, The Astrophysical Journal, 933, 205

    work page 2022

  66. [66]

    V ., Khodachenko, M

    Zaqarashvili, T. V ., Khodachenko, M. L., & Rucker, H. O. 2011, Astron. Astro- phys., 529, A82 Zel’dovich, Y . B. & Barenblatt, G. 1957, Appl. Math. Mech., 21, 718 Zel’dovich, Y . B. & Kompaneets, A. 1950, in Collection of Papers dedicated to the 70th birthday of A.F.Ioffe, ed. P. I. Lukirskii (Moscow: Izd. Akad. Nauk. USSR), 61–72 Zel’dovich, Y . B. & Ra...

    work page 2011

  67. [67]

    2011, Journal of Computational Physics, 230, 1035

    Ziegler, U. 2011, Journal of Computational Physics, 230, 1035

    work page 2011

  68. [68]

    Zweibel, E. G. 1989, Astrophys. J., 340, 550

    work page 1989

  69. [69]

    Zweibel, E. G. 1994, in NATO Advanced Science Institutes (ASI) Series C, V ol. 422, NATO Advanced Science Institutes (ASI) Series C, ed. D. Lynden-Bell, 73

    work page 1994

  70. [70]

    Zweibel, E. G. & Brandenburg, A. 1997, Astrophys. J., 478, 563

    work page 1997

  71. [71]

    Zweibel, E. G. & Yamada, M. 2009, Ann. Rev. Astron. Astrophys., 47, 291 Article number, page 16

    work page 2009