Monoenergetic Particles Research Articles

Multiwavelength Galactic center gamma-ray observations explained by a unified cosmic-ray dynamics model

Context. High-energy (HE) and very high-energy (VHE) gamma-ray observations from the Galactic center (GC) detected extended emission correlated with the morphology of the central molecular zone (CMZ). Emission in both bands is expected to be produced by hadronic interaction between cosmic rays (CRs) and ambient gas. Aims. We examine if our three previously proposed scenarios for the CR sources and dynamics, which are consistent with the VHE gamma-ray data (1–100 TeV), also match the HE gamma-ray observations (10–300 GeV). Additionally, we analyze the effect of the isotropic Galactic CR “sea” inside the CMZ. Methods. We generated synthetic gamma-ray maps considering a simplified isotropic diffusion, but more realistic dynamics with two diffusion zones (in and out of the CMZ) and polar advection, for mono-energetic particles of 3 TeV. Additionally, we considered two gas distributions for the CMZ (with and without an inner cavity), and CR populations injected from the clusters of young massive stars (the Arches Cluster, the Quintuplet Cluster, and the nuclear star cluster), plus the supernova Sgr A East. Results. Only the combination of more realistic CR dynamics, the CMZ with an inner cavity, CR injection from all proposed sources, and a CR sea similar to that observed in the Solar System reproduced the current HE and VHE gamma-ray detection from the CMZ and was consistent with the observed gamma-rays from Sagittarius A*. Conclusions. The HE and VHE gamma-rays observations of the GC can be reproduced by a unified model for the CRs.

Laser Proton Acceleration from a Near-Critical Imploding Gas Target.

The interaction between relativistic intense laser pulses and near-critical-density targets has been sought after in order to increase the efficiency of laser-plasma energy coupling, particularly for laser-driven proton acceleration. To achieve the density regime for high-repetition-rate applications, one elusive approach is to use gas targets, provided that stringent target density profile requirements are met. These include reaching the critical plasma density while maintaining micron-scale density gradients. In this Letter, we present a novel scheme for achieving the necessary requirements using optical laser pulses to transversely shape the target and create a colliding shock wave in both planar and cylindrical geometries. Utilizing this approach, we experimentally demonstrated stable proton acceleration and achieved up to 5MeV in a monoenergetic distribution and particle numbers above 10^{8} Sr^{-1} MeV^{-1} using a 1.5J energy on-target laser pulse. The Letter also reports for the first time an extend series of 200 consecutive shots that demonstrates the robustness of the approach and its maturity for applications. These results open the door for future work in controlling gas targets and optimizing the acceleration process for more energetic multipetawatt laser systems.

Energy-dependent implementation of secondary electron emission models in continuum kinetic sheath simulations

Abstract The plasma-material interactions present in multiple fusion and propulsion concepts between the flow of plasma through a channel and a material wall drive the emission of secondary electrons. This emission is capable of altering the fundamental structure of the sheath region, significantly changing the expected particle fluxes to the wall. The emission spectrum is separated into two major energy regimes, a peak of elastically backscattered primary electrons at the incoming energy, and cold secondary electrons inelastically emitted directly from the material. The ability of continuum kinetic simulations to accurately represent the secondary electron emission is limited by relevant models being formulated in terms of monoenergetic particle interactions which cannot be applied directly to the discrete distribution function. As a result, rigorous implementation of energy-dependent physics is often neglected in favor of simplified, constant models. We present here a novel implementation of semi-empirical models in the boundary of continuum kinetic simulations which allows the full range of this emission to be accurately captured in physically-relevant regimes.

Modeling the Galactic center gamma-ray emission with more realistic cosmic-ray dynamics

