Search NASA⌕ Search

SEARCH · Search NASA

Results for “density 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 685 records · Page 38

Infrared and X-Ray Spectroscopy of the Kes 75 Supernova Remnant Shell: Characterizing the Dust and Gas Properties

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component with a temperature of approximately 1.5 keY, or with two thermal components with temperatures of 1.5 and 0.2 keY. Previous studies suggest that the hot component may originate from reverse-shocked supernova (SN) ejecta. However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, and Fe, as expected from SN ejecta, or for the IR spectral signatures characteristic of confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR emission. The X-ray and IR emission in the shell are spatially correlated, suggesting that the dust particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows that the dust is heated to a temperature of approximately 140 K by a relatively dense, hot plasma that also gives rise to the hot X-my emission component. The density inferred from the IR emission is significantly higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-my emitting gas. The total mass of the warm dust component is at least 1.3 x 10(exp -2) x solar mass, assuming no significant dust destruction has occurred in the shell. The IR data also reveal the presence of an additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. Our IR analysis therefore provides an independent verification of the cooler component of the X-ray emission. The complementary analyses of the X-ray and IR emission provide quantitative estimates of density and filling factors of the clumpy medium swept up by the SNR.

Temim, Tea↗

INFRARED AND X-RAY SPECTROSCOPY OF THE Kes 75 SUPERNOVA REMNANT SHELL: CHARACTERIZING THE DUST AND GAS PROPERTIES

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component with a temperature of {approx}1.5 keV, or with two thermal components with temperatures of 1.5 and 0.2 keV. Previous studies suggest that the hot component may originate from reverse-shocked supernova (SN) ejecta. However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, and Fe, as expected from SN ejecta, or for the IR spectral signatures characteristic of confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR emission. The X-ray and IR emission in the shell are spatially correlated, suggesting that the dust particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows that the dust is heated to a temperature of {approx}140 K by a relatively dense, hot plasma that also gives rise to the hot X-ray emission component. The density inferred from the IR emission is significantly higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-ray emitting gas. The total mass of the warm dust component is at least 1.3 Multiplication-Sign 10{sup -2} M{sub Sun }, assuming no significant dust destruction has occurred in the shell. The IR data also reveal the presence of an additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. Our IR analysis therefore provides an independent verification of the cooler component of the X-ray emission. The complementary analyses of the X-ray and IR emission provide quantitative estimates of density and filling factors of the clumpy medium swept up by the SNR.

79 ASTRONOMY AND ASTROPHYSICS↗

On the flux and the energy spectrum of interstellar ions in the solar system

The flux density of ions created by ionization of interstellar neutral particles in the solar system and picked up by the solar wind is calculated as a function of the neutral particles. For atomic hydrogen the flux density is estimated to exceed 10,000/sq cm/sec over the distance range from a few to nearly 100 AU. The velocity space distribution of the interstellar ions is calculated under the assumption of no significant energy diffusion but with inclusion of adiabatic effects as well as a possible strong pitch angle diffusion. The energy spectrum is highly nonthermal and much broader than that of the solar wind ions; interstellar protons are easily distinguishable from solar wind protons by their location in velocity space. If charge exchange is an important contributor to the ionization of hydrogen, the observed local intensity of interstellar protons should exhibit time variations correlated with the density changes of the solar wind stream structure.

Vasyliunas, V. M.↗

Multiple Asteroid Systems: Dimensions and Thermal Properties from Spitzer Space Telescope and Ground-based Observations

