Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical Cherenkov radiation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Frequency-domain calculation of Smith–Purcell radiation for metallic and dielectric gratings

The intensity of Smith–Purcell radiation from metallic and dielectric gratings (silicon, silica) is compared in a frequency-domain simulation. Additionally, the numerical model is discussed and verified with the Frank–Tamm formula for Cherenkov radiation. For 30 keV electrons, rectangular dielectric gratings are less efficient than their metallic counterparts, by an order of magnitude for silicon, and two orders of magnitude for silica. For all gratings studied, radiation intensity oscillates with grating tooth height due to electromagnetic resonances in the grating. 3D and 2D numerical models are compared.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spontaneous emission near the electron plasma frequency in a plasma with a runaway electron tail

Spontaneous emission of radiation with frequencies near the electron plasma frequency is studied for a plasma which consists of both thermal and runaway electrons. It is found that a substantial enhancement of the spontaneous radiation intensity can occur in this frequency regime via a Cherenkov resonance with the runaway electrons. Numerical analysis indicates that, for reasonable estimates of densities and energies, the plasma-frequency radiation can attain levels greater than the peak thermal emission at the second gyroharmonic.

Freund, H. P.↗

Cherenkov and Transition Radiation as Low-Energy Background Sources in SuperCDMS Detectors

The Weakly Interacting Massive Particle (WIMP) has historically been a prime candidate for dark matter due to its elegant compatibility with the Minimally Supersymmetric Standard Model (MSSM). Recent dark matter experiments have ruled out much of the GeV-TeV mass range predicted by the MSSM, however, and the newest generation of direct detection experiments, such as SuperCDMS SNOLAB, have begunto explore sub-GeV dark matter. As experiments push towards lower masses, sensitivity to eV-scale electron recoil events has become increasingly important. A variety of unexplained excesses at energy deposits of 1-100 eV have been found inmany such low-mass experiments, across different detection techniques and at different excess rates. These low-energy events are not thought to be dark matter, and they must be understood and effectively removed to further improve detector energy resolution. This thesis will outline the simulation of one possible source of low-energy excess at SuperCDMS: optical photons produced from charges in uniform motion, specifically from Cherenkov radiation (CR), transition radiation (TR), and the intermediate hybrid transition-Cherenkov radiation (HR). In the latter case, the standard formulae for HR contain divergences at certain angles in transparent media, which are an artifact of their derivation and become problematic to simulate. To rsolve this, we present a novel approach to normalize the divergent HR peaks, which allows HR to transition smoothly between TR and CR in a numerical simulation. We also add TR and CR as a physics process to the SuperCDMS Monte Carlo simulation package SuperSim, based on the Geant4 simulation toolkit. To verify the physics of our addition, we show some test simulations in comparison to theoretical predictions of optical radiation intensity. We also use the simulation to make some preliminary predictions on the contribution of TR and CR to the low energy background, with the expectation that more thorough analyses will be conducted in the future.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A potential probe of Fundamental Interactions using Multi-Modal Cosmic Ray events

The Horizon-T experiment in Tien Shan is based on the idea of measuring the time at which EAS disc passes the observation level with nanosecond accuracy. The detector system consists of ten charged particles registration points located at distances of up to several hundred meters from each other. The points are equipped with detectors based on registration of Cherenkov radiation in glass and registration of scintillation light in polystyrene. The detectors register the arrival times of charged particles at the observation level with a resolution of ~2 ns, as required to study the spatial and temporal characteristics of the EAS and the structure of the multi-modal events specifically. Over the period of the Horizon-T data taking since 2017 to present, a large number of multi-modal events were detected. The data has presented numerous challenges that show the direction towards the further development of the detector system and of the analysis methods and techniques that could be applied to these multi-modal events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Axions and SN 1987A: Axion trapping

If an axion of mass between about 10(exp -3) and 10 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become trapped and radiated from an axion sphere. The trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a diffusion approximation for axion-energy transport. The axion opacity due to inverse nucleon-nucleon, axion bremsstrahlung is computed; and then the numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande II (KII) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The estimate of the axion mass is confirmed above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration - the most sensitive barometer of axion cooling - it is concluded that for an axion mass greater than about 3 eV axion emission would not have had a significant effect on the neutrino bursts detected by KII and IMB. It is strongly suggested that an axion with mass in the interval 10(exp -3) to 3 eV is excluded by the observation of neutrinos from SN 1987A.

Burrows, Adam↗

Axions and SN 1987A: Axion trapping

If an axion of mass between about 10(exp -3) eV and 1 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass less than about 10(exp -2) eV, axions produced deep inside the neutron star simply stream out; in a previous paper this case has been addressed. Remarkably, for an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become 'trapped' and radiated from an axion sphere. In this paper the trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a leakage approximation scheme for axion-energy transport. The axion opacity is computed due to inverse nucleon-nucleon, axion bremsstrahlung, and numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande 2 (K2) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The earlier estimate is confirmed and refined of the axion mass above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration--the most sensitive barometer of axion cooling--it is concluded that for an axion mass of greater than about 0.3 eV, axion emission would not have had a significant effect on the neutrino bursts detected by K2 and IMB. The present work, together with the previous work, strongly suggests that an axion with mass in the interval 10(exp -3) eV to 0.3 eV is excluded by SN 1987A.

Burrows, Adam↗

Characterization of the miniPlanacon XPM85112-S-R2D2 MCP-PMT with custom modified backend electronics

We report the results of the measurements of three pieces of the new Photonis miniPLANACON microchannel-plate photomultipliers (MCP-PMTs) intended for use in the demanding environment of the Large Hadron Collider (LHC) beamline as a part of the AFP Time-of-Flight detector. These photomultipliers were modified in cooperation with the manufacturer by using a custom backend and were subjected to numerous tests, with the focus on the rate capability and crosstalk behaviour. We determined that the two of them with a lower MCP resistance are able to operate without significant saturation at an anode current density of 1 μA/cm 2 . These two are, therefore, suitable for the intended use and are currently installed as part of the AFP detector packages.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