Search NASA⌕ Search

SEARCH · Search NASA

Results for “SOURCES”

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.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A helium-burning white dwarf binary as a supersoft X-ray source

Abstract Type Ia supernovae are cosmic distance indicators 1,2 , and the main source of iron in the Universe 3,4 , but their formation paths are still debated. Several dozen supersoft X-ray sources, in which a white dwarf accretes hydrogen-rich matter from a non-degenerate donor star, have been observed 5 and suggested as Type Ia supernovae progenitors 6–9 . However, observational evidence for hydrogen, which is expected to be stripped off the donor star during the supernova explosion 10 , is lacking. Helium-accreting white dwarfs, which would circumvent this problem, have been predicted for more than 30 years (refs. 7,11,12 ), including their appearance as supersoft X-ray sources, but have so far escaped detection. Here we report a supersoft X-ray source with an accretion disk whose optical spectrum is completely dominated by helium, suggesting that the donor star is hydrogen-free. We interpret the luminous and supersoft X-rays as resulting from helium burning near the surface of the accreting white dwarf. The properties of our system provide evidence for extended pathways towards Chandrasekhar-mass explosions based on helium accretion, in particular for stable burning in white dwarfs at lower accretion rates than expected so far. This may allow us to recover the population of the sub-energetic so-called Type Iax supernovae, up to 30% of all Type Ia supernovae 13 , within this scenario.

Science & Technology - Other Topics↗
DOE OSTI2020

Source of 12 C in Calvin–Benson cycle intermediates and isoprene emitted from plant leaves fed with 13 CO 2

Feeding 14 CO 2 was crucial to uncovering the path of carbon in photosynthesis. Feeding 13 CO 2 to photosynthesizing leaves emitting isoprene has been used to develop hypotheses about the sources of carbon for the methylerythritol 4-phosphate pathway, which makes the precursors for terpene synthesis in chloroplasts and bacteria. Both photosynthesis and isoprene studies found that products label very quickly (<10 min) up to 80–90% but the last 10–20% of labeling requires hours indicating a source of 12 C during photosynthesis and isoprene emission. Furthermore, studies with isoprene showed that the proportion of slow label could vary significantly. This was interpreted as a variable contribution of carbon from sources other than the Calvin–Benson cycle (CBC) feeding the methylerythritol 4-phosphate pathway. Here, we measured the degree of label in isoprene and photosynthetic metabolites 20 min after beginning to feed 13 CO 2 . Isoprene labeling was the same as labeling of photosynthesis intermediates. High temperature reduced the label in isoprene and photosynthesis intermediates by the same amount indicating no role for alternative carbon sources for isoprene. A model assuming glucose, fructose, and/or sucrose reenters the CBC as ribulose 5-phosphate through a cytosolic shunt involving glucose 6-phosphate dehydrogenase was consistent with the observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The miniJPAS survey: AGN and host galaxy coevolution of X-ray-selected sources

Studies indicate strong evidence of a scaling relation in the local Universe between the supermassive black hole mass (MBH) and the stellar mass of their host galaxies (M * ). They even show similar histories across cosmic times of their differential terms: the star formation rate (SFR) and black hole accretion rate (BHAR). However, a clear picture of this coevolution is far from being understood. We selected an X-ray sample of active galactic nuclei (AGN) up to z = 2.5 in the miniJPAS footprint. Their X-ray to infrared spectral energy distributions (SEDs) have been modeled with the CIGALE code, constraining the emission to 68 bands, from which 54 are the narrow filters from the miniJPAS survey. For a final sample of 308 galaxies, we derived their physical properties, such as their M * , SFR, star formation history (SFH), and the luminosity produced by the accretion process of the central BH (L AGN ). For a subsample of 113 sources, we also fit their optical spectra to obtain the gas velocity dispersion from the broad emission lines and estimated the M BH . We calculated the BHAR in physical units depending on two radiative efficiency regimes. We find that the Eddington ratios (λ Edd ) and its popular proxy (L X /M * ) have a difference of 0.6 dex, on average, and a KS test indicates that they come from different distributions. Our sources exhibit a considerable scatter on the MBH - M * scaling relation, which can explain the difference between λ Edd and its proxy. We also modeled three evolution scenarios for each source to recover the integral properties at z = 0. Using the SFR and BHAR, we show a notable diminution in the scattering between MBH - M * . For the last scenario, we considered the SFH and a simple energy budget for the AGN accretion, and we retrieved a relation similar to the calibrations known for the local Universe. Our study covers ~1 deg 2 in the sky and is sensitive to biases in luminosity. Nevertheless, we show that, for bright sources, the link between the differential values (SFR and BHAR) and their decoupling based on an energy limit is the key that leads to the local M BH - M * scaling relation. In the future, we plan to extend this methodology to a thousand degrees of the sky using JPAS with an X-ray selection from eROSITA, to obtain an unbiased distribution of BHAR and Eddington ratios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On Alfvénic turbulence of solar wind streams observed by Solar Orbiter during March 2022 perihelion and their source regions