We collected mid-IR spectra from 5.2 to 38 microns using the Spitzer Space Telescope Infrared Spectrograph of 28 asteroids representative of all established types of binary groups. Photometric light curves were also obtained for 14 of them during the Spitzer observations to provide the context of the observations and reliable estimates of their absolute magnitudes. The extracted mid-IR spectra were analyzed using a modified standard thermal model (STM) and a thermophysical model (TPM) that takes into account the shape and geometry of the large primary at the time of the Spitzer observation. We derived a reliable estimate of the size, albedo, and beaming factor for each of these asteroids, representing three main taxonomic groups: C, S, and X. For large (volume-equivalent system diameter Deq > 130 km) binary asteroids, the TPM analysis indicates a low thermal inertia (Lambda < or = approx.100 J/1/2 s/K/sq m2) and their emissivity spectra display strong mineral features, implying that they are covered with a thick layer of thermally insulating regolith. The smaller (surface-equivalent system diameter Deff < 17 km) asteroids also show some emission lines of minerals, but they are significantly weaker, consistent with regoliths with coarser grains, than those of the large binary asteroids. The average bulk densities of these multiple asteroids vary from 0.7-1.7 g/cu cm (P-, C-type) to approx. 2 g/cu cm (S-type). The highest density is estimated for the M-type (22) Kalliope (3.2 +/- 0.9 g/cu cm). The spectral energy distributions (SEDs) and emissivity spectra, made available as a supplement document, could help to constrain the surface compositions of these asteroids.

Marchis, F.↗

Route to high-energy density polymeric nitrogen t-N via He-N compounds

Polymeric nitrogen, stabilized by compressing pure molecular nitrogen, has yet to be recovered to ambient conditions, precluding its application as a high-energy density material. Here we suggest a route for synthesis of a tetragonal polymeric nitrogen, denoted t-N, via He-N compounds at high pressures. Using first-principles calculations with structure searching, we predict a class of nitrides with stoichiometry HeN 4 that are energetically stable (relative to a mixture of solid He and N 2 ) above 8.5 GPa. At high pressure, HeN 4 comprises a polymeric channel-like nitrogen framework filled with linearly arranged helium atoms. The nitrogen framework persists to ambient pressure on decompression after removal of helium, forming pure polymeric nitrogen, t-N. t-N is dynamically and mechanically stable at ambient pressure with an estimated energy density of ~11.31 kJ/g, marking it out as a remarkable high-energy density material. As a result, this expands the known polymeric forms of nitrogen and indicates a route to its synthesis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

BeZnCdSe quantum-well ridge-waveguide laser diodes under low threshold room-temperature continuous-wave operation

Low threshold current ridge-waveguide BeZnCdSe quantum-well laser diodes (LDs) have been developed by completely etching away the top p-type BeMgZnSe/ZnSe:N short-period superlattice cladding layer, which can suppress the leakage current that flows laterally outside of the electrode. The waveguide LDs are covered with a thick SiO{sub 2} layer and planarized with chemical-mechanical polishing and a reactive ion etching process. Room-temperature lasing under continuous-wave condition is achieved with the laser cavity formed by the cleaved waveguide facets coated with high-reflectivity dielectric films. For a 4 μm-wide green LD lasing around a wavelength of 535 nm, threshold current and voltage of 7.07 mA and 7.89 V are achieved for a cavity length of 300 μm, and the internal differential quantum efficiency, internal absorption loss, gain constant, and nominal transparency current density are estimated to be 27%, 4.09 cm{sup −1}, 29.92 (cm × μm)/kA and 6.35 kA/(cm{sup 2 }× μm), respectively. This compact device can realize a significantly improved performance with much lower threshold power consumption, which would benefit the potential application for ZnSe-based green LDs as light sources in full-color display and projector devices installed in consumer products such as pocket projectors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Monolithically integrated InGaAs/GaAs/AlGaAs quantum well laser grown by MOCVD on exact Ge/Si(001) substrate

We report on realization of the InGaAs/GaAs/AlGaAs quantum well laser grown by metallorganic chemical vapor deposition on a virtual Ge-on-Si(001) substrate. The Ge buffer layer has been grown on a nominal Si(001) substrate by solid-source molecular beam epitaxy. Such Ge buffer possessed rather good crystalline quality and smooth surface and so provided the subsequent growth of the high-quality A{sub 3}B{sub 5} laser structure. The laser operation has been demonstrated under electrical pumping at 77 K in the continuous wave mode and at room temperature in the pulsed mode. The emission wavelengths of 941 nm and 992 nm have been obtained at 77 K and 300 K, respectively. The corresponding threshold current densities were estimated as 463 A/cm{sup 2} at 77 K and 5.5 kA/cm{sup 2} at 300 K.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Atomic scale observation of phase separation and formation of silicon clusters in Hf higk-{kappa} silicates

