SEARCH · Engineering Papers

Results for “Source and propagation modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1981

On the transport and propagation of relativistic electrons in galaxies

Analytic solutions to the time-dependent equation describing the transport of relativistic electrons that are being convected and diffused, and may be either gaining or losing energy, are described. The calculations yield a body of solutions from which specific ones can be chosen for the modelling of synchrotron radiation in the Galaxy and in other galaxies. The solutions allow for (1) spatially varying convection velocity; (2) an energy and spatially dependent diffusion coefficient; (3) adiabatic acceleration effects; (4) a term describing spatially dependent synchrotron loss and/or Fermi acceleration; (5) an injection source for the electrons that may be impulsive, steady-state, or variable-injection; (6) a spatial injection or spatially more diffuse source; (7) arbitrary energy dependence of the injection source; and (8) spatial and energy boundaries where the spatial or energy form of the convection, diffusion, and energy supply terms in the basic equation can change character.

Lerche, I.↗

NASA NTRS1985

A General-applications Direct Global Matrix Algorithm for Rapid Seismo-acoustic Wavefield Computations

A new matrix method for rapid wave propagation modeling in generalized stratified media, which has recently been applied to numerical simulations in diverse areas of underwater acoustics, solid earth seismology, and nondestructive ultrasonic scattering is explained and illustrated. A portion of recent efforts jointly undertaken at NATOSACLANT and NORDA Numerical Modeling groups in developing, implementing, and testing a new fast general-applications wave propagation algorithm, SAFARI, formulated at SACLANT is summarized. The present general-applications SAFARI program uses a Direct Global Matrix Approach to multilayer Green's function calculation. A rapid and unconditionally stable solution is readily obtained via simple Gaussian ellimination on the resulting sparsely banded block system, precisely analogous to that arising in the Finite Element Method. The resulting gains in accuracy and computational speed allow consideration of much larger multilayered air/ocean/Earth/engineering material media models, for many more source-receiver configurations than previously possible. The validity and versatility of the SAFARI-DGM method is demonstrated by reviewing three practical examples of engineering interest, drawn from ocean acoustics, engineering seismology and ultrasonic scattering.

Schmidt, H.↗

NASA NTRS1994

Response of the thermosphere and ionosphere to geomagnetic storms

Four numerical simulations have been performed, at equinox, using a coupled thermosphere-ionosphere model, to illustrate the response of the upper atmosphere to geomagnetic storms. The storms are characterized by an increase in magnetospheric energy input at high latitude for a 12-hour period; each storm commences at a different universal time (UT). The initial response at high latitude is that Joule heating raises the temperature of the upper thermosphere and ion drag drives high-velocity neutral winds. The heat source drives a global wind surge, from both polar regions, which propagates to low latitudes and into the opposite hemisphere. The surge has the character of a large-scale gravity wave with a phase speed of about 600 m/s. Behind the surge a global circulation of magnitude 100 m/s is established at middle latitudes, indicating that the wave and the onset of global circulation are manifestations of the same phenomena. A dominant feature of the response is the penetration of the surge into the opposite hemisphere where it drives poleward winds for a few hours. The global wind surge has a preference for the night sector and for the longitude of the magnetic pole and therefore depends on the UT start time of the storm. A second phase of the meridional circulation develops after the wave interaction but is also restricted, in this case by the buildup of zonal winds via the Coriolis interaction. Conservation of angular momentum may limit the buildup of zonal wind in extreme cases. The divergent wind field drives upwelling and composition change on both height and pressure surfaces. The composition bulge responds to both the background and the storm-induced horizontal winds; it does not simply rotate with Earth. During the storm the disturbance wind modulates the location of the bulge; during the recovery the background winds induce a diurnal variation in its position. Equatorward winds in sunlight produce positive ionospheric changes during the main driving phase of the storm. Negative ionospheric phases are caused by increases of molecular nitrogen in regions of sunlight, the strength of which depends on longitude and the local time of the sector during the storm input. Regions of positive phase in the ionosphere persist in the recovery period due to decreases in mean molecular mass in regions of previous downwelling. Ion density changes, expressed as a ratio of disturbed to quiet values, exhibit a diurnal variation that is driven by the location of the composition bulge; this variation explains the ac component of the local time variation of the observed negative storm phase.

