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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

EVOLUTION OF THE COSMIC MICROWAVE BACKGROUND POWER SPECTRUM ACROSS WILKINSON MICROWAVE ANISOTROPY PROBE DATA RELEASES: A NONPARAMETRIC ANALYSIS

Using a nonparametric function estimation methodology, we present a comparative analysis of the Wilkinson Microwave Anisotropy Probe (WMAP) 1-, 3-, 5-, and 7-year data releases for the cosmic microwave background (CMB) angular power spectrum with respect to the following key questions. (1) How well is the power spectrum determined by the data alone? (2) How well is the {Lambda}CDM model supported by a model-independent, data-driven analysis? (3) What are the realistic uncertainties on peak/dip locations and heights? Our results show that the height of the power spectrum is well determined by data alone for multipole l approximately less than 546 (1-year), 667 (3-year), 804 (5-year), and 842 (7-year data). We show that parametric fits based on the {Lambda}CDM model are remarkably close to our nonparametric fits in l-regions where data are sufficiently precise. In contrast, the power spectrum for an H{Lambda}CDM model is progressively pushed away from our nonparametric fit as data quality improves with successive data realizations, suggesting incompatibility of this particular cosmological model with respect to the WMAP data sets. We present uncertainties on peak/dip locations and heights at the 95% (2{sigma}) level of confidence and show how these uncertainties translate into hyperbolic 'bands' on the acoustic scale (l{sub A} ) and peak shift ({phi}{sub m}) parameters. Based on the confidence set for the 7-year data, we argue that the low-l upturn in the CMB power spectrum cannot be ruled out at any confidence level in excess of about 10% ( Almost-Equal-To 0.12{sigma}). Additional outcomes of this work are a numerical formulation for minimization of a noise-weighted risk function subject to monotonicity constraints, a prescription for obtaining nonparametric fits that are closer to cosmological expectations on smoothness, and a method for sampling cosmologically meaningful power spectrum variations from the confidence set of a nonparametric fit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A SELF-SIMILAR EXPANSION MODEL FOR USE IN SOLAR WIND TRANSIENT PROPAGATION STUDIES

Since the advent of wide-angle imaging of the inner heliosphere, a plethora of techniques have been developed to investigate the three-dimensional structure and kinematics of solar wind transients, such as coronal mass ejections, from their signatures in single- and multi-spacecraft imaging observations. These techniques, which range from the highly complex and computationally intensive to methods based on simple curve fitting, all have their inherent advantages and limitations. In the analysis of single-spacecraft imaging observations, much use has been made of the fixed {phi} fitting (FPF) and harmonic mean fitting (HMF) techniques, in which the solar wind transient is considered to be a radially propagating point source (fixed {phi}, FP, model) and a radially expanding circle anchored at Sun centre (harmonic mean, HM, model), respectively. Initially, we compare the radial speeds and propagation directions derived from application of the FPF and HMF techniques to a large set of STEREO/Heliospheric Imager (HI) observations. As the geometries on which these two techniques are founded constitute extreme descriptions of solar wind transients in terms of their extent along the line of sight, we describe a single-spacecraft fitting technique based on a more generalized model for which the FP and HM geometries form the limiting cases. In addition to providing estimates of a transient's speed and propagation direction, the self-similar expansion fitting (SSEF) technique provides, in theory, the capability to estimate the transient's angular extent in the plane orthogonal to the field of view. Using the HI observations, and also by performing a Monte Carlo simulation, we assess the potential of the SSEF technique.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

MASS/RADIUS CONSTRAINTS ON THE QUIESCENT NEUTRON STAR IN M13 USING HYDROGEN AND HELIUM ATMOSPHERES

The mass and radius of the neutron star (NS) in low-mass X-ray binaries can be obtained by fitting the X-ray spectrum of the NS in quiescence, and the mass and radius constrains the properties of dense matter in NS cores. A critical ingredient for spectral fits is the composition of the NS atmosphere: hydrogen atmospheres are assumed in most prior work, but helium atmospheres are possible if the donor star is a helium white dwarf. Here we perform spectral fits to XMM-Newton, Chandra, and ROSAT data of a quiescent NS in the globular cluster M13. This NS has the smallest inferred radius from previous spectral fitting. Assuming an atmosphere composed of hydrogen, we find a significantly larger radius, more consistent with those from other quiescent NSs. With a helium atmosphere (an equally acceptable fit), we find even larger values for the radius.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