Hafnium silicate films were fabricated by RF reactive magnetron sputtering technique. Fine microstructural analyses of the films were performed by means of high-resolution transmission electron microscopy and atom probe tomography. A thermal treatment of as-grown homogeneous films leads to a phase separation process. The formation of SiO{sub 2} and HfO{sub 2} phases as well as pure Si one was revealed. This latter was found to be amorphous Si nanoclusters, distributed uniformly in the film volume. Their mean diameter and density were estimated to be about 2.8 nm and (2.9 {+-} 0.4) x 10{sup 17} Si-ncs/cm{sup 3}, respectively. The mechanism of the decomposition process was proposed. The obtained results pave the way for future microelectronic and photonic applications of Hf-based high-{kappa} dielectrics with embedded Si nanoclusters.

36 MATERIALS SCIENCE↗

Long-term Variability of H{sub 2}CO Masers in Star-forming Regions

We present results of a multi-epoch monitoring program on variability of 6 cm formaldehyde (H{sub 2}CO) masers in the massive star-forming region NGC 7538 IRS 1 from 2008 to 2015, conducted with the Green Bank Telescope, the Westerbork Radio Telescope , and the Very Large Array. We found that the similar variability behaviors of the two formaldehyde maser velocity components in NGC 7538 IRS 1 (which was pointed out by Araya and collaborators in 2007) have continued. The possibility that the variability is caused by changes in the maser amplification path in regions with similar morphology and kinematics is discussed. We also observed 12.2 GHz methanol and 22.2 GHz water masers toward NGC 7538 IRS 1. The brightest maser components of CH{sub 3}OH and H{sub 2}O species show a decrease in flux density as a function of time. The brightest H{sub 2}CO maser component also shows a decrease in flux density and has a similar LSR velocity to the brightest H{sub 2}O and 12.2 GHz CH{sub 3}OH masers. The line parameters of radio recombination lines and the 20.17 and 20.97 GHz CH{sub 3}OH transitions in NGC 7538 IRS 1 are also reported. In addition, we observed five other 6 cm formaldehyde maser regions. We found no evidence of significant variability of the 6 cm masers in these regions with respect to previous observations, the only possible exception being the maser in G29.96−0.02. All six sources were also observed in the H{sub 2}{sup 13}CO isotopologue transition of the 6 cm H{sub 2}CO line; H{sub 2}{sup 13}CO absorption was detected in five of the sources. Estimated column density ratios [H{sub 2}{sup 12}CO]/[H{sub 2}{sup 13}CO] are reported.

79 ASTRONOMY AND ASTROPHYSICS↗

An Application of Multivariate Statistical Analysis for Query-Driven Visualization

Driven by the ability to generate ever-larger, increasingly complex data, there is an urgent need in the scientific community for scalable analysis methods that can rapidly identify salient trends in scientific data. Query-Driven Visualization (QDV) strategies are among the small subset of techniques that can address both large and highly complex datasets. This paper extends the utility of QDV strategies with a statistics-based framework that integrates non-parametric distribution estimation techniques with a new segmentation strategy to visually identify statistically significant trends and features within the solution space of a query. In this framework, query distribution estimates help users to interactively explore their query's solution and visually identify the regions where the combined behavior of constrained variables is most important, statistically, to their inquiry. Our new segmentation strategy extends the distribution estimation analysis by visually conveying the individual importance of each variable to these regions of high statistical significance. We demonstrate the analysis benefits these two strategies provide and show how they may be used to facilitate the refinement of constraints over variables expressed in a user's query. We apply our method to datasets from two different scientific domains to demonstrate its broad applicability.

97 MATHEMATICS AND COMPUTING↗

The resolved outflow from 3C 48

