Uncertainties in polarimetric 3D reconstructions of coronal mass ejections

This work is aimed at quantifying the uncertainties in the 3D reconstruction of the location of coronal mass ejections (CMEs) obtained with the polarization ratio technique. The method takes advantage of the different distributions along the line of sight (LOS) of total (tB) and polarized (pB) brightnesses to estimate the average location of the emitting plasma. To this end, we assumed two simple electron density distributions along the LOS (a constant density and Gaussian density profiles) for a plasma blob and synthesized the expected tB and pB for different distances $z$ of the blob from the plane of the sky (POS) and different projected altitudes $\rho$. Reconstructed locations of the blob along the LOS were thus compared with the real ones, allowing a precise determination of uncertainties in the method. Independently of the analytical density profile, when the blob is centered at a small distance from the POS (i.e. for limb CMEs) the distance from the POS starts to be significantly overestimated. Polarization ratio technique provides the LOS position of the center of mass of what we call folded density distribution, given by reflecting and summing in front of the POS the fraction of density profile located behind that plane. On the other hand, when the blob is far from the POS, but with very small projected altitudes (i.e. for halo CMEs, $\rho < 1.4$ R$_\odot$), the inferred distance from that plane is significantly underestimated. Better determination of the real blob position along the LOS is given for intermediate locations, and in particular when the blob is centered at an angle of $20^\circ$ from the POS. These result have important consequences not only for future 3D reconstruction of CMEs with polarization ratio technique, but also for the design of future coronagraphs aimed at providing a continuous monitoring of halo-CMEs for space weather prediction purposes

High-performance blob-based iterative three-dimensional reconstruction in electron tomography using multi-GPUs

<p>Abstract</p> <p>Background</p> <p>Three-dimensional (3D) reconstruction in electron tomography (ET) has emerged as a leading technique to elucidate the molecular structures of complex biological specimens. Blob-based iterative methods are advantageous reconstruction methods for 3D reconstruction in ET, but demand huge computational costs. Multiple graphic processing units (multi-GPUs) offer an affordable platform to meet these demands. However, a synchronous communication scheme between multi-GPUs leads to idle GPU time, and a weighted matrix involved in iterative methods cannot be loaded into GPUs especially for large images due to the limited available memory of GPUs.</p> <p>Results</p> <p>In this paper we propose a multilevel parallel strategy combined with an asynchronous communication scheme and a blob-ELLR data structure to efficiently perform blob-based iterative reconstructions on multi-GPUs. The asynchronous communication scheme is used to minimize the idle GPU time so as to asynchronously overlap communications with computations. The blob-ELLR data structure only needs nearly 1/16 of the storage space in comparison with ELLPACK-R (ELLR) data structure and yields significant acceleration.</p> <p>Conclusions</p> <p>Experimental results indicate that the multilevel parallel scheme combined with the asynchronous communication scheme and the blob-ELLR data structure allows efficient implementations of 3D reconstruction in ET on multi-GPUs.</p

Demonstration of the event identification capabilities of the NEXT-White detector