EMITTING ELECTRONS SPECTRA AND ACCELERATION PROCESSES IN THE JET OF Mrk 421: FROM THE LOW STATE TO THE GIANT FLARE STATE

We investigate the electron energy distributions (EEDs) and the acceleration processes in the jet of Mrk 421 through fitting the spectral energy distributions (SEDs) in different active states in the frame of a one-zone synchrotron self-Compton model. After assuming two possible EEDs formed in different acceleration models: the shock-accelerated power law with exponential cut-off (PLC) EED and the stochastic-turbulence-accelerated log-parabolic (LP) EED, we fit the observed SEDs of Mrk 421 in both low and giant flare states using the Markov Chain Monte Carlo method which constrains the model parameters in a more efficient way. The results from our calculations indicate that (1) the PLC and LP models give comparably good fits for the SED in the low state, but the variations of model parameters from low state to flaring can be reasonably explained only in the case of the PLC in the low state; and (2) the LP model gives better fits compared to the PLC model for the SED in the flare state, and the intra-day/night variability observed at GeV-TeV bands can be accommodated only in the LP model. The giant flare may be attributed to the stochastic turbulence re-acceleration of the shock-accelerated electrons in the low state. Therefore, we may conclude that shock acceleration is dominant in the low state, while stochastic turbulence acceleration is dominant in the flare state. Moreover, our result shows that the extrapolated TeV spectra from the best-fit SEDs from optical through GeV with the two EEDs are different. It should be considered with caution when such extrapolated TeV spectra are used to constrain extragalactic background light models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE DUST PROPERTIES OF z {approx} 3 MIPS-LBGs FROM PHOTOCHEMICAL MODELS

The stacked spectral energy distribution (SED) 24 {mu}m Lyman break galaxies (MIPS-LBGs) detected by the Multiband Imaging Photometer for Spitzer (MIPS) is fitted by means of the spectrophotometric model GRASIL with an ''educated'' fitting approach which benefits from the results of chemical evolution models. The star formation rate-age-metallicity degeneracies of SED modeling are broken by using star formation history (SFH) and chemical enrichment history suggested by chemical models. The dust mass, dust abundance, and chemical pattern of elements locked in the dust component are also directly provided by chemical models. Using our new ''fitting'' approach, we derive the total mass M{sub tot}, stellar mass M{sub *}, gas mass M{sub g} , dust mass M{sub d} , age, and star formation rate (SFR) of the stacked MIPS-LBG in a self-consistent way. Our estimate of M{sub *} = 8 Multiplication-Sign 10{sup 10} of the stacked MIPS-LBG agrees with other works based on UV-optical SED fitting. We suggest that the MIPS-LBGs at z {approx} 3 are young (0.3-0.6 Gyr), massive (M{sub tot} {approx} 10{sup 11} M{sub Sun }), dusty (M{sub d} {approx} 10{sup 8} M{sub Sun }), and metal-rich (Z {approx} Z{sub Sun }) progenitors of elliptical galaxies undergoing a strong burst of star formation (SFR {approx} 200 M{sub Sun} yr{sup -1}). Our estimate of M{sub d} = 7 Multiplication-Sign 10{sup 7} M{sub Sun} of the stacked MIPS-LBG is about a factor of eight lower than the estimated value based on single temperature graybody fitting, suggesting that self-consistent SED models are needed to estimate dust mass. By comparing with Milky Way molecular cloud and dust properties, we suggest that denser and dustier environments and flatter dust size distribution are likely in high-redshift massive star-forming galaxies. These dust properties, as well as the different types of SFHs, can cause different SED shapes between high-redshift star-forming ellipticals and local starburst templates. This discrepancy of SED shapes could in turn explain the non-detection at submillimeter wavelengths of IR luminous (L{sub IR} Succeeds-Above-Single-Line-Equals-Sign 10{sup 12} L{sub Sun }) MIPS-LBGs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Properties of polycyclic aromatic hydrocarbons in the northwest photon dominated region OF NGC 7023. I. PAH size, charge, composition, and structure distribution

