Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FIT”

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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

FERMI-LAT OBSERVATIONS AND A BROADBAND STUDY OF SUPERNOVA REMNANT CTB 109

CTB 109 (G109.1-1.0) is a Galactic supernova remnant (SNR) with a hemispherical shell morphology in X-rays and in the radio band. In this work, we report the detection of {gamma}-ray emission coincident with CTB 109, using 37 months of data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. We study the broadband characteristics of the remnant using a model that includes hydrodynamics, efficient cosmic-ray (CR) acceleration, nonthermal emission, and a self-consistent calculation of the X-ray thermal emission. We find that the observations can be successfully fit with two distinct parameter sets, one where the {gamma}-ray emission is produced primarily by leptons accelerated at the SNR forward shock and the other where {gamma}-rays produced by forward shock accelerated CR ions dominate the high-energy emission. Consideration of thermal X-ray emission introduces a novel element to the broadband fitting process, and while it does not rule out either the leptonic or the hadronic scenarios, it constrains the parameter sets required by the model to fit the observations. Moreover, the model that best fits the thermal and nonthermal emission observations is an intermediate case, where both radiation from accelerated electrons and hadrons contribute almost equally to the {gamma}-ray flux observed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE {gamma}-RAY SPECTRUM OF GEMINGA AND THE INVERSE COMPTON MODEL OF PULSAR HIGH-ENERGY EMISSION

We reanalyze the Fermi spectra of the Geminga and Vela pulsars. We find that the spectrum of Geminga above the break is well approximated by a simple power law without the exponential cutoff, making Geminga's spectrum similar to that of Crab. Vela's broadband {gamma}-ray spectrum is equally well fit with both the exponential cutoff and the double power-law shapes. In the broadband double power-law fits, for a typical Fermi spectrum of a bright {gamma}-ray pulsar, most of the errors accumulate due to the arbitrary parameterization of the spectral roll-off. In addition, a power law with an exponential cutoff gives an acceptable fit for the underlying double power-law spectrum for a very broad range of parameters, making such fitting procedures insensitive to the underlying Fermi photon spectrum. Our results have important implications for the mechanism of pulsar high-energy emission. A number of observed properties of {gamma}-ray pulsars-i.e., the broken power-law spectra without exponential cutoffs and stretching in the case of Crab beyond the maximal curvature limit, spectral breaks close to or exceeding the maximal breaks due to curvature emission, patterns of the relative intensities of the leading and trailing pulses in the Crab repeated in the X-ray and {gamma}-ray regions, presence of profile peaks at lower energies aligned with {gamma}-ray peaks-all point to the inverse Compton origin of the high-energy emission from majority of pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE MILKY WAY TOMOGRAPHY WITH SLOAN DIGITAL SKY SURVEY. IV. DISSECTING DUST