Fuller-Rowell, T. J.↗

NASA NTRS1995

Energization in regions of CIRs unconnected to shocks are probably not the result of cross-field transport

Corotating energetic ion populations (CEIPs) associated with the forward and reverse shocks of corotating interaction regions (CIRs) are observed in CIRs at places where models say are magnetically unconnected to either shock. Such disconnections between CEIPs and shocks are common and have been documented with data from Pioneers 10 and 11 and confirmed with data from Ulysses. They pose a problem for models that account for these CEIPs in terms of ion energization at the shocks followed by ion propagation along field lines. Two possible resolutions to this problem have been suggested: diffusion of the ions across field lines and extension of the ion energization process to regions beyond the shock waves. Here we quantitatively examine the first of these possibilities. We give the Green's function solution to the convection-diffusion equation applied to idealized CIR geometry, with a source at the reverse shock -- the main producer of CEIPs. Two kinds of diffusion are considered: resonant diffusion and stochastic field line diffusion. We find that for resonant diffusion the computed ratio is many orders of magnitudes below the observed ratio. For stochastic field line diffusion, the computed ratio approximately equals the observed ratio if a diffusion coefficient appropriate to the free solar wind is used. It is several orders of magnitude below the observed ratio, however, if a diffusion coefficient appropriate to CIRs is used. We conclude that cross-field diffusion probably does not account for the presence of energetic ions in regions of CIRs that are magnetically unconnected to its shock waves. We suggest that the alternative possibility -- the energetic ions in regions magnetically unconnected to shocks result from an acceleration process that is independent of shocks -- be pursued to the point where quantitative tests can be performed.

Intriligator, Devrie S.↗

NASA NTRS2006

Comparison of Orbital Parameters for GEO Debris Predicted by LEGEND and Observed by MODEST: Can Sources of Orbital Debris be Identified?

Since 2002 the National Aeronautics and Space Administration (NASA) has carried out an optical survey of the debris environment in the geosynchronous Earth-orbit (GEO) region with the Michigan Orbital Debris Survey Telescope (MODEST) in Chile. The survey coverage has been similar for 4 of the 5 years allowing us to follow the orbital evolution of Correlated Targets (CTs), both controlled and un-controlled objects, and Un-Correlated Targets (UCTs). Under gravitational perturbations the distributions of uncontrolled objects, both CTs and UCTs, in GEO orbits will evolve in predictable patterns, particularly evident in the inclination and right ascension of the ascending node (RAAN) distributions. There are several clusters (others have used a "cloud" nomenclature) in observed distributions that show evolution from year to year in their inclination and ascending node elements. However, when MODEST is in survey mode (field-of-view approx.1.3deg) it provides only short 5-8 minute orbital arcs which can only be fit under the assumption of a circular orbit approximation (ACO) to determine the orbital parameters. These ACO elements are useful only in a statistical sense as dedicated observing runs would be required to obtain sufficient orbital coverage to determine a set of accurate orbital elements and then to follow their evolution. Identification of the source(s) for these "clusters of UCTs" would be advantageous to the overall definition of the GEO orbital debris environment. This paper will set out to determine if the ACO elements can be used to in a statistical sense to identify the source of the "clustering of UCTs" roughly centered on an inclination of 12deg and a RAAN of 345deg. The breakup of the Titan 3C-4 transtage on February 21, 1992 has been modeled using NASA s LEGEND (LEO-to-GEO Environment Debris) code to generate a GEO debris cloud. Breakup fragments are created based on the NASA Standard Breakup Model (including fragment size, area-to-mass (A/M), and delta-V distributions). Once fragments are created, they are propagated forward in time with a subroutine GEOPROP. Perturbations included in GEOPROP are those due to solar/lunar gravity, radiation pressure, and major geopotential terms. The question to be addressed: are the UCTs detected by MODEST in this inclination/RAAN region related to the Titan 3C-4 breakup? Discussion will include the observational biases in attempting to detect a specific, uncontrolled target during given observing session. These restrictions include: (1) the length of the observing session which is 8 hours or less at any given date or declination; (2) the assumption of ACO elements for detected object when the breakup model predicts debris with non-zero eccentricities; (3) the size and illumination or brightness of the debris predicted by the model and the telescope/sky limiting magnitude.