We investigate the properties of the high-velocity outflow driven by the young radio jet of 3C 48, a compact-steep-spectrum source. We use the Space Telescope Imaging Spectrograph on board the Hubble Space Telecope to obtain (1) low-resolution UV and optical spectra and (2) multi-slit medium-resolution spectra of the ionized outflow. With supporting data from ground-based spectrographs, we are able to accurately measure the ratios of diagnostic emission lines such as [O III] λ5007, [O III] λ3727, [N II] λ6548, Hα, Hβ, [Ne V] λ3425, and [Ne III] λ3869. We fit the observed emission-line ratios using a range of ionization models, powered by active galactic nucleus (AGN) radiation and shocks, produced by the MAPPINGS code. We have determined that AGN radiation is likely the dominant ionization source. The outflow's density is estimated to be in the range n = 10{sup 3}-10{sup 4} cm{sup –3}, the mass is ∼6 × 10{sup 6} M {sub ☉}, and the metallicity is likely equal to or higher than solar. Compared with the typical outflows associated with more evolved radio jets, this young outflow is denser, less massive, and more metal rich. Multi-slit observations allow us to construct a two-dimensional velocity map of the outflow that shows a wide range of velocities with distinct velocity components, suggesting a wide-angle clumpy outflow.

79 ASTRONOMY AND ASTROPHYSICS↗

Hot Gas in the Wolf–Rayet Nebula NGC 3199

The Wolf–Rayet (WR) nebula NGC 3199 has been suggested to be a bow shock around its central star, WR 18, which is presumably a runaway star, because optical images of the nebula show a dominating arc of emission southwest of the star. We present the XMM-Newton detection of extended X-ray emission from NGC 3199, unveiling the powerful effect of the fast wind from WR 18. The X-ray emission is brighter in the region southeast of the star and an analysis of the spectral properties of the X-ray emission reveals abundance variations: (i) regions close to the optical arc present nitrogen-rich gas enhanced by the stellar wind from WR 18 and (ii) gas at the eastern region exhibits abundances close to those reported for the nebular abundances derived from optical studies, which is a signature of an efficient mixing of the nebular material with the stellar wind. The dominant plasma temperature and electron density are estimated to be T ≈ 1.2 × 10{sup 6} K and n {sub e} = 0.3 cm{sup −3} with an X-ray luminosity in the 0.3–3.0 keV energy range of L {sub X} = 2.6 × 10{sup 34} erg s{sup −1}. Combined with information derived from Herschel and the recent Gaia first data release, we conclude that WR 18 is not a runaway star and that the formation, chemical variations, and the shape of NGC 3199 depend on the initial configuration of the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

GIANT MOLECULAR CLOUD EVOLUTIONS IN THE NEARBY SPIRAL GALAXY M33

We present a giant molecular cloud (GMC) catalog of M33, containing 71 GMCs in total, based on wide-field and high-sensitivity CO(J = 3-2) observations with a spatial resolution of 100 pc using the ASTE 10 m telescope. Employing archival optical data, we identify 75 young stellar groups (YSGs) from the excess of the surface stellar density, and estimate their ages by comparing with stellar evolution models. A spatial comparison among the GMCs, YSGs, and H II regions enable us to classify GMCs into four categories: Type A, showing no sign of massive star formation (SF); Type B, being associated only with H II regions; Type C, with both H II regions and <10 Myr old YSGs; and Type D, with both H II regions and 10-30 Myr YSGs. Out of 65 GMCs (discarding those at the edges of the observed fields), 1 (1%), 13 (20%), 29 (45%), and 22 (34%) are Types A, B, C, and D, respectively. We interpret these categories as stages in a GMC evolutionary sequence. Assuming that the timescale for each evolutionary stage is proportional to the number of GMCs, the lifetime of a GMC with a mass >10{sup 5} M{sub Sun} is estimated to be 20-40 Myr. In addition, we find that the dense gas fraction as traced by the CO(J = 3-2)/CO(J = 1-0) ratio is enhanced around SF regions. This confirms a scenario where dense gas is preferentially formed around previously generated stars, and will be the fuel for the next stellar generation. In this way, massive SF gradually propagates in a GMC until gas is exhausted.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of Computed Tomography Techniques for Material Identification in Low Level and Intermediate Level Waste - 16648