It has been recently accepted that the standard classification of the solar wind solely according to flow speed is outdated, and particular interest has been devoted to the study of the origin and evolution of so-called Alfvénic slow solar wind streams and to what extent such streams resemble or differ from fast wind. In March 2022, Solar Orbiter completed its first nominal phase perihelion passage. During this interval, it observed several Alfvénic streams, allowing for characterization of fluctuations in three slow wind intervals (AS1-AS3) and comparison with a fast wind stream (F) at almost the same heliocentric distance. This work makes use of Solar Orbiter plasma parameters from the Solar Wind Analyzer (SWA) and magnetic field measurements from the magnetometer (MAG). The magnetic connectivity to the solar sources of selected solar wind intervals was reconstructed using a ballistic extrapolation based on measured solar wind speed down to the (spherical) source surface at 2.5 R s below which a potential field extrapolation was used to map back to the Sun. The source regions were identified using SDO/AIA observations. A spectral analysis of in situ measured magnetic field and velocity fluctuations was performed to characterize correlations, Alfvénicity, normalized cross-helicity, and residual energy in the frequency domain as well as intermittency of the fluctuations and spectral energy transfer rate estimated via mixed third-order moments. A machine learning technique was used to separate proton core, proton beam, and alpha particles and to study v − b correlations for the different ion populations in order to evaluate the role played by each population in determining the Alfvénic content of solar wind fluctuations. The comparison between fast wind and Alfvénic slow wind intervals highlights the differences between the two solar wind regimes: The fast wind is characterized by larger amplitude fluctuations, and magnetic and velocity fluctuations are closer to equipartition of energy. In fact the Alfvénic slow wind streams appear to be on a spectrum of wind types, with AS1, originating from open field lines neighboring active regions and displaying similarities with the fast wind in terms of fluctuation amplitude and turbulence characteristics, but not with respect to the alpha particles and proton beams. The other two slow streams differed both in their sources as well as plasma characteristics, with AS2 coming from the expansion of a narrow coronal hole corridor and AS3 from a region straddling a pseudostreamer. The latter displayed the coldest and highest density but the slowest stream with the smallest fluctuation amplitude and greatest magnetic energy excess. It also showed the largest scatter in proton beam speeds and the greatest difference in speed between proton beam and alpha particles. This study shows how the old fast–slow solar wind dichotomy, already called into question by the observations of slower Alfvénic solar wind streams, should further be refined, as the Alfvénic slow wind, originating in different solar wind regions, show significant differences in density, temperature, and proton and alpha-particle properties in the inner heliosphere. The observations presented here provide the starting point for a better understanding of the origin and evolution of different solar wind streams as well as the evolving turbulence contained within.

magnetohydrodynamics (MHD)↗
DOE OSTI2020

Simple, compact, high-resolution monochromatic x-ray source for characterization of x-ray calorimeter arrays

X-ray calorimeters routinely achieve very high spectral resolution, typically a few eV full width at half maximum (FWHM). Measurements of calorimeter line shapes are usually dominated by the natural linewidth of most laboratory calibration sources. This compounds the data acquisition time necessary to statistically sample the instrumental line broadening and can add systematic uncertainty if the intrinsic line shape of the source is not well known. To address these issues, we have built a simple, compact monochromatic x-ray source using channel cut crystals. A commercial x-ray tube illuminates a pair of channel cut crystals that are aligned in a dispersive configuration to select the Kα1 line of the x-ray tube anode material. The entire device, including the x-ray tube, can be easily hand-carried by one person and may be positioned manually or using a mechanical translation stage. The output monochromatic beam provides a collimated image of the anode spot with magnification of unity in the dispersion direction (typically 100 μm–200 μm for the x-ray tubes used here) and is unfocused in the cross-dispersion direction so that the source image in the detector plane appears as a line. We measured output count rates as high as 10 count/s/pixel for the Hitomi soft x-ray spectrometer, which had 819 μm square pixels. We implemented different monochromator designs for energies of 5.4 keV (one design) and 8.0 keV (two designs), which have effective theoretical FWHM energy resolution of 0.125 eV, 0.197 eV, and 0.086 eV, respectively; these are well-suited for optimal calibration measurements of state-of-the art x-ray calorimeters. We measured an upper limit for the energy resolution of our Cr Kα1 monochromator of 0.7 eV FWHM at 5.4 keV, consistent with the theoretical prediction of 0.125 eV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Long duration x-ray source development for x-ray diffraction at the National Ignition Facility