Polycyclic aromatic hydrocarbon (PAH) emission in the Spitzer Infrared Spectrograph spectral map of the northwest photon dominated region (PDR) in NGC 7023 was analyzed exclusively using PAH spectra from the NASA Ames PAH IR Spectroscopic Database (www.astrochem.org/pahdb). The 5-15 μm spectrum at each pixel is fitted using a non-negative-least-squares fitting approach. The fits are of good quality, allowing decomposition of the PAH emission into four subclasses: size, charge, composition, and hydrogen adjacency (structure). Maps tracing PAH subclass distributions across the region paint a coherent astrophysical picture. Once past some 20 seconds of arc from HD 200775, the emission is dominated by the more stable, large, symmetric, compact PAH cations with smaller, neutral PAHs taking over along the lines-of-sight toward the more distant molecular cloud. The boundary between the PDR and the denser cloud material shows up as a distinct discontinuity in the breakdown maps. Noteworthy is the requirement for PANH cations to fit the bulk of the 6.2 and 11.0 μm features and the indication of PAH photo-dehydrogenation and fragmentation close to HD 200775. Decomposition of the spectral maps into 'principal' subclass template spectra provides additional insight into the behavior of each subclass. However, the general applicability of this computationally more efficient approach is presently undetermined. This is the first time the spectra of individual PAHs are exclusively used to fit the 5-15 μm region and analyze the spatial behavior of the aromatic infrared bands, providing fundamental, new information about astronomical PAH subpopulations including their dependence on, and response to, changes in local conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Modeling the near-ultraviolet band of GK stars. III. Dependence on abundance pattern

We extend the grid of non-LTE (NLTE) models presented in Paper II to explore variations in abundance pattern in two ways: (1) the adoption of the Asplund et al. (GASS10) abundances, (2) for stars of metallicity, [M/H], of –0.5, the adoption of a non-solar enhancement of α-elements by +0.3 dex. Moreover, our grid of synthetic spectral energy distributions (SEDs) is interpolated to a finer numerical resolution in both T {sub eff} (ΔT {sub eff} = 25 K) and log g (Δlog g = 0.25). We compare the values of T {sub eff} and log g inferred from fitting LTE and NLTE SEDs to observed SEDs throughout the entire visible band, and in an ad hoc 'blue' band. We compare our spectrophotometrically derived T {sub eff} values to a variety of T {sub eff} calibrations, including more empirical ones, drawn from the literature. For stars of solar metallicity, we find that the adoption of the GASS10 abundances lowers the inferred T {sub eff} value by 25-50 K for late-type giants, and NLTE models computed with the GASS10 abundances give T {sub eff} results that are marginally in better agreement with other T {sub eff} calibrations. For stars of [M/H] = –0.5 there is marginal evidence that adoption of α-enhancement further lowers the derived T {sub eff} value by 50 K. Stellar parameters inferred from fitting NLTE models to SEDs are more dependent than LTE models on the wavelength region being fitted, and we find that the effect depends on how heavily line blanketed the fitting region is, whether the fitting region is to the blue of the Wien peak of the star's SED, or both.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

ESTIMATING THE PROPERTIES OF HARD X-RAY SOLAR FLARES BY CONSTRAINING MODEL PARAMETERS

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23. We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non-Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, E{sub c} , a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that E{sub c} lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 {+-} 4 keV with 68% probability. Using the Bayesian/MCMC approach, we also estimate for the first time probability density functions for the total number of flare-accelerated electrons and the energy they carry for each flare studied. For the 2002 July 23 event, these probability density functions are asymmetric with long tails orders of magnitude higher than the most probable value, caused by the poorly constrained value of the low-energy cutoff. The most probable electron power is estimated at 10{sup 28.1} erg s{sup -1}, with a 68% credible interval estimated at 10{sup 28.1}-10{sup 29.0} erg s{sup -1}, and a 95% credible interval estimated at 10{sup 28.0}-10{sup 30.2} erg s{sup -1}. For the 2005 January 19 flare spectrum, the probability density functions for the total number of flare-accelerated electrons and their energy are much more symmetric and narrow: the most probable electron power is estimated at 10{sup 27.66{+-}0.01} erg s{sup -1} (68% credible intervals). However, in this case the uncertainty due to systematic sources of error is estimated to dominate the uncertainty due to statistical sources of error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Non-local thermodynamic equilibrium 1.5D modeling of red giant stars