Context. Very-high-energy gamma-ray observations of the Galactic center (GC) show extended emission that is strongly correlated with the morphology of the central molecular zone (CMZ). The best explanation for that emission is a hadronic interaction between cosmic rays (CRs) and ambient gas, where a CR central and continuous source accelerates protons up to 1 PeV (“PeVatron”). However, current models assume very simplistic CR dynamics. Aims. Our goal is to verify if more realistic CR dynamics for the GC environment are consistent with current gamma-ray observations, and whether they could be constrained by upcoming observations with the Cherenkov Telescope Array (CTA). Methods. We generated synthetic gamma-ray maps using a CR transport model with spherical injection, different diffusion regimes (in and out of the CMZ), polar advection, and mono-energetic particles of 1 PeV, and including different CR populations injected from the Arches, Quintuplet, and nuclear clusters of young massive stars, plus supernova Sgr A East. We adopted two different 3D gas distributions consistent with the observed gas column density, either with or without an inner cavity. Results. In order to reproduce the existing observations detected by the High Energy Stereoscopic System (HESS), a ring-like gas distribution, with its mass set by the standard Galactic CO-to-H2 conversion factor, and CR acceleration from all relevant sources are required. For a conversion factor one order of magnitude lower, injection rates that are ten times higher are needed. We show that CTA will be able to differentiate between models with different CR dynamics, proton sources, and CMZ morphologies, owing to its unprecedented sensitivity and angular resolution. Conclusions. More realistic CR dynamics suggest that the CMZ has a large inner cavity and that the GC PeVatron is a composite CR population accelerated by the Arches, Quintuplet, and nuclear star clusters, and Sgr A East.

A personalized Monte Carlo study of tumor and critical organ doses for Trans-arterial Radioembolization patients.

Trans-arterial Radioembolization (TARE) is an intra-arterial treatment method for liver malignancies. In this procedure, the therapeutic tumor dose is significant for predicting the treatment effectiveness while the dose absorbed in an organ at risk provides an understanding of its tolerance to radiation. This study proposes a Monte Carlo (MC) approach for determining absorbed organ doses for patients undergoing TARE treatment. The technique is based on the use of a voxel-based partial body model generated for each patient from his/her anatomical image data to represent the critical body structures more realistically. These structures are first segmented from image slices to create an image block which is then incorporated into a radiation transport package (MCNP6.2) to perform MC simulations. When used along with the parameters specific to a patient's treatment, such as lung-shunt factor (LSF), tumor-to-normal liver ratio (TLR), fractional uptakes, and administered activity, this approach allowed more accurate simulation of radiation interactions and hence provided absorbed doses specific to a TARE patient. The MC method also calculated the absorbed doses in organs or tissues that were close to target tissues for which the MIRD formalism makes no predictions. MIRD calculations were found to overestimate the absorbed doses by as much as 11% in lungs, 5% in liver, and 20% in tumor volumes. This raises concerns about the treatment's efficacy when estimating the correct activity to be administered to a patient. When each patient simulation was repeated with a 90Y source spectrum to reflect the distribution of varying beta energies, the liver and the lungs were observed to receive relatively lower doses than those obtained with monoenergetic beta particles. Thus, it can be stated that the approach adopted in this study offers a more precise model of the patient's critical tissues and serves as a personalized dosimetric tool for TARE treatment planning.

Energetic particle dynamics in a simplified model of a solar wind magnetic switchback

Context.Recent spacecraft observations in the inner heliosphere have revealed the presence of local Alfvénic reversals of the magnetic field, while the field magnitude remains almost constant. These are called magnetic switchbacks (SBs) and are very common in the plasma environment close to the Sun explored by the Parker Solar Probe satellite.Aims.A simple numerical model of a magnetic field reversal with constant magnitude is used in order to explore the influence of SBs on the propagation of energetic particles within a range of energy typical of solar energetic particles.Methods.We model the reversal as a region of space of adjustable size bounded by two rotational discontinuities. By means of test particle simulations, beams of mono-energetic particles can be injected upstream of the SB with various initial pitch- and gyro-phase angles. In each simulation, the particle energy may also be changed.Results.Particle dynamics is highly affected by the ratio between the particle gyroradius and the size of the SB, with multiple pitch-angle scatterings occurring when the particle gyroradius is of the order of the SB size. Further, particle motion is extremely sensitive to the initial conditions, implying a transition to chaos; for some parameters of the system, a large share of particles is reflected backwards upstream as they interact with the SB. These results could have a profound impact on our understanding of solar energetic particle transport in the inner heliosphere, and therefore possible comparisons with in situ spacecraft data are discussed.