We present the results of experiments to produce a 10 ns-long, quasi-monochromatic x-ray source. This effort is needed to support time-resolved x-ray diffraction (XRDt) measurements of phase transitions during laser-driven dynamic compression experiments at the National Ignition Facility. To record XRDt of phase transitions as they occur, we use high-speed (∼1 ns) gated hybrid CMOS detectors, which record multiple frames of data over a timescale of a few to tens of ns. Consequently, to make effective use of these imagers, XRDt needs the x-ray source to be narrow in energy and uniform in time as long as the sensors are active. The x-ray source is produced by a laser irradiated Ge foil. Our results indicate that the x-ray source lasts during the whole duration of the main laser pulse. Both time-resolved and time-integrated spectral data indicate that the line emission is dominated by the He-α complex over higher energy emission lines. Time-integrated spectra agree well with a one-dimensional Cartesian simulation using HYDRA that predicts a conversion efficiency of 0.56% when the incident intensity is 2 × 1015 W/cm2 on a Ge backlighter.

Werellapatha, K. (ORCID:0000000293154054)↗
DOE OSTI2021

Strong effect of scandium source purity on chemical and electronic properties of epitaxial Sc x Al 1–x N/GaN heterostructures

Epitaxial multilayer heterostructures of Sc x Al 1–x N/GaN with Sc contents x = 0.11–0.45 are found to exhibit significant differences in structural quality, chemical impurity levels, and electronic properties depending on the starting Sc source impurity levels. A higher purity source leads to a 2–3 orders of magnitude reduction in the carbon, oxygen, and fluorine unintentional doping densities in MBE-grown Sc x Al 1–x N/GaN multilayers. Electrical measurements of Sc x Al 1–x N/n + GaN single heterostructure barriers show a 5–7 orders of magnitude reduction in the electrical leakage for films grown with a higher purity Sc source at most Sc contents. The measured chemical and electrical properties of epitaxial Sc x Al 1–x N highlight the importance of the starting Sc source material purity for epitaxial device applications that need these highly piezoelectric and/or ferroelectric transition-metal nitride alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A mathematical framework for ejecta cloud dynamics with application to source models and piezoelectric mass measurements

We present a mathematical framework for describing the dynamical evolution of an ejecta cloud generated by a generic ejecta source model. We consider a piezoelectric sensor fielded in the path of an ejecta cloud, for experimental configurations in which the ejecta are created at a singly shocked planar surface and fly ballistically through vacuum to the stationary sensor. To do so, we introduce the concept of a time- and velocity-dependent ejecta “areal mass function.” We derive expressions for the analytic (“true”) accumulated ejecta areal mass at the sensor and the measured (“inferred”) value obtained via the standard method for analyzing piezoelectric voltages. In this way, we derive an exact expression and upper bound for the error imposed upon a piezoelectric ejecta mass measurement (in a perfect system) by the assumption of instantaneous creation, which is commonly required for momentum diagnostic analyses. This error term is zero for truly instantaneous source models; otherwise, the standard piezoelectric analysis is guaranteed to overestimate the true mass. When combined with a piezoelectric dataset, this framework provides a unique solution for the ejecta particle velocity distribution, subject to the assumptions inherent in the data analysis. The framework also leads to strong boundary conditions that any ejecta source model must satisfy in order to be consistent with apparently global properties of piezoelectric measurements from a wide range of experiments. We demonstrate this methodology by applying it to the Richtmyer–Meshkov instability+self-similar velocity distribution ejecta source model currently under development at Los Alamos National Laboratory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spectral characterization of an SPDC source with a fast broadband spectrometer

Knowing the properties of single photons produced in a Spontaneous Parametric Down-Conversion (SPDC) source can be crucial for specific applications and uses. In particular, the spectral properties are of key relevance. Here, we investigate a commercial SPDC source using our fast broadband spectrometer. Our analysis is a valid method for other SPDC sources, as well as other single-photon generation techniques, thus providing a good example of how to use this spectrometer design. We calibrate the spectrometer using known lines of the argon emission spectrum. We show that the two down-converted photons from the SPDC source have different spectral properties depending on the pump power, and under which condition we measured spectrally similar down-converted photons. Finally, we were able to reconstruct and investigate the spectral information for the pump photon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

