High-energy Data Research Articles

A Very High Energy γ-Ray Survey toward the Cygnus Region of the Galaxy

Abstract We present results from deep observations toward the Cygnus region using 300 hr of very high energy (VHE) γ-ray data taken with the VERITAS Cerenkov telescope array and over 7 yr of high-energy γ-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse γ-ray emission in the northern sky, the Cygnus region provides a promising area to probe the origins of cosmic rays. We report the identification of a potential Fermi-LAT counterpart to VER J2031+415 (TeV J2032+4130) and resolve the extended VHE source VER J2019+368 into two source candidates (VER J2018+367* and VER J2020+368*) and characterize their energy spectra. The Fermi-LAT morphology of 3FGL J2021.0+4031e (the Gamma Cygni supernova remnant) was examined, and a region of enhanced emission coincident with VER J2019+407 was identified and jointly fit with the VERITAS data. By modeling 3FGL J2015.6+3709 as two sources, one located at the location of the pulsar wind nebula CTB 87 and one at the quasar QSO J2015+371, a continuous spectrum from 1 GeV to 10 TeV was extracted for VER J2016+371 (CTB 87). An additional 71 locations coincident with Fermi-LAT sources and other potential objects of interest were tested for VHE γ-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula flux. We interpret these observations in a multiwavelength context and present the most detailed γ-ray view of the region to date.

Deep Learning and Its Application to LHC Physics

Machine learning has played an important role in the analysis of high-energy physics data for decades. The emergence of deep learning in 2012 allowed for machine learning tools which could adeptly handle higher-dimensional and more complex problems than previously feasible. This review is aimed at the reader who is familiar with high-energy physics but not machine learning. The connections between machine learning and high-energy physics data analysis are explored, followed by an introduction to the core concepts of neural networks, examples of the key results demonstrating the power of deep learning for analysis of LHC data, and discussion of future prospects and concerns.

Coexistence of ferroelectric phases and electric field induced structural transformation in sodium potassium bismuth titanate ceramics

Lead-free sodium potassium bismuth titanate (1−x)[Na0.5Bi0.5TiO3]-x[K0.5Bi0.5TiO3] (NKBTx) ceramics have been believed to be a rhombohedral crystal system with R3c symmetry up to x = 0.16 mole fraction, right now. For the first time, on the contrary, we have disclosed the existence of a low-symmetry (MA)-type monoclinic phase (69%) along with rhombohedral R3c (31%) in the NKBT0.10 sample from Rietveld studies on conventional and high-energy synchrotron powder diffraction data in the unpoled state. The peak splitting and reduced FWHM of {110}pc and {111}pc reflections after poling indicate the field-induced structural transformation. Rietveld analysis reveals that the monoclinic (Cc) phase is partially transformed to rhombohedral (R3c) at 50 kV/cm. Two distinct dielectric anomalies are noticed in the dielectric curve at 120 and 310 °C corresponding to freezing temperature (Tf) and glass-like transition temperature (Tm), respectively. The dielectric diffusivity (γ) estimated from the modified Curie-Weiss law reveals that the sample has a relaxor behavior. Structural inhomogeneity in the nanoscale region due to cation disorder in the A-site causes relaxing nature. A well-saturated square-like ferroelectric hysteresis loop is obtained with a remnant polarization and a coercive field of 35.8 μC/cm2 and 47 kV/cm, respectively. A piezoelectric constant (d33) of 85 pC/N is obtained for the sample which is somewhat higher than that of pure Na0.5Bi0.5TiO3 ceramics.

High-energy nuclear data uncertainties propagated to MYRRHA safety parameters

