Engineering topics

Pinto, Philip A.

Publications and source records attributed to Pinto, Philip A..

10 records

DOE OSTI2021

The ABACUS cosmological N -body code

We present abacus, a fast and accurate cosmological N-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near- and far-field forces, reducing the former to direct 1/r 2 summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of 10 -5 . We express the simulation time-step as an event-driven ‘pipeline’, incorporating asynchronous events such as completion of co-processor work, input/output, and network communication. abacus has been used to produce the largest suite of N-body simulations to date, the abacussummit suite of 60 trillion particles, incorporating on-the-fly halo finding. abacus enables the production of mock catalogues of the volume and resolution required by the coming generation of cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1994

SN 1993J: A Type IIb supernova

The evolution of the bright Type II supernova discovered last year in M81, SN 1993J, is consistent with that expected for the explosion of a star which on the main sequence had a mass of 13-16 Solar Mass but which, owing to mass exchange with a binary companion (a intially approximately 3-5 AU, depending upon the actual presupernova radius and the masses of the two stars) lost almost all of its hydrogen-rich envelope during late helium burning. At the time of explosion, the helium core mass was 4.0 +/- 0.5 Solar Mass and the hydrogen envelope, 0.20 +/- 0.05 Solar Mass. The envelope was helium and nitrogen-rich (carbon-deficient) and the radius of the star, 4 +/- 1 x 10(exp 13) cm. The luminosity of the presupernova star was 3 + 1 x 10(exp 38) ergs/s, with the companion star contributing an additional approximately 10(exp 38) ergs/s. The star may have been a pulsating variable at the time of the explosion. For an explosion energy near 10(exp 51) ergs (KE at infinity) and an assumed distance of 3.3 Mpc, a mass of Ni-56 in the range 0.07 +/- 0.01 Solar Mass was produced and ejected. This presciption gives a light curve which compares favorably with the bolomatric observations. Color photometry is more restrictive and requires a model in which the hydrogen-envelope mass is low and the mixing of hydrogen inward has been small, but in which appreciable Ni-56 has been mixed outward into the helium and heavy-element core. It is possible to obtain good agreement with B and V light curves during the first 50 days, but later photometry, especially in bands other than B and V, will require a non-local thermo-dynamic equilibrium (LTE) spectral calculation for comparison. Based upon our model, we predict a flux of approximately 10(exp -5)(3.3 Mpc/D)(exp 2) photons/sq cm/s in the 847 keV line of CO-56 at peak during 1993 August. It may be easier to detect the Computonized continuum which peaks at a few times 10(exp -4) photons /s/sq cm/MeV at 40 keV a few months after the explosion (though neither of these signals were, or should have been, detected by the Compton Gamma-Ray observatory (CGRO). The presupernova star was filing its Roche lobe at the time of the explosion and thus its envelope was highly deformed (about 3:2). The companion star is presently embedded in the supernova, but should becopme visable at age 3 yr (perhaps earlier in the ultraviolet) when the supernova has faded below 10(exp 38) ergs/s. Indeed, if 'kicks' have not played an important role, it is still bound to the neutron star.

Woosley, S. E.↗

NASA NTRS1994

Evidence for Ni-56 yields Co-56 yields Fe-56 decay in type Ia supernovae

In the prevailing picture of Type Ia supernovae (SN Ia), their explosive burning produces Ni-56, and the radioactive decay chain Ni-56 yields Co-56 yields Fe-56 powers the subsequent emission. We test a central feature of this theory by measuring the relative strengths of a (Co III) emission feature near 5900 A and a (Fe III) emission feature near 4700 A. We measure 38 spectra from 13 SN Ia ranging from 48 to 310 days after maximum light. When we compare the observations with a simple multilevel calculation, we find that the observed Fe/Co flux ratio evolves as expected when the Fe-56/Co-56 abundance ratio follows from Ni-56 yields Co-56 yields Fe-56 decay. From this agreement, we conclude that the cobalt and iron atoms we observe through SN Ia emission lines are produced by the radioactive decay of Ni-56, just as predicted by a wide range of models for SN Ia explosions.

Kuchner, Marc J.↗

NASA NTRS1993

Airborne spectrophotometry of SN 1987A from 1.7 to 12.6 microns - Time history of the dust continuum and line emission

Spectrophotometric observations of SN 1987A from the Kuiper Airborne Observatory are presented for five epochs at 60, 260, 415, 615, and 775 days after the explosion. The low-resolution (lambda/Delta lambda = 50-100) spectra of SN 1987A are combined with data from other wavelengths to model the continuum, subtract the continuum from the spectra to determine line strengths and reveal molecular bands, separate the atomic continuum radiation from the dust continuum, and derive constraints on the grain temperatures and optical depths. A scenario for the evolution of SN 1987A and that of the ejecta from which it arises is obtained on the basis of the analysis of the continuum emission.

Wooden, Diane H.↗

NASA NTRS1993

Spectrum formation in supernovae - Numerical techniques

The study combines several novel techniques for spectrum simulation in the Eddington computer program which solves the comoving frame equation of transfer coupled with the statistical and radiative equilibrium equations. One of these is a generalization of the accelerated lambda iteration (ALI) scheme to include an approximate frequency-derivative operator. This greatly enhances the convergence rate of ALI in optically thick, high-velocity shear flows. Another is a partial linearization technique which is capable of efficiently solving a very large number of rate equations on a moderately sized computer. An expansion opacity and emissivity approximation is derived which makes it possible to determine the effect on the transfer and statistical equilibrium of a very large number of lines not explicitly represented in the frequency grid and additionally to treat line-blanketing from species not explicitly included in the rate equations. The utility of these techniques is illustrated with models of two supernovae.

Eastman, Ronald G.↗

NASA NTRS1992

A distance-independent calibration of the luminosity of type Ia supernovae and the Hubble constant

The absolute magnitude of SNe Ia at maximum is calibrated here using radioactive decay models for the light curve and a minimum of assumptions. The absolute magnitude parameter space is studied using explosion models and a range of rise times, and absolute B magnitudes at maximum are used to derive a range of the H0 and the distance to the Virgo Cluster from SNe Ia. Rigorous limits for H0 of 45 and 105 km/s/Mpc are derived.

Leibundgut, Bruno↗

NASA NTRS1991

Spectra of two very old supernovae - SN 1986J and SN 1980K

Spectra are presented for SN 1986J in NGC 891, observed 4 and 7 years after the explosion, and SN 1980K, observed 9 years after maximum light. A narrow-line and a broad-line component are noted in SN 1986J; these are respectively attributed to a circumstellar shell and the actual stellar interior (without hydrogen). SN 1980K continues to emit in very broad lines. Comparison with earlier observations suggests that this SN halted its exponential flux decline in the early 1980s, and is currently emitting at a constant rate. In both SNs studies, the identity of the energy source which sustains emission remains uncertain.

Leibundgut, Bruno↗

NASA NTRS1989

Hard emission at late times from SN 1987A

A model for the explosion that has been successful in predicting and explaining the evolution of SN 1987A during its first one and one-half years is used to calculate the future photometric evolution of the supernova in the UV-optical-IR, X-ray, and gamma-ray bands using Monte Carlo techniques. Special attention is given to the contribution from radioactive isotopes other than Co-56, notably Co-57, Ti-44, and Na-22, and to the possible appearance, probably within the next year, of X-rays from an accreting magnetic neutron star at the center of the supernova. The signature of a radio pulsar like the Crab is also considered. In both cases the time history of the bolometric light curve may be a more sensitive diagnostic than the existence of either pulses or hard emission, but will become confused in the near future by contributions from the rare radioactivities. A source as faint as the pulsed emission of the Crab pulsar will never be discernable in the bolometric light curve. Pulsed hard emission is best sought near 30 keV near the end of 1989 (and thereafter), although more sensitive instrumentation than hitherto employed may be necessary.

Woosley, S. E.↗

NASA NTRS1988

SN 1987A - Out on the tail

The expected light curve of Supernova 1987A is calculated for the next 10 yr based solely upon energy input from radioactive decay of Co-56 and Co-57. During the third year, energy input from Co-57 begins to dominate. How the bolometric (ultraviolet-optical-infrared) light curve deviates from that given by 100 percent deposition of the gamma-rays is a measure of the mixing that has gone on in the supernova. Current measurements favor a mixed model. The radioactive light curve decreases beneath 10 to the 39th ergs/s at an age of about 900 days, by which time the neutron star remnant should have made its presence known.

Pinto, Philip A.↗

NASA NTRS1988

The theory of gamma-ray emergence in supernova 1987A

It was anticipated that the decay of Co-56 to Fe-56 in SN1987A would give rise to detectable gamma-ray emission at 847 and 1238 keV with a peak flux of 0.001 photon/sq cm/s about one year after the explosion. Both these lines were detected in August, with a strength within a factor of two of 0.001 photon/sq cm/s, but about six months earlier than predicted. It is shown here that the early emergence of gamma-rays can be accounted for in a 'mixed' model in which an approximately isotropic process destroys chemical segregation with respect to radial mass coordinate and velocity.

Pinto, Philip A.↗