An implicit solution for Asay foil trajectories generated by separable, sustained-production ejecta source models

We present a simple implicit solution for the time-dependent trajectory of a thin Asay foil ejecta diagnostic for the general case where the impinging ejecta cloud is generated by a source function characterized by an arbitrary (sustained) time dependence and a time-independent (stationary) particle velocity distribution. In the limit that the source function time dependence becomes a delta function, this solution—which is amenable to rapid numerical calculations of arbitrary accuracy—exactly recovers a previously published solution for the special case of instantaneous ejecta production. We also derive simple expressions for the free-surface arrival (catch-up) time as well as the true ejecta areal mass accumulation on the accelerating foil and place bounds on the level of error incurred when applying instant-production mass solutions to a sustained-production trajectory. We demonstrate these solutions with example calculations for hypothetical source functions spanning a wide range of ejecta production durations, velocity distributions, and temporal behaviors. These calculations demonstrate how the foil trajectory is often insensitive to the temporal dependence of the source function, instead being dominated by the velocity distribution. We quantify this insensitivity using a “compatibility score” metric. Under certain conditions, one may capitalize upon this insensitivity to obtain a good approximation of the second integral of the velocity distribution from the observed foil trajectory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗
DOE OSTI2026

Titanium alloy response sensitivity to variations in spectral reconstructions of National Ignition Facility xenon line-emission x-ray sources

Thermomechanical shock experiments on the National Ignition Facility (NIF) aim to study high strain rate dynamic material response. In such experiments, the NIF laser is used to generate high fluence x-ray emission sources, which irradiate material samples of interest. Under sufficiently high x-ray energy deposition, thermomechanical impulses are generated in the materials. While it is known that the characteristics of x-ray generated impulses vary as a function of incident x-ray spectra, it remains unclear how spectral assumptions and uncertainties in NIF spectral reconstructions affect our interpretation of impulsive loading. Here, in this paper, we simulate the response of a standard titanium alloy baseline sample to synthetic analytically derived and measured NIF xenon line-emission x-ray sources with a radiation hydrodynamics code. We vary the source spectral characteristics based on different source reconstruction techniques to understand the resulting variation in baseline sample response and compare the simulated response with experimental results. We find that the response is highly sensitive to assumptions made about the spectral contents and that knowledge of spectral uncertainties bounds our understanding of the resulting material response. The results of this effort help to extend our ability to use baseline material samples to extract quantitative properties from x-ray experiments on the NIF.

Alloys↗
DOE OSTI2020

Homogeneous broadening in the spectrum of a broadband erbium-doped fibre light source

Interference of wave trains separated out in an inhomogeneous emission spectrum of a superluminescent erbium-doped fibre source has been studied by scanning the wavelength of a spectrum analyser. The results demonstrate that, as the spectral resolution of the spectrum analyser decreases to a level corresponding to a sharp drop in interference fringe visibility relative to its original level, near 100 %, fringe visibility becomes dependent on the wavelength in the inhomogeneous spectrum of the source. The reason for this is that, at some characteristic resolution δλ (δλ > 0.5 nm in our case), the coherence time of wave trains is not limited by the spectral resolution of the spectrum analyser and is related to the inhomogeneous spectrum of the source, whose limiting width can be thought of as the homogeneous linewidth in the spectrum of the source. This feature can be used for evaluating homogeneous line broadening. (paper)

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Sound propagation in realistic interactive 3D scenes with parameterized sources using deep neural operators

We address the challenge of acoustic simulations in three-dimensional (3D) virtual rooms with parametric source positions, which have applications in virtual/augmented reality, game audio, and spatial computing. The wave equation can fully describe wave phenomena such as diffraction and interference. However, conventional numerical discretization methods are computationally expensive when simulating hundreds of source and receiver positions, making simulations with parametric source positions impractical. To overcome this limitation, we propose using deep operator networks to approximate linear wave-equation operators. This enables the rapid prediction of sound propagation in realistic 3D acoustic scenes with parametric source positions, achieving millisecond-scale computations. By learning a compact surrogate model, we avoid the offline calculation and storage of impulse responses for all relevant source/listener pairs. Our experiments, including various complex scene geometries, show good agreement with reference solutions, with root mean squared errors ranging from 0.02 to 0.10 Pa. Notably, our method signifies a paradigm shift as—to our knowledge—no prior machine learning approach has achieved precise predictions of complete wave fields within realistic domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The Source Permeator System and Tritium Transport in the TEX PbLi Loop

