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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Implementation and application of the peak scaling method for temperature measurement in the laser heated diamond anvil cell

In this work, a new design for a double-sided high-pressure diamond anvil cell laser heating set-up is described. The prototype is deployed at beamline 12.2.2 of the Advanced Light Source at Lawrence Berkeley National Lab. Our compact design features shortened mechanical lever arms, which results in more stable imaging optics, and thus more user friendly and more reliable temperature measurements based on pyrometry. A modification of the peak scaling method was implemented for pyrometry, including an iterative method to determine the absolute peak temperature, thus allowing for quasi-real time temperature mapping of the actual hotspot within a laser-heated diamond anvil cell without any assumptions on shape, size, and symmetry of the hotspot and without any assumptions to the relationship between fitted temperature and peak temperature. This is important since we show that the relationship between peak temperature and temperature obtained by fitting the Planck function against the thermal emission spectrum averaged over the entire hotspot is not constant but depends on variable fitting parameters (in particular, the size and position of the fitting window). The accuracy of the method is confirmed through measuring melting points of metal wires at ambient pressure. Having absolute temperature maps in real time allows for more differentiated analyses of laser heating experiments. Finally, we present such an example of the pressure variations within a heated hotspot of AgI at a loaded base pressure of 3.8 GPa.

42 ENGINEERING↗
DOE OSTI2019

Direct diabatization and analytic representation of coupled potential energy surfaces and couplings for the reactive quenching of the excited 2 Σ + state of OH by molecular hydrogen

In this work, we have employed extended multiconfiguration quasidegenerate perturbation theory, fourfold-way diabatic molecular orbitals, and configurational uniformity to develop a global three-state diabatic representation of the potential energy surfaces and their couplings for the electronically nonadiabatic reaction OH* + H 2 → H 2 O + H, where * denotes electronic excitation to the A 2 Σ + state. To achieve sign consistency of the computed diabatic couplings, we developed a graphics processing unit-accelerated algorithm called the cluster-growing algorithm. Having obtained consistent signs of the diabatic couplings, we fit the diabatic matrix elements (which consist of the diabatic potentials and the diabatic couplings) to analytic representations. Adiabatic potential energy surfaces are generated by diagonalizing the 3 × 3 diabatic potential energy matrix. The comparisons between the fitted and computed diabatic matrix elements and between the originally computed adiabatic potential energy surfaces and those generated from the fits indicate that the current fit is accurate enough for dynamical studies, and it may be used for quantal or semiclassical dynamics calculations.

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

Direct measurement of storage and loss behavior in AFM force–distance experiments using the modified Fourier transformation

Force–distance curve experiments are commonly performed in atomic force microscopy (AFM) to obtain the viscoelastic characteristics of materials, such as the storage and loss moduli or compliances. The classic methods used to obtain these characteristics consist of fitting a viscoelastic material model to the experimentally obtained AFM data. Here, we demonstrate a new method that utilizes the modified discrete Fourier transform to approximate the storage and loss behavior of a material directly from the data, without the need for a fit. Additionally, one may still fit a model to the resulting storage and loss behavior if a parameterized description of the material is desired. In contrast to fitting the data to a model chosen a priori, departing from a model-free description of the material's frequency behavior guides the selection of the model, such that the user may choose the one that is most appropriate for the particular material under study. To this end, we also include modified Fourier domain descriptions of commonly used viscoelastic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

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

Revisiting the Anisotropic Complex Refractive Indices of Sodium Nitrate for Interpretation of the Reflectance Spectrum of Pressed Pellets