Propagation of high-energy (above 20 MeV) nuclear data uncertainties on the safety related neutronic responses in accelerator driven systems has been assessed. The total core power and production of radionuclides contributing to radiation source terms were focused on. The article features a method based on the Monte Carlo sampling of random nuclear data files from the covariance matrices generated from the sets of reaction cross sections obtained with model calculations of high-energy particle interactions with matter or picked up from already existing nuclear data libraries. It has been demonstrated that nuclear data uncertainties do not need to be propagated through particle transport calculations to obtain uncertainties on the responses. This advantage allowed to investigate the convergence of the sample average to the best estimate. The number of random nuclear data file sets needed to obtain reliable uncertainty on the total core power is around 300 that results in the uncertainty of 14%. The uncertainties on the concentrations of nuclides most important for the safety assessment that are accumulated in lead–bismuth eutectic during irradiation, range from 5 to 60%. Concentrations of some nuclides exemplified by Tritium converge much slower than neutron multiplicities so that several thousands of samples are needed to ensure reliable uncertainty estimates.

Initializing type-2 residual stresses in crystal plasticity finite element simulations utilizing high-energy diffraction microscopy data

Crystal plasticity finite element models have advanced in recent years by increasing their fidelity via incorporating physically based deformation mechanisms and more realistic microstructures. In order to better understand their effectiveness, researchers have looked at comparing their results with that of experiments at the grain level. In most of these efforts, the importance of incorporating initial residual stresses, representing the true grain morphologies, and imposing correct boundary conditions, which may have a significant effect on the results, is often overlooked. This work utilizes an available dataset of high-energy X-ray diffraction microscopy experiments conducted on a titanium alloy, Ti-7Al, providing complete grain averaged elastic strain tensors to investigate these three issues in crystal plasticity models. In this study, a method to initialize grain-level residual stresses is formulated and its effect on the correlation of results between simulations and experiments, is evaluated. The effect of grain morphology on the simulation results is evaluated by comparing results from simulations using exact 3D grain morphology with that of simulations using tessellated grain structures. In addition, the importance of applying physically realistic boundary conditions, obtained directly from experiments is also investigated. The findings of this work indicate that initial residual stresses and physically realistic boundary conditions play a key role in the results obtained from CPFE simulations and it is imperative that more importance is given to them. On the other hand, using exact grain morphology is not of much benefit while comparing data on the grain averaged scale as long as the microstructure contains information about the position of the grains, as is the case with tessellations, although grain morphology plays an important role in predicting intragranular strain localization.

Analyticity constraints for hadron amplitudes: Going high to heal low-energy issues

Analyticity constitutes a rigid constraint on hadron scattering amplitudes. This property is used to relate models in different energy regimes. Using meson photoproduction as a benchmark, we show how to test contemporary low-energy models directly against high-energy data. This method pinpoints deficiencies of the models and treads a path to further improvement. The implementation of this technique enables one to produce more stable and reliable partial waves for future use in hadron spectroscopy and new physics searches.

Exploring the origin of broad-band emissions of Mrk 501 with a two-zone model

Abstract We propose a two-zone synchrotron self-Compton (SSC) model, including an inner gamma-ray emitting region with spherical shape and a conical radio emitting region located at the extended jet, to alleviate the long-standing “bulk Lorentz factor crisis” in blazars. In this model, the spectral energy distributions (SEDs) of blazars are produced by considering the gamma-ray emitting region inverse Compton scattering of both the synchrotron photons itself and the ambient photons from the radio emitting region. Applying the model to Mrk 501, we obtain that the radio emitting region has a comoving length of ∼0.15 pc and is located at sub-parsec scale from the central engine by modeling the radio data; the flux of the Compton scattering of the ambient photons is so low that it can be neglected safely. The characteristic hard gamma-ray spectrum can be explained by the superposition of two SSC processes, and the model can approximately explain the very high energy (VHE) data. The insights into the spectral shape and the inter-band correlations under the flaring state will provide us with a diagnostic for the bulk Lorentz factor of radio emitting region, where the low and upper limits of 8 and 15 are preferred, and for the two-zone SSC model itself. In addition, our two-zone SSC model shows that the gamma-ray emitting region creates flare on the timescale of merely a few hours, and the long time outbursts more likely originate from the extended radio emitting region.