In this work, permeation is investigated for the introduction of hydrogen isotopes into PbLi for the Tritium Extraction eXperiment (TEX). TEX is a forced-convection PbLi loop under construction at Idaho National Laboratory (INL) that will test the vacuum permeator (VP) method of tritium extraction from PbLi. The source permeator (SP) delivers atomic hydrogen (H, D, and T) from a gas-phase reservoir into the PbLi via a permeable dense metal membrane. A modular and a fixed source permeator system are presented. In the modular design, PbLi flows through the inside of a tubular membrane and gas-phase hydrogen is introduced on the outside of the membrane. Atomic hydrogen permeates radially inward through the membrane into the PbLi. In the fixed design, PbLi flows into an expansion chamber with closed-ended tubular membranes inserted. Gas-phase hydrogen is introduced on the inside of the closed-ended membranes and atomic hydrogen permeates radially outward into the flowing PbLi. Hydrogen-transport models based on steady-state mass transport through PbLi and permeation through the metal membrane were developed to assess the operation of the source permeator relative to experimental variables and to allow understanding of uncertain parameter effects, such as PbLi hydrogen-transport properties and the effective hydrogen permeability of the vacuum permeator. This modeling effort considers iron as the source permeator material and vanadium as the vacuum permeator material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantifying and Visualizing Uncertainty for Source Localisation in Electrocardiographic Imaging

Electrocardiographic imaging (ECGI) presents a clinical opportunity to noninvasively understand the sources of arrhythmias for individual patients. To help increase the effectiveness of ECGI, we provide new ways to visualise associated measurement and modelling errors. Here, we study source localisation uncertainty in two steps: First, we perform Monte Carlo simulations of a simple inverse ECGI source localisation model with error sampling to understand the variations in ECGI solutions. Second, we present multiple visualisation techniques, including confidence maps, level-sets, and topology-based visualisations, to better understand uncertainty in source localization. Our approach offers a new way to study uncertainty in the ECGI pipeline.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

A Proof of the Asymptotic Variance of Path Length Estimators for Single-Collision Monte Carlo Source Iteration in the Thick Diffusion Limit

Here, we prove a theorem relating the variance of path length estimators for single-collision Monte Carlo source iteration to a parameter that becomes infinitesimally small in an important physical regime arising in radiative transfer. In our usage, “single-collision Monte Carlo source iteration” refers to Monte Carlo Boltzmann transport methods in which each Monte Carlo particle history includes no more than a single collision, and the physics of multiple scattering is modeled by lagging the scattering source term and iterating until this term converges. Our theorem can be used to construct variance reduction techniques which improve the order of the estimator variance. This enables calculations that would otherwise require impractically large sample sizes to achieve practical estimator uncertainties. We believe this is the first postulation of a theorem relating estimator variance to a limiting case parameter for single-collision Monte Carlo source iteration, and the first proof of such a theorem. We illustrate the theorem’s value with an example in which the authors of a transport method used the theorem to design a variance reduction technique that improved the uncertainty of their solution by a factor of about 500 for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material.

Mathematics and Computing↗
DOE OSTI2021

The Highest Energy HAWC Sources are Likely Leptonic and Powered by Pulsars

The HAWC Collaboration has observed gamma rays at energies above 56 TeV from a collection of nine sources. It has been suggested that this emission could be hadronic in nature, requiring that these systems accelerate cosmic-ray protons or nuclei up to PeV-scale energies. In this paper, we instead show that the spectra of these objects favor a leptonic (inverse Compton) origin for their emission. More specifically, the gamma-ray emission from these objects can be straightforwardly accommodated within a model in which $\sim \mathcal{O}(10\%)$ of the host pulsar's spindown power is transferred into the acceleration of electrons and positrons with a power-law spectrum that extends to several hundred TeV or higher. The spectral break that is observed among these sources is naturally explained within the context of this simple model, and occurs at the energy where the timescale for energy losses matches the age of the pulsar. In contrast, this spectral feature cannot be straightforwardly accommodated in hadronic scenarios. Furthermore, hadronic models predict that these sources should produce more emission at GeV-scale energies than is observed. In light of these considerations, we conclude that HAWC's highest energy sources should be interpreted as TeV halos or pulsar wind nebulae, which produce their emission through inverse Compton scattering, and are powered by the rotational kinetic energy of their host pulsar.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