Reflectance spectroscopy is notoriously confounding in that the measured spectral response is highly dependent upon the morphology of the sample. Fortunately, all such perturbations are neatly encoded by the complex refractive index of the sample. Herein, we seek to quantitatively model the measured infrared reflectance spectrum of a specularly flat pressed pellet sample of the birefringent compound, sodium nitrate. Single crystals of sodium nitrate were synthesized via a slow evaporation process and spectroscopically analyzed using polarization-dependent infrared single-angle reflectance spectroscopy. The anisotropic complex refractive index was measured from 7500 to 300 cm-1 (1.33 to 33.33 µm). The deduced anisotropic optical constants were found to be consistent with those previously reported. Once measured and validated, the optical constants were applied to model the pressed pellet reflectance spectrum. It was evident that an average of the anisotropic refractive indices was insufficient to account for the measured pellet reflectance. In order to account for contributions of all possible microcrystalline orientations within the pellet, the Python package PYELLI was used to calculate a basis set of orientation-dependent reflection spectra spanning the distinct ? and ? Euler rotations of the uniaxial crystal. When the population of orientations was allowed to vary freely in a spectral fit analysis, the fit-deduced orientations were tightly clustered along f = 45º, hinting at residual anisotropy in the pressed pellet sample. Conversely, an equally valid spectral fit (with marginally worse fit metric) was obtained when the population was constrained to an isotropic distribution of orientations. Subsequent non-zero cross-polarization reflectance measurements likewise suggested anisotropy in the pellet. However, both grazing incidence wide-angle x-ray scattering and scanning electron microscopy measurements revealed that the microcrystal orientations at the surface of the pressed pellet sample were isotropically distributed (and that the average crystallite size was larger than ?/10). Application of the measured complex refractive indices for modeling the reflectance spectrum of the pressed pellet, and rectification of these seemingly contradictory observations will be discussed.

Wilhelm, Michael J.↗
DOE OSTI2013

THE VIRUS-P EXPLORATION OF NEARBY GALAXIES (VENGA): SURVEY DESIGN, DATA PROCESSING, AND SPECTRAL ANALYSIS METHODS

We present the survey design, data reduction, and spectral fitting pipeline for the VIRUS-P Exploration of Nearby Galaxies (VENGA). VENGA is an integral field spectroscopic survey, which maps the disks of 30 nearby spiral galaxies. Targets span a wide range in Hubble type, star formation activity, morphology, and inclination. The VENGA data cubes have 5.''6 FWHM spatial resolution, {approx}5 A FWHM spectral resolution, sample the 3600 A-6800 A range, and cover large areas typically sampling galaxies out to {approx}0.7R{sub 25}. These data cubes can be used to produce two-dimensional maps of the star formation rate, dust extinction, electron density, stellar population parameters, the kinematics and chemical abundances of both stars and ionized gas, and other physical quantities derived from the fitting of the stellar spectrum and the measurement of nebular emission lines. To exemplify our methods and the quality of the data, we present the VENGA data cube on the face-on Sc galaxy NGC 628 (a.k.a. M 74). The VENGA observations of NGC 628 are described, as well as the construction of the data cube, our spectral fitting method, and the fitting of the stellar and ionized gas velocity fields. We also propose a new method to measure the inclination of nearly face-on systems based on the matching of the stellar and gas rotation curves using asymmetric drift corrections. VENGA will measure relevant physical parameters across different environments within these galaxies, allowing a series of studies on star formation, structure assembly, stellar populations, chemical evolution, galactic feedback, nuclear activity, and the properties of the interstellar medium in massive disk galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Local variability in the orbit of Saturn's F ring

We present an analysis of the orbit of Saturn's F ring using images recorded by the Imaging Science Subsystem of the Cassini spacecraft. A total of 9805 observations have been made from 10 image sequences obtained between 2006 November 23 and 2009 July 28. Each sequence of up to 240 images spans a single orbit of the F ring, allowing 10 independent high-precision estimates of the ring orbit to be made over this ∼3 year period. The ring has been modeled as an inclined uniformly precessing ellipse. The results show a variability in the orbital elements with, for example, the semi-major axis scattered between 140211.2 ± 0.1 km and 140232.9 ± 0.4 km and the fitted periapses locked to the value obtained from a combined fit using the entire three-year span of observations. We show that the observed scatter between the individual estimates of the ring orbit reflect the differing past histories of the particular segments of ring being fitted and that the values are scattered within the limits expected from a single gravitational encounter with the nearby moon, Prometheus. In the combined fit, the scatter averages out to reveal a small systematic bias with respect to the results of Bosh et al. and Albers et al. We believe this is a consequence of the proximity of Prometheus to the ring in the image sequences chosen for this analysis. Finally, we note a close empirical commensurability between the apsidal precession rate, ϖ-dot , of the F ring and the synodic period between Prometheus and the F ring, such that n {sub Prom} – n {sub Fring} ≈ 2 ϖ-dot {sub Fring}, where n {sub Prom} and n {sub Fring} are the mean motions, and discuss its implications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

