Comparison of unfolding methods for the inference of runaway electron energy distribution from γ-ray spectroscopic measurements

2021 ◽  
Vol 16 (12) ◽  
pp. C12005
Author(s):  
E. Panontin ◽  
A. Dal Molin ◽  
M. Nocente ◽  
G. Croci ◽  
J. Eriksson ◽  
...  
Keyword(s):  
Iterative Method ◽  
Energy Distribution ◽  
Expectation Maximization ◽  
Maximum Energy ◽  
Runaway Electrons ◽  
X Ray ◽  
Maximum Likelihood Expectation Maximization ◽  
The Mean ◽  
Value Decomposition ◽  
Robust To Noise

Abstract Unfolding techniques are employed to reconstruct the 1D energy distribution of runaway electrons from Bremsstrahlung hard X-ray spectrum emitted during plasma disruptions in tokamaks. Here we compare four inversion methods: truncated singular value decomposition, which is a linear algebra technique, maximum likelihood expectation maximization, which is an iterative method, and Tikhonov regularization applied to χ 2 and Poisson statistics, which are two minimization approaches. The reconstruction fidelity and the capability of estimating cumulative statistics, such as the mean and maximum energy, have been assessed on both synthetic and experimental spectra. The effect of measurements limitations, such as the low energy cut and few number of counts, on the final reconstruction has also been studied. We find that the iterative method performs best as it better describes the statistics of the experimental data and is more robust to noise in the recorded spectrum.

scholarly journals Monte Carlo Simulation for Polychromatic X-Ray Fluorescence Computed Tomography with Sheet-Beam Geometry

2017 ◽  
Vol 2017 ◽  
pp. 1-10 ◽  
Author(s):  
Shanghai Jiang ◽  
Peng He ◽  
Luzhen Deng ◽  
Mianyi Chen ◽  
Biao Wei
Keyword(s):  
Computed Tomography ◽  
Monte Carlo Simulation ◽  
Monte Carlo ◽  
Expectation Maximization ◽  
Huge Amount ◽  
X Ray ◽  
Beam Geometry ◽  
Weight Concentration ◽  
Imaging Model ◽  
Maximum Likelihood Expectation Maximization

X-ray fluorescence computed tomography (XFCT) based on sheet beam can save a huge amount of time to obtain a whole set of projections using synchrotron. However, it is clearly unpractical for most biomedical research laboratories. In this paper, polychromatic X-ray fluorescence computed tomography with sheet-beam geometry is tested by Monte Carlo simulation. First, two phantoms (A and B) filled with PMMA are used to simulate imaging process through GEANT 4. Phantom A contains several GNP-loaded regions with the same size (10 mm) in height and diameter but different Au weight concentration ranging from 0.3% to 1.8%. Phantom B contains twelve GNP-loaded regions with the same Au weight concentration (1.6%) but different diameter ranging from 1 mm to 9 mm. Second, discretized presentation of imaging model is established to reconstruct more accurate XFCT images. Third, XFCT images of phantoms A and B are reconstructed by filter back-projection (FBP) and maximum likelihood expectation maximization (MLEM) with and without correction, respectively. Contrast-to-noise ratio (CNR) is calculated to evaluate all the reconstructed images. Our results show that it is feasible for sheet-beam XFCT system based on polychromatic X-ray source and the discretized imaging model can be used to reconstruct more accurate images.

scholarly journals Singular value decomposition analysis of back projection operator of maximum likelihood expectation maximization PET image reconstruction

2018 ◽  
Vol 52 (3) ◽  
pp. 337-345
Author(s):  
Vencel Somai ◽  
David Legrady ◽  
Gabor Tolnai
Keyword(s):  
Monte Carlo ◽  
Maximum Likelihood ◽  
Singular Value Decomposition ◽  
Expectation Maximization ◽  
Physical Modelling ◽  
Singular Value ◽  
Processing Unit ◽  
Back Projection ◽  
Maximum Likelihood Expectation Maximization ◽  
Value Decomposition

Abstract Background In emission tomography maximum likelihood expectation maximization reconstruction technique has replaced the analytical approaches in several applications. The most important drawback of this iterative method is its linear rate of convergence and the corresponding computational burden. Therefore, simplifications are usually required in the Monte Carlo simulation of the back projection step. In order to overcome these problems, a reconstruction code has been developed with graphical processing unit based Monte Carlo engine which enabled full physical modelling in the back projection. Materials and methods Code performance was evaluated with simulations on two geometries. One is a sophisticated scanner geometry which consists of a dodecagon with inscribed circle radius of 8.7 cm, packed on each side with an array of 39 × 81 LYSO detector pixels of 1.17 mm sided squares, similar to a Mediso nanoScan PET/CT scanner. The other, simplified geometry contains a 38,4mm long interval as a voxel space, detector pixels are assigned in two parallel sections each containing 81 crystals of a size 1.17×1.17 mm. Results We have demonstrated that full Monte Carlo modelling in the back projection step leads to material dependent inhomogeneities in the reconstructed image. The reasons behind this apparently anomalous behaviour was analysed in the simplified system by means of singular value decomposition and explained by different speed of convergence. Conclusions To still take advantage of the higher noise stability of the full physical modelling, a new filtering technique is proposed for convergence acceleration. Some theoretical considerations for the practical implementation and for further development are also presented.