Barker, Edwin S.↗

NASA NTRS1981

Directionality effects in the transfer of X-rays from a magnetized atmosphere: Beam pulse shape

A formalism is presented for radiation transfer in two normal polarization modes in finite and semiinfinite plane parallel uniform atmospheres with a magnetic field perpendicular to the surface and arbitrary propagation angles. This method is based on the coupled integral equations of transfer, including emission, absorption, and scattering. Calculations are performed for atmosphere parameters typical of X-ray pulsars. The directionality of the escaping radiation is investigated for several cases, varying the input distributions. Theoretical pencil beam profiles and X-ray pulse shapes are obtained assuming the radiation is emitted from the polar caps of spinning neutron stars. Implications for realistic models of accreting magnetized X-ray sources are briefly discussed.

Meszaros, P.↗

NASA NTRS2019

Forecasting the Ambient Solar Wind with Numerical Models. I. On the Implementation of an Operational Framework

The ambient solar wind conditions in interplanetary space and in the near-Earth environment are determined by activity on the Sun. Steady solar wind streams modulate the propagation behavior of interplanetary coronal mass ejections and are themselves an important driver of recurrent geomagnetic storm activity. The knowledge of the ambient solar wind flows and fields is thus an essential component of successful space weather forecasting. Here, we present an implementation of an operational framework for operating, validating, and optimizing models of the ambient solar wind flow on the example of Carrington Rotation 2077. We reconstruct the global topology of the coronal magnetic field using the potential field source surface model (PFSS) and the Schatten current sheet model (SCS) and discuss three empirical relationships for specifying the solar wind conditions near the Sun, namely the Wang–Sheeley (WS) model, the distance from the coronal hole boundary model (DCHB), and the Wang-Sheeley–Arge (WSA) model. By adding uncertainty in the latitude about the sub-Earth point, we select an ensemble of initial conditions and map the solutions to Earth by the Heliospheric Upwind eXtrapolation (HUX) model. We assess the forecasting performance from a continuous variable validation and find that the WSA model most accurately predicts the solar wind speed time series (RMSE≈83 km/s). We note that the process of ensemble forecasting slightly improves the forecasting performance of all solar wind models investigated. We conclude that the implemented framework is well suited for studying the relationship between coronal magnetic fields and the properties of the ambient solar wind flow in the near-Earth environment.

Reiss, Martin A.↗

NASA NTRSDate not supplied

Towards an Optimal Estimation Retrieval of Cirrus Cloud Optical and Microphysical Properties Using Hyperspectral Shortwave Instruments and A Fast Radiative Transfer Algorithm

Cirrus cloud retrieval products (here, cloud optical depth, effective particle size, and cloud top height) are important inputs into numerical weather and climate models. Uncertainties in such retrieval products, as a matter of course, propagate downstream, impacting model calculations. Improvements in high-quality global cirrus cloud optical and microphysical data products from satellite observations are needed to understand and reduce retrieval uncertainties. Hyperspectral instruments produce high-resolution and information-dense radiance spectra, thus offering the opportunity to reduce uncertainties in retrieval products. An optimal estimation-based cirrus cloud optical and microphysical product retrieval is under development for the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and the Earth Surface Mineral Dust Source Investigation (EMIT) instruments, and the forthcoming Climate Absolute Radiance and Refractivity Observatory Pathfinder (CLARREO-Pathfinder). The spectral coverage of all three instruments includes the ultraviolet, visible, and near-infrared. In this study a very fast radiative transfer model, the Principal Component-based Radiative Transfer Model in the solar spectral region (PCRTM-Solar), is used for forward modeling computations. This will reduce the computational burden in the forward radiance and Jacobian calculations as the cost function is minimized, enabling the entire AVIRIS, EMIT, and CLARREO-Pathfinder spectral range to be used in the optimal estimation framework. Using the entire spectrum will maximize the information content, resulting in a more robust and more accurate retrieval. As a first step, the retrieval is being designed for single layer ice clouds over open ocean water. Preliminary results will be shown.

Jeffrey Mast↗

NASA NTRS1974