A spatial characterization of the Sagittarius dwarf galaxy tidal tails

We measure the spatial density of F turnoff stars in the Sagittarius dwarf tidal stream, from Sloan Digital Sky Survey data, using statistical photometric parallax. We find a set of continuous, consistent parameters that describe the leading Sgr stream's position, direction, and width for 15 stripes in the north Galactic cap, and three stripes in the south Galactic cap. We produce a catalog of stars that has the density characteristics of the dominant leading Sgr tidal stream that can be compared with simulations. We find that the width of the leading (north) tidal tail is consistent with recent triaxial and axisymmetric halo model simulations. The density along the stream is roughly consistent with common disruption models in the north, but possibly not in the south. We explore the possibility that one or more of the dominant Sgr streams has been misidentified, and that one or more of the ''bifurcated'' pieces is the real Sgr tidal tail, but we do not reach definite conclusions. If two dwarf progenitors are assumed, fits to the planes of the dominant and ''bifurcated'' tidal tails favor an association of the Sgr dwarf spheroidal galaxy with the dominant southern stream and the ''bifurcated'' stream in the north. In the north Galactic cap, the best fit Hernquist density profile for the smooth component of the stellar halo is oblate, with a flattening parameter q = 0.53, and a scale length of r {sub 0} = 6.73. The southern data for both the tidal debris and the smooth component of the stellar halo do not match the model fits to the north, although the stellar halo is still overwhelmingly oblate. Finally, we verify that we can reproduce the parameter fits on the asynchronous MilkyWay@home volunteer computing platform.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

CONSTRAINTS ON THE SYNCHROTRON SHOCK MODEL FOR THE FERMI GRB 090820A OBSERVED BY GAMMA-RAY BURST MONITOR

Discerning the radiative dissipation mechanism for prompt emission in gamma-ray bursts (GRBs) requires detailed spectroscopic modeling that straddles the {nu}F{sub {nu}} peak in the 100 keV-1 MeV range. Historically, empirical fits such as the popular Band function have been employed with considerable success in interpreting the observations. While extrapolations of the Band parameters can provide some physical insight into the emission mechanisms responsible for GRBs, these inferences do not provide a unique way of discerning between models. By fitting physical models directly, this degeneracy can be broken, eliminating the need for empirical functions; our analysis here offers a first step in this direction. One of the oldest, and leading, theoretical ideas for the production of the prompt signal is the synchrotron shock model. Here we explore the applicability of this model to a bright Fermi gamma-ray burst monitor (GBM) burst with a simple temporal structure, GRB 090820A. Our investigation implements, for the first time, thermal and non-thermal synchrotron emissivities in the RMFIT forward-folding spectral analysis software often used in GBM burst studies. We find that these synchrotron emissivities, together with a blackbody shape, provide at least as good a match to the data as the Band GRB spectral fitting function. This success is achieved in both time-integrated and time-resolved spectral fits.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

MULTI-FREQUENCY STUDIES OF RADIO RELICS IN THE GALAXY CLUSTERS A4038, A1664, AND A786

We present a multi-frequency study of radio relics associated with the galaxy clusters A4038, A1664, and A786. Radio images, integrated spectra, spectral index maps, and fits to the integrated spectra in the framework of the adiabatic compression model are presented. Images of the relic in A4038 at 150, 240, and 606 MHz with the Giant Meterwave Radio Telescope have revealed extended ultra-steep spectrum ({alpha} {approx} -1.8 to -2.7) emission of extent 210 Multiplication-Sign 80 kpc{sup 2}. The model of passively evolving radio lobes compressed by a shock fits the integrated spectrum best. The relic with a circular morphology at the outskirts of the cluster A1664 has an integrated spectral index of {approx} - 1.10 {+-} 0.06 and is best fit by the model of radio lobes lurking for {approx}4 Multiplication-Sign 10{sup 7} yr. The relic near A786 has a curved spectrum and is best fit by a model of radio lobes lurking for {approx}3 Multiplication-Sign 10{sup 7} yr. At 4.7 GHz, a compact radio source, possibly the progenitor of the A786 relic, is detected near the center of the radio relic. The A786 radio relic is thus likely a lurking radio galaxy rather than a site of cosmological shock as has been considered in earlier studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