scholarly journals Multi GPU parallelization of maximum likelihood expectation maximization method for digital rock tomography data

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Jaya Prakash ◽  
Umang Agarwal ◽  
Phaneendra K. Yalavarthy
Keyword(s):  
Maximum Likelihood ◽  
Expectation Maximization ◽  
Graphics Processing Units ◽  
Rock Physics ◽  
Micro Computed Tomography ◽  
X Ray ◽  
Iterative Schemes ◽  
Digital Rock ◽  
Maximum Likelihood Expectation Maximization ◽  
Fdk Algorithm

AbstractDigital rock is an emerging area of rock physics, which involves scanning reservoir rocks using X-ray micro computed tomography (XCT) scanners and using it for various petrophysical computations and evaluations. The acquired micro CT projections are used to reconstruct the X-ray attenuation maps of the rock. The image reconstruction problem can be solved by utilization of analytical (such as Feldkamp–Davis–Kress (FDK) algorithm) or iterative methods. Analytical schemes are typically computationally more efficient and hence preferred for large datasets such as digital rocks. Iterative schemes like maximum likelihood expectation maximization (MLEM) are known to generate accurate image representation over analytical scheme in limited data (and/or noisy) situations, however iterative schemes are computationally expensive. In this work, we have parallelized the forward and inverse operators used in the MLEM algorithm on multiple graphics processing units (multi-GPU) platforms. The multi-GPU implementation involves dividing the rock volumes and detector geometry into smaller modules (along with overlap regions). Each of the module was passed onto different GPU to enable computation of forward and inverse operations. We observed an acceleration of $$\sim 30$$ ∼ 30 times using our multi-GPU approach compared to the multi-core CPU implementation. Further multi-GPU based MLEM obtained superior reconstruction compared to traditional FDK algorithm.

scholarly journals The spectral energy distribution of an X-ray pair halo from a gamma-ray source with a power-law continuum

2021 ◽  
Vol 2145 (1) ◽  
pp. 012013
Author(s):  
A Eungwanichayapant ◽  
W Luangtip
Keyword(s):  
Energy Distribution ◽  
Power Law ◽  
Gamma Ray ◽  
Energy Levels ◽  
Spectral Energy Distribution ◽  
High Energy ◽  
Maximum Energy ◽  
Spectral Energy ◽  
Power Indices ◽  
X Ray

Abstract Interactions between Very High Energy (VHE) gamma-rays from Active Galactic Nuclei (AGNs) and infrared photons from the Extragalactic Background Light (EBL) can start electromagnetic cascades. If the extragalactic magnetic field near a host galaxy is strong enough (∼1 µG), the cascades would develop isotropically around the AGN. As a result, the electron/positron pairs created along the development of the cascades would create an X-ray halo via synchrotron radiation process. It is believed that the VHE gamma-ray spectra from the AGNs could be approximated by a power-law model which is truncated at high energy end (i.e. maximum energy). In this work we studied the X-ray Spectral Energy Distribution (SED) of the halo generated from the AGN spectra with different power indices and maximum energy levels. The results showed that the SEDs were slightly higher and broader, as they were obtaining higher flux if the power indices were lower. On the other hand, the SEDs were sensitive to the maximum energy levels between 100-300 TeV. More flux could be obtained from the higher maximum energy. However, we found that the SED becomes insensitive to the varied parameters when the maximum energy and the power index are > 500 TeV and < 1.5, respectively.

Isomorphous replacement in electron diffraction of wet crystals of rat hemoglobin

1983 ◽  
Vol 41 ◽  
pp. 446-447
Author(s):  
William F. Tivol ◽  
Murray Vernon King ◽  
D. F. Parsons
Keyword(s):  
Electron Diffraction ◽  
Heavy Atom ◽  
Isomorphous Substitution ◽  
X Ray Diffraction ◽  
Salt Solutions ◽  
X Ray ◽  
Isomorphous Replacement ◽  
Structure Factors ◽  
Diffraction Structure ◽  
The Mean

Feasibility of isomorphous substitution in electron diffraction is supported by a calculation of the mean alteration of the electron-diffraction structure factors for hemoglobin crystals caused by substituting two mercury atoms per molecule, following Green, Ingram & Perutz, but with allowance for the proportionality of f to Z3/4 for electron diffraction. This yields a mean net change in F of 12.5%, as contrasted with 22.8% for x-ray diffraction.Use of the hydration chamber in electron diffraction opens prospects for examining many proteins that yield only very thin crystals not suitable for x-ray diffraction. Examination in the wet state avoids treatments that could cause translocation of the heavy-atom labels or distortion of the crystal. Combined with low-fluence techniques, it enables study of the protein in a state as close to native as possible.We have undertaken a study of crystals of rat hemoglobin by electron diffraction in the wet state. Rat hemoglobin offers a certain advantage for hydration-chamber work over other hemoglobins in that it can be crystallized from distilled water instead of salt solutions.