Spectra for two-dimensional (2D) stars in the 1.5D approximation are created from synthetic spectra of one-dimensional (1D) non-local thermodynamic equilibrium (NLTE) spherical model atmospheres produced by the PHOENIX code. The 1.5D stars have the spatially averaged Rayleigh-Jeans flux of a K3-4 III star while varying the temperature difference between the two 1D component models (ΔT {sub 1.5D}) and the relative surface area covered. Synthetic observable quantities from the 1.5D stars are fitted with quantities from NLTE and local thermodynamic equilibrium (LTE) 1D models to assess the errors in inferred T {sub eff} values from assuming horizontal homogeneity and LTE. Five different quantities are fit to determine the T {sub eff} of the 1.5D stars: UBVRI photometric colors, absolute surface flux spectral energy distributions (SEDs), relative SEDs, continuum normalized spectra, and TiO band profiles. In all cases except the TiO band profiles, the inferred T {sub eff} value increases with increasing ΔT {sub 1.5D}. In all cases, the inferred T {sub eff} value from fitting 1D LTE quantities is higher than from fitting 1D NLTE quantities and is approximately constant as a function of ΔT {sub 1.5D} within each case. The difference between LTE and NLTE for the TiO bands is caused indirectly by the NLTE temperature structure of the upper atmosphere, as the bands are computed in LTE. We conclude that the difference between T {sub eff} values derived from NLTE and LTE modeling is relatively insensitive to the degree of the horizontal inhomogeneity of the star being modeled and largely depends on the observable quantity being fit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Untangling the nature of spatial variations of cold dust properties in star forming galaxies

We investigate the far-infrared (IR) dust emission for 20 local star forming galaxies from the Key Insights on Nearby Galaxies: A Far-IR Survey with Herschel (KINGFISH) sample. We model the far-IR/submillimeter spectral energy distribution (SED) using images from Spitzer Space Telescope and Herschel Space Observatory. We calculate the cold dust temperature (T{sub c} ) and emissivity (β) on a pixel by pixel basis (where each pixel ranges from 0.1 to 3 kpc{sup 2}) using a two-temperature modified blackbody fitting routine. Our fitting method allows us to investigate the resolved nature of temperature and emissivity variations by modeling from the galaxy centers to the outskirts (physical scales of ∼15-50 kpc, depending on the size of the galaxy). We fit each SED in two ways: (1) fit T{sub c} and β simultaneously, (2) hold β constant and fit T{sub c} . We compare T{sub c} and β with star formation rates (calculated from L{sub Hα} and L{sub 24μm}), the luminosity of the old stellar population (traced through L{sub 3.6μm}), and the dust mass surface density (traced by 500 μm luminosity, L{sub 500}). We find a significant trend between SFR/L{sub 500} and T{sub c} , implying that the flux of hard UV photons relative to the amount of dust is significantly contributing to the heating of the cold, or diffuse, dust component. We also see a trend between L{sub 3.6}/L{sub 500} and β, indicating that the old stellar population contributes to the heating at far-IR/submillimeter wavelengths. Finally, we find that when β is held constant, T{sub c} exhibits a strongly decreasing radial trend, illustrating that the shape of the far-IR SED is changing radially through a galaxy, thus confirming on a sample almost double in size the trends observed in Galametz et al.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Star formation and black hole growth at z ≅ 4.8