Supervised classification of electric power transmission line nominal voltage from high-resolution aerial imagery

For many researchers, government agencies, and emergency responders, access to the geospatial data of US electric power infrastructure is invaluable for analysis, planning, and disaster recovery. Historically, however, access to high quality geospatial energy data has been limited to few agencies because of commercial licenses restrictions, and those resources which are widely accessible have been of poor quality, particularly with respect to reliability. Recent efforts to develop a highly reliable and publicly accessible alternative to the existing datasets were met with numerous challenges – not the least of which was filling the gaps in power transmission line voltage ratings. To address the line voltage rating problem, we developed and tested a basic methodology that fuses knowledge and techniques from power systems, geography, and machine learning domains. Specifically, we identified predictors of nominal voltage that could be extracted from aerial imagery and developed a tree-based classifier to classify nominal line voltage ratings. Overall, we found that line support height, support span, and conductor spacing are the best predictors of voltage ratings, and that the classifier built with these predictors had a reliable predictive accuracy (that is, within one voltage class for four out of the five classes sampled). We applied our approach to a study area in Minnesota.

Growth Kinetics and Microstructure of PbTe Films Produced on Si and BaF2 Substrates by a Modified Hot-Wall Method

Lead telluride films have been grown on Si (100) and BaF2 (100) substrates by a modified hot-wall method using a graphite reaction chamber. According to X-ray diffraction, X-ray microanalysis, and scanning electron microscopy characterization results, the average growth rate of PbTe films having compositions within the homogeneity range of lead telluride increases with increasing lead vapor partial pressure and decreases with increasing tellurium vapor partial pressure, independent of the nature of the substrate. The rate of PbTe film growth has been shown to be maximal in the initial stage of the process and decrease monotonically over time, independent of the nature of the substrate. Independent of the growth time, the average growth rate of the PbTe films on the Si (100) substrates is considerably higher than that on the BaF2 (100) substrates. Reflection high-energy electron diffraction data indicate that the texture of the PbTe films on Si (100) corresponds to the substrate orientation and that the misorientation angle of the mosaic blocks does not exceed 20°. On the BaF2 (100) substrates, we observe epitaxial PbTe film growth with the orientation relationship (100), [011] PbTe ║ (100), [011] BaF2.

Extragalactic diffuse γ-rays from dark matter annihilation: revised prediction and full modelling uncertainties

Recent high-energy data from Fermi-LAT on the diffuse γ-ray background have been used to set among the best constraints on annihilating TeV cold dark matter candidates. In order to assess the robustness of these limits, we revisit and update the calculation of the isotropic extragalactic γ-ray intensity from dark matter annihilation. The emission from halos with masses ⩾ 1010 M⊙ provides a robust lower bound on the predicted intensity. The intensity including smaller halos whose properties are extrapolated from their higher mass counterparts is typically 5 times higher, and boost from subhalos yields an additional factor ∼ 1.5. We also rank the uncertainties from all ingredients and provide a detailed error budget for them. Overall, our fiducial intensity is a factor 5 lower than the one derived by the Fermi-LAT collaboration in their latest analysis. This indicates that the limits set on extragalactic dark matter annihilations could be relaxed by the same factor. We also calculate the expected intensity for self-interacting dark matter in massive halos and find the emission reduced by a factor 3 compared to the collisionless counterpart. The next release of the CLUMPY code will provide all the tools necessary to reproduce and ease future improvements of this prediction.

Single-Scan Dual-Energy CT Using Primary Modulation.