Measures required to assure long term safe storage of Low-level and intermediate-level short-lived (LIL/SL) waste generated by the nuclear-power industry in France include the characterization of the waste and identification of materials within the waste. Digital Radiography (DR) techniques using x-ray at energies in the 6 MeV range have been shown to be effective for validation of the integrity of the concrete containers, and it was considered that the use of this technique together with the addition of Computed Tomography (CT) techniques could effectively identify materials within the waste. This paper describes the initial evaluation of two DR imaging techniques and two CT techniques using x-ray energy spectra up to 9 MeV to inspect a surrogate waste package, which comprised an outer concrete cylindrical container of 1400 mm diameter with 148 mm wall thickness, and an overall volume of 2 meters cubed, within which a 200 liter drum containing representative heterogeneous inactive waste, was placed in the center and the annulus between the 200 liter drum and the outer cylindrical package was filled with cement grout. The surrogate waste within the 200 liter drum comprised three separate layers of high density, medium density and low density mixed materials. The equipment used was a linear accelerator producing high energy spectra at either 6 MeV or 9 MeV and a high-precision CT manipulator with X, Y, Z and rotation motions, which were common to all techniques investigated. Two different digital detectors were used in the DR and CT performance evaluations, a Digital Detector Array (DDA) comprising a columnar Thallium-doped Cesium-Iodide scintillator screen in intimate contact with a 2 k by 2 k array of amorphous silicon detectors, 200 μm x 200 μm in area, and a Linear Diode Array (LDA) comprising a line of Cadmium-Tungstate columns, 1 cm thick and 400 μm x 400 μm in area. The work comprised evaluations of linear accelerator spectra and the attenuation effects of the waste container and contents, image acquisition performance of each detector, and CT techniques using area detector and linear detector. The results presented in this paper include mathematical evaluation of the two CT techniques, and analysis of measurement data, digital images, CT reconstruction images, and estimated material densities. The data was compared with recorded details and photographs of the actual configuration of materials within the surrogate waste. The Digital Detector Array can generate CT inspection volumes with considerable resolution and detail for the low-density and medium-density sections of the surrogate waste, however, the high-density sections were not imaged with sufficient detail and the images contained a large number of artifacts mainly due to the detectors lack of high energy photon stopping power. Materials identification with the Digital Detector Array, while possible for certain materials, for materials in the steel to plastic density range the attenuation values of the voxels in the reconstruction are too similar and difficult to resolve from the artifacts present. The Linear Diode Array can generate CT inspection volumes of regions of interest of the low-density and medium-density sections with considerable resolution and contrast, and potentially acceptable levels of resolution and contrast of the high-density section of the waste. While a few artifacts were present in the reconstructions, the majority are not due to the fundamental properties of the detector and are considered to be correctable. Materials identification with the LDA inspection data appear to be viable, with steel, copper, and higher-density metals being readily be distinguished from other materials on the basis of voxel value. Future work is needed to obtain a better definition of optimum spatial resolution and inspection speed, and to develop a system of 'calibration rods' which would consist of a series of 50 mm diameter cylinders of 7 materials; Plastic, Teflon, Concrete, Aluminum, Copper, Tungsten and Lead. The cylinders will be distributed at a certain radius in the surrogate waste and will be used as reference points in the evaluation and demonstration of the contrastive performance for the low, medium and high-density cases. The ability to identify the materials within a sealed waste container, by nondestructive examination, enables the verification of waste acceptance criteria without the need to open the waste containers. This eliminates the need for complex and expensive equipment to open and remotely handle the intermediate level waste, avoids potential hazards and risk of personnel dose uptake, and does not increase the volume of waste caused by repackaging of materials after examination and the decontamination or disposal of equipment used.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain-Corrected Satellite Sea Surface Salinity Using the GMOA S2s Forecast System

The ENSO phenomenon has a significant global socio-economic impact and has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help improve the initialization of the thermocline, while satellite SST aids in constraining surface heat-fluxes, leading to improved coupled system sub-seasonal to seasonal forecasts. However, few studies have focused on improving the near-surface density and mixing through satellite sea surface salinity (SSS) assimilation. The few ocean models that assimilate satellite SSS, bias correct to normalize towards the near-surface Argo data for expediency. This assumption is likely inadequate in rainy regions, where buoyant water forms a fresh surface lens. In previous work, we showed that adjusting SSS to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 improves the near-surface density and mixed layer depth, leading to deeper thermocline and improved NINO3.4 SST forecasts. We now utilize the Soil Moisture and Ocean Salinity rain-corrected (SMOS_RC) SSS, available in SMOS-CATDS products, to represent Sb more accurately at the first model layer (e.g., 5 m). Rather than a diffusivity model as RIM, SMOS_RC uses a statistical correction dependent on Integrated Multi-satellitE Retrievals for GPM (IMERG) rain rates, established on observed SMOS SSS decreases related to Sb in the presence of rain (Supply et al., 2020). For all experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is, and a separate reanalysis assimilates SMOS_RC. We assess the impact on near-surface and subsurface dynamics by validating against observations and explore how SSS assimilation (SMOS vs SMOS_RC) impacts ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (S2S-v3, Molod et al. 2020). We show that improved estimates of density and near-surface mixing led to more accurate coupled air/sea interaction and better ENSO forecasts. The increased SSS, resulting from the removal of the instantaneous rain effect, modifies the ocean state by enhancing mixing and deepening the thermocline.