International Cosmic Ray Conference, 13th, University of Denver, Denver, Colo., August 17-30, 1973, Proceedings. Volumes 1 & 2

Topics discussed include diffuse gamma rays, gamma ray sources (in particular, the Crab Nebula and Crab Pulsar), X-ray sources, isotopic and nuclear composition, spectra of nuclei, searches for antinuclei, studies of VH and VVH nuclei, measurements of electrons and positrons, the spectra of high-energy electrons, cosmic gamma rays, cosmic ray sources, interplanetary radial gradients, the theory and observation of electron modulation, short-period cosmic ray intensity variations, atmospheric and coupling effects on neutron monitors and muon telescopes, secondary particles and photons in the atmosphere, harmonics of diurnal variations, geomagnetic effects and cutoffs, planetary and interplanetary effects, solar particle access observations, long-term cosmic ray intensity modulation by solar activity and interplanetary effects, Forbush decreases, modulation and anisotropy of low-energy particles, the theory and models of solar particle propagation, solar flare particle composition, and relativistic solar particles. Individual items are announced in this issue.

Source record↗

NASA NTRS2003

The Consequences of Alfven Waves and Parallel Potential Drops in the Auroral Zone

The goal of this research is to examine the causes of field-aligned plasma acceleration in the auroral zone using satellite data and numerical simulations. A primary question to be addressed is what causes the field-aligned acceleration of electrons (leading to precipitation) and ions (leading to upwelling ions) in the auroral zone. Data from the Fast Auroral SnapshoT (FAST) and Polar satellites is used when the two satellites are in approximate magnetic conjunction and are in the auroral region. FAST is at relatively low altitudes and samples plasma in the midst of the auroral acceleration region while Polar is at much higher altitudes and can measure plasmas and waves propagating towards the Earth. Polar can determine the sources of energy streaming earthward from the magnetotail, either in the form of field-aligned currents, electromagnetic waves or kinetic particle energy, that ultimately leads to the acceleration of plasma in the auroral zone. After identifying and examining several events, numerical simulations are run that bridges the spatial region between the two satellites. The code is a one-dimensional, long system length particle in cell simulation that has been developed to model the auroral region. A main goal of this research project is to include Alfven waves in the simulation to examine how these waves can accelerate plasma in the auroral zone.

Schriver, David↗

NASA NTRS2006

A non-PRE double-peaked burst with oscillations: burning front propagation and stalling

Non-photospheric-radius-expansion (non-PRE) double-peaked bursts may be explained in terms of spreading (and temporary stalling) of thermonuclear flames from a rotational pole on the neutron star surface, as we argued in a previous study. Here we analyze Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) data of such a burst from the low mass X-ray binary (LMXB) system 4U 1636-536, and show that our model (with ignition at high latitudes) can qualitatively explain the observed burst profile, and spectral evolution. Moreover, the evolution of the source radius inferred from the data shows a strong signature of temporary stalling of the burning front, which is an essential ingredient of our model. This implies that an understanding of thermonuclear flame spreading on neutron stars can be achieved by a simultaneous study of the evolution of intensity and spectrum of these bursts. We also report the discovery of millisecond period brightness oscillations from this burst, which is the first such observation from a non-PRE double-peaked burst. Our model can explain the corresponding oscillation amplitude during the first (weaker) peak, and the absence of oscillations during the second peak. We discuss how observations of oscillations during non-PRE double-peaked bursts provide an additional t 001 for understanding thermonuclear flame spreading successfully.

Bhattacharyya, Sudip↗

NASA NTRS2020

GRACE gravitational measurements of tsunamis after the 2004, 2010, and 2011 great earthquakes

The 2004 Sumatra, 2010 Maule, and 2011 Tohoku great earthquakes triggered tsunamis as large as a few decimeters over a few 100 km in the open ocean. The transient ocean mass redistribution propagating as tsunamis changed the Earth’s gravity field enough to perturb the GRACE satellites’ orbits at ~ 500 km above the surface. The on-board microwave ranging system detected inter-satellite acceleration anomalies of up to 1.0–4.0 nm/s2. There is good agreement between GRACE measurements and tsunami models for the three events. Complementarily to buoys, ocean bottom pressure sounders, and satellite altimeters, GRACE is sensitive to the long-wavelength spatial scale of tsunamis and provides an independent source of information for assessing alternate early earthquake and tsunami models. Our study demonstrates an innovative way of applying GRACE and GRACE Follow-On data to detect transient geophysical mass changes which cannot be observed by the conventional monthly Level-2 and mascon solutions.