Compared with conventional computed tomography (CT), dual-energy CT (DECT) provides better material differentiation but requires projection data acquired with two different effective x-ray spectra, limiting DECT applications to specialized scanners. We propose a hardware-based method, known as PM-DECT, which utilizes primary beam modulation to enable single-scan DECT on a conventional CT scanner. PM-DECT inserts an attenuation sheet with a spatially varying pattern-primary beam modulator-between the x-ray source and imaged object. During a CT scan, the modulator selectively hardens the x-ray beam, thereby increasing the average photon energy at specific detector pixel locations. Thus, PM-DECT simultaneously acquires high and low energy data at each projection angle. From the sparse projection data, high and low energy CT images are jointly reconstructed and simultaneously decomposed into basis materials via an iterative CT reconstruction algorithm with gradient weighting and an improved version of similarity based regularization. Studies on Catphan 600 and anthropomorphic head phantoms demonstrate that PM-DECT retains a high level of spatial resolution compared with conventional CT scans. Electron density values calculated from decomposed images indicate a limited error of 1.12% for PM-DECT. Comparison against a two-scan DECT technique shows that PM-DECT's image reconstruction from sparse data sets contributes only 0.66% error. By granting the opportunity for high-quality single-scan DECT on conventional CT scanners via limited hardware modification, PM-DECT has the potential to liberate DECT from specialized scanners, extending clinical availability, and implementation.

Multi-epoch VLBA Imaging of 20 New TeV Blazars: Apparent Jet Speeds

Abstract We present 88 multi-epoch Very Long Baseline Array (VLBA) images (most at an observing frequency of 8 GHz) of 20 TeV blazars, all of the high-frequency-peaked BL Lac (HBL) class, that have not been previously studied at multiple epochs on the parsec scale. From these 20 sources, we analyze the apparent speeds of 43 jet components that are all detected at four or more epochs. As has been found for other TeV HBLs, the apparent speeds of these components are relatively slow. About two-thirds of the components have an apparent speed that is consistent (within 2σ) with no motion, and some of these components may be stationary patterns whose apparent speed does not relate to the underlying bulk flow speed. In addition, a superluminal tail to the apparent speed distribution of the TeV HBLs is detected for the first time, with eight components in seven sources having a 2σ lower limit on the apparent speed exceeding . We combine the data from these 20 sources with an additional 18 sources from the literature to analyze the complete apparent speed distribution of all 38 TeV HBLs that have been studied with very long baseline interferometry at multiple epochs. The highest 2σ apparent speed lower limit considering all sources is . This suggests that bulk Lorentz factors of up to about 4, but probably not much higher, exist in the parsec-scale radio-emitting regions of these sources, consistent with estimates obtained in the radio by other means such as brightness temperatures. This can be reconciled with the high Lorentz factors estimated from the high-energy data if the jet has velocity structures consisting of different emission regions with different Lorentz factors. In particular, we analyze the current apparent speed data for the TeV HBLs in the context of a model with a fast central spine and a slower outer layer.

Time evolution of gamma rays from supernova remnants

We present a systematic phenomenological study focused on the time evolution of the non-thermal radiation - from radio waves to gamma rays - emitted by typical supernova remnants via hadronic and leptonic mechanisms, for two classes of progenitors: thermonuclear and core-collapse. To this aim, we develop a numerical tool designed to model the evolution of the cosmic-ray spectrum inside a supernova remnant, and compute the associated multi-wavelength emission. We demonstrate the potential of this tool in the context of future population studies based on large collection of high-energy gamma-ray data. We discuss and explore the relevant parameter space involved in the problem, and focus in particular on their impact on the maximum energy of accelerated particles, in order to study the effectiveness and duration of the PeVatron phase. We outline the crucial role of the ambient medium through which the shock propagates during the remnant evolution. In particular, we point out the role of dense clumps in creating a significant hardening in the hadronic gamma-ray spectrum.

Enhanced source location privacy mechanism for WSNs