Long-period modulated superstructures in Pd-12.5 at.%Ce and Pd-15 at.%Ce alloys

1990 ◽  
Vol 48 (4) ◽  
pp. 164-165
Author(s):  
Noriyuki Kuwano ◽  
Masaru Itakura ◽  
Kensuke Oki
Keyword(s):  
Electron Microscopy ◽  
Mean Value ◽  
Heavy Elements ◽  
Arc Furnace ◽  
X Ray ◽  
Long Period ◽  
Ultra High Purity ◽  
Heat Treated ◽  
Atomic Scattering ◽  
The Mean

Pd-Ce alloys exhibit various anomalies in physical properties due to mixed valences of Ce, and the anomalies are thought to be strongly related with the crystal structures. Since Pd and Ce are both heavy elements, relative magnitudes of (fcc-fpd) are so small compared with <f> that superlattice reflections, even if any, sometimes cannot be detected in conventional x-ray powder patterns, where fee and fpd are atomic scattering factors of Ce and Pd, and <f> the mean value in the crystal. However, superlattices in Pd-Ce alloys can be analyzed by electron microscopy, thanks to the high detectability of electron diffraction. In this work, we investigated modulated superstructures in alloys with 12.5 and 15.0 at.%Ce.Ingots of Pd-Ce alloys were prepared in an arc furnace under atmosphere of ultra high purity argon. The disc specimens cut out from the ingots were heat-treated in vacuum and electrothinned to electron transparency by a jet method.

Localization of calcium and zinc in the sperm storage tubules of chicken, quail and turkey using X-ray microanalysis

Reproduction ◽  
2000 ◽  
pp. 331-336 ◽  
Author(s):  
L Holm ◽  
H Ekwall ◽  
GJ Wishart ◽  
Y Ridderstrale
Keyword(s):  
Calcium Concentration ◽  
Sperm Storage ◽  
X Ray ◽  
Elemental Concentrations ◽  
Low Concentrations ◽  
Intracellular Content ◽  
Wet Weight ◽  
The Mean ◽  
Sperm Storage Tubules ◽  
Scanning Electron

Sperm storage tubules from the utero-vaginal junction of chickens, quails and turkeys were analysed for calcium and zinc using X-ray microanalysis of ultra-rapidly frozen tissue in a scanning electron microscope. This technique enabled the tubular fluid surrounding the stored spermatozoa and the intracellular content of the cells of the sperm storage tubules to be analysed separately and, by using standards with known concentrations, their elemental concentrations were estimated. The mean (+/- SEM) concentration of calcium in the tubular fluid from chickens, quails and turkeys was 17 +/- 3, 19 +/- 3 and 17 +/- 4 mmol kg(-1) wet weight, respectively. The intracellular calcium concentration of the cells of the tubules did not differ significantly from these values and was also similar in the mucosal epithelial cells of the utero-vaginal junction. Zinc was localized in the cells of turkey sperm storage tubules and tubular fluid, but at low concentrations. No zinc could be detected in corresponding structures from chickens and quails. The concentration of calcium in the tubular fluid is within the range known to inhibit the motility of spermatozoa, supporting this function for calcium during storage. Zinc is known to depress turkey sperm metabolism and it may also be involved in inducing quiescence of spermatozoa during storage in this species.

Use of the electro-separation method for improvement of the utility value of winter rapeseeds

2013 ◽  
Vol 27 (4) ◽  
pp. 419-424 ◽  
Author(s):  
S.J. Kovalyshyn ◽  
O.P. Shvets ◽  
S. Grundas ◽  
J. Tys
Keyword(s):  
Corona Discharge ◽  
Mechanical Damage ◽  
Edible Oil ◽  
Separation Method ◽  
Utility Value ◽  
X Ray ◽  
Winter Rapeseed ◽  
The Mean ◽  
Discharge Field ◽  
Prototype Apparatus

Abstract The paper presents the results of a study of the use of electro-separation methods for improvement of the utility value of 5 winter rapeseed cultivars. The process of electro-separation of rapeseed was conducted on a prototype apparatus built at the Laboratory of Application of Electro-technologies in Agriculture, Lviv National Agriculture University. The process facilitated separation of damaged, low quality seeds from the sowing material. The initial mean level of mechanically damaged seeds in the winter rapeseed cultivars studied varied within the range of 15.8-20.1%. Verification of the amount of seeds with mechanical damage was performed on X-ray images of seeds acquired by means of a digital X-ray apparatus. In the course of analysis of the X-ray images, it was noted that the mean level of mechanical damage to the seeds after the electro-separation was in the range of 2.1-3.8%. The application of the method of separation of rapeseeds in the corona discharge field yielded a significant reduction of the level of seeds with mechanical damage. The application of the method in practice may effectively contribute to improvement of the utility value of sowing material or seed material for production of edible oil.

Sign in / Sign up

Export Citation Format

Share Document