Khosro Ghobadi-Far↗

NASA NTRS1996

High-Resolution X-Ray Imaging of Colliding Radio-Jet Galaxies

We received ROSAT data for four program objects:3C31,3C278,3C449,and NGC1044. The first three sources were observed with the ROSAT HRI instrument. Our plan was to use the HRI to image the hot gas distribution in a few pairs of strongly disturbed interacting elliptical galaxies which are also strong radio sources having a bent-jet source morphology. The PSPC was used for NGC1044 in order to obtain a flux measurement to use in planning future High Resolution Imager (HRI) observations of that source. Though we never requested such HRI observations of NGC1044, others have used those archival PSPC data from our project for other research projects and analyses. The goal of the program was to elucidate the detailed distribution of hot gas into which the jets flow. The X-ray data were consequently analyzed in conjunction with existing VLA radio maps, optical broad-band and H-alpha Charge Couple device (CCD) images, and optical kinematic data to constrain models for the propagation of ballistic jets in interacting galaxies. We were able to test and validate the claimed causal connection between tidal interaction, the presence of gas, and the onset of activity in galaxies. The full multi-wavelength multi-observatory analyses described here are still on-going and will be published in the future. Because of the relevance of this research to on-going work in the field of active galaxies, the grant was used to support travel to several scientific meetings where our x-ray analysis, numerical modeling, and related radio results were presented and discussed.

Born, Kirk D.↗

NASA NTRS1998

Current Limitations on VLBI Accuracy

The contribution of VLBI to geophysics and geodesy arises from its ability to measure distances between stations in a network and to determine the orientation of stations in the network as well as the orientation of the network with respect to the external reference frame of extragalactic radio objects. Integrating nearly two decades of observations provides useful information about station positions and velocities and the orientation of the Earth, but the complications of the real world and the limitations of observing, modeling and analysis prevent recovery of all effects. Of the factors that limit the accuracy of seemingly straightforward geodetic parameters, the neutral propagation medium has been subject to the greatest scrutiny, but the treatment of the mapping function, the wet component and spatial/temporal inhomogeneities is still improving. These affect both the terrestrial scale and consistency over time. The modeling of non-secular site motions (tides and loading) has increased in sophistication, but there are some differences between the models and the observations. VLBI antennas are massive objects, so their behavior is quite unlike GPS monuments, but antenna deformations add some (generally) unmodeled signal. Radio sources used in geodetic VLBI observations are selected for strength and (relative) absence of structure, but apparent changes in position can leak into geodetic parameters. A linear rate of change of baseline or site parameters is the simplest model and its error improves with time span. However, in most cases the VLBI data distribution is insufficient to look for real non-linear behavior that might affect the average rate. A few sites have multiple VLBI antennas, and some show small differences in rate. VLBI intrinsically measures relative positions and velocities, but individual site positions and velocities are generally more useful. The creation of the VLBI terrestrial reference frame, which transforms relative information into individual results, is an empirical process that has intrinsic errors. While UT1 is uniquely measured by VLBI, the geographical distribution and availability of VLBI stations, especially in the southern hemisphere, and the consistency of the VLBI terrestrial reference frame may limit the accuracy of Earth orientation measurements. The effects of particular error sources on geodetic and geophysical parameters derived from VLBI data will be illustrated.

Ma, Chopo↗

NASA NTRS1974

Propagation of sound through a real jet flowfield

The sound propagation of harmonic disturbances through a real jet flowfield, which contains mean flow gradients in pressure and velocities, has been studied. A finite-difference approximation for the equations which govern the acoustic disturbance is obtained for a subsonic axisymmetric jet in terms of acoustic pressure and acoustic energy flux. Experimental results are used to describe the mean flowfield of a model jet. The equations are solved numerically using a Newton-type iterative scheme. The directivity at various distances from a point source located at two jet diameters downstream of the jet exit on its center line compared qualitatively well with the numerical result of Schubert (1972) and with the experimental result of Grande (1966).

