In modeling the charged alpha particle transport in hot-spot plasmas of inertial confinement fusion, the energy-losing rate is a major concern in the Monte Carlo simulations of alpha particle transport of a radiative-hydrodynamic code. However, the traditionally used energy stopping-power only describes the averaged energy-losing rate of the incident charged particles, whereas the variance of the energy exchange with the background particles is generally ignored. In this paper, the variance of charged particle collisions is studied by both analytical derivation and Monte Carlo simulations. An expression of the divergence of the charged particle energy-losing rate is given for the first time, which can be directly used for practical estimations. It indicates that when the areal density of the target particles along the incident particle path length is low, the divergence of the lost energy would be much larger than the average value, and the traditionally used energy stopping-power would be no longer sufficient to describe the charged particle Coulomb collisions. It helps to obtain a more comprehensive understanding about the charged particle transport in plasmas.
Proton radiography is a widely used experimental method to diagnose the electric and magnetic (EM) fields in high-energy-density plasmas. In proton radiography, the probe protons are typically assumed to be deflected only by the EM fields, whereas the Coulomb scattering caused by the charged particles in the target plasmas is generally ignored. However, at high plasma densities, the presence of Coulomb scattering could reduce the proton flux perturbations recorded on the detector and influence the inversion of the EM fields from experiments. In this paper, a theoretical model is developed for the first time to describe the proton flux distribution on the detector when the EM field deflections and Coulomb scattering coexist in deflecting the probe proton trajectories. Our theory indicates that the Coulomb scattering could decrease the signal contrast of the probed EM fields, which is determined not only by the strengths of the EM field deflections and Coulomb scattering but also by the spatial gradient of the EM fields. Monte Carlo simulations are also conducted to validate our theoretical model. It would be helpful to interpret the proton radiography experiments quantitatively.
In the alpha particle transport in ICF hotspot, previous models focus mainly on how the incident particles lose their energy but lost sight of how the target particles will respond to this lost energy. In this paper, we developed a novel single-scattering model based on the Monte Carlo method, which abandons the stopping-power and models every single-scattering event in the alpha particle life. It enables to describe both the energy stopping of the incident alpha particle and the target particles response to the collisions. With this model, it shows that the target DT-ions at the ICF hotspot boundary will be non-Maxwellian distributed after colliding with the high-energy alpha particles, which refers to a much higher fusion reactivity compared with a Maxwellian one. At the same time, this model gives a longer and dispersed alpha particle range in hotspot plasmas and suggests that the traditionally used stopping power models would overestimate the stopping ability of the target particles.
When simulating a charged particle trajectory in plasmas with Monte Carlo methods, the traditional single-scattering model is accurate but computationally expensive. The traditional multiple-scattering model reduces the computation cost but assumes an exponential distribution of the scattering angle, which overestimates the cumulative scattering angle. In this paper, we develop a new Monte Carlo method for the charged particle transport simulations. The scattering angle is naturally sampled from a modified screened Rutherford scattering cross section, in which a greatly amplified minimum deflection angle is used. With this model, the simulations give the same scattering angle distribution with the single-scattering model but only take a comparable computation cost with the multiple-scattering model when the small scattering angle and Markovian approximations can be fulfilled. This model would be especially useful for the study of inertial-confinement-fusion, which is sensitive to the alpha-particle energy deposition.
The magnetic fields generated in plasmas have extensive influences on many processes of the inertial confinement fusion and the astrophysics. Therefore, the quantitative diagnosis of the magnetic field is quite essential. Proton radiography is a widely used experimental technique to diagnose the electric field or magnetic field in high-energy-density plasma. The effective explanation of the results of proton radiography depends on the reliability and availability of the inversion method. Traditional inversion methods can only provide one- or two-dimensional structure of the self-generated magnetic field. In this study, it is found that there is an Abel transformation relationship between the deflection velocity and the magnetic field with column symmetry, which allows us to reconstruct the three-dimensional structure of the magnetic field for the first time. We theoretically deduce the process of reconstructing the cylindrical magnetic field through proton radiography with the Abel inversion algorithm. The feasibility of this method is verified by numerical simulation as well. Based on this inversion method, we reanalyze the proton radiography experimental results of Li et al. (2016 Nat. Commun. 7 13081) on the self-generated magnetic field of plasma jets. The inversion results show that the maximum magnetic field intensity is about 1.9 times the traditional inversion results. We discuss a new proton radiography inversion method for the existence of magnetic fields with cylindrical symmetry in thiswork, which will contributes to an intensive understanding of the self-generated electromagnetic field and its spatiotemporal evolution related to the laser fusion and the laboratory astrophysics.
In proton radiography, degeneracy of electric and magnetic fields in deflecting the probe protons can prevent full interpretation of proton flux perturbations in the detection plane. In this paper, theoretical analyses and numerical simulations suggest that the contributions of the electric and magnetic fields can be separately obtained by analyzing the difference between the flux distributions of two discriminated proton energies in a single shot of proton radiography. To eliminate the influence of field evolution on the separation, a strategy is proposed in which slow field evolution is assumed or an approximate estimate of field growth is made. This could help achieve a clearer understanding of the radiographic process and allow further quantitative analysis.