The complexity of DNA damage by radiation follows a Gamma distribution: insights from the Microdosimetric Gamma Model.

DNA damage is the main predictor of response to radiation therapy for cancer. Its Q8 quantification and characterization are paramount for treatment optimization, particularly in advanced modalities such as proton and alpha-targeted therapy. We present a novel approach called the Microdosimetric Gamma Model (MGM) to address this important issue. The MGM uses the theory of microdosimetry, specifically the mean energy imparted to small sites, as a predictor of DNA damage properties. MGM provides the number of DNA damage sites and their complexity, which were determined using Monte Carlo simulations with the TOPAS-nBio toolkit for monoenergetic protons and alpha particles. Complexity was used together with a illustrative and simplistic repair model to depict the differences between high and low LET radiations. DNA damage complexity distributions were were found to follow a Gamma distribution for all monoenergetic particles studied. The MGM functions allowed to predict number of DNA damage sites and their complexity for particles not simulated with microdosimetric measurements (yF) in the range of those studied. Compared to current methods, MGM allows for the characterization of DNA damage induced by beams composed of multi-energy components distributed over any time configuration and spatial distribution. The output can be plugged into ad hoc repair models that can predict cell killing, protein recruitment at repair sites, chromosome aberrations, and other biological effects, as opposed to current models solely focusing on cell survival. These features are particularly important in targeted alpha-therapy, for which biological effects remain largely uncertain. The MGM provides a flexible framework to study the energy, time, and spatial aspects of ionizing radiation and offers an excellent tool for studying and optimizing the biological effects of these radiotherapy modalities.

Intra-brain vascular models within the ICRP mesh-type adult reference phantoms for applications to internal dosimetry

Objective. Phantoms of the International Commission on Radiological Protection provide a framework for standardized dosimetry. The modeling of internal blood vessels—essential to tracking circulating blood cells exposed during external beam radiotherapy and to account for radiopharmaceutical decays while still in blood circulation—is, however, limited to the major inter-organ arteries and veins. Intra-organ blood is accounted for only through the assignment of a homogeneous mixture of parenchyma and blood [single-region (SR) organs]. Our goal was to develop explicit dual-region (DR) models of intra-organ blood vasculature of the adult male brain (AMB) and adult female brain (AFB). Approach. A total of 4000 vessels were created amongst 26 vascular trees. The AMB and AFB models were then tetrahedralized for coupling to the PHITS radiation transport code. Absorbed fractions were computed for monoenergetic alpha particles, electrons, positrons, and photons for both decay sites within the blood vessels and for tissues outside these vessels. Radionuclide S-values were computed for 22 and 10 radionuclides commonly employed in radiopharmaceutical therapy and nuclear medicine diagnostic imaging, respectively. Main results. For radionuclide decays, values of S(brain tissue ← brain blood) assessed in the traditional manner (SR) were higher than those computed using our DR models by factors of 1.92, 1.49, and 1.57 for therapeutic alpha-emitters, beta-emitters, and Auger electron-emitters, respectively in the AFB and by factors of 1.65, 1.37, and 1.42 for these same radionuclide categories in the AMB. Corresponding ratios of SR and DR values of S(brain tissue ← brain blood) were 1.34 (AFB) and 1.26 (AMB) for four SPECT radionuclides, and were 1.32 (AFB) and 1.24 (AMB) for six common PET radionuclides. Significance. The methodology employed in this study can be explored in other organs of the body for proper accounting of blood self-dose for that fraction of the radiopharmaceutical still in general circulation.

Stopping Power of Matter for a Beam of Monoenergetic Alpha Particles

Stopping Power of Matter for a Beam of Monoenergetic Alpha Particles

Development of a voxel S-value database for patient internal radiation dosimetry