PSR J2030+3641: RADIO DISCOVERY AND GAMMA-RAY STUDY OF A MIDDLE-AGED PULSAR IN THE NOW IDENTIFIED FERMI-LAT SOURCE 1FGL J2030.0+3641

In a radio search with the Green Bank Telescope of three unidentified low Galactic latitude Fermi Large Area Telescope (LAT) sources, we have discovered the middle-aged pulsar J2030+3641 associated with 1FGL J2030.0+3641 (2FGL J2030.0+3640). Following the detection of gamma-ray pulsations using a radio ephemeris, we have obtained a phase-coherent timing solution based on gamma-ray and radio pulse arrival times which spans the entire Fermi mission. With a rotation period of 0.2 s, a spin-down luminosity of 3 Multiplication-Sign 10{sup 34} erg s{sup -1}, and a characteristic age of 0.5 Myr, PSR J2030+3641 is a middle-aged neutron star with spin parameters similar to those of the exceedingly gamma-ray-bright and radio-undetected Geminga. Its gamma-ray flux is 1% that of Geminga, primarily because of its much larger distance, as suggested by the large integrated column density of free electrons, DM = 246 pc cm{sup -3}. We fit the gamma-ray light curve, along with limited radio polarimetric constraints, to four geometrical models of magnetospheric emission, and while none of the fits have high significance some are encouraging and suggest that further refinements of these models may be worthwhile. We argue that not many more non-millisecond radio pulsars may be detected along the Galactic plane that are responsible for LAT sources, but that modified methods to search for gamma-ray pulsations should be productive-PSR J2030+3641 would have been found blindly in gamma rays if only {approx}> 0.8 GeV photons had been considered, owing to its relatively flat spectrum and location in a region of high soft background.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

MODELING THE NEAR-UV BAND OF GK STARS. II. NON-LTE MODELS

We present a grid of atmospheric models and synthetic spectral energy distributions (SEDs) for late-type dwarfs and giants of solar and 1/3 solar metallicity with many opacity sources computed in self-consistent non-local thermodynamic equilibrium (NLTE), and compare them to the LTE grid of Short and Hauschildt (Paper I). We describe, for the first time, how the NLTE treatment affects the thermal equilibrium of the atmospheric structure (T({tau}) relation) and the SED as a finely sampled function of T{sub eff}, log g, and [A/H] among solar metallicity and mildly metal-poor red giants. We compare the computed SEDs to the library of observed spectrophotometry described in Paper I across the entire visible band, and in the blue and red regions of the spectrum separately. We find that for the giants of both metallicities, the NLTE models yield best-fit T{sub eff} values that are 30-90 K lower than those provided by LTE models, while providing greater consistency between log g values, and, for Arcturus, T{sub eff} values, fitted separately to the blue and red spectral regions. There is marginal evidence that NLTE models give more consistent best-fit T{sub eff} values between the red and blue bands for earlier spectral classes among the solar metallicity GK giants than they do for the later classes, but no model fits the blue-band spectrum well for any class. For the two dwarf spectral classes that we are able to study, the effect of NLTE on derived parameters is less significant. We compare our derived T{sub eff} values to several other spectroscopic and photometric T{sub eff} calibrations for red giants, including one that is less model dependent based on the infrared flux method (IRFM). We find that the NLTE models provide slightly better agreement to the IRFM calibration among the warmer stars in our sample, while giving approximately the same level of agreement for the cooler stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

OBSERVATION AND SPECTRAL MEASUREMENTS OF THE CRAB NEBULA WITH MILAGRO