In experiments searching for neutrinoless double-beta decay, the possibility of identifying the two emitted electrons is a powerful tool in rejecting background events and therefore improving the overall sensitivity of the experiment. In this paper we present the first measurement of the efficiency of a cut based on the different event signatures of double and single electron tracks, using the data of the NEXT-White detector, the first detector of the NEXT experiment operating underground. Using a 228Th calibration source to produce signal-like and background-like events with energies near 1.6 MeV, a signal efficiency of 71.6 ± 1.5 stat± 0.3 sys% for a background acceptance of 20.6 ± 0.4 stat± 0.3 sys% is found, in good agreement with Monte Carlo simulations. An extrapolation to the energy region of the neutrinoless double beta decay by means of Monte Carlo simulations is also carried out, and the results obtained show an improvement in background rejection over those obtained at lower energies. [Figure not available: see fulltext.

Hydrodynamic simulations with the Godunov SPH

We present results based on an implementation of the Godunov Smoothed Particle Hydrodynamics (GSPH), originally developed by Inutsuka (2002), in the GADGET-3 hydrodynamic code. We first review the derivation of the GSPH discretization of the equations of moment and energy conservation, starting from the convolution of these equations with the interpolating kernel. The two most important aspects of the numerical implementation of these equations are (a) the appearance of fluid velocity and pressure obtained from the solution of the Riemann problem between each pair of particles, and (b the absence of an artificial viscosity term. We carry out three different controlled hydrodynamical three-dimensional tests, namely the Sod shock tube, the development of Kelvin-Helmholtz instabilities in a shear flow test, and the "blob" test describing the evolution of a cold cloud moving against a hot wind. The results of our tests confirm and extend in a number of aspects those recently obtained by Cha (2010): (i) GSPH provides a much improved description of contact discontinuities, with respect to SPH, thus avoiding the appearance of spurious pressure forces; (ii) GSPH is able to follow the development of gas-dynamical instabilities, such as the Kevin--Helmholtz and the Rayleigh-Taylor ones; (iii) as a result, GSPH describes the development of curl structures in the shear-flow test and the dissolution of the cold cloud in the "blob" test. We also discuss in detail the effect on the performances of GSPH of changing different aspects of its implementation. The results of our tests demonstrate that GSPH is in fact a highly promising hydrodynamic scheme, also to be coupled to an N-body solver, for astrophysical and cosmological applications. [abridged]Comment: 19 pages, 13 figures, MNRAS accepted, high resolution version can be obtained at http://adlibitum.oats.inaf.it/borgani/html/papers/gsph_hydrosim.pd

High-ISO long-exposure image denoising based on quantitative blob characterization

Blob detection and image denoising are fundamental, sometimes related tasks in computer vision. In this paper, we present a computational method to quantitatively measure blob characteristics using normalized unilateral second-order Gaussian kernels. This method suppresses non-blob structures while yielding a quantitative measurement of the position, prominence and scale of blobs, which can facilitate the tasks of blob reconstruction and blob reduction. Subsequently, we propose a denoising scheme to address high-ISO long-exposure noise, which sometimes spatially shows a blob appearance, employing a blob reduction procedure as a cheap preprocessing for conventional denoising methods. We apply the proposed denoising methods to real-world noisy images as well as standard images that are corrupted by real noise. The experimental results demonstrate the superiority of the proposed methods over state-of-the-art denoising methods

Radiogenic backgrounds in the NEXT double beta decay experiment

Natural radioactivity represents one of the main backgrounds in the search for neutrinoless double beta decay. Within the NEXT physics program, the radioactivity- induced backgrounds are measured with the NEXT-White detector. Data from 37.9 days of low-background operations at the Laboratorio Subterráneo de Canfranc with xenon depleted in 136Xe are analyzed to derive a total background rate of (0.84±0.02) mHz above 1000 keV. The comparison of data samples with and without the use of the radon abatement system demonstrates that the contribution of airborne-Rn is negligible. A radiogenic background model is built upon the extensive radiopurity screening campaign conducted by the NEXT collaboration. A spectral fit to this model yields the specific contributions of 60Co, 40K, 214Bi and 208Tl to the total background rate, as well as their location in the detector volumes. The results are used to evaluate the impact of the radiogenic backgrounds in the double beta decay analyses, after the application of topological cuts that reduce the total rate to (0.25±0.01) mHz. Based on the best-fit background model, the NEXT-White median sensitivity to the two-neutrino double beta decay is found to be 3.5σ after 1 year of data taking. The background measurement in a Qββ±100 keV energy window validates the best-fit background model also for the neutrinoless double beta decay search with NEXT-100. Only one event is found, while the model expectation is (0.75±0.12) events. [Figure not available: see fulltext.]