Source location privacy in wireless sensor networks is an important security issue when a wireless sensor network is used in monitoring valuable assets or the source is a sensitive object. However, the open nature of a sensor network makes relatively easy for an adversary to detect message flows and trace back the message hop-by-hop to it source by moving in the reverse direction of the flows. Many schemes have been proposed recently to provide source location privacy but all performs poorly, because either with insufficient safety period or high energy consumption or data delivery delay. In this paper, after analysing privacy weakness of the cyclic entrapment method (CEM), we propose an enhanced scheme by introducing two mechanisms to provide strong location privacy against skilled adversary. The proposed strategy is principally based on false paths instead of loops and ensures that the arrival angle formed between the real and the real path is minimal. We evaluate our scheme based on two criteria: safety period, and ratio delivery. The simulation results show that our scheme improve safety period while keeping packets delivery ratio similar as in CEM method.

Origin of the High-energy Neutrino Flux at IceCube

Abstract We discuss the spectrum of the different components in the astrophysical neutrino flux reaching the Earth, and the possible contribution of each component to the high-energy IceCube data. We show that the diffuse flux from cosmic ray (CR) interactions with gas in our galaxy implies just two events among the 54-event sample. We argue that the neutrino flux from CR interactions in the intergalactic (intracluster) space depends critically on the transport parameter δ describing the energy dependence in the diffusion coefficient of galactic CRs. Our analysis motivates a neutrino spectrum with a drop at PeV energies that fits the data well, including the non-observation of the Glashow resonance at 6.3 PeV. We also show that a CR flux described by an unbroken power law may produce a neutrino flux with interesting spectral features (bumps and breaks) related to changes in the CR composition.

Consistent searches for SMEFT effects in non-resonant dijet events

We investigate the bounds which can be placed on generic new-physics contributions to dijet production at the LHC using the framework of the Standard Model Effective Field Theory, deriving the first consistently-treated EFT bounds from non-resonant high-energy data. We recast an analysis searching for quark compositeness, equivalent to treating the SM with one higher-dimensional operator as a complete UV model. In order to reach consistent, model-independent EFT conclusions, it is necessary to truncate the EFT effects consistently at order 1/Λ2 and to include the possibility of multiple operators simultaneously contributing to the observables, neither of which has been done in previous searches of this nature. Furthermore, it is important to give consistent error estimates for the theoretical predictions of the signal model, particularly in the region of phase space where the probed energy is approaching the cutoff scale of the EFT. There are two linear combinations of operators which contribute to dijet production in the SMEFT with distinct angular behavior; we identify those linear combinations and determine the ability of LHC searches to constrain them simultaneously. Consistently treating the EFT generically leads to weakened bounds on new-physics parameters. These constraints will be a useful input to future global analyses in the SMEFT framework, and the techniques used here to consistently search for EFT effects are directly applicable to other off-resonance signals.

A Feasibility Study of Low-Dose Single-Scan Dual-Energy Cone-Beam CT in Many-View Under-Sampling Framework.

A single-scan dual-energy low-dose cone-beam CT (CBCT) imaging technique that exploits a multi-slit filter is proposed in this paper. The multi-slit filter installed between the x-ray source and the scanned object is reciprocated during a scan. The x-ray beams through the slits would generate relatively low-energy x-ray projection data, while the filtered beams would make high-energy projection data. An iterative image reconstruction algorithm that uses an adaptive-steepest-descent method to minimize image total-variation under the constraint of data fidelity was applied to reconstructing the image from the low-energy projection data. Since the high-energy projection data suffer from a substantially high noise level due to the beam filtration, we have developed a new algorithm that exploits the joint sparsity between the low- and high-energy CT images for image reconstruction of the high-energy CT image. The proposed image reconstruction algorithm uses a gradient magnitude image (GMI) of the low-energy CT image by regularizing the difference of GMIs of the low- and high-energy CT images to be minimized. The feasibility of the proposed technique has been demonstrated by the use of various phantoms in the experimental CBCT setup. Furthermore, based on the proposed dual-energy imaging, a material differentiation was performed and its potential utility has been shown. The proposed imaging technique produced promising results for its potential application to a low-dose single-scan dual-energy CBCT.

Comprehensive Analysis of the Geoeffective Solar Event of 21 June 2015: Effects on the Magnetosphere, Plasmasphere, and Ionosphere Systems