Liu, C. H.↗

NASA NTRS1993

SeaTrack: Ground station orbit prediction and planning software for sea-viewing satellites

An orbit prediction software package (Sea Track) was designed to assist High Resolution Picture Transmission (HRPT) stations in the acquisition of direct broadcast data from sea-viewing spacecraft. Such spacecraft will be common in the near future, with the launch of the Sea viewing Wide Field-of-view Sensor (SeaWiFS) in 1994, along with the continued Advanced Very High Resolution Radiometer (AVHRR) series on NOAA platforms. The Brouwer-Lyddane model was chosen for orbit prediction because it meets the needs of HRPT tracking accuracies, provided orbital elements can be obtained frequently (up to within 1 week). Sea Track requires elements from the U.S. Space Command (NORAD Two-Line Elements) for the satellite's initial position. Updated Two-Line Elements are routinely available from many electronic sources (some are listed in the Appendix). Sea Track is a menu-driven program that allows users to alter input and output formats. The propagation period is entered by a start date and end date with times in either Greenwich Mean Time (GMT) or local time. Antenna pointing information is provided in tabular form and includes azimuth/elevation pointing angles, sub-satellite longitude/latitude, acquisition of signal (AOS), loss of signal (LOS), pass orbit number, and other pertinent pointing information. One version of Sea Track (non-graphical) allows operation under DOS (for IBM-compatible personal computers) and UNIX (for Sun and Silicon Graphics workstations). A second, graphical, version displays orbit tracks, and azimuth-elevation for IBM-compatible PC's, but requires a VGA card and Microsoft FORTRAN.

Lambert, Kenneth S.↗

NASA NTRS2006

Understanding X-Ray Source Motions in a Solar Flare Loop

RHESSI images of a solar flare on 2002 November 28 showed a 3-6 keV hard X-ray source that was initially located at the flare loop top, split and propagated to the foot points of the loop during the flare rise phase, and then propagated back up to the loop top during the declining phase of the flare (Sai, Holman, & Dennis 2006). Higher energy X-ray sources were located lower in the legs of the loop during this period of source evolution, with X-rays above 25 keV seen only at the foot points. Sui, Holman, & Dennis suggested that this spatial evolution reflected the evolution of the spectral index and low-energy cutoff to the distribution of accelerated electrons in the flare. We construct a model flare loop and electron distribution injected at the top of this loop to reproduce the source evolution of the November 28 flare. We determine the constraints on the loop model and the evolution of the accelerated electron distribution. We also study the implications of the model for energy deposition into the loop plasma, and the integrated and imaged X-ray spectra. This work is supported in part by the RHESSI Project and the NASA Guest Investigator Program.

Holman, Gordon D.↗

NASA NTRS2012

The Extragalactic Background Light and the Gamma-ray Opacity of the Universe

The extragalactic background light (EBL) is one of the fundamental observational quantities in cosmology. All energy releases from resolved and unresolved extragalactic sources, and the light from any truly diffuse background, excluding the cosmic microwave background (CMB), contribute to its intensity and spectral energy distribution. It therefore plays a crucial role in cosmological tests for the formation and evolution of stellar objects and galaxies, and for setting limits on exotic energy releases in the universe. The EBL also plays an important role in the propagation of very high energy gamma-rays which are attenuated en route to Earth by pair producing gamma-gamma interactions with the EBL and CMB. The EBL affects the spectrum of the sources, predominantly blazars, in the approx 10 GeV to 10 TeV energy regime. Knowledge of the EBL intensity and spectrum will allow the determination of the intrinsic blazar spectrum in a crucial energy regime that can be used to test particle acceleration mechanisms and VHE gamma-ray production models. Conversely, knowledge of the intrinsic gamma-ray spectrum and the detection of blazars at increasingly higher redshifts will set strong limits on the EBL and its evolution. This paper reviews the latest developments in the determination of the EBL and its impact on the current understanding of the origin and production mechanisms of gamma-rays in blazars, and on energy releases in the universe. The review concludes with a summary and future directions in Cherenkov Telescope Array techniques and in infrared ground-based and space observatories that will greatly improve our knowledge of the EBL and the origin and production of very high energy gamma-rays.

Dwek, Eli↗