The Crab Nebula was detected with the Milagro experiment at a statistical significance of 17 standard deviations over the lifetime of the experiment. The experiment was sensitive to approximately 100 GeV-100 TeV gamma-ray air showers by observing the particle footprint reaching the ground. The fraction of detectors recording signals from photons at the ground is a suitable proxy for the energy of the primary particle and has been used to measure the photon energy spectrum of the Crab Nebula between {approx}1 and {approx}100 TeV. The TeV emission is believed to be caused by inverse-Compton upscattering of ambient photons by an energetic electron population. The location of a TeV steepening or cutoff in the energy spectrum reveals important details about the underlying electron population. We describe the experiment and the technique for distinguishing gamma-ray events from the much more-abundant hadronic events. We describe the calculation of the significance of the excess from the Crab and how the energy spectrum is fitted. The differential photon energy spectrum, including the statistical errors from the fit, obtained using a simple power-law hypothesis for data between 2005 September and 2008 March is (6.5 {+-} 0.4) Multiplication-Sign 10{sup -14}(E/10 TeV){sup -3.1{+-}0.1}(cm{sup 2} s TeV ){sup -1} between {approx}1 TeV and {approx}100 TeV. Allowing for a possible exponential cutoff, the photon energy spectrum is fitted as (2.5{sup +0.7}{sub -0.4}) Multiplication-Sign 10{sup -12}(E/3 TeV){sup -2.5{+-}0.4}exp (- E/32{sup +39}{sub -18} TeV) (cm{sup 2} s TeV){sup -1}. The results are subject to an {approx}30% systematic uncertainty in the overall flux and an {approx}0.1 systematic uncertainty in the power-law indices quoted. Uncertainty in the overall energy scale has been absorbed into these errors. Fixing the spectral index to values that have been measured below 1 TeV by IACT experiments (2.4-2.6), the fit to the Milagro data suggests that Crab exhibits a spectral steepening or cutoff between about 20-40 TeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

TRANSIT TIMING OBSERVATIONS FROM KEPLER. IV. CONFIRMATION OF FOUR MULTIPLE-PLANET SYSTEMS BY SIMPLE PHYSICAL MODELS

Eighty planetary systems of two or more planets are known to orbit stars other than the Sun. For most, the data can be sufficiently explained by non-interacting Keplerian orbits, so the dynamical interactions of these systems have not been observed. Here we present four sets of light curves from the Kepler spacecraft, each which of shows multiple planets transiting the same star. Departure of the timing of these transits from strict periodicity indicates that the planets are perturbing each other: the observed timing variations match the forcing frequency of the other planet. This confirms that these objects are in the same system. Next we limit their masses to the planetary regime by requiring the system remain stable for astronomical timescales. Finally, we report dynamical fits to the transit times, yielding possible values for the planets' masses and eccentricities. As the timespan of timing data increases, dynamical fits may allow detailed constraints on the systems' architectures, even in cases for which high-precision Doppler follow-up is impractical.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

AN EXAMINATION OF THE X-RAY SOURCES IN THE GLOBULAR CLUSTER NGC 6652

We observed the globular cluster NGC 6652 with Chandra for 47.5 ks, detecting six known X-ray sources, as well as five previously undetected X-ray sources. Source A (XB 1832-330) is a well-known bright low-mass X-ray binary (LXMB). The second brightest source, B, has a spectrum that fits well to either a power-law model ({Gamma} {approx} 1.3) or an absorbed hot gas emission model (kT {approx} 34 keV). Its unabsorbed 0.5-10 keV luminosity (L{sub X} = 1.6{sup +0.1}{sub -0.1} Multiplication-Sign 10{sup 34} erg s{sup -1}) is suggestive of a neutron star primary; however, Source B exhibits unusual variability for an LMXB, varying by over an order of magnitude on timescales of {approx}100 s. Source C's spectrum contains a strong low-energy component below {approx}1 keV. Its spectrum is well fit to a simplified magnetic cataclysmic variable (CV) model, thus the soft component may be explained by a hot polar cap of a magnetic CV. Source D has an average L{sub X} (0.5-10 keV) {approx}9 Multiplication-Sign 10{sup 32} erg s{sup -1} and its spectrum is well fit to a neutron star atmosphere model. This is indicative of a quiescent neutron star LXMB, suggesting Source D may be the third known LMXB in NGC 6652. Source E has L{sub X} (0.5-10 keV) {approx}3 Multiplication-Sign 10{sup 32} erg s{sup -1}, while Source F has L{sub X} (0.5-10 keV) {approx}1 Multiplication-Sign 10{sup 32} erg s{sup -1}. Their relatively hard X-ray spectra are well-fit by power-law or plasma emission models. Five newly detected fainter sources have luminosities between 1 and 5 Multiplication-Sign 10{sup 31} erg s{sup -1}. NGC 6652 has an unusually flat X-ray luminosity function compared to other globular clusters, which may be connected to its extremely high central density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

