Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance estimation”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Modification of Eccentric Gaze-Holding

Clear vision and accurate localization of objects in the environment are prerequisites for reliable performance of motor tasks. Space flight confronts the crewmember with a stimulus rearrangement that requires adaptation to function effectively with the new requirements of altered spatial orientation and motor coordination. Adaptation and motor learning driven by the effects of cerebellar disorders may share some of the same demands that face our astronauts. One measure of spatial localization shared by the astronauts and those suffering from cerebellar disorders that is easily quantified, and for which a neurobiological substrate has been identified, is the control of the angle of gaze (the "line of sight"). The disturbances of gaze control that have been documented to occur in astronauts and cosmonauts, both in-flight and postflight, can be directly related to changes in the extrinsic gravitational environment and intrinsic proprioceptive mechanisms thus, lending themselves to description by simple non-linear statistical models. Because of the necessity of developing robust normal response populations and normative populations against which abnormal responses can be evaluated, the basic models can be formulated using normal, non-astronaut test subjects and subsequently extended using centrifugation techniques to alter the gravitational and proprioceptive environment of these subjects. Further tests and extensions of the models can be made by studying abnormalities of gaze control in patients with cerebellar disease. A series of investigations were conducted in which a total of 62 subjects were tested to: (1) Define eccentric gaze-holding parameters in a normative population, and (2) explore the effects of linear acceleration on gaze-holding parameters. For these studies gaze-holding was evaluated with the subjects seated upright (the normative values), rolled 45 degrees to both the left and right, or pitched back 30 and 90 degrees. In a separate study the further effects of acceleration on gaze stability was examined during centrifugation (+2 G (sub x) and +2 G (sub z) using a total of 23 subjects. In all of our investigations eccentric gaze-holding was established by having the subjects acquire an eccentric target (+/-30 degrees horizontal, +/- 15 degrees vertical) that was flashed for 750 msec in an otherwise dark room. Subjects were instructed to hold gaze on the remembered position of the flashed target for 20 sec. Immediately following the 20 sec period, subjects were cued to return to the remembered center position and to hold gaze there for an additional 20 sec. Following this 20 sec period the center target was briefly flashed and the subject made any corrective eye movement back to the true center position. Conventionally, the ability to hold eccentric gaze is estimated by fitting the natural log of centripetal eye drifts by linear regression and calculating the time constant (G) of these slow phases of "gaze-evoked nystagmus". However, because our normative subjects sometimes showed essentially no drift (tau (sub c) = m), statistical estimation and inference on the effect of target direction was performed on values of the decay constant theta = 1/(tau (sub c)) which we found was well modeled by a gamma distribution. Subjects showed substantial variance of their eye drifts, which were centrifugal in approximately 20 % of cases, and > 40% for down gaze. Using the ensuing estimated gamma distributions, we were able to conclude that rightward and leftward gaze holding were not significantly different, but that upward gaze holding was significantly worse than downward (p<0.05). We also concluded that vertical gaze holding was significantly worse than horizontal (p<0.05). In the case of left and right roll, we found that both had a similar improvement to horizontal gaze holding (p<0.05), but didn't have a significant effect on vertical gaze holding. For pitch tilts, both tilt angles significantly decreased gaze-holding ility in all directions (p<0.05). Finally, we found that hyper-g centrifugation significantly decreased gaze holding ability in the vertical plane. The main findings of this study are as follows: (1) vertical gaze-holding is less stable than horizontal, (2) gaze-holding to upward targets is less stable than to downward targets, (3) tilt affects gaze holding, and (4) hyper-g affects gaze holding. This difference between horizontal and vertical gaze-holding may be ascribed to separate components of the velocity-to-position neural integrator for eye movements, and to differences in orbital mechanics. The differences between upward and downward gaze-holding may be ascribed to an inherent vertical imbalance in the vestibular system. Because whole body tilt and hyper-g affects gaze-holding, it is implied that the otolith organs have direct connections to the neural integrator and further studies of astronaut gaze-holding are warranted. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied to normals who maybe exposed to environments which may modify the central integrator and require monitoring, and to evaluate patients with gaze-evoked nystagmus by comparing to the above established normative criteria.

Reschke, M. F.↗
DOE OSTI2017

Coupling sample paths to the thermodynamic limit in Monte Carlo estimators with applications to gene expression

Many biochemical systems appearing in applications have a multiscale structure so that they converge to piecewise deterministic Markov processes in a thermodynamic limit. The statistics of the piecewise deterministic process can be obtained much more efficiently than those of the exact process. We explore the possibility of coupling sample paths of the exact model to the piecewise deterministic process in order to reduce the variance of their difference. We then apply this coupling to reduce the computational complexity of a Monte Carlo estimator. Motivated by the rigorous results in , we show how this method can be applied to realistic biological models with nontrivial scalings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2018

Variability in Galactic Cosmic Radiation- Induced DNA Damage Response in Inbred Mice Is Modulated by Genetics

In radiation biology, the ability to predict cancer risk associated with exposure to low doses of high-LET (Linear Energy Transfer) ionizing radiation remains a challenge. Epidemiological methods lack the sensitivity and power to provide detailed risk estimates for cancer and ignore individual sensitivity. We have hypothesized that DNA repair capacity is the primary factor differentiating peoples radiation sensitivity. We previously showed in immortalized human cell lines that characterizing the dose and time dependence of p53-binding protein 1 (53BP1) foci formation in the nucleus following X-rays exposure is sufficient to predict DNA repair response to any other LET in the same cell line. We now tested this hypothesis across a population of mice with different genetic background. Fibroblast cells were extracted and cultivated from 76 individual mice from 15 different strains and exposed to HZE (high (H) atomic number (Z) and energy (E) galactic cosmic ray particles) particles and X-rays. Individual radiation sensitivities were investigated by high throughput measurement of DNA repair kinetics that evaluated 53bp1 foci numbers as a surrogate for DNA double-strand breaks at various times post-irradiation. Instead of just counting foci which can be hard to distinguish for high-LET or high doses, we also took into account the track structure of high-LET particles to compute the remaining number of unrepaired tracks as a function of time post-irradiation. As expected, the percentage of unrepaired track over a 48 hours follow-up period increased with LET. In addition, repair rate was modulated by genetics, with animals from the same strain showing small variance while large rate differences were observed between strains. Radiation strain sensitivity ranking was estimated based on repair rates from exposure to each LET evaluated in this work, and ranking for high-LET correlated better with ranking from high dose of X-ray, not low dose. At the in-vivo level, drops in T-cells and B-cells number measured 24 hours after 0.1 Gy (Gray) X-ray exposure, correlated with slower DNA repair kinetic in fibroblast cells of the same strains of mice. At the genomic level, mouse genome wide association (GWA) analysis identified seven significant genetic loci on chromosomes 2, 3, 7, 10, 11, 13 and 19 with different significance depending on the LET. Interestingly, for the two highest LET, a common locus on Chromosome 10 was identified with high enrichment for DNA repair associated genes.Overall, this work suggests that repair kinetics of primary skin fibroblasts is a good surrogate marker for in-vivo radiation sensitivities in other tissues and that this response is modulated by genetics. Our study also confirms that DNA repair kinetics following high doses of X-ray can be used to predict radiation sensitivity to high-LET.

High-LET Radiation↗
DOE OSTI2021

Accurate Simulation of Both Sensitivity and Variability for Amazonian Photosynthesis: Is It Too Much to Ask?

Abstract Estimates of Amazon rainforest gross primary productivity (GPP) differ by a factor of 2 across a suite of three statistical and 18 process models. This wide spread contributes uncertainty to predictions of future climate. We compare the mean and variance of GPP from these models to that of GPP at six eddy covariance (EC) towers. Only one model's mean GPP across all sites falls within a 99% confidence interval for EC GPP, and only one model matches EC variance. The strength of model response to climate drivers is related to model ability to match the seasonal pattern of the EC GPP. Models with stronger seasonal swings in GPP have stronger responses to rain, light, and temperature than does EC GPP. The model to data comparison illustrates a trade‐off inherent to deterministic models between accurate simulation of a mean (average) and accurate responsiveness to drivers. The trade‐off exists because all deterministic models simplify processes and lack at least some consequential driver or interaction. If a model's sensitivities to included drivers and their interactions are accurate, then deterministically predicted outcomes have less variability than is realistic. If a GPP model has stronger responses to climate drivers than found in data, model predictions may match the observed variance and seasonal pattern but are likely to overpredict GPP response to climate change. High or realistic variability of model estimates relative to reference data indicate that the model is hypersensitive to one or more drivers.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Efficient computation of parameter confidence intervals

An important step in system identification of aircraft is the estimation of stability and control derivatives from flight data along with an assessment of parameter accuracy. When the maximum likelihood estimation technique is used, parameter accuracy is commonly assessed by the Cramer-Rao lower bound. It is known, however, that in some cases the lower bound can be substantially different from the parameter variance. Under these circumstances the Cramer-Rao bounds may be misleading as an accuracy measure. This paper discusses the confidence interval estimation problem based on likelihood ratios, which offers a more general estimate of the error bounds. Four approaches are considered for computing confidence intervals of maximum likelihood parameter estimates. Each approach is applied to real flight data and compared.

Murphy, Patrick C.↗
DOE OSTI2019

Mass variance from archival X-ray properties of Dark Energy Survey Year-1 galaxy clusters

ABSTRACT Using archival X-ray observations and a lognormal population model, we estimate constraints on the intrinsic scatter in halo mass at fixed optical richness for a galaxy cluster sample identified in Dark Energy Survey Year-One (DES-Y1) data with the redMaPPer algorithm. We examine the scaling behaviour of X-ray temperatures, TX, with optical richness, λRM, for clusters in the redshift range 0.2 < z < 0.7. X-ray temperatures are obtained from Chandra and XMM observations for 58 and 110 redMaPPer systems, respectively. Despite non-uniform sky coverage, the TX measurements are $\gt 50{{\ \rm per\ cent}}$ complete for clusters with λRM > 130. Regression analysis on the two samples produces consistent posterior scaling parameters, from which we derive a combined constraint on the residual scatter, $\sigma _{\ln T \, |\, \lambda }= 0.275 \pm 0.019$. Joined with constraints for TX scaling with halo mass from the Weighing the Giants program and richness–temperature covariance estimates from the LoCuSS sample, we derive the richness-conditioned scatter in mass, $\sigma _{\ln M \, |\, \lambda }= 0.30 \pm 0.04\, _{({\rm stat})} \pm 0.09\, _{({\rm sys})}$, at an optical richness of approximately 100. Uncertainties in external parameters, particularly the slope and variance of the TX–mass relation and the covariance of TX and λRM at fixed mass, dominate the systematic error. The $95{{\ \rm per\ cent}}$ confidence region from joint sample analysis is relatively broad, $\sigma _{\ln M \, |\, \lambda }\in [0.14, \, 0.55]$, or a factor 10 in variance.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Spatial Statistical Data Fusion (SSDF)

As remote sensing for scientific purposes has transitioned from an experimental technology to an operational one, the selection of instruments has become more coordinated, so that the scientific community can exploit complementary measurements. However, tech nological and scientific heterogeneity across devices means that the statistical characteristics of the data they collect are different. The challenge addressed here is how to combine heterogeneous remote sensing data sets in a way that yields optimal statistical estimates of the underlying geophysical field, and provides rigorous uncertainty measures for those estimates. Different remote sensing data sets may have different spatial resolutions, different measurement error biases and variances, and other disparate characteristics. A state-of-the-art spatial statistical model was used to relate the true, but not directly observed, geophysical field to noisy, spatial aggregates observed by remote sensing instruments. The spatial covariances of the true field and the covariances of the true field with the observations were modeled. The observations are spatial averages of the true field values, over pixels, with different measurement noise superimposed. A kriging framework is used to infer optimal (minimum mean squared error and unbiased) estimates of the true field at point locations from pixel-level, noisy observations. A key feature of the spatial statistical model is the spatial mixed effects model that underlies it. The approach models the spatial covariance function of the underlying field using linear combinations of basis functions of fixed size. Approaches based on kriging require the inversion of very large spatial covariance matrices, and this is usually done by making simplifying assumptions about spatial covariance structure that simply do not hold for geophysical variables. In contrast, this method does not require these assumptions, and is also computationally much faster. This method is fundamentally different than other approaches to data fusion for remote sensing data because it is inferential rather than merely descriptive. All approaches combine data in a way that minimizes some specified loss function. Most of these are more or less ad hoc criteria based on what looks good to the eye, or some criteria that relate only to the data at hand.

Braverman, Amy J.↗
DOE OSTI2016

MODELING GD-1 GAPS IN A MILKY WAY POTENTIAL

The GD-1 star stream is currently the best available for identifying density fluctuations, “gaps,” along its length as a test of the LCDM prediction of large numbers of dark matter sub-halos orbiting in the halo. Density variations of some form are present, since the variance of the density along the stream is three times that expected from the empirically estimated variation in the filtered mean star counts. The density variations are characterized with filters that approximate the shape of sub-halo, gravitationally induced stream gaps. The filters locate gaps and measure their amplitude, leading to a measurement of the distribution of gap widths. To gain an understanding of the factors influencing the gap width distribution, a suite of collisionless n-body simulations for a GD-1-like orbit in a Milky-Way-like potential provides a dynamically realistic statistical prediction of the gap distribution. The simulations show that every location in the stream has been disturbed to some degree by a sub-halo. The small gaps found via the filtering are largely noise. Larger gaps, those longer than 1 kpc, or 10° for GD-1, are the source of the excess variance. The suite of stream simulations shows that sub-halos at the predicted inner halo abundance or possibly somewhat higher can produce the required large-scale density variations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Validation of a technique for estimating outgoing longwave radiation from HIRS radiance observations

Simultaneous observations by the Earth Radiation Budget Experiment (ERBE) scanning radiometer and the High-Resolution Infrared Sounder (HIRS) on board the NOAA-9 spacecraft have been used to validate a multispectral technique for estimating the outgoing longwave radiation (OLR) from the earth-atmosphere system. Results from approximately 100 000 collocated observations show that the HIRS technique provides instantaneous OLR estimates that agree with the ERBE observations just as well as different ERBE scanners agree with each other--about 5 W m(exp -2) rms. Although there are differences between the HIRS and ERBE estimates that depend upon the scene type and time of day, the HIRS technique explained more than 99% of the variance of the ERBE observations for both day and night observations. The results suggest that the HIRS OLR technique is a suitable replacement for the Advanced Very High Resolution Radiometer technique now used by the National Oceanic and Atmospheric Administration for operational estimates of the OLR.

Ellingson, Robert G.↗
DOE OSTI2019

Mass Variance from Archival X-ray Properties of Dark Energy Survey Year-1 Galaxy Clusters

Using archival X-ray observations and a log-normal population model, we estimate constraints on the intrinsic scatter in halo mass at fixed optical richness for a galaxy cluster sample identified in Dark Energy Survey Year-One (DES-Y1) data with the redMaPPer algorithm. We examine the scaling behavior of X-ray temperatures, $T_X$, with optical richness, $\lambda_{RM}$, for clusters in the redshift range $0.2 50\%$ complete for clusters with $\lambda_{RM} > 130$. Regression analysis on the two samples produces consistent posterior scaling parameters, from which we derive a combined constraint on the residual scatter, $\sigma_{\ln Tx | \lambda} = 0.275 \pm 0.019$. Joined with constraints for $T_X$ scaling with halo mass from the Weighing the Giants program and richness--temperature covariance estimates from the LoCuSS sample, we derive the richness-conditioned scatter in mass, $\sigma_{\ln M | \lambda} = 0.30 \pm 0.04\, _{({\rm stat})} \pm 0.09\, _{({\rm sys})}$, at an optical richness of approximately 70. Uncertainties in external parameters, particularly the slope and variance of the $T_X$--mass relation and the covariance of $T_X$ and $\lambda_{RM}$ at fixed mass, dominate the systematic error. The $95\%$ confidence region from joint sample analysis is relatively broad, $\sigma_{\ln M | \lambda} \in [0.14, \, 0.55]$, or a factor ten in variance.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Adaptive Error Estimation in Linearized Ocean General Circulation Models

Data assimilation methods are routinely used in oceanography. The statistics of the model and measurement errors need to be specified a priori. This study addresses the problem of estimating model and measurement error statistics from observations. We start by testing innovation based methods of adaptive error estimation with low-dimensional models in the North Pacific (5-60 deg N, 132-252 deg E) to TOPEX/POSEIDON (TIP) sea level anomaly data, acoustic tomography data from the ATOC project, and the MIT General Circulation Model (GCM). A reduced state linear model that describes large scale internal (baroclinic) error dynamics is used. The methods are shown to be sensitive to the initial guess for the error statistics and the type of observations. A new off-line approach is developed, the covariance matching approach (CMA), where covariance matrices of model-data residuals are "matched" to their theoretical expectations using familiar least squares methods. This method uses observations directly instead of the innovations sequence and is shown to be related to the MT method and the method of Fu et al. (1993). Twin experiments using the same linearized MIT GCM suggest that altimetric data are ill-suited to the estimation of internal GCM errors, but that such estimates can in theory be obtained using acoustic data. The CMA is then applied to T/P sea level anomaly data and a linearization of a global GFDL GCM which uses two vertical modes. We show that the CMA method can be used with a global model and a global data set, and that the estimates of the error statistics are robust. We show that the fraction of the GCM-T/P residual variance explained by the model error is larger than that derived in Fukumori et al.(1999) with the method of Fu et al.(1993). Most of the model error is explained by the barotropic mode. However, we find that impact of the change in the error statistics on the data assimilation estimates is very small. This is explained by the large representation error, i.e. the dominance of the mesoscale eddies in the T/P signal, which are not part of the 21 by 1" GCM. Therefore, the impact of the observations on the assimilation is very small even after the adjustment of the error statistics. This work demonstrates that simult&neous estimation of the model and measurement error statistics for data assimilation with global ocean data sets and linearized GCMs is possible. However, the error covariance estimation problem is in general highly underdetermined, much more so than the state estimation problem. In other words there exist a very large number of statistical models that can be made consistent with the available data. Therefore, methods for obtaining quantitative error estimates, powerful though they may be, cannot replace physical insight. Used in the right context, as a tool for guiding the choice of a small number of model error parameters, covariance matching can be a useful addition to the repertory of tools available to oceanographers.

Chechelnitsky, Michael Y.↗
DOE OSTI2014

ITAR: A modified TAR method to determine depth dose distribution for an ophthalmic device that performs kilovoltage x-ray pencil-beam stereotaxy

Purpose: New technology has been developed to treat age-related macular degeneration (AMD) using 100 kVp pencil-beams that enter the patient through the radio-resistant sclera with a depth of interest between 1.6 and 2.6 cm. Measurement of reference and relative dose in a kilovoltage x-ray beam with a 0.42 cm diameter field size and a 15 cm source to axis distance (SAD) is a challenge that is not fully addressed in current guidelines to medical physicists. AAPM's TG-61 gives dosimetry recommendations for low and medium energy x-rays, but not all of them are feasible to follow for this modality. Methods: An investigation was conducted to select appropriate equipment for the application. PTW's Type 34013 Soft X-Ray Chamber (Freiburg, Germany) and CIRS's Plastic Water LR (Norfolk, VA) were found to be the best available options. Attenuation curves were measured with minimal scatter contribution and thus called Low Scatter Tissue Air Ratio (LSTAR). A scatter conversion coefficient (C{sub scat}) was derived through Monte Carlo radiation transport simulation using MCNPX (LANL, Los Alamos, NM) to quantify the difference between a traditional TAR curve and the LSTAR curve. A material conversion coefficient (C{sub mat}) was determined through experimentation to evaluate the difference in attenuation properties between water and Plastic Water LR. Validity of performing direct dosimetry measurements with a source to detector distance other than the treatment distance, and therefore a different field size due to a fixed collimator, was explored. A method—Integrated Tissue Air Ratio (ITAR)—has been developed that isolates each of the three main radiological effects (distance from source, attenuation, and scatter) during measurement, and integrates them to determine the dose rate to the macula during treatment. Results: LSTAR curves were determined to be field size independent within the range explored, indicating that direct dosimetry measurements may be performed with a source to detector distance of 20 cm even though the SAD is 15 cm during treatment. C{sub scat} varied from 1.102 to 1.106 within the range of depths of interest. The experimental variance among repeated measurements of C{sub mat} was larger than depth dependence, so C{sub mat} was estimated as1.019 for all depths of interest. Conclusions: Equipment selection, measurement techniques, and formalism for the determination of dose rate to the macula during stereotaxy for AMD have been determined and are strongly recommended by the authors of this paper to be used by clinical medical physicists.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1991

On the sea-state bias of the Geosat altimeter

The sea-state bias in a satellite altimeter's range measurement is caused by the influence of ocean waves on the radar return pulse; it results in an estimate of sea level that is too low according to some function of the wave height. This bias is here estimated for Geosat by correlating collinear differences of altimetric sea-surface heights with collinear differences of significant wave heights (H1/3). Corrections for satellite orbit error are estimated simultaneously with the sea-state bias. Based on twenty 17-day repeat cycles of the Geosat Exact Repeat Mission, the solution for the sea-state bias is 2.6 + or - 0.2 percent of H1/3. The least-squares residuals, however, show a correlation with wind speed U, so the traditional model of the bias has been supplemented with a second term: H1/3 + alpha-2H1/3U. This second term produces a small, but statistically significant, reduction in variance of the residuals. Both systematic and random errors in H1/3 and U tend to bias the estimates of alpha-1 and alpha-2, which complicates comparisons of the results with ground-based measurements of the sea-state bias.

Ray, Richard D.↗
DOE OSTI2022

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗
DOE OSTI2026

A Bayesian inferencing framework for ultrasound wave speed measurements in metal additive manufacturing

Process-related changes during metal additive manufacturing introduce microstructural variability in the material properties of printed parts, directly affecting component reliability. Accurate estimation of these property variations with part performance are essential for quality assurance. Ultrasound testing offers a non-destructive means to estimate mechanical properties and detect defects; however, conventional analysis methods often neglect the influence of microstructural variability, limiting their effectiveness. Here, this research presents a Bayesian inference technique for quantifying wave speed uncertainty from ultrasound measurements of metal additive manufactured parts. By integrating prior ultrasound data with a Bayesian model, the proposed approach generates posterior density estimates of wave speed that systematically account for manufacturing-induced variability and uncertainty. The novelty of this research lies in applying a Bayesian framework to analyze experimental ultrasound measurements within the context of metal additive manufacturing variability. The method enhances the accuracy of wave speed estimation by 64%, defect position by 50% and increases confidence associated with wave speed variance by 30% across different porosity levels, thereby providing a robust foundation for improved decision-making and increased reliability in additively manufactured components.

Additive manufacturing↗
DOE OSTI2015

Diagnosing multiplicative error with lensing magnification of type Ia supernovae

Weak lensing causes spatially coherent fluctuations in the flux of Type Ia supernovae (SNe Ia). This lensing magnification allows for weak lensing measurement independent of cosmic shear. It is free of the shape measurement errors associated with cosmic shear and can therefore be used to diagnose and calibrate multiplicative error. Although this lensing magnification is difficult to accurately measure in auto correlation, its cross correlation with cosmic shear and galaxy distribution in an overlapping area can be measured to a significantly higher accuracy. Therefore, these cross correlations can put useful constraints on multiplicative error, and the obtained constraint is free of cosmic variance in the weak lensing field. We present two methods implementing this idea and estimate their performances. We find that, with ∼1 million SNe Ia that can be achieved with the proposed D2k survey with the LSST telescope, a multiplicative error of ∼0.5% for source galaxies at z{sub s}∼1 can be detected and a larger multiplicative error can be corrected to the level of 0.5%. It is therefore a promising approach to control the multiplicative error to the sub-percent level required for stage IV projects. The combination of the two methods even has the potential to diagnose and calibrate galaxy intrinsic alignment, which is another major systematic error in cosmic shear cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A general reconstruction of the recent expansion history of the universe

Distance measurements are currently the most powerful tool to study the expansion history of the universe without specifying its matter content nor any theory of gravitation. Assuming only an isotropic, homogeneous and flat universe, in this work we introduce a model-independent method to reconstruct directly the deceleration function via a piecewise function. Including a penalty factor, we are able to vary continuously the complexity of the deceleration function from a linear case to an arbitrary (n+1)-knots spline interpolation. We carry out a Monte Carlo (MC) analysis to determine the best penalty factor, evaluating the bias-variance trade-off, given the uncertainties of the SDSS-II and SNLS supernova combined sample (JLA), compilations of baryon acoustic oscillation (BAO) and H(z) data. The bias-variance analysis is done for three fiducial models with different features in the deceleration curve. We perform the MC analysis generating mock catalogs and computing their best-fit. For each fiducial model, we test different reconstructions using, in each case, more than 10{sup 4} catalogs in a total of about 5× 10{sup 5}. This investigation proved to be essential in determining the best reconstruction to study these data. We show that, evaluating a single fiducial model, the conclusions about the bias-variance ratio are misleading. We determine the reconstruction method in which the bias represents at most 10% of the total uncertainty. In all statistical analyses, we fit the coefficients of the deceleration function along with four nuisance parameters of the supernova astrophysical model. For the full sample, we also fit H{sub 0} and the sound horizon r{sub s}(z{sub d}) at the drag redshift. The bias-variance trade-off analysis shows that, apart from the deceleration function, all other estimators are unbiased. Finally, we apply the Ensemble Sampler Markov Chain Monte Carlo (ESMCMC) method to explore the posterior of the deceleration function up to redshift 1.3 (using only JLA) and 2.3 (JLA+BAO+H(z)). We obtain that the standard cosmological model agrees within 3σ level with the reconstructed results in the whole studied redshift intervals. Since our method is calibrated to minimize the bias, the error bars of the reconstructed functions are a good approximation for the total uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