REVIEW 4 major objections 3 minor 49 references
Kinetic simulations of fusion ignition with hot-spot ablator mix
T0 review · 4 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Localized carbon ablator mix is worse for fusion ignition than the same mass spread uniformly: the spike cools, contracts, and becomes an alpha-particle barrier that also inflates inferred ion temperatures.
desk verdict Plausible mechanism for localized ablator mix severity, but the 1D slab geometry overstates the alpha-barrier component and the quantitative claim needs a geometry check. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the localized carbon mix spike, a fully ionized carbon population embedded in a deuterium-tritium hot spot. The argument is carried by a truncated kinetic ion model in which each species' velocity distribution is written as an isotropic part $f_0$ plus a first-order anisotropic part $f_1$ (Eq. 1); this makes two-dimensional kinetic simulation feasible while still capturing non-Maxwellian tails. The load-bearing mechanism is radiative cooling and contraction: carbon's bremsstrahlung losses cool the spike, it compresses to several times the surrounding density, the $n_e^2$ scaling of radiation accelerates the cooling, and the contracted spike becomes an $\alpha$-particle stopper, a heat sink, and a driver of fluid flows. A second mechanism is a boundary-layer kinetic effect in which long-mean-free-path fast fuel ions are absorbed at the spike's temperature gradient, depleting the distribution tail and lowering the fusion rate.
What would settle it
Run the same two-dimensional kinetic setup with the carbon jet replaced by a uniform distribution of the same 1.9% carbon mass; the paper's claim predicts the uniform case still self-heats toward ignition while the localized case does not. If the uniform case burns no better, or if a seeded carbon spike of the same mass produces no line-of-sight neutron-temperature difference, the central claim would be contradicted.
Extended reading notes
Core claim
On its own terms, the paper's central discovery is that localized mix is more detrimental than uniform mix because of its radiative cooling and contraction. A fully ionized carbon spike totalling 1.9% of the hot-spot mass, initialized as a 5-micrometer Gaussian, cools the surrounding electrons and then the ions through increased bremsstrahlung; the spike contracts to several times the hot-spot density, and since radiative losses scale as $n_e^2$, the contraction feeds back into faster cooling. The dense cool spike then has high $\alpha$ stopping power, so fast $\alpha$ particles born elsewhere in the hot spot slow down inside it and deposit their energy where it is radiated away; the spike therefore acts like an internal shell, partially separating the two sides of the hot spot and reducing the effective areal density below the threshold value. A separate kinetic effect arises because fast fuel ions have long mean free paths and are absorbed at the spike's temperature gradient, depleting the tail of the deuterium-tritium distribution and cutting the fusion reactivity by a few percent, up to about 10% near the edges. Finally, the cooling-driven inflow of fuel broadens the fusion neutron spectrum: synthetic spectra fitted with a thermal model give 4464 eV along one line of sight and 4526 eV along another, versus 4026 eV when fluid motion is neglected, an inflation comparable to the unexplained experimental discrepancy.
Load-bearing premise
The simulation depends on the truncated kinetic expansion of Eq. (1), which is valid only while the fluid velocity is small compared to the thermal velocity; the two-dimensional neutron-broadening result is analyzed at t = 40 ps, when the fluid kinetic-to-total energy ratio has reached about 0.1, near the edge of that assumption.
Editorial extensions
If this is right
- Ignition thresholds depend on the spatial structure of ablator mix, not just its mass: a carbon load that is tolerable when spread evenly can extinguish ignition when concentrated in a jet or spike.
- Localized mix lowers the effective hot-spot areal density by acting as an internal alpha-particle barrier, so a target whose total areal density meets the ignition criterion may still fail if a mix spike is present.
- Inferred ion temperatures from neutron spectra can be inflated by several hundred electronvolts and vary with line of sight when mix-driven flows are present, complicating comparisons between neutron and X-ray diagnostics.
- In a hot spot with multiple mix jets, the regions between jets are expected to have fusion reactivity reduced by roughly 10% relative to a Maxwellian plasma, so burn-rate estimates from fluid codes will run high.
- The fusion burn shape and the temperature or X-ray emission shape should differ in a mixed hot spot: the burn profile is flattened relative to the temperature profile, which could be used to distinguish localized from uniform mix.
Reading between the lines
- This suggests a testable diagnostic strategy: because the neutron-temperature excess and its line-of-sight variation track the radial inflow into the spike, measuring those variations on shots with known mix asymmetries could quantify mix localization independent of X-ray imaging.
- By the same radiative-contraction logic, higher-Z ablator dopants such as tungsten should produce a comparable internal barrier at smaller mass fractions, since bremsstrahlung scales with the square of the ion charge; this could be tested by rerunning the same setup with the carbon charge replaced.
- If the kinetic tail depletion around a mix spike is as large as the two-dimensional run suggests, radiation-hydrodynamic codes that assume Maxwellian reactivity will systematically overestimate yield in mixed implosions; a simple local correction for mix edges could close much of that gap without full kinetic simulation.
- The contraction timescale of tens of picoseconds implies mix damage worsens as the hot spot evolves, so time-resolved neutron or X-ray measurements over a 40 ps window should show growing broadening; if the broadening saturates instead, the feedback loop in the model may be weaker than assumed.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper applies a truncated ion Vlasov-Fokker-Planck model (Eqs. (1)-(15)) to the stagnated ICF hot-spot with carbon ablator mix. In 1D Cartesian geometry, it compares a Gaussian carbon spike carrying 1.9% of the hot-spot mass with the same carbon mass uniformly distributed, and finds that the localized mix cools radiatively, contracts, increases alpha stopping and radiative losses, and acts as an alpha barrier, so ignition is suppressed in the localized case but not in the uniform-mix case. A 2D jet simulation is used to estimate fluid-flow-induced broadening of the DT neutron spectrum, giving inferred ion temperatures about 10% higher than the thermal value and a line-of-sight asymmetry. Kinetic corrections to fusion reactivity near the mix boundaries are post-processed using Eq. (16), yielding reductions from several percent near the mix up to about 20% at the hot-spot edge.
Significance. If the conclusions hold, they are important for ICF mix physics: localized ablator mix is a plausible experimental morphology, and the predicted neutron-spectrum broadening and line-of-sight asymmetry offer falsifiable diagnostic signatures. The paper is transparent about several model assumptions, such as the optically thin radiation treatment and the estimated reabsorption optical depth of about 0.2, and it specifies numerical parameters sufficiently for reproduction. The explicit synthetic neutron-spectrum calculation and the direct comparison of two mix morphologies are strengths. The paper does not provide machine-checked proofs or public code, but the equations are stated in enough detail that the central simulations are in principle reproducible.
major comments (4)
- [§2 and §5, Figs. 1 and 4] The headline comparison between localized and uniform mix is made only in one-dimensional Cartesian geometry, where the "localized" carbon is an infinite slab perpendicular to x. In this geometry every alpha particle crossing the hot-spot must pass through the mix region, which maximizes the barrier contribution to the claimed severity. The two-dimensional jet in Fig. 4 is still infinite in the third direction and is not compared against a uniform-mix control. Consequently the quantitative claim that "localized mix is more detrimental than the uniform mix" (abstract and §5) is not yet established for a finite three-dimensional Rayleigh-Taylor spike; alpha particles could travel around such a spike. The radiative-cooling/contraction mechanism may survive, but the comparison that supports the headline conclusion needs either a 2D uniform-mix control or a 3D simulation, or the claim must be explicitly restricted to slab/ridge geometries.
- [§2 and §3 (Eqs. (6), (16))] The kinetic reduction of fusion reactivity is not part of the simulated dynamics: the fusion source Fa in Eq. (6) uses the Maxwellian reactivity of ref. [35], while the kinetic rate from Eq. (16) is computed as a post-processing diagnostic from the stored distributions. Thus the abstract's statement that "a purely kinetic effect reduces fusion reactivity by several percent" is a diagnostic estimate, not a self-consistent prediction of yield or ignition threshold. The simulations should either feed the kinetic rate back into the fusion source term, or the text should clearly label the reactivity reduction as a one-way diagnostic and avoid drawing conclusions about ignition-threshold changes from it.
- [§4 and §5] No convergence study or error estimate is reported for any of the quantitative outputs. The runs use 288 cells and a 3000-point velocity grid with a 20 fs time-step in one dimension and equivalent parameters in two dimensions, but there is no demonstration that the 10-20% reactivity reductions, the 4.0-4.5 keV neutron temperatures, or the localized-vs-uniform ignition contrast are converged with respect to spatial resolution, velocity resolution, or time-step. Because several of the headline numbers are of order ten percent, a resolution study, or at least a two-resolution comparison, is needed before those numbers can be taken at face value.
- [§2 and §5 (Fig. 4b)] The neutron-spectrum analysis is performed at t=40 ps, the time at which Fig. 4b shows the deuterium-tritium fluid kinetic energy reaching about 10% of its total energy. Using the relation between this ratio and the fluid velocity, this corresponds to |u| on the order of half the thermal velocity, which is not deeply inside the regime where retaining only f0 and f1 in Eq. (1) is clearly accurate, and the flow contribution to the synthetic spectra is exactly the quantity being measured. The authors should check the sensitivity of the inferred temperatures to higher-order terms, for example by testing a more conservative case or by reporting the velocity-to-thermal-speed ratio explicitly, and should discuss whether the expansion in Eq. (1) remains controlled at the time used for the main diagnostic claim.
minor comments (3)
- [Title page] The title in the manuscript text reads "ablat or mix" and should read "ablator mix."
- [Eq. (15)] Please state the units of the numerical coefficient 1.69e-32 so that the expression can be evaluated without guessing; the text gives electron temperature in eV and densities in cm^-3 but not the coefficient's units.
- [Fig. 4 caption] The color-bar label "EK/U" is not defined in the caption; the text defines the ratio as mn|u|^2/(2U), so the caption should define EK explicitly.
Circularity Check
No significant circularity: all headline results are emergent simulation outputs; the kinetic-reactivity and neutron-broadening claims are explicitly post-processed diagnostics.
full rationale
Walking the derivation chain: the localized-vs-uniform comparison is a simulation output generated from explicitly stated initial conditions (Gaussian carbon spike vs uniform carbon of equal mass) and the evolution equations (1)-(15); the claimed radiative cooling, contraction, increased alpha stopping power, and neutron-spectrum broadening are emergent, not encoded in the inputs. The fusion reactivity reduction is not used as a simulation input—the paper states that 'The fusion term Fa in equation (6) used the Maxwellian fusion reactivity formula [35]' and the kinetic reactivity is post-processed via Eq. (16) from the simulated distribution, so the claimed reduction is a diagnostic of the simulated distribution, not a fitted quantity. Similarly, the neutron-broadening temperature overestimate follows by comparing synthetic spectra with and without fluid velocity using Eq. (17), and the statement that the magnitude 'is similar to the known discrepancy at the National Ignition Facility [3]' is an ex post comparison, not an inverse fit. The only overlapping-author citation is [33] (Tzoufras, Bell, Norreys, Tsung) for the standard truncated Fokker-Planck and Krook collision forms; it is not used as an unverified premise or uniqueness theorem, and the model equations are described in the text. The 1D/2D Cartesian geometry limits extrapolation to finite 3D mix spikes, but that is an external-validity concern, not circularity.
Assumptions & free parameters
free parameters (4)
- Hot-spot initial conditions =
T=4.9 keV, density=100 g/cm^3, areal density=0.35 g/cm^2
- Carbon mix profile =
1.9% of hot-spot mass, Gaussian waist 5 um (1D and 2D azimuth), 20 um radial in 2D
- Flux limiter fraction =
0.05
- Coulomb logarithm =
~3
assumptions (7)
- domain assumption Diffusive approximation truncating the distribution after first-order anisotropy (Eq. 1)
- domain assumption Krook form for the anisotropic collision operator C1a (Eq. 13)
- domain assumption Electron fluid model with quasi-neutrality and electrostatic approximation
- domain assumption Optically thin bremsstrahlung radiation with open radiation boundaries
- domain assumption Bremsstrahlung emission formula from Atzeni and Meyer-ter Vehn (Eq. 15)
- standard math Maxwellian fusion reactivity (Bosch-Hale) used for evolving the distribution
- domain assumption Carbon is fully ionized
Cite this review
Pith. "Pith review of Kinetic simulations of fusion ignition with hot-spot ablator mix." pith.science (2026). https://pith.science/paper/PQGC62SC
@misc{pith2026190805577,
author = {Pith},
title = {Pith review of: Kinetic simulations of fusion ignition with hot-spot ablator mix},
year = {2026},
howpublished = {\url{https://pith.science/paper/PQGC62SC}},
note = {Machine review of arXiv:1908.05577}
}
read the original abstract
Inertial confinement fusion fuel suffers increased X-ray radiation losses when carbon from the capsule ablator mixes into the hot-spot. Here we present one and two-dimensional ion Vlasov-Fokker-Planck simulations that resolve hot-spot self heating in the presence of a localized spike of carbon mix, totalling 1.9% of the hot-spot mass. The mix region cools and contracts over tens of picoseconds, increasing its alpha particle stopping power and radiative losses. This makes a localized mix region more severe than an equal amount of uniformly distributed mix. There is also a purely kinetic effect that reduces fusion reactivity by several percent, since faster ions in the tail of the distribution are absorbed by the mix region. Radiative cooling and contraction of the spike induces fluid motion, causing neutron spectrum broadening. This artificially increases the inferred experimental ion temperatures and gives line of sight variations.
Figures
Reference graph
Works this paper leans on
-
[35]
A. R. Bell, R. G. Evans, and D. J. Nicholas, Phys. Rev. Lett. 46, 243 (1981)
work page 1981
-
[1]
is no longer valid. Inserting equation ( 1) into the non-relativistic Vlasov- Fokker-Planck equation and averaging over angles leads to the kinetic equations for each ion species [32, 33], ∂f0a ∂t + v 3 ∇.f 1a + qa 3mav2 ∂ ∂v (v2E.f 1a) = C0a + Fa, (6) ∂f1a ∂t + v∇f0a + qa ma E ∂f0a ∂v = C1a. (7) In addition to these ion equations, the electrons are model...
- [2]
-
[3]
is substituted in to the Fokker-Planck form given by equation (2) in ref. [8], it yields the velocity dependent Krook form [30, 33], C1a = − ∑ b Γ abnb v3 (f1a − fMa) , (13) fMa = −uDT ∂f0a ∂v . (14) Although this is an approximation to the full Fokker- Planck form, it is most accurate for the supra-thermal ions responsible for the fusion reactivity reduc...
-
[4]
S. Atzeni and J. Meyer-ter Vehn, The Physics of Inertial Fusion (Oxford University Press, Oxford, 2004)
work page 2004
-
[5]
is small compared to the second term. This rules out mod- elling of the implosion (as previously achieved with one- dimensional ion-kinetic codes [31]) but does allow effi- cient two-dimensional modelling of the stagnated hot- spot. As such, simulations will be initialized in the time of peak fuel compression. This approximation also breaks down during the ...
-
[6]
used the Maxwellian fusion reactivity formula [35], giving the vol- umetric reaction rate r as a function of temperature. There may be kinetic corrections to this formula, so the more accurate reactivity will be post-processed by inte- grating the full hydrogen distribution function from the simulation. The hydrogen ion density was reduced and the alpha p...
-
[7]
This increases the stopping power for alpha par- ticles. Furthermore, the energy deposited in the mix re- gion is more efficiently radiated away and wasted, since the radiative rate scales as n2 e. Once it cools and contracts, the mix region acts as an effective heat sink and barrier for the alpha particles, partially separating the two sides of the hot-spot...
Show all 49 references
-
[8]
The fusion rate is reduced in Knudsen layers at the edge of the hot-spot and near the mix region
The volumetric fusion reaction rate was calcu- lated by numerically integrating the fusion cross section across the distribution function of the fuel ions in three velocity space dimensions, r(x) = ∫ d3v′ d3v fD(x, v′)fT(x, v)σ(vr)vr, (16) where vr = |v − v′| is the relative s...
-
[9]
Therefore the inflow of hot fuel into the mix region and dense fuel shell is enough to significantly broaden the fusion neutron distribution
to the spectra yields T = 4026 eV in the control case, T = 4464 eV along the x direction and T = 4526 eV in the y direction. Therefore the inflow of hot fuel into the mix region and dense fuel shell is enough to significantly broaden the fusion neutron distribution. The magnitud...
2014
-
[10]
Le Pape, L
S. Le Pape, L. F. Berzak Hopkins, L. Divol, A. Pak, E. L. Dewald, S. Bhandarkar, L. R. Bennedetti, T. Bunn, J. Biener et al ., Phys. Rev. Lett. 120, 245003 (2018)
2018
-
[11]
Kagan, D
G. Kagan, D. Svyatskiy, H. G. Rinderknecht, M. J. Rosenberg, A. B. Zylstra, C.-K. Huang and C. J. McDevitt, Phys. Rev. Lett. 115, 105002 (2015)
2015
-
[12]
O. A. Hurricane, D. A. Callahan, D. T. Casey, P. M. Cel- liers, C. Cerjan, E. L. Dewald, T. R Dittrich, T. D¨ oppner, D. E. Hinkel et al. , Nature 506, 343 (2014)
2014
-
[13]
H. Sio, J. A. Frenje, A. Le, S. Atzeni, T. J. T. Kwan, M. Gatu Johnson, G. Kagan, C. Stoeckl, C. K. Li et al. , Phys. Rev. Lett. 122, 035001 (2019)
2019
-
[14]
Sherlock and S
M. Sherlock and S. Rose, High Energy Density Physics 5, 27 (2009)
2009
-
[15]
The boundary conditions are open for the radiation and periodic for the plasma
is only valid for weakly coupled plasma and not the dense fuel shell, the model imposes a smooth exponential cut-off of W for regions of plasma with temperature below 1 .5 keV. The boundary conditions are open for the radiation and periodic for the plasma. The plasma is assumed...
-
[16]
Michta, F
D. Michta, F. Graziani, T. Luu, and J. Pruet, Phys. of Plasmas 17, 012707 (2010)
2010
-
[17]
The calculation was performed to find the broadening for lines of sight in the x and y directions [Fig
assumes that the neutrons are not scattered by collisions as they es- cape. The calculation was performed to find the broadening for lines of sight in the x and y directions [Fig. 5]. The numerical integration of equation ( 17) for the simulation data at t = 40 ps was compared ...
-
[18]
B. E. Peigney, O. Larroche and V. Tikhonchuk, J. of Comp. Phys. 278, 416 (2014)
2014
-
[19]
Molvig, N
K. Molvig, N. M. Hoffman, B. J. Albright, E. M. Nelson and R. B. Webster, Phys. Rev. Lett. 109, 095001 (2012)
2012
-
[20]
B. I. Cohen, A. M. Dimits, G. B. Zimmerman and S. C. Wilks, Phys. of Plasmas 21, 122701 (2014)
2014
-
[21]
Albright, K
B. Albright, K. Molvig, C.-K. Huang, A. Simakov, E. Dodd, N. Hoffman, G. Kagan and P. Schmit, Phys. of Plasmas 20, 122705 (2013)
2013
-
[22]
Inglebert, B
A. Inglebert, B. Canaud and O. Larroche, Europhysics Letters 107, 65003 (2014)
2014
-
[23]
T. Ma, P. K. Patel, N. Izumi, P. T. Springer, M. H. Key, L. J. Atherton, M. A. Barrios, L. R. Benedetti, R. Bionta et al. , Phys. of Plasmas 24, 056311 (2017)
2017
-
[24]
Larroche, H
O. Larroche, H. G. Rinderknecht and M. J. Rosenberg, Phys. Rev. E 98, 031201(R) (2018)
2018
-
[25]
Larroche, Phys
O. Larroche, Phys. of Plasmas 19, 122706 (2012)
2012
-
[26]
T. Ma, P. K. Patel, N. Izumi, P. T. Springer, M. H. Key, L. J. Atherton, L. R. Benedetti, D. K. Bradley, D. A. Callahan et al. , Phys. Rev. Lett. 111, 085004 (2013)
2013
-
[27]
S. P. Regan, R. Epstein, B. A. Hammel, L. J. Suter, H. A. Scott, M. A. Barrios, D. K. Bradley, D. A. Callahan, C. 8 Cerjan et al. , Phys. Rev. Lett. 111, 045001 (2013)
2013
-
[28]
T. Ma, O. A. Hurricane, D. A. Callahan, M. A. Barrios, D. T. Casey, E. L. Dewald, T. R. Dittrich, T. D¨ oppner, S. W. Haan et al. , Phys. Rev. Lett. 114, 145004 (2015)
2015
-
[29]
H. G. Rinderknecht, H. Sio, C. K. Li, A. B. Zylstra, M. J. Rosenberg, P. Amendt, J. Delettrez, C. Bellei, J. A. Frenje et al. , Phys. Rev. Lett. 112, 135001 (2014)
2014
-
[30]
H. G. Rinderknecht, H. Sio, C. K. Li, N. Hoffman, A. B. Zylstra, M. J. Rosenberg, J. A. Frenje, M. Gatu Johnson, F. H. S´ eguinet al. , Phys. of Plasmas 21, 056311 (2014)
2014
-
[31]
B. A. Hammel, H. A. Scott, S. P. Regan, C. Cerjan, D. S. Clark, M. J. Edwards, R. Epstein, S. H. Glenzer, S. W. Haan et al. , Phys. of Plasmas 18, 056310 (2011)
2011
-
[32]
V. A. Smalyuk, R. E. Tipton, J. E. Pino, D. T. Casey, G. P. Grim, B. A. Remington, D. P. Rowley, S. V. Weber, M. Barrios et al. , Phys. Rev. Lett. 112, 025002 (2014)
2014
-
[33]
N. B. Meezan, M. J. Edwards, O. A. Hurricane, P. K. Patel, D. A. Callahan, W. W Hsing, R. P. J. Town, F. Albert, P. A. Amendt et al. , Plasma Phys. and Cont. Fusion 59, 014021 (2016)
2016
-
[34]
D. S. Clark, C. R. Weber, J. L. Milovich, J. D. Salmonson, A. L. Kritcher, S. W. Haan, B. A. Ham- mel, D. E. Hinkel, O. A. Hurricane et al. , Phys. of Plas- mas 23, 056302 (2016)
2016
-
[36]
R. J. Henchen, M. Sherlock, W. Rozmus, J. Katz, D. Cao, J. P. Palastro, and D. H. Froula, Phys. Rev. Lett. 121, 125001 (2018)
2018
-
[37]
Y. T. Lee and R. More, Phys. of Fluids 27, 1273 (1984)
1984
-
[38]
M. N. Rosenbluth, W. M. MacDonald, and D. L. Judd, Phys. Rev. 107, 1 (1957)
1957
-
[39]
M. J. Keskinen, Phys. Rev. Lett. 103, 055001 (2009)
2009
-
[40]
W. T. Taitano, A. N. Simakov, L. Chacon, and B. Keenan, Phys. of Plasmas 25, 056310 (2018)
2018
-
[41]
Thomas, M
A. Thomas, M. Tzoufras, A. Robinson, R. Kingham, C. Ridgers, M. Sherlock, and A. Bell, J. of Comp. Phys. 231, 1051 (2012)
2012
-
[42]
Tzoufras, A
M. Tzoufras, A. Bell, P. Norreys, and F. Tsung, J. of Comp. Phys. 230, 6475 (2011)
2011
-
[43]
Chang and G
J. Chang and G. Cooper, J. of Comp. Phys. 6, 1 (1970)
1970
-
[44]
H. S. Bosch and G. Hale, Nuclear fusion 32, 611 (1992)
1992
-
[45]
C. R. Weber, D. S. Clark, A. W. Cook, D. C. Eder, S. W. Haan, B. A. Hammel, D. E. Hinkel, O. S. Jones, M. M. Marinak et al. , Phys. of Plasmas 22, 032702 (2015)
2015
-
[46]
Gatu Johnson, J
M. Gatu Johnson, J. P. Knauer, C. J. Cerjan, M. J. Eckart, G. P. Grim, E. P. Hartouni, R. Hatarik, J. D. Kilkenny, D. H. Munro et al. , Phys. Rev. E 94, 021202(R) (2016)
2016
-
[47]
K. M. Woo, R. Betti, D. Shvarts, A. Bose, D. Patel, R. Yan, P. Y. Chang, O. M. Mannion, R. Epstein et al. , Physics of Plasmas 25, 052704 (2018)
2018
-
[48]
Murphy, Phys
T. Murphy, Phys. of Plasmas 21, 072701 (2014)
2014
-
[49]
Brysk, Plasma Physics 15, 611 (1973)
H. Brysk, Plasma Physics 15, 611 (1973)
1973
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.