HERSCHEL PACS OBSERVATIONS AND MODELING OF DEBRIS DISKS IN THE TUCANA-HOROLOGIUM ASSOCIATION

We present Herschel PACS photometry of 17 B- to M-type stars in the 30 Myr old Tucana-Horologium Association. This work is part of the Herschel Open Time Key Programme 'Gas in Protoplanetary Systems'. 6 of the 17 targets were found to have infrared excesses significantly greater than the expected stellar IR fluxes, including a previously unknown disk around HD30051. These six debris disks were fitted with single-temperature blackbody models to estimate the temperatures and abundances of the dust in the systems. For the five stars that show excess emission in the Herschel PACS photometry and also have Spitzer IRS spectra, we fit the data with models of optically thin debris disks with realistic grain properties in order to better estimate the disk parameters. The model is determined by a set of six parameters: surface density index, grain size distribution index, minimum and maximum grain sizes, and the inner and outer radii of the disk. The best-fitting parameters give us constraints on the geometry of the dust in these systems, as well as lower limits to the total dust masses. The HD105 disk was further constrained by fitting marginally resolved PACS 70 {mu}m imaging.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

OBSERVING SIMULATED PROTOSTARS WITH OUTFLOWS: HOW ACCURATE ARE PROTOSTELLAR PROPERTIES INFERRED FROM SEDs?

The properties of unresolved protostars and their local environment are frequently inferred from spectral energy distributions (SEDs) using radiative transfer modeling. In this paper, we use synthetic observations of realistic star formation simulations to evaluate the accuracy of properties inferred from fitting model SEDs to observations. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud including the effects of protostellar outflows. To obtain the dust temperature distribution and SEDs of the forming protostars, we post-process the simulations using HYPERION, a state-of-the-art Monte Carlo radiative transfer code. We find that the ORION and HYPERION dust temperatures typically agree within a factor of two. We compare synthetic SEDs of embedded protostars for a range of evolutionary times, simulation resolutions, aperture sizes, and viewing angles. We demonstrate that complex, asymmetric gas morphology leads to a variety of classifications for individual objects as a function of viewing angle. We derive best-fit source parameters for each SED through comparison with a pre-computed grid of radiative transfer models. While the SED models correctly identify the evolutionary stage of the synthetic sources as embedded protostars, we show that the disk and stellar parameters can be very discrepant from the simulated values, which is expected since the disk and central source are obscured by the protostellar envelope. Parameters such as the stellar accretion rate, stellar mass, and disk mass show better agreement, but can still deviate significantly, and the agreement may in some cases be artificially good due to the limited range of parameters in the set of model SEDs. Lack of correlation between the model and simulation properties in many individual instances cautions against overinterpreting properties inferred from SEDs for unresolved protostellar sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

TYPE Ia SUPERNOVA PROPERTIES AS A FUNCTION OF THE DISTANCE TO THE HOST GALAXY IN THE SDSS-II SN SURVEY

We use Type Ia supernovae (SNe Ia) discovered by the Sloan Digital Sky Survey-II SN Survey to search for dependencies between SN Ia properties and the projected distance to the host-galaxy center, using the distance as a proxy for local galaxy properties (local star formation rate, local metallicity, etc.). The sample consists of almost 200 spectroscopically or photometrically confirmed SNe Ia at redshifts below 0.25. The sample is split into two groups depending on the morphology of the host galaxy. We fit light curves using both MLCS2K2 and SALT2, and determine color (A{sub V} , c) and light-curve shape ({Delta}, x{sub 1}) parameters for each SN Ia, as well as its residual in the Hubble diagram. We then correlate these parameters with both the physical and the normalized distances to the center of the host galaxy and look for trends in the mean values and scatters of these parameters with increasing distance. The most significant (at the 4{sigma} level) finding is that the average fitted A{sub V} from MLCS2K2 and c from SALT2 decrease with the projected distance for SNe Ia in spiral galaxies. We also find indications that supernovae (SNe) in elliptical galaxies tend to have narrower light curves if they explode at larger distances, although this may be due to selection effects in our sample. We do not find strong correlations between the residuals of the distance moduli with respect to the Hubble flow and the galactocentric distances, which indicates a limited correlation between SN magnitudes after standardization and local host metallicity.

79 ASTRONOMY AND ASTROPHYSICS↗