A full-halo coronal mass ejection (CME) left the Sun on 21 June 2015 from active region (AR) NOAA 12371. It encountered Earth on 22 June 2015 and generated a strong geomagnetic storm whose minimum Dst value was −204 nT. The CME was associated with an M2-class flare observed at 01:42 UT, located near disk center (N12 E16). Using satellite data from solar, heliospheric, and magnetospheric missions and ground-based instruments, we performed a comprehensive Sun-to-Earth analysis. In particular, we analyzed the active region evolution using ground-based and satellite instruments (Big Bear Solar Observatory (BBSO), Interface Region Imaging Spectrograph (IRIS), Hinode, Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO), Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), covering H $\upalpha$ , EUV, UV, and X-ray data); the AR magnetograms, using data from SDO/Helioseismic and Magnetic Imager (HMI); the high-energy particle data, using the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) instrument; and the Rome neutron monitor measurements to assess the effects of the interplanetary perturbation on cosmic-ray intensity. We also evaluated the 1 – 8 A soft X-ray data and the ${\sim}\, 1$ MHz type III radio burst time-integrated intensity (or fluence) of the flare in order to predict the associated solar energetic particle (SEP) event using the model developed by Laurenza et al. (Space Weather 7(4), 2009). In addition, using ground-based observations from lower to higher latitudes (International Real-time Magnetic Observatory Network (INTERMAGNET) and European Quasi-Meridional Magnetometer Array (EMMA)), we reconstructed the ionospheric current system associated with the geomagnetic sudden impulse (SI). Furthermore, Super Dual Auroral Radar Network (SuperDARN) measurements were used to image the global ionospheric polar convection during the SI and during the principal phases of the geomagnetic storm. In addition, to investigate the influence of the disturbed electric field on the low-latitude ionosphere induced by geomagnetic storms, we focused on the morphology of the crests of the equatorial ionospheric anomaly by the simultaneous use of the Global Navigation Satellite System (GNSS) receivers, ionosondes, and Langmuir probes onboard the Swarm constellation satellites. Moreover, we investigated the dynamics of the plasmasphere during the different phases of the geomagnetic storm by examining the time evolution of the radial profiles of the equatorial plasma mass density derived from field line resonances detected at the EMMA network ( $1.5 < \mathrm{L} < 6.5$ ). Finally, we present the general features of the geomagnetic response to the CME by applying innovative data analysis tools that allow us to investigate the time variation of ground-based observations of the Earth’s magnetic field during the associated geomagnetic storm.

Structure fluctuation in Gd- and Mg-substituted BaTiO3 with cubic structure

Gd- and Mg-substituted barium titanate crystals (Ba1−xGdx)(Ti1−x/2Mgx/2)O3 with the Gd composition range of x > 0.1 show no ferroelectric phase transition. To discuss atomic thermal fluctuations in such unusual barium titanate crystals with cubic symmetry, accurate crystal structure analyses were performed by the maximum entropy method (MEM)/Rietveld method with the split-atom model and anharmonic temperature factor model using high-energy synchrotron radiation X-ray diffraction data. The off-centered positions of Gd ions at the Ba site in the probability density function were clearly observed at a low temperature but smeared in the high-temperature region by the over-barrier motion of Gd ions.

The Orbit of the Gamma-Ray Binary 1FGL J1018.6−5856

Abstract Gamma-ray binaries are a small subclass of the high mass X-ray binary population that exhibit emission across the whole electromagnetic spectrum. We present the radial velocities of 1FGL J1018.6−5856 based on the observations obtained with the Southern African Large Telescope. We combine our measurements with those published in the literature to get a broad phase coverage. The mass function obtained supports a neutron star compact object, although a black hole mass is possible for the very low inclination angles. The improved phase coverage allows constraints to be placed on the orbital eccentricity (e = 0.31 ± 0.16), which agrees with the estimates from the high-energy data.