We report Herschel/SPIRE, Spitzer and Wide-field Infrared Survey Explorer observations of 44 z ≅ 4.8 optically selected active galactic nuclei (AGNs). This flux-limited sample contains the highest mass black holes (BHs) at this redshift. Ten of the objects were detected by Herschel and five show emission that is not clearly associated with the AGNs. The star formation (SF) luminosity (L{sub SF}) obtained by fitting the spectral energy distribution (SED) with standard SF templates, taking into account AGN contribution, is in the range 10{sup 46.62}-10{sup 47.21} erg s{sup –1} corresponding to SF rates of 1090-4240 M{sub ☉} yr{sup –1}. Fitting with very luminous submillimeter galaxy SEDs gives SF rates that are smaller by 0.05 dex when using all bands and 0.1 dex when ignoring the 250 μm band. A 40 K graybody fits to only the 500 μm fluxes reduce L{sub SF} by about a factor of two. A stacking analysis of 29 undetected sources gives significant signals in all three bands. A SF template fit indicates L{sub SF} = 10{sup 46.19-46.23} erg s{sup –1} depending on the assumed AGN contribution. A 40 K fit to the stacked 500 μm flux gives L{sub SF} = 10{sup 45.95} erg s{sup –1}. The mean BH mass (M{sub BH}) and AGN luminosity (L{sub AGN}) of the detected sources are significantly higher than those of the undetected ones. The spectral differences are seen all the way from UV to far infrared wavelengths. The mean optical-UV spectra are similar to those predicted for thin accretion disks around BHs with similar masses and accretion rates. We suggest two alternative explanations to the correlation of L{sub SF}, L{sub AGN} and M{sub BH}, one involving no AGN feedback and the second involving moderate feedback that affects, but does not totally quench, SF in three-quarters of the sources. We compare our L{sub SF} and L{sub AGN} to lower redshift samples and show a new correlation between L{sub SF} and M{sub BH}. We also examine several rather speculative ideas about the host galaxy properties including the possibility that the detected sources are above the SF mass sequence (MS) at z ≅ 4.8, perhaps in mergers, and most of the undetected sources are on the MS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The light curve of SN 1987A revisited: constraining production masses of radioactive nuclides

We revisit the evidence for the contribution of the long-lived radioactive nuclides {sup 44}Ti, {sup 55}Fe, {sup 56}Co, {sup 57}Co, and {sup 60}Co to the UVOIR light curve of SN 1987A. We show that the V-band luminosity constitutes a roughly constant fraction of the bolometric luminosity between 900 and 1900 days, and we obtain an approximate bolometric light curve out to 4334 days by scaling the late time V-band data by a constant factor where no bolometric light curve data is available. Considering the five most relevant decay chains starting at {sup 44}Ti, {sup 55}Co, {sup 56}Ni, {sup 57}Ni, and {sup 60}Co, we perform a least squares fit to the constructed composite bolometric light curve. For the nickel isotopes, we obtain best fit values of M({sup 56}Ni) = (7.1 ± 0.3) × 10{sup –2} M {sub ☉} and M({sup 57}Ni) = (4.1 ± 1.8) × 10{sup –3} M {sub ☉}. Our best fit {sup 44}Ti mass is M({sup 44}Ti) = (0.55 ± 0.17) × 10{sup –4} M {sub ☉}, which is in disagreement with the much higher (3.1 ± 0.8) × 10{sup –4} M {sub ☉} recently derived from INTEGRAL observations. The associated uncertainties far exceed the best fit values for {sup 55}Co and {sup 60}Co and, as a result, we only give upper limits on the production masses of M({sup 55}Co) < 7.2 × 10{sup –3} M {sub ☉} and M({sup 60}Co) < 1.7 × 10{sup –4} M {sub ☉}. Furthermore, we find that the leptonic channels in the decay of {sup 57}Co (internal conversion and Auger electrons) are a significant contribution and constitute up to 15.5% of the total luminosity. Consideration of the kinetic energy of these electrons is essential in lowering our best fit nickel isotope production ratio to [{sup 57}Ni/{sup 56}Ni] = 2.5 ± 1.1, which is still somewhat high but is in agreement with gamma-ray observations and model predictions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Kinematic modeling of the Milky Way using the RAVE and GCS stellar surveys