PurposePersonalized dosimetry with high accuracy drew great attention in clinical practices. Voxel S-value (VSV) convolution has been proposed to speed up absorbed dose calculations. However, the VSV method is efficient for personalized internal radiation dosimetry only when there are pre-calculated VSVs of the radioisotope. In this work, we propose a new method for VSV calculation based on the developed mono-energetic particle VSV database of γ, β, α, and X-ray for any radioisotopes. MethodsMono-energetic VSV database for γ, β, α, and X-ray was calculated using Monte Carlo methods. Radiation dose was first calculated based on mono-energetic VSVs for [F-18]-FDG in 10 patients. The estimated doses were compared with the values obtained from direct Monte Carlo simulation for validation of the proposed method. The number of VSVs used in calculation was optimized based on the estimated dose accuracy and computation time. ResultsThe generated VSVs showed a great consistency with the results calculated using direct Monte Carlo simulation. For [F-18]-FDG, the proposed VSV method with number of VSV of 9 shows the best relative average organ absorbed dose uncertainty of 3.25% while the calculation time was reduced by 99% and 97% compared to the Monte Carlo simulation and traditional multiple VSV methods, respectively. ConclusionsIn this work, we provided a method to generate the VSV kernels for any radioisotope based on the pre-calculated mono-energetic VSV database and significantly reduced the time cost for the multiple VSVs dosimetry approach. A software was developed to generate VSV kernels for any radioisotope in 19 mediums.

Stopping Power of Matter for a Beam of Monoenergetic Alpha Particles

The results of using the statistics of multiple scattering to describe the dependence of the stopping power S of matter on the energy E0 of an alpha particle beam are presented. It is shown that the application of a new technique based on taking into account the dependence of the charge state of the beam ions on the ratio of the ion velocity to the minimum velocity of the substance electrons makes it possible to calculate S adequately to the experimental results in a wide range of particle energies E0.

Measurement of the insensitive surface layer thickness of a PIN photodiode based on alpha-particle spectrometry

In this work, the insensitive layer thickness of a PIN photodiode (SFH206K - Osram) has been measured by varying the incident angle of a collimated monoenergetic alpha particle beam. This technique is based on the fact that, except for the intrinsic uncertainties of the emission and straggling of alpha particles, their trajectories at different incident angles are distinct, with smaller ones corresponding to the incidence perpendicular to the detector surface. The result obtained (711 ± 23) nm, less than 1% of the intrinsic layer thickness, besides validating the employed method, demonstrates that the diode herein investigated is suitable for high resolution charged particle spectrometry.

RADIATION FRICTION OF RELATIVISTIC CHARGED PARTICLES MOVING IN A PERIODIC FIELD

The motion of relativistic charged particles beam in an external periodic field is considered, taking into account the influence of incoherent fields produced by particles on this motion. On the basis of the dynamics of individual particles motion under the action of the pair interaction forces each of them we derived the coefficient of friction. The expression for the friction force, which describes the average change in the momentum of charged particles per unit time, in the case of motion of an initially monoenergetic particle beam, is obtained.

A mesh-based model of liver vasculature: implications for improved radiation dosimetry to liver parenchyma for radiopharmaceuticals

PurposeTo develop a model of the internal vasculature of the adult liver and demonstrate its application to the differentiation of radiopharmaceutical decay sites within liver parenchyma from those within organ blood.MethodComputer-generated models of hepatic arterial (HA), hepatic venous (HV), and hepatic portal venous (HPV) vascular trees were algorithmically created within individual lobes of the ICRP adult female and male livers (AFL/AML). For each iteration of the algorithm, pressure, blood flow, and vessel radii within each tree were updated as each new vessel was created and connected to a viable bifurcation site. The vascular networks created inside the AFL/AML were then tetrahedralized for coupling to the PHITS radiation transport code. Specific absorbed fractions (SAF) were computed for monoenergetic alpha particles, electrons, positrons, and photons. Dual-region liver models of the AFL/AML were proposed, and particle-specific SAF values were computed assuming radionuclide decays in blood within two locations: (1) sites within explicitly modeled hepatic vessels, and (2) sites within the hepatic blood pool residing outside these vessels to include the capillaries and blood sinuses. S values for 22 and 10 radionuclides commonly used in radiopharmaceutical therapy and imaging, respectively, were computed using the dual-region liver models and compared to those obtained in the existing single-region liver model.ResultsLiver models with virtual vasculatures of ~ 6000 non-intersecting straight cylinders representing the HA, HPV, and HV circulations were created for the ICRP reference. For alpha emitters and for beta and auger-electron emitters, S values using the single-region models were approximately 11% (AML) to 14% (AFL) and 11% (AML) to 13% (AFL) higher than the S values obtained using the dual-region models, respectively.ConclusionsThe methodology employed in this study has shown improvements in organ parenchymal dosimetry through explicit consideration of blood self-dose for alpha particles (all energies) and for electrons at energies below ~ 100 keV.