We use Sloan Digital Sky Survey (SDSS) photometry of 73 million stars to simultaneously constrain best-fit main-sequence stellar spectral energy distribution (SED) and amount of dust extinction along the line of sight toward each star. Using a subsample of 23 million stars with Two Micron All Sky Survey (2MASS) photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, R{sub V} , can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes (|b| < 30 Degree-Sign ). Our results are in good agreement with the extinction normalization given by the Schlegel et al. (SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with recent study by Schlafly et al. The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses disfavor the reddening law of O'Donnell, but support the models by Fitzpatrick and Cardelli et al. For the latter, we find a ratio of the total to selective absorption to be R{sub V} = 3.0 {+-} 0.1(random){+-}0.1 (systematic) over most of the high-latitude sky. At low Galactic latitudes (|b| < 5 Degree-Sign ), we demonstrate that the SFD map cannot be reliably used to correct for extinction because most stars are embedded in dust, rather than behind it, as is the case at high Galactic latitudes. We analyze three-dimensional maps of the best-fit R{sub V} and find that R{sub V} = 3.1 cannot be ruled out in any of the 10 SEGUE stripes at a precision level of {approx}0.1-0.2. Our best estimate for the intrinsic scatter of R{sub V} in the regions probed by SEGUE stripes is {approx}0.2. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed 'dusty parallax relation', which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format. These data can be used for studies of stellar number density distribution, the distribution of dust properties, for selecting sources whose SED differs from SEDs for high-latitude main-sequence stars, and for estimating distances to dust clouds and, in turn, to molecular gas clouds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

REGULATION OF BLACK HOLE WINDS AND JETS ACROSS THE MASS SCALE

We present a study of the mechanical power generated by both winds and jets across the black hole mass scale. We begin with the study of ionized X-ray winds and present a uniform analysis using Chandra grating spectra. The high-quality grating spectra facilitate the characterization of the outflow velocity, ionization, and column density of the absorbing gas. We find that the kinetic power of the winds, derived from these observed quantities, scales with increasing bolometric luminosity as log (L {sub wind,42}/C{sub v} ) = (1.58 {+-} 0.07)log (L {sub Bol,42}) - (3.19 {+-} 0.19). This suggests that supermassive black holes may be more efficient than stellar-mass black holes in launching winds, per unit filling factor, C{sub v} . If the black hole binary (BHB) and active galactic nucleus (AGN) samples are fit individually, the slopes flatten to {alpha}{sup BHB} = 0.91 {+-} 0.31 and {alpha}{sup AGN} = 0.63 {+-} 0.30 (formally consistent within errors). The broad fit and individual fits both characterize the data fairly well, and the possibility of common slopes may point to common driving mechanisms across the mass scale. For comparison, we examine jet production, estimating jet power based on the energy required to inflate local bubbles. The jet relation is log (L {sub Jet,42}) = (1.18 {+-} 0.24)log (L {sub Bondi,42}) - (0.96 {+-} 0.43). The energetics of the bubble associated with Cygnus X-1 are particularly difficult to determine, and the bubble could be a background supernova remnant. If we exclude Cygnus X-1 from our fits, then the jets follow a relation consistent with the winds, but with a higher intercept, log (L {sub Jet,42}) = (1.34 {+-} 0.50)log (L {sub Bondi,42}) - (0.80 {+-} 0.82). The formal consistency in the wind and jet scaling relations, when assuming that L {sub Bol} and L {sub Bondi} are both proxies for mass accretion rate, suggests that a common launching mechanism may drive both flows; magnetic processes, such as magnetohydrodynamics and magnetocentrifugal forces, are viable possibilities. We also examine winds that are moving at especially high velocities, v > 0.01c. These ultra-fast outflows tend to resemble the jets more than the winds in terms of outflow power, indicating that we may be observing a regime in which winds become jets. A transition at approximately L {sub Bol} Almost-Equal-To 10{sup -2} L {sub Edd} is apparent when outflow power is plotted versus Eddington fraction. At low Eddington fractions, the jet power is dominant, and at high Eddington fractions, the wind power is dominant. This study allows for the total power from black hole accretion, both mechanical and radiative, to be characterized in a simple manner and suggests possible connections between winds and jets. X-ray wind data and jet cavity data will enable stronger tests.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

TIMING AND SPECTRAL PROPERTIES OF Be/X-RAY PULSAR EXO 2030+375 DURING A TYPE I OUTBURST

We present results from a study of broadband timing and spectral properties of EXO 2030+375 using a Suzaku observation. Pulsations with a period of 41.41 s and strong energy-dependent pulse profiles were clearly detected up to 100 keV. Narrow dips are seen in the profiles up to {approx}70 keV. The presence of prominent dips at several phases in the profiles up to such high energy ranges was not seen before. At higher energies, these dips gradually disappeared and the profile appeared to be single-peaked. The 1.0-200.0 keV broadband spectrum is found to be well described by a partial covering high-energy cutoff power-law model. Several low-energy emission lines are also detected in the pulsar spectrum. We fitted the spectrum using neutral as well as partially ionized absorbers along with the above continuum model yielding similar parameter values. The partial covering with a partially ionized absorber resulted in a marginally better fit. The spectral fitting did not require any cyclotron feature in the best-fit model. To investigate the changes in spectral parameters at dips, we carried out pulse-phase-resolved spectroscopy. During the dips, the value of the additional column density was estimated to be high compared to other pulse phases. While using a partially ionized absorber, the value of the ionization parameter is also higher at the dips. This may be the reason for the presence of dips up to higher energies. No other spectral parameters show any systematic variation with pulse phases of the pulsar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

MODELING THE NUCLEAR INFRARED SPECTRAL ENERGY DISTRIBUTION OF TYPE II ACTIVE GALACTIC NUCLEI

We present results from model fitting to the spectral energy distribution (SED) of a homogeneous sample of Seyfert II galaxies drawn from the 12 {mu}m Galaxy Sample. Imaging and nuclear flux measurements are presented in an accompanying paper. Here we add Spitzer/IRS observations to further constrain the SEDs after careful subtraction of a starburst component. We use the library of CLUMPY torus models from Nenkova et al. and also test the two-phase models recently produced by Stalevski et al. We find that photometric and spectroscopic observations in the mid-IR ({lambda} {approx}> 5 {mu}m) are crucial to properly constrain the best-fit torus models. About half of our sources show clear near-IR excess of their SEDs above the best-fit models. This problem can be less severe when using the Stalevski et al. models. The nature of this emission is not clear since best-fitted blackbody temperatures are very high ({approx}1700-2500 K) and the Type II classification of our sources would correspond to a small probability to peer directly into the hottest regions of the torus. Crucially, the derived torus parameters are quite robust when using CLUMPY models, independently of whether or not the sources require an additional blackbody component. Our findings suggest that tori are characterized by N{sub 0}{approx}>5, {sigma} {approx}> 40, {tau} {approx}< 25, Angle i {approx}> 40 Degree-Sign , Y {approx}< 50, and A {sup los} {sub v} {approx} 100-300, where N{sub 0} is the number of clouds in the equatorial plane of the torus, {sigma} is the characteristic opening angle of the cloud distribution, {tau} is the opacity of a single cloud, Angle i is the line-of-sight orientation of the torus, Y is the ratio of the inner to the outer radii, and A {sup los} {sub v} is the total opacity along the line of sight. From these, we can determine typical torus sizes and masses of 0.1-5.0 pc and 10{sup 4}-10{sup 6} M {sub Sun }, respectively. We find tentative evidence that those nuclei with detected hidden broad-line regions are characterized by lower levels of extinction than those without one. Finally, we find no correlation between the torus properties and the presence of circumnuclear or more global star formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

REVISITING THE SCALING RELATIONS OF BLACK HOLE MASSES AND HOST GALAXY PROPERTIES

New kinematic data and modeling efforts in the past few years have substantially expanded and revised dynamical measurements of black hole masses (M {sub .}) at the centers of nearby galaxies. Here we compile an updated sample of 72 black holes and their host galaxies, and present revised scaling relations between M {sub .} and stellar velocity dispersion ({sigma}), V-band luminosity (L), and bulge stellar mass (M {sub bulge}), for different galaxy subsamples. Our best-fitting power-law relations for the full galaxy sample are log{sub 10}(M {sub .}) = 8.32 + 5.64log{sub 10}({sigma}/200 km s{sup -1}), log{sub 10}(M {sub .}) = 9.23 + 1.11log{sub 10}(L/10{sup 11} L {sub Sun }), and log{sub 10}(M {sub .}) = 8.46 + 1.05log{sub 10}(M {sub bulge}/10{sup 11} M {sub Sun }). A log-quadratic fit to the M {sub .}-{sigma} relation with an additional term of {beta}{sub 2} [log{sub 10}({sigma}/200 km s{sup -1})]{sup 2} gives {beta}{sub 2} = 1.68 {+-} 1.82 and does not decrease the intrinsic scatter in M {sub .}. Including 92 additional upper limits on M {sub .} does not change the slope of the M {sub .}-{sigma} relation. When the early- and late-type galaxies are fit separately, we obtain similar slopes of 5.20 and 5.06 for the M {sub .}-{sigma} relation but significantly different intercepts-M {sub .} in early-type galaxies are about two times higher than in late types at a given sigma. Within early-type galaxies, our fits to M {sub .}({sigma}) give M {sub .} that is about two times higher in galaxies with central core profiles than those with central power-law profiles. Our M {sub .}-L and M {sub .}-M {sub bulge} relations for early-type galaxies are similar to those from earlier compilations, and core and power-law galaxies yield similar L- and M {sub bulge}-based predictions for M {sub .}. When the conventional quadrature method is used to determine the intrinsic scatter in M {sub .}, our data set shows weak evidence for increased scatter at M {sub bulge} < 10{sup 11} M {sub Sun} or L{sub V} < 10{sup 10.3} L {sub Sun }, while the scatter stays constant for 10{sup 11} < M {sub bulge} < 10{sup 12.3} M {sub Sun} and 10{sup 10.3} < L{sub V} < 10{sup 11.5} L {sub Sun }. A Bayesian analysis indicates that a larger sample of M {sub .} measurements would be needed to detect any statistically significant trend in the scatter with galaxy properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE TEMPERATURES OF RED SUPERGIANTS

We present a re-appraisal of the temperatures of red supergiants (RSGs) using their optical and near-infrared spectral energy distributions (SEDs). We have obtained data of a sample of RSGs in the Magellanic Clouds using VLT+XSHOOTER, and we fit MARCS model atmospheres to different regions of the spectra, deriving effective temperatures for each star from (1) the TiO bands, (2) line-free continuum regions of the SEDs, and (3) the integrated fluxes. We show that the temperatures derived from fits to the TiO bands are systematically lower than the other two methods by several hundred kelvin. The TiO fits also dramatically overpredict the flux in the near-IR, and imply extinctions which are anomalously low compared to neighboring stars. In contrast, the SED temperatures provide good fits to the fluxes at all wavelengths other than the TiO bands, are in agreement with the temperatures from the flux integration method, and imply extinctions consistent with nearby stars. After considering a number of ways to reconcile this discrepancy, we conclude that three-dimensional effects (i.e., granulation) are the most likely cause, as they affect the temperature structure in the upper layers where the TiO lines form. The continuum, however, which forms at much deeper layers, is apparently more robust to such effects. We therefore conclude that RSG temperatures are much warmer than previously thought. We discuss the implications of this result for stellar evolution and supernova progenitors, and provide relations to determine the bolometric luminosities of RSGs from single-band photometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

GAMMA-RAY BURST AFTERGLOW LIGHT CURVES FROM A LORENTZ-BOOSTED SIMULATION FRAME AND THE SHAPE OF THE JET BREAK

The early stages of decelerating gamma-ray burst (GRB) afterglow jets have been notoriously difficult to resolve numerically using two-dimensional hydrodynamical simulations even at very high resolution, due to the extreme thinness of the blast wave and high outflow Lorentz factors. However, these resolution issues can be avoided by performing the simulations in a boosted frame, which makes it possible to calculate afterglow light curves from numerically computed flows in sufficient detail to accurately quantify the shape of the jet break and the post-break steepening of the light curve. Here, we study afterglow jet breaks for jets with opening angles of 0.05, 0.1, and 0.2 radians decelerating in a surrounding medium of constant density, observed at various angles ranging from on-axis to the edge of the jet. A single set of scale-invariant functions describing the time evolution of afterglow synchrotron spectral break frequencies and peak flux, depending only on jet opening angle and observer angle, are all that is needed to reconstruct light curves for arbitrary explosion energy, circumburst density and synchrotron particle distribution power law slope p. These functions are presented in the paper. Their time evolutions change directly following the jet break, although an earlier reported temporary post-break steepening of the cooling break is found to have been resolution-induced. We compare synthetic light curves to fit functions using sharp power law breaks as well as smooth power law transitions. We confirm our earlier finding that the measured jet break time is very sensitive to the angle of the observer and can be postponed significantly. We find that the difference in temporal indices across the jet break is larger than theoretically anticipated and is about -(0.5 + 0.5p) below the cooling break and about -(0.25 + 0.5p) above the cooling break, both leading to post-break slopes of roughly about 0.25 - 1.3p, although different observer angles, jet opening angles and heuristic descriptions of the break introduce a wide range of temporal indices. Nevertheless, the post-break slope from our constant density interstellar medium simulations is sufficiently steep to be hard to reconcile with post-break slopes measured for the Swift sample, suggesting that Swift GRBs mostly do not explode in a homogeneous medium or that the jet breaks are hidden from view by additional physics such as prolonged energy injection or viewing angle effects. A comparison between different smooth power law fit functions shows that although smooth power law transitions of the type introduced by Harrison et al. often provide better fits, smooth power law transitions of the type introduced by Beuermann et al. or even sharp power law fits are easier to interpret in terms of the underlying model. Light curves and spectral break and peak flux evolution functions will be made publicly available online at http://cosmo.nyu.edu/afterglowlibrary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

EVIDENCE FOR A PHOTOSPHERIC COMPONENT IN THE PROMPT EMISSION OF THE SHORT GRB 120323A AND ITS EFFECTS ON THE GRB HARDNESS-LUMINOSITY RELATION

The short GRB 120323A had the highest flux ever detected with the Gamma-Ray Burst Monitor on board the Fermi Gamma-Ray Space Telescope. Here we study its remarkable spectral properties and their evolution using two spectral models: (1) a single emission component scenario, where the spectrum is modeled by the empirical Band function (a broken power law), and (2) a two-component scenario, where thermal (a Planck-like function) emission is observed simultaneously with a non-thermal component (a Band function). We find that the latter model fits the integrated burst spectrum significantly better than the former, and that their respective spectral parameters are dramatically different: when fit with a Band function only, the E{sub peak} of the event is unusually soft for a short gamma-ray burst (GRB; 70 keV compared to an average of 300 keV), while adding a thermal component leads to more typical short GRB values (E{sub peak} {approx} 300 keV). Our time-resolved spectral analysis produces similar results. We argue here that the two-component model is the preferred interpretation for GRB 120323A based on (1) the values and evolution of the Band function parameters of the two component scenario, which are more typical for a short GRB, and (2) the appearance in the data of a significant hardness-intensity correlation, commonly found in GRBs, when we employee two-component model fits; the correlation is non-existent in the Band-only fits. GRB 110721A, a long burst with an intense photospheric emission, exhibits the exact same behavior. We conclude that GRB 120323A has a strong photospheric emission contribution, observed for the first time in a short GRB. Magnetic dissipation models are difficult to reconcile with these results, which instead favor photospheric thermal emission and fast cooling synchrotron radiation from internal shocks. Finally, we derive a possibly universal hardness-luminosity relation in the source frame using a larger set of GRBs (L{sub i}{sup Band} = (1.59 {+-} 0.84) x 10{sup 50} (E{sub peak,i}){sup 1.33{+-}0.07} erg s{sup -1}), which could be used as a possible redshift estimator for cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE GAMMA-RAY BURST HUBBLE DIAGRAM AND ITS IMPLICATIONS FOR COSMOLOGY

In this paper, we continue to build support for the proposal to use gamma-ray bursts (GRBs) as standard candles in constructing the Hubble diagram at redshifts beyond the current reach of Type Ia supernova observations. We confirm that correlations among certain spectral and light-curve features can indeed be used as luminosity indicators, and demonstrate from the most up-to-date GRB sample appropriate for this work that the {Lambda}CDM model optimized with these data is characterized by parameter values consistent with those in the concordance model. Specifically, we find that ({Omega}{sub m},{Omega}{sub {Lambda}}){approx}(0.25{sub -0.06}{sup +0.05}, 0.75{sub -0.05}{sup +0.06}), which are consistent, to within 1{sigma}, with (0.29, 0.71) obtained from the 9 yr Wilkinson Microwave Anisotropy Probe data. We also carry out a comparative analysis between {Lambda}CDM and the R{sub h} = ct universe and find that the optimal {Lambda}CDM model fits the GRB Hubble diagram with a reduced {chi}{sup 2}{sub dof}{approx}2.26, whereas the fit using R{sub h} = ct results in a {chi}{sup 2}{sub dof}{approx}2.14. In both cases, about 20% of the events lie at least 2{sigma} away from the best-fit curves, suggesting that either some contamination by non-standard GRB luminosities is unavoidable or that the errors and intrinsic scatter associated with the data are being underestimated. With these optimized fits, we use three statistical tools-the Akaike information criterion, the Kullback information criterion, and the Bayes information criterion-to show that, based on the GRB Hubble diagram, the likelihood of R{sub h} = ct being closer to the correct model is {approx}85%-96%, compared to {approx}4%-15% for {Lambda}CDM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

COMPLETE INFRARED SPECTRAL ENERGY DISTRIBUTIONS OF MILLIMETER DETECTED QUASARS AT z > 5

We present Herschel far-infrared (FIR) photometry of 11 quasars at redshift z > 5 that have previously been detected at 1.2 mm. We perform full spectral energy distribution (SED) fits over the wavelength range {lambda}{sub rest} {approx} 0.1-400 {mu}m for those objects with good Herschel detections. These fits reveal the need for an additional FIR component besides the emission from a dusty active galactic nucleus (AGN)-powered torus. This additional FIR component has temperatures of T{sub FIR} {approx} 40-60 K with luminosities of L{sub 8-1000{mu}m} {approx} 10{sup 13} L{sub Sun} (accounting for 25%-60% of the bolometric FIR luminosity). If the FIR dust emission is due to star formation it would suggest star formation rates in excess of 1000 solar masses per year. We show that at long wavelengths ({lambda}{sub rest} {approx}> 50 {mu}m) the contribution of the AGN-powered torus emission is negligible. This explains how previous FIR studies of high-redshift quasars that relied on single-component fits to (ground-based) observations at {lambda}{sub obs} {approx}> 350 {mu}m reached T{sub FIR} and L{sub FIR} values similar to our complete SED fits. Stacking the Herschel data of four individually undetected sources reveals a significant average signal in the PACS bands but not in SPIRE. The average SED of sources with individual Herschel detections shows a striking surplus in near- and mid-infrared (MIR) emission when compared to common AGN templates. The comparison between two average SEDs (sources with and without individual Herschel detections) matched in the UV/optical indicates that for these objects the strength of the MIR emission may correlate with the strength of the FIR emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

EVIDENCE FOR ACCRETION RATE CHANGE DURING TYPE I X-RAY BURSTS

The standard approach for time-resolved X-ray spectral analysis of thermonuclear bursts involves subtraction of the pre-burst emission as background. This approach implicitly assumes that the persistent flux remains constant throughout the burst. We reanalyzed 332 photospheric radius expansion bursts observed from 40 sources by the Rossi X-Ray Timing Explorer, introducing a multiplicative factor f{sub a} to the persistent emission contribution in our spectral fits. We found that for the majority of spectra the best-fit value of f{sub a} is significantly greater than 1, suggesting that the persistent emission typically increases during a burst. Elevated f{sub a} values were not found solely during the radius expansion interval of the burst, but were also measured in the cooling tail. The modified model results in a lower average value of the {chi}{sup 2} fit statistic, indicating superior spectral fits, but not yet to the level of formal statistical consistency for all the spectra. We interpret the elevated f{sub a} values as an increase of the mass accretion rate onto the neutron star during the burst, likely arising from the effects of Poynting-Robertson drag on the disk material. We measured an inverse correlation of f{sub a} with the persistent flux, consistent with theoretical models of the disk response. We suggest that this modified approach may provide more accurate burst spectral parameters, as well as offering a probe of the accretion disk structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

MEASUREMENT OF SPIN-ORBIT MISALIGNMENT AND NODAL PRECESSION FOR THE PLANET AROUND PRE-MAIN-SEQUENCE STAR PTFO 8-8695 FROM GRAVITY DARKENING

PTFO 8-8695b represents the first transiting exoplanet candidate orbiting a pre-main-sequence star (van Eyken et al. 2012, ApJ, 755, 42). We find that the unusual lightcurve shapes of PTFO 8-8695 can be explained by transits of a planet across an oblate, gravity-darkened stellar disk. We develop a theoretical framework for understanding precession of a planetary orbit's ascending node for the case when the stellar rotational angular momentum and the planetary orbital angular momentum are comparable in magnitude. We then implement those ideas to simultaneously and self-consistently fit two separate lightcurves observed in 2009 December and 2010 December. Our two self-consistent fits yield M{sub p} = 3.0 M{sub Jup} and M{sub p} = 3.6 M{sub Jup} for assumed stellar masses of M{sub *} = 0.34 M{sub Sun} and M{sub *} = 0.44 M{sub Sun} respectively. The two fits have precession periods of 293 days and 581 days. These mass determinations (consistent with previous upper limits) along with the strength of the gravity-darkened precessing model together validate PTFO 8-8695b as just the second hot Jupiter known to orbit an M-dwarf. Our fits show a high degree of spin-orbit misalignment in the PTFO 8-8695 system: 69 Degree-Sign {+-} 2 Degree-Sign or 73. Degree-Sign 1 {+-} 0. Degree-Sign 5, in the two cases. The large misalignment is consistent with the hypothesis that planets become hot Jupiters with random orbital plane alignments early in a system's lifetime. We predict that as a result of the highly misaligned, precessing system, the transits should disappear for months at a time over the course of the system's precession period. The precessing, gravity-darkened model also predicts other observable effects: changing orbit inclination that could be detected by radial velocity observations, changing stellar inclination that would manifest as varying vsin i, changing projected spin-orbit alignment that could be seen by the Rossiter-McLaughlin effect, changing transit shapes over the course of the precession, and differing lightcurves as a function of wavelength. Our measured planet radii of 1.64 R{sub Jup} and 1.68 R{sub Jup} in each case are consistent with a young, hydrogen-dominated planet that results from a ''hot-start'' formation mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

REAL OR INTERLOPER? THE REDSHIFT LIKELIHOODS OF z > 8 GALAXIES IN THE HUDF12

In the absence of spectra, the technique of fitting model galaxy template spectra to observed photometric fluxes has become the workhorse method for determining the redshifts and other properties for high-z galaxy candidates. In this paper, we present an analysis of the most recent and possibly most distant galaxies (z {approx} 8-12) discovered in the Hubble Ultra Deep Field (HUDF) using a more robust method of redshift estimation based on Markov Chain Monte Carlo (MCMC) fitting, in contrast to the ''best fit'' models obtained using simpler {chi}{sup 2} minimization techniques. The advantage of MCMC fitting is the ability to accurately estimate the probability density function of the redshift for each object as well as any input model parameters. This makes it possible to derive accurate, credible intervals by properly marginalizing over all other input model parameters. We apply our method to 13 recently identified sources in the HUDF and show that, despite claims based on {chi}{sup 2} minimization, none of these sources can be securely ruled out as low redshift interlopers (z < 4) due to the low signal-to-noise of currently available observations. There is an average probability of 21% that these sources are low redshift interlopers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE AGES OF 55 GLOBULAR CLUSTERS AS DETERMINED USING AN IMPROVED ΔV{sup HB}{sub TO} METHOD ALONG WITH COLOR-MAGNITUDE DIAGRAM CONSTRAINTS, AND THEIR IMPLICATIONS FOR BROADER ISSUES

Ages have been derived for 55 globular clusters (GCs) for which Hubble Space Telescope Advanced Camera for Surveys photometry is publicly available. For most of them, the assumed distances are based on fits of theoretical zero-age horizontal-branch (ZAHB) loci to the lower bound of the observed distributions of HB stars, assuming reddenings from empirical dust maps and metallicities from the latest spectroscopic analyses. The age of the isochrone that provides the best fit to the stars in the vicinity of the turnoff (TO) is taken to be the best estimate of the cluster age. The morphology of isochrones between the TO and the beginning part of the subgiant branch (SGB) is shown to be nearly independent of age and chemical abundances. For well-defined color-magnitude diagrams (CMDs), the error bar arising just from the 'fitting' of ZAHBs and isochrones is ≈ ± 0.25 Gyr, while that associated with distance and chemical abundance uncertainties is ∼ ± 1.5-2 Gyr. The oldest GCs in our sample are predicted to have ages of ≈13.0 Gyr (subject to the aforementioned uncertainties). However, the main focus of this investigation is on relative GC ages. In conflict with recent findings based on the relative main-sequence fitting method, which have been studied in some detail and reconciled with our results, ages are found to vary from mean values of ≈12.5 Gyr at [Fe/H] ∼< – 1.7 to ≈11 Gyr at [Fe/H] ∼> –1. At intermediate metallicities, the age-metallicity relation (AMR) appears to be bifurcated: one branch apparently contains clusters with disk-like kinematics, whereas the other branch, which is displaced to lower [Fe/H] values by ≈0.6 dex at a fixed age, is populated by clusters with halo-type orbits. The dispersion in age about each component of the AMR is ∼ ± 0.5 Gyr. There is no apparent dependence of age on Galactocentric distance (R{sub G}) nor is there a clear correlation of HB type with age. As previously discovered in the case of M3 and M13, subtle variations have been found in the slope of the SGB in the CMDs of other metal-poor ([Fe/H] ∼< – 1.5) GCs. They have been tentatively attributed to cluster-to-cluster differences in the abundance of helium. Curiously, GCs that have relatively steep 'M13-like' SGBs tend to be massive systems, located at small R{sub G}, that show the strongest evidence of in situ formation of multiple stellar populations. The clusters in the other group are typically low-mass systems (with 2-3 exceptions, including M3) that, at the present time, should not be able to retain the matter lost by mass-losing stars due either to the development of GC winds or to ram-pressure stripping by the halo interstellar medium. The apparent separation of the two groups in terms of their present-day gas retention properties is difficult to understand if all GCs were initially ∼20 times their current masses. The lowest-mass systems, in particular, may have never been massive enough to retain enough gas to produce a significant population of second-generation stars. In this case, the observed light element abundance variations, which are characteristic of all GCs, were presumably present in the gas out of which the observed cluster stars formed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

MODELING THE NON-RECYCLED FERMI GAMMA-RAY PULSAR POPULATION

We use Fermi Gamma-ray Space Telescope detections and upper limits on non-recycled pulsars obtained from the Large Area Telescope (LAT) to constrain how the gamma-ray luminosity L{sub γ} depends on the period P and the period derivative P-dot . We use a Bayesian analysis to calculate a best-fit luminosity law, or dependence of L{sub γ} on P and P-dot , including different methods for modeling the beaming factor. An outer gap (OG) magnetosphere geometry provides the best-fit model, which is L{sub γ}∝P{sup -a} P-dot {sup b} where a = 1.36 ± 0.03 and b = 0.44 ± 0.02, similar to but not identical to the commonly assumed L{sub γ}∝√( E-dot )∝P{sup -1.5} P-dot {sup 0.5}. Given upper limits on gamma-ray fluxes of currently known radio pulsars and using the OG model, we find that about 92% of the radio-detected pulsars have gamma-ray beams that intersect our line of sight. By modeling the misalignment of radio and gamma-ray beams of these pulsars, we find an average gamma-ray beaming solid angle of about 3.7π for the OG model, assuming a uniform beam. Using LAT-measured diffuse fluxes, we place a 2σ upper limit on the average braking index and a 2σ lower limit on the average surface magnetic field strength of the pulsar population of 3.8 and 3.2 × 10{sup 10} G, respectively. We then predict the number of non-recycled pulsars detectable by the LAT based on our population model. Using the 2 yr sensitivity, we find that the LAT is capable of detecting emission from about 380 non-recycled pulsars, including 150 currently identified radio pulsars. Using the expected 5 yr sensitivity, about 620 non-recycled pulsars are detectable, including about 220 currently identified radio pulsars. We note that these predictions significantly depend on our model assumptions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

CHARACTERIZING THE ORBITAL AND DYNAMICAL STATE OF THE HD 82943 PLANETARY SYSTEM WITH KECK RADIAL VELOCITY DATA

We present an updated analysis of radial velocity data of the HD 82943 planetary system based on 10 yr of measurements obtained with the Keck telescope. Previous studies have shown that the HD 82943 system has two planets that are likely in 2:1 mean-motion resonance (MMR), with orbital periods about 220 and 440 days. However, alternative fits that are qualitatively different have also been suggested, with two planets in a 1:1 resonance or three planets in a Laplace 4:2:1 resonance. Here we use χ{sup 2} minimization combined with a parameter grid search to investigate the orbital parameters and dynamical states of the qualitatively different types of fits, and we compare the results to those obtained with the differential evolution Markov chain Monte Carlo method. Our results support the coplanar 2:1 MMR configuration for the HD 82943 system, and show no evidence for either the 1:1 or three-planet Laplace resonance fits. The inclination of the system with respect to the sky plane is well constrained at 20{sup +4.9}{sub -5.5} degrees, and the system contains two planets with masses of about 4.78 M{sub J} and 4.80 M{sub J} (where M{sub J} is the mass of Jupiter) and orbital periods of about 219 and 442 days for the inner and outer planet, respectively. The best fit is dynamically stable with both eccentricity-type resonant angles θ{sub 1} and θ{sub 2} librating around 0°.

79 ASTRONOMY AND ASTROPHYSICS↗