Veronica Ruiz Xomchuk↗

Minimum entropy filtering for a single output non-Gaussian stochastic system using state transformation

This paper presents a novel filter design for the single-output stochastic non-linear systems subjected to non-Gaussian noises and the proposed assumptions. Based on a state transformation, the unmeasurable states of the systems can be estimated where non-linear terms in the systems have been eliminated. It has been shown that the estimation error is linearly dynamical regarding to the presented vector-valued filter gain which can be optimised by minimising the entropy-based performance criterion. In addition, the convergence of the presented algorithm is analysed in mean-square sense and a numerical example is given to verify the effectiveness of the presented filtering algorithm. Meanwhile, the extended Kalman filter, unscented particle filter and minimum entropy filter are given for the comparisons of the filtering performance. Following the presented framework, some extensions of the presented filtering algorithm are discussed to indicate the flexibility of the filter design. The contribution of this paper can be summarised as establishing a novel minimum entropy filtering framework which consists of model transformation, entropy optimisation and convergence analysis.

42 ENGINEERING↗

Estimation of distributions via multilevel Monte Carlo with stratified sampling

We design and implement a novel algorithm for computing a multilevel Monte Carlo (MLMC) estimator of the joint cumulative distribution function (CDF) of a vector-valued quantity of interest in problems with random input parameters and initial conditions. Our approach combines MLMC with stratified sampling of the input sample space by replacing standard Monte Carlo at each level with stratified Monte Carlo initialized with proportionally allocated samples. We show that the resulting stratified MLMC (sMLMC) algorithm is more efficient than its standard MLMC counterpart due to the additional variance reduction provided by the stratification of the random parameter's domain, especially at the coarsest levels. Additional computational cost savings are obtained by smoothing the indicator function with a Gaussian kernel, which proves to be an efficient and robust alternative to recently developed polynomial-based techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determinants of the reliability of ultrasound tomography sound speed estimates as a surrogate for volumetric breast density

Purpose: High breast density, as measured by mammography, is associated with increased breast cancer risk, but standard methods of assessment have limitations including 2D representation of breast tissue, distortion due to breast compression, and use of ionizing radiation. Ultrasound tomography (UST) is a novel imaging method that averts these limitations and uses sound speed measures rather than x-ray imaging to estimate breast density. The authors evaluated the reproducibility of measures of speed of sound and changes in this parameter using UST. Methods: One experienced and five newly trained raters measured sound speed in serial UST scans for 22 women (two scans per person) to assess inter-rater reliability. Intrarater reliability was assessed for four raters. A random effects model was used to calculate the percent variation in sound speed and change in sound speed attributable to subject, scan, rater, and repeat reads. The authors estimated the intraclass correlation coefficients (ICCs) for these measures based on data from the authors’ experienced rater. Results: Median (range) time between baseline and follow-up UST scans was five (1–13) months. Contributions of factors to sound speed variance were differences between subjects (86.0%), baseline versus follow-up scans (7.5%), inter-rater evaluations (1.1%), and intrarater reproducibility (∼0%). When evaluating change in sound speed between scans, 2.7% and ∼0% of variation were attributed to inter- and intrarater variation, respectively. For the experienced rater’s repeat reads, agreement for sound speed was excellent (ICC = 93.4%) and for change in sound speed substantial (ICC = 70.4%), indicating very good reproducibility of these measures. Conclusions: UST provided highly reproducible sound speed measurements, which reflect breast density, suggesting that UST has utility in sensitively assessing change in density.

60 APPLIED LIFE SCIENCES↗