Preliminary benchmarks and analysis of boundary conditions in a trenched microstructured silicon radiation detector

Microstructured neutron detectors have the benefit of enhanced neutron detection efficiency as compared to planar devices, achieved by etching 6LiF-filled trenches on the top surface of a silicon PIN diode. This sensor geometry results in a complex electric field distribution and depletion characteristics within the diode under reverse bias. For the first time on record, the effects of a fixed oxide charge on the microstructured device depletion characteristics and mobile carrier transport is investigated. Prototype detectors were fabricated with non-conformal surface doping. Capacitance voltage and current voltage measurements were performed for these prototypes and compared with COMSOL Multiphysics simulations. A spectral response from an 241Am alpha particle source was acquired and analyzed. It was found that monoenergetic alpha particles produce three prominent peaks in the pulse height spectrum output by the device. The peaks were confirmed by simulations to correlate with dead layers and incident trajectories into the microstructure. It was also found that significant differences in pulse rise time result, corresponding with events arriving in a low-field region in the fins and a high-field region in the bulk. Geant4 was utilized for radiation transport, interaction modeling, and benchmarking the spectral data. The results of this simulation work provide confidence in the ability to attain and benchmark electrical characteristics and spectral data for semiconductor radiation detectors employing complex microstructures.

Galactic center gamma-ray production by cosmic rays from stellar winds and Sgr A East

Context. The High Energy Stereoscopic System, the Major Atmospheric Gamma-ray Imaging Cherenkov Telescope, and the Very Energetic Radiation Imaging Telescope Array System have observed diffuse gamma-ray emission strongly correlated with the central molecular zone in the Galactic center. The most accepted scenario to generate this emission is via a hadronic interaction between cosmic rays (CRs) and ambient gas, where CRs are accelerated from a central and continuous source of 1 PeV protons. Aims. We explore the influence of the three-dimensional (3D) shape of the central molecular zone on the indirect observation of the CR energy density via gamma-ray detection. Methods. We simulated synthetic gamma-ray maps using a CR diffusion model with spherical injection, one isotropic diffusion coefficient, no advection, and mono-energetic particles of 1 PeV. Also, we used two different 3D gas distributions considering the observed gas column density, both with and without an inner cavity. Results. We find that when using a persistent CR source, a disk-like gas distribution is needed to reproduce the existing CR indirect observations. This is in agreement with the continuous gas distribution implied by some dynamical models and studies based on the comparison of emission and absorption molecular lines. However, it contradicts several models of the central molecular zone, which imply that this structure has a significant inner cavity. This tension can be reconciled by an additional, impulsive CR injection. Conclusions. If the central molecular zone has a cavity, a composite CR population, coming from the stellar winds of the Wolf-Rayet stars in the central 0.5 pc and the supernova Sgr A East, produces a good match to the observed gamma-ray morphology in the Galactic center.

The relation between microdosimetry and induction of direct damage to DNA by alpha particles