We investigate the kinematic parameters of the Milky Way disk using the Radial Velocity Experiment (RAVE) and Geneva-Copenhagen Survey (GCS) stellar surveys. We do this by fitting a kinematic model to the data and taking the selection function of the data into account. For stars in the GCS we use all phase-space coordinates, but for RAVE stars we use only (ℓ, b, v {sub los}). Using the Markov Chain Monte Carlo technique, we investigate the full posterior distributions of the parameters given the data. We investigate the age-velocity dispersion relation for the three kinematic components (σ {sub R}, σ{sub φ}, σ {sub z}), the radial dependence of the velocity dispersions, the solar peculiar motion (U {sub ☉}, V {sub ☉}, W {sub ☉}), the circular speed Θ{sub 0} at the Sun, and the fall of mean azimuthal motion with height above the midplane. We confirm that the Besançon-style Gaussian model accurately fits the GCS data but fails to match the details of the more spatially extended RAVE survey. In particular, the Shu distribution function (DF) handles noncircular orbits more accurately and provides a better fit to the kinematic data. The Gaussian DF not only fits the data poorly but systematically underestimates the fall of velocity dispersion with radius. The radial scale length of the velocity dispersion profile of the thick disk was found to be smaller than that of the thin disk. We find that correlations exist between a number of parameters, which highlights the importance of doing joint fits. The large size of the RAVE survey allows us to get precise values for most parameters. However, large systematic uncertainties remain, especially in V {sub ☉} and Θ{sub 0}. We find that, for an extended sample of stars, Θ{sub 0} is underestimated by as much as 10% if the vertical dependence of the mean azimuthal motion is neglected. Using a simple model for vertical dependence of kinematics, we find that it is possible to match the Sgr A* proper motion without any need for V {sub ☉} being larger than that estimated locally by surveys like GCS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A Hubble space telescope search for a sub-Earth-sized exoplanet in the GJ 436 system

The detection of small planets orbiting nearby stars is an important step toward the identification of Earth twins. In previous work using the Spitzer Space Telescope, we found evidence to support at least one sub-Earth-sized exoplanet orbiting the nearby mid-M dwarf star GJ 436. As a follow up, here we used the Hubble Space Telescope (HST) to investigate the existence of one of these candidate planets, UCF-1.01, by searching for two transit signals as it passed in front of its host star. Interpretation of the data hinges critically on correctly modeling and removing the Wide Field Camera 3 (WFC3) instrument systematics from the light curves. Building on previous HST work, we demonstrate that WFC3 analyses need to explore the use of a quadratic function to fit a visit-long time-dependent systematic. This is important for establishing absolute transit and eclipse depths in the white light curves of all transiting systems. The work presented here exemplifies this point by putatively detecting the primary transit of UCF-1.01 with the use of a linear trend. However, using a quadratic trend, we achieve a better fit to the white light curves and a reduced transit depth that is inconsistent with previous Spitzer measurements. Furthermore, quadratic trends with or without a transit model component produce comparable fits to the available data. Using extant WFC3 transit light curves for GJ 436b, we further validate the quadratic model component by achieving photon-limited model fit residuals and consistent transit depths over multiple epochs. We conclude that, when we fit for a quadratic trend, our new data contradict the prediction of a sub-Earth-sized planet orbiting GJ 436 with the size, period, and ephemeris posited from the Spitzer data by a margin of 3.1σ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

CONSTRAINING THE MILKY WAY'S HOT GAS HALO WITH O VII AND O VIII EMISSION LINES

The Milky Way hosts a hot (≈2 × 10{sup 6} K), diffuse, gaseous halo based on detections of z = 0 O VII and O VIII absorption lines in quasar spectra and emission lines in blank-sky spectra. Here we improve constraints on the structure of the hot gas halo by fitting a radial model to a much larger sample of O VII and O VIII emission line measurements from XMM-Newton/EPIC-MOS spectra compared to previous studies (≈650 sightlines). We assume a modified β-model for the halo density distribution and a constant-density Local Bubble from which we calculate emission to compare with the observations. We find an acceptable fit to the O VIII emission line observations with χ{sub red}{sup 2} (dof) = 1.08 (644) for best-fit parameters of n{sub o}r{sub c}{sup 3β}=1.35±0.24 cm{sup –3} kpc{sup 3β} and β = 0.50 ± 0.03 for the hot gas halo and negligible Local Bubble contribution. The O VII observations yield an unacceptable χ{sub red}{sup 2} (dof) = 4.69 (645) for similar best-fit parameters, which is likely due to temperature or density variations in the Local Bubble. The O VIII fitting results imply hot gas masses of M(<50 kpc) = 3.8{sub −0.3}{sup +0.3}×10{sup 9} M{sub ⊙} and M(<250 kpc) = 4.3{sub −0.8}{sup +0.9}×10{sup 10} M{sub ⊙}, accounting for ≲50% of the Milky Way's missing baryons. We also explore our results in the context of optical depth effects in the halo gas, the halo gas cooling properties, temperature and entropy gradients in the halo gas, and the gas metallicity distribution. The combination of absorption and emission line analyses implies a sub-solar gas metallicity that decreases with radius, but that also must be ≥0.3 Z {sub ☉} to be consistent with the pulsar dispersion measure toward the Large Magellanic Cloud.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Interstellar dust properties from a survey of X-ray halos

Interstellar dust grains produce X-ray halos around bright sources due to small-angle X-ray scattering. Numerous studies have examined these halos, but no systematic study has yet tested the available halo data against the large number of well-defined dust models in circulation. We have therefore obtained the largest sample to date of X-ray dust halos from XMM-Newton and Chandra, and fitted them with 14 commonly used dust grain models, including comparisons with the optical extinction, A{sub V}, where available in the literature. Our main conclusions are summarized as follows. (1) Comparing A{sub V} with N{sub H} values measured via X-ray spectral fits, we find a ratio of A{sub V}/N{sub H} (10{sup 21} cm{sup −2}) = 0.48 ± 0.06, in agreement with previous work. (2) Out of 35 halos, 27 could be fit by one or more grain models, with the most successful models having maximum grain radius a{sub max}<0.4 μm and fewer large grains than the less successful models. This suggests that the diffuse ISM does not contain a signicant presence of grains with a{sub max}>0.5 μm. (3) Most halos were best fit assuming a single dust cloud dominated the scattering, rather than smoothly distributed dust along the sightline. (4) Eight sources could not be fit with the models considered here, most of which were along distant (d>5 kpc) sight lines through the Galactic thin disk. (5) Some sight lines had halos with observed X-ray scattering optical depth τ{sub sca}/A{sub V} that were signicantly different than expected. This may result from an inhomogeneous dust distribution across the halo extraction area.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

VLBI FOR GRAVITY PROBE B. IV. A NEW ASTROMETRIC ANALYSIS TECHNIQUE AND A COMPARISON WITH RESULTS FROM OTHER TECHNIQUES

When very long baseline interferometry (VLBI) observations are used to determine the position or motion of a radio source relative to reference sources nearby on the sky, the astrometric information is usually obtained via (1) phase-referenced maps or (2) parametric model fits to measured fringe phases or multiband delays. In this paper, we describe a 'merged' analysis technique which combines some of the most important advantages of these other two approaches. In particular, our merged technique combines the superior model-correction capabilities of parametric model fits with the ability of phase-referenced maps to yield astrometric measurements of sources that are too weak to be used in parametric model fits. We compare the results from this merged technique with the results from phase-referenced maps and from parametric model fits in the analysis of astrometric VLBI observations of the radio-bright star IM Pegasi (HR 8703) and the radio source B2252+172 nearby on the sky. In these studies we use central-core components of radio sources 3C 454.3 and B2250+194 as our positional references. We obtain astrometric results for IM Peg with our merged technique even when the source is too weak to be used in parametric model fits, and we find that our merged technique yields astrometric results superior to the phase-referenced mapping technique. We used our merged technique to estimate the proper motion and other astrometric parameters of IM Peg in support of the NASA/Stanford Gravity Probe B mission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Analysis of the theoretical bias in dark matter direct detection

Fitting the model ''A'' to dark matter direct detection data, when the model that underlies the data is ''B'', introduces a theoretical bias in the fit. We perform a quantitative study of the theoretical bias in dark matter direct detection, with a focus on assumptions regarding the dark matter interactions, and velocity distribution. We address this problem within the effective theory of isoscalar dark matter-nucleon interactions mediated by a heavy spin-1 or spin-0 particle. We analyze 24 benchmark points in the parameter space of the theory, using frequentist and Bayesian statistical methods. First, we simulate the data of future direct detection experiments assuming a momentum/velocity dependent dark matter-nucleon interaction, and an anisotropic dark matter velocity distribution. Then, we fit a constant scattering cross section, and an isotropic Maxwell-Boltzmann velocity distribution to the simulated data, thereby introducing a bias in the analysis. The best fit values of the dark matter particle mass differ from their benchmark values up to 2 standard deviations. The best fit values of the dark matter-nucleon coupling constant differ from their benchmark values up to several standard deviations. We conclude that common assumptions in dark matter direct detection are a source of potentially significant bias.

79 ASTRONOMY AND ASTROPHYSICS↗