In radiopharmaceutical treatments α-particles are employed to treat tumor cells. However, the mechanism that drives the biological effect induced is not well known. Being ionizing radiation, α-particles can affect biological organisms by producing damage to the DNA, either directly or indirectly. Following the principle that microdosimetry theory accounts for the stochastic way in which radiation deposits energy in sub-cellular sized volumes via physical collisions, we postulate that microdosimetry represents a reasonable framework to characterize the statistical nature of direct damage induction by α-particles to DNA. We used the TOPAS-nBio Monte Carlo package to simulate direct damage produced by monoenergetic alpha particles to different DNA structures. In separate simulations, we obtained the frequency-mean lineal energy () and dose-mean lineal energy () of microdosimetric distributions sampled with spherical sites of different sizes. The total number of DNA strand breaks, double strand breaks (DSBs) and complex strand breaks per track were quantified and presented as a function of either or The probability of interaction between a track and the DNA depends on how the base pairs are compacted. To characterize this variability on compactness, spherical sites of different size were used to match these probabilities of interaction, correlating the size-dependent specific energy () with the damage induced. The total number of DNA strand breaks per track was found to linearly correlate with and when using what we defined an effective volume as microdosimetric site, while the yield of DSB per unit dose linearly correlated with or being larger for compacted than for unfolded DNA structures. The yield of complex breaks per unit dose exhibited a quadratic behavior with respect to and a greater difference among DNA compactness levels. Microdosimetric quantities correlate with the direct damage imparted on DNA.

Polymer electrets and their applications

AbstractPolymer electrets have revealed great potential application in electromechanical devices because of the low weight, large quasi‐piezoelectric sensitivity, and excellent flexibility. For an electret, a permanent and macroscopic electric field exists on the surface, principally led by a macroscopic electrostatic charge on the surface or a net orientation of polar groups inside the object. Here, progress in the development of polymer electrets is reviewed. After a brief retrospect of the research courses and those typical polymer electrets that are classified into fluorine polymer and nonfluorine polymer, we present a survey on the charging methods, including corona, soft X‐ray, contact, thermal and monoenergetic particle beams. The latest representative applications (i.e., power harvesting, sensors, field effect transistors, and biomedicine) based on polymer electrets are also summarized. Finally, we complete this review with a discussion on perspectives and challenges in this field.

Magnetic confinement of effectively unmagnetized plasma particles

A purely magnetic applied field may provide plasma confinement under conditions where the bulk of the plasma is effectively free of the applied magnetic field. The applied magnetic field surrounds the bulk of the plasma, and plasma particles that are incident on the applied magnetic field can be reflected back into the effectively unmagnetized region of plasma. The concept belongs to a class of magnetic plasma confinement approaches studied long ago, for which some experimental results indicated that classical (collision-based) cross-magnetic-field transport may occur. However, multiple magnetic coils are required to be immersed within the confined plasma, and rapid plasma loss may occur if material structures are present, which pass through the plasma (e.g., to hold the immersed coils in place). In the work reported, the concept is studied in combination with magnetic plasma expulsion [R. E. Phillips and C. A. Ordonez, Phys. Plasmas 25, 012508 (2018)], which would be employed to keep plasma away from material structures that pass through the plasma. A planar model is used for the study. A classical trajectory Monte Carlo simulation is carried out on particles that are independently incident on the applied magnetic field. With monoenergetic incident particles, the results indicate that the applied magnetic field can reflect all independently incident particles in certain regions of parameter space. Prospects for achieving three-dimensional magnetic confinement of an effectively unmagnetized plasma with a Maxwellian velocity distribution are discussed.

Saturn-ring proton backlighters for the National Ignition Facility.

Proton radiography is a well-established technique for measuring electromagnetic fields in high-energy-density plasmas. Fusion reactions producing monoenergetic particles, such as D3He, are commonly used as a source, produced by a capsule implosion. Using smaller capsules for radiography applications is advantageous as the source size decreases, but on the National Ignition Facility (NIF), this can introduce complications from increasing blow-by light, since the phase plate focal spot size is much larger than the capsules. We report a demonstration of backlighter targets where a "Saturn" ring is placed around the capsule to block this light. The nuclear performance of the backlighters is unperturbed by the addition of a ring. We also test a ring with an equatorial cutout, which severely affects the proton emission and is not viable for radiography applications. These results demonstrate the general viability of Saturn ring backlighter targets for use on the NIF.