Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large deviations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measurement of the $t\bar{t} t\bar{t}$ production cross section in pp collisions at $ \sqrt{s} $ = 13 TeV with the ATLAS detector

A measurement of four-top-quark production using proton-proton collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS detector at the Large Hadron Collider corresponding to an integrated luminosity of 139 fb -1 is presented. Events are selected if they contain a single lepton (electron or muon) or an opposite-sign lepton pair, in association with multiple jets. The events are categorised according to the number of jets and how likely these are to contain b-hadrons. A multivariate technique is then used to discriminate between signal and background events. The measured four-top-quark production cross section is found to be 26$^{+17}_{-15}$ fb, with a corresponding observed (expected) significance of 1.9 (1.0) standard deviations over the background-only hypothesis. The result is combined with the previous measurement performed by the ATLAS Collaboration in the multilepton final state. The combined four-top-quark production cross section is measured to be 24$^{+7}_{-6}$ fb, with a corresponding observed (expected) signal significance of 4.7 (2.6) standard deviations over the background-only predictions. It is consistent within 2.0 standard deviations with the Standard Model expectation of 12.0 ± 2.4 fb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Molecular Vibration Induced Plasmon Decay

Noble metal nanoparticles, when interacting with an external electric field, give rise to a phenomenon called surface plasmon resonance characterized by collective valence electron oscillations, and because of this unique feature these systems are employed for a large and heterogeneous number of applications. To improve their performance, it is necessary to develop a deep knowledge of the main factors affecting the plasmon lifetime. In order to answer this question, in this work a linear silver chain is investigated as a simplified model for silver nanorods. Through a non-adiabatic molecular dynamics approach the role of nuclear dynamics on transverse plasmon lifetime is investigated. A strong dependence of plasmon dynamics on the specific nature of nuclear motions is found: nuclear motions along the chain do not affect the transverse plasmon lifetime while motions causing a deviation from linearity of the wire have an important impact on the plasmon dynamics causing its decay. As the vibrational energy increases, the decay becomes faster because of an accelerated loss of symmetry and shows a weak Landau-like mechanism. The unveiled molecular nature of the plasmon decay on a linear wire can be representative of possible decay mechanisms taking place in larger systems. Finally, these results improve the knowledge of plasmon dynamics, and can be helpful for an efficient design of more performing materials.

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

Role of Bottlebrush Additives on the Structure of Block Copolymers in the Bulk and Thin Films

Blending block copolymers (BCP) with additives is a useful approach for controlling BCP morphology and properties. In athermal systems, blends of BCPs with polymer additives having very high molecular (M n ) mass generally result in macrophase separation. Bottlebrush polymers, which consist of a linear backbone and grafted side chains, present an interesting alternative where the overall M n of the system can be very large but the low M n side chains may drive miscibility with the BCP. Here, in this study, a bottlebrush with a polynorbornene backbone and polystyrene (PS) side chains is blended with PS-b-poly(methyl methacrylate) (PS-b-PMMA) of varying M n , and the resulting morphologies are examined in both the bulk and thin films. Two different M n of PS-b-PMMA were used in the bulk study, and the analysis of small-angle X-ray scattering data shows that the blends were miscible and lamellar at all concentrations. This deviates from reference series of both low and high M n linear polymer additives, which either showed morphological transitions from lamellae to cylinders (low M n ) or were immiscible at all mass fractions studied (high M n ). The relative molecular mass of the side chain (N SC ) and the corresponding component in the BCP (N A ) dictate the distribution of the bottlebrush throughout the BCP, analogous to BCP/linear blends or grafted nanoparticles in a homopolymer matrix. The studies on thin films show a thickness dependence for bottlebrush mass fractions at or above 0.17, a behavior which may be driven by conformational changes of the bottlebrush upon confinement.

36 MATERIALS SCIENCE↗
DOE OSTI2013

INTERPRETING THE ROLE OF THE MAGNETIC FIELD FROM DUST POLARIZATION MAPS

Dust polarization observations from the Submillimeter Array (SMA) and the Caltech Submillimeter Observatory (CSO) are analyzed with the goal of providing a general tool to interpret the role of the magnetic field in molecular clouds. Magnetic field and dust emission gradient orientations are observed to show distinct patterns and features. The angle δ between these two orientations can be interpreted as a magnetic field alignment deviation, assuming the emission gradient orientation to coincide with the density gradient orientation in the magnetohydrodynamics force equation. In SMA high-resolution (collapsing) cores, additional symmetry properties in δ can reveal accretion and outflow zones. All these observational findings suggest the angle δ to be a relevant quantity that can assess the role of the magnetic field. Indeed, when comparing this angle with the (projection-free) magnetic field significance Σ{sub B} (introduced by Koch and coworkers in 2012), it is demonstrated that |δ| yields an approximation to the change in Σ{sub B}. Thus, changes in the magnetic field alignment deviation δ trace changes in the role of the magnetic field. The angle δ is observationally straightforward to determine, providing a tool to distinguish between zones of minor or significant magnetic field impact. This is exemplified by the CSO M+0.25 + 0.01, Mon R2, CO+0.02 – 0.02, M–0.02 – 0.07 sources and by the SMA high-resolution data from W51 e2, W51 North, Orion BN/KL and g5.89. Additional CSO sources are analyzed, providing further support of this result. Finally, based on the different features found in our sample of 31 sources in total, covering sizes from large-scale complexes to collapsing cores, a schematic evolutionary scenario is proposed. Here, the significance of the magnetic field evolves both with position and scale, and can be assessed with the angle δ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE CO-TO-H{sub 2} CONVERSION FACTOR AND DUST-TO-GAS RATIO ON KILOPARSEC SCALES IN NEARBY GALAXIES

We present ∼kiloparsec spatial resolution maps of the CO-to-H{sub 2} conversion factor (α{sub CO}) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for α{sub CO} and the DGR by assuming that the DGR is approximately constant on kiloparsec scales. With this assumption, we can combine maps of dust mass surface density, CO-integrated intensity, and H I column density to solve for both α{sub CO} and the DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high-resolution far-IR maps from the Herschel key program KINGFISH, {sup 12}CO J = (2-1) maps from the IRAM 30 m large program HERACLES, and H I 21 cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our α{sub CO} results on the more typically used {sup 12}CO J = (1-0) scale and show using literature measurements that variations in the line ratio do not affect our results. In total, we derive 782 individual solutions for α{sub CO} and the DGR. On average, α{sub CO} = 3.1 M{sub ☉} pc{sup –2} (K km s{sup –1}){sup –1} for our sample with a standard deviation of 0.3 dex. Within galaxies, we observe a generally flat profile of α{sub CO} as a function of galactocentric radius. However, most galaxies exhibit a lower α{sub CO} value in the central kiloparsec—a factor of ∼2 below the galaxy mean, on average. In some cases, the central α{sub CO} value can be factors of 5-10 below the standard Milky Way (MW) value of α{sub CO,{sub MW}} = 4.4 M{sub ☉} pc{sup –2} (K km s{sup –1}){sup –1}. While for α{sub CO} we find only weak correlations with metallicity, the DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate α{sub CO} for studies of nearby galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Ages of the Thin Disk, Thick Disk, and the Halo from Nearby White Dwarfs

We present a detailed analysis of the white dwarf luminosity functions derived from the local 40 pc sample and the deep proper motion catalog of Munn et al. Many previous studies have ignored the contribution of thick disk white dwarfs to the Galactic disk luminosity function, which results in an erroneous age measurement. We demonstrate that the ratio of thick/thin disk white dwarfs is roughly 20% in the local sample. Simultaneously fitting for both disk components, we derive ages of 6.8–7.0 Gyr for the thin disk and 8.7 ± 0.1 Gyr for the thick disk from the local 40 pc sample. Similarly, we derive ages of 7.4–8.2 Gyr for the thin disk and 9.5–9.9 Gyr for the thick disk from the deep proper motion catalog, which shows no evidence of a deviation from a constant star formation rate in the past 2.5 Gyr. We constrain the time difference between the onset of star formation in the thin disk and the thick disk to be 1.6{sub −0.4}{sup +0.3} Gyr. The faint end of the luminosity function for the halo white dwarfs is less constrained, resulting in an age estimate of 12.5{sub −3.4}{sup +1.4} Gyr for the Galactic inner halo. This is the first time that ages for all three major components of the Galaxy have been obtained from a sample of field white dwarfs that is large enough to contain significant numbers of disk and halo objects. The resultant ages agree reasonably well with the age estimates for the oldest open and globular clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Evaluation of observed and modelled aerosol lifetimes using radioactive tracers of opportunity and an ensemble of 19 global models

Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 ( 137 Cs) and xenon-133 ( 133 Xe) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. 137 Cs size distribution measurements taken close to the power plant suggested that accumulation-mode (AM) sulfate aerosols were the main carriers of cesium. Hence, 137 Cs can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas 133 Xe behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of 137 Cs that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and 133 Xe emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled 137 Cs and 133 Xe concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime τ e , calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95 % confidence interval 13.1–15.7 days). The equivalent modelled τ e lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.4 ± 2.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer ( 133 Xe) concentrations in the Arctic as well but to a smaller extent than for the aerosol ( 137 Cs) tracer. As a result, this indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multi-century dynamics of the climate and carbon cycle under both high and net negative emissions scenarios

Abstract. Future climate projections from Earth system models (ESMs) typically focus on the timescale of this century. We use a set of five ESMs and one Earth system model of intermediate complexity (EMIC) to explore the dynamics of the Earth's climate and carbon cycles under contrasting emissions trajectories beyond this century to the year 2300. The trajectories include a very-high-emissions, unmitigated fossil-fuel-driven scenario, as well as a mitigation scenario that diverges from the first scenario after 2040 and features an “overshoot”, followed by a decrease in atmospheric CO2 concentrations by means of large net negative CO2 emissions. In both scenarios and for all models considered here, the terrestrial system switches from being a net sink to either a neutral state or a net source of carbon, though for different reasons and centered in different geographic regions, depending on both the model and the scenario. The ocean carbon system remains a sink, albeit weakened by carbon cycle feedbacks, in all models under the high-emissions scenario and switches from sink to source in the overshoot scenario. The global mean temperature anomaly is generally proportional to cumulative carbon emissions, with a deviation from proportionality in the overshoot scenario that is governed by the zero emissions commitment. Additionally, 23rd century warming continues after the cessation of carbon emissions in several models in the high-emissions scenario and in one model in the overshoot scenario. While ocean carbon cycle responses qualitatively agree in both globally integrated and zonal mean dynamics in both scenarios, the land models qualitatively disagree in zonal mean dynamics, in the relative roles of vegetation and soil in driving C fluxes, in the response of the sink to CO2, and in the timing of the sink–source transition, particularly in the high-emissions scenario. The lack of agreement among land models on the mechanisms and geographic patterns of carbon cycle feedbacks, alongside the potential for lagged physical climate dynamics to cause warming long after CO2 concentrations have stabilized, points to the possibility of surprises in the climate system beyond the 21st century time horizon, even under relatively mitigated global warming scenarios, which should be taken into consideration when setting global climate policy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Statistical Techniques For Real-time Anomaly Detection Using Spark Over Multi-source VMware Performance Data

Anomaly detection refers to the identi cation of an irregular or unusual pat- tern which deviates from what is standard, normal, or expected. Such deviated patterns typically correspond to samples of interest and are assigned different labels in different domains, such as outliers, anomalies, exceptions, or malware. Detecting anomalies in fast, voluminous streams of data is a formidable chal- lenge. This paper presents a novel, generic, real-time distributed anomaly detection framework for heterogeneous streaming data where anomalies appear as a group. We have developed a distributed statistical approach to build a model and later use it to detect anomaly. As a case study, we investigate group anomaly de- tection for a VMware-based cloud data center, which maintains a large number of virtual machines (VMs). We have built our framework using Apache Spark to get higher throughput and lower data processing time on streaming data. We have developed a window-based statistical anomaly detection technique to detect anomalies that appear sporadically. We then relaxed this constraint with higher accuracy by implementing a cluster-based technique to detect sporadic and continuous anomalies. We conclude that our cluster-based technique out- performs other statistical techniques with higher accuracy and lower processing time.

Apache Spark↗
DOE OSTI2017

An experimental investigation of gas fuel injection with X-ray radiography

In this paper, an outward-opening compressed natural gas, direct injection fuel injector has been studied with single-shot x-ray radiography. Three dimensional simulations have also been performed to compliment the x-ray data. Argon was used as a surrogate gas for experimental and safety reasons. This technique allows the acquisition of a quantitative mapping of the ensemble-average and standard deviation of the projected density throughout the injection event. Two dimensional, ensemble average and standard deviation data are presented to investigate the quasi-steady-state behavior of the jet. Upstream of the stagnation zone, minimal shot-to-shot variation is observed. Downstream of the stagnation zone, bulk mixing is observed as the jet transitions to a subsonic turbulent jet. From the time averaged data, individual slices at all downstream locations are extracted and an Abel inversion was performed to compute the radial density distribution, which was interpolated to create three dimensional visualizations. The Abel reconstructions reveal that upstream of the stagnation zone, the gas forms an annulus with high argon density and large density gradients. Inside this annulus, a recirculation region with low argon density exists. Downstream, the jet transitions to a fully turbulent jet with Gaussian argon density distributions. This experimental data is intended to serve as a quantitative benchmark for simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Combined Recipe for Clinical Target Volume and Planning Target Volume Margins

Purpose: To develop a combined recipe for clinical target volume (CTV) and planning target volume (PTV) margins. Methods and Materials: A widely accepted PTV margin recipe is M{sub geo} = aΣ{sub geo} + bσ{sub geo}, with Σ{sub geo} and σ{sub geo} standard deviations (SDs) representing systematic and random geometric uncertainties, respectively. On the basis of histopathology data of breast and lung tumors, we suggest describing the distribution of microscopic islets around the gross tumor volume (GTV) by a half-Gaussian with SD Σ{sub micro}, yielding as possible CTV margin recipe: M{sub micro} = ƒ(N{sub i}) × Σ{sub micro}, with N{sub i} the average number of microscopic islets per patient. To determine ƒ(N{sub i}), a computer model was developed that simulated radiation therapy of a spherical GTV with isotropic distribution of microscopic disease in a large group of virtual patients. The minimal margin that yielded D{sub min} <95% in maximally 10% of patients was calculated for various Σ{sub micro} and N{sub i}. Because Σ{sub micro} is independent of Σ{sub geo}, we propose they should be added quadratically, yielding for a combined GTV-to-PTV margin recipe: M{sub GTV-PTV} = √([aΣ{sub geo}]{sup 2} + [ƒ(N{sub i})Σ{sub micro}]{sup 2}) + bσ{sub geo}. This was validated by the computer model through numerous simultaneous simulations of microscopic and geometric uncertainties. Results: The margin factor ƒ(N{sub i}) in a relevant range of Σ{sub micro} and N{sub i} can be given by: ƒ(N{sub i}) = 1.4 + 0.8log(N{sub i}). Filling in the other factors found in our simulations (a = 2.1 and b = 0.8) yields for the combined recipe: M{sub GTV-PTV} = √((2.1Σ{sub geo}){sup 2} + ([1.4 + 0.8log(N{sub i})] × Σ{sub micro}){sup 2}) + 0.8σ{sub geo}. The average margin difference between the simultaneous simulations and the above recipe was 0.2 ± 0.8 mm (1 SD). Calculating M{sub geo} and M{sub micro} separately and adding them linearly overestimated PTVs by on average 5 mm. Margin recipes based on tumor control probability (TCP) instead of D{sub min} criteria yielded similar results. Conclusions: A general recipe for GTV-to-PTV margins is proposed, which shows that CTV and PTV margins should be added in quadrature instead of linearly.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Quantifying the performance of in vivo portal dosimetry in detecting four types of treatment parameter variations

Purpose: To quantify the ability of electronic portal imaging device (EPID) dosimetry used during treatment (in vivo) in detecting variations that can occur in the course of patient treatment. Methods: Images of transmitted radiation from in vivo EPID measurements were converted to a 2D planar dose at isocenter and compared to the treatment planning dose using a prototype software system. Using the treatment planning system (TPS), four different types of variability were modeled: overall dose scaling, shifting the positions of the multileaf collimator (MLC) leaves, shifting of the patient position, and changes in the patient body contour. The gamma pass rate was calculated for the modified and unmodified plans and used to construct a receiver operator characteristic (ROC) curve to assess the detectability of the different parameter variations. The detectability is given by the area under the ROC curve (AUC). The TPS was also used to calculate the impact of the variations on the target dose–volume histogram. Results: Nine intensity modulation radiation therapy plans were measured for four different anatomical sites consisting of 70 separate fields. Results show that in vivo EPID dosimetry was most sensitive to variations in the machine output, AUC = 0.70 − 0.94, changes in patient body habitus, AUC = 0.67 − 0.88, and systematic shifts in the MLC bank positions, AUC = 0.59 − 0.82. These deviations are expected to have a relatively small clinical impact [planning target volume (PTV) D{sub 99} change <7%]. Larger variations have even higher detectability. Displacements in the patient’s position and random variations in MLC leaf positions were not readily detectable, AUC < 0.64. The D{sub 99} of the PTV changed by up to 57% for the patient position shifts considered here. Conclusions: In vivo EPID dosimetry is able to detect relatively small variations in overall dose, systematic shifts of the MLC’s, and changes in the patient habitus. Shifts in the patient’s position which can introduce large changes in the target dose coverage were not readily detected.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Searches for large-scale anisotropy in the arrival directions of cosmic rays detected above energy of 10{sup 19} eV at the Pierre Auger observatory and the telescope array

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above 10{sup 19} eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above 10{sup 19} eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗
DOE OSTI2011

Gogny HFB prediction of nuclear structure properties

Large scale mean field calculations from proton to neutron drip lines have been performed using the Hartree-Fock-Bogoliubov method based on the Gogny nucleon-nucleon effective interaction. This extensive study has shown the ability of the method to reproduce bulk nuclear structure data available experimentally. This includes nuclear masses, radii, matter densities, deformations, moment of inertia as well as collective mode (low energy and giant resonances). In particular, the first mass table based on a Gogny-Hartree-Fock-Bogolyubov calculation including an explicit and coherent account of all the quadrupole correlation energies is presented. The rms deviation with respect to essentially all the available mass data is 798 keV. Nearly 8000 nuclei have been studied under the axial symmetry hypothesis and going beyond the mean-field approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Laser deposition of large-area thin films

A new method for fabricating large-area thin films of uniform thickness on a rotating substrate is proposed. Its distinctive features are (i) the presence of a diaphragm, partially transmitting the evaporated material, between the target and substrate and (ii) the translatory motion of the rotating substrate with respect to the target at a certain velocity. The method proposed makes it possible to obtain thin films of uniform thickness on substrates with sizes limited by only the deposition chamber size. The method is experimentally verified by depositing thin CuO films on silicon substrates placed over the radius of a disk 300 mm in diameter. The deviation of the film thickness from the average value does not exceed {+-}3% throughout the entire radius, which confirms good prospects of this method for microelectronics, optical industry, and other modern technologies. (laser technology)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Residual foreground contamination in the WMAP data and bias in non-Gaussianity estimation

We analyze whether there is any residual foreground contamination in the cleaned WMAP 7 years data for the differential assemblies, Q, V and W. We calculate the correlation between the foreground map, from which long wavelength correlations have been subtracted, and the foreground reduced map for each differential assembly after applying the Galaxy and point sources masks. We find positive correlations for all the differential assemblies, with high statistical significance. For Q and V, we find that a large fraction of the contamination comes from pixels where the foreground maps have positive values larger than three times the rms values. These findings imply the presence of residual contamination from Galactic emissions and unresolved point sources. We redo the analysis after masking the extended point sources cataloque of Scodeller et al. [7] and find a drop in the correlation and corresponding significance values. To quantify the effect of the residual contamination on the search for primordial non-Gaussianity in the CMB we add estimated contaminant fraction to simulated Gaussian CMB maps and calculate the characteristic non-Gaussian deviation shapes of Minkowski Functionals that arise due to the contamination. We find remarkable agreement of these deviation shapes with those measured from WMAP data, which imply that a major fraction of the observed non-Gaussian deviation comes from residual foreground contamination. We also compute non-Gaussian deviations of Minkowski Functionals after applying the point sources mask of Scodeller et al. and find a decrease in the overall amplitudes of the deviations which is consistent with a decrease in the level of contamination.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-J-131: Reproducibility of Positioning Error Due to Temporarily Indwelled Urethral Catheter for Urethra-Sparing Prostate IMRT

Purpose: The purpose of this study was to prospectively assess the reproducibility of positioning errors due to temporarily indwelled catheter in urethra-sparing image-guided (IG) IMRT. Methods: Ten patients received urethra-sparing prostate IG-IMRT with implanted fiducials. After the first CT scan was performed in supine position, 6-Fr catheter was indwelled into urethra, and the second CT images were taken for planning. While the PTV received 80 Gy, 5% dose reduction was applied for the urethral PRV along the catheter. Additional CT scans were also performed at 5th and 30th fraction. Positions of interests (POIs) were set on posterior edge of prostate at beam isocenter level (POI1) and cranial and caudal edge of prostatic urethra on the post-indwelled CT images. POIs were copied into the pre-indwelled, 5th and 30th fraction’s CT images after fiducial matching on these CT images. The deviation of each POI between pre- and post-indwelled CT and the reproducibility of prostate displacement due to catheter were evaluated. Results: The deviation of POI1 caused by the indwelled catheter to the directions of RL/AP/SI (mm) was 0.20±0.27/−0.64±2.43/1.02±2.31, respectively, and the absolute distances (mm) were 3.15±1.41. The deviation tends to be larger if closer to the caudal edge of prostate. Compared with the pre-indwelled CT scan, a median displacement of all POIs (mm) were 0.3±0.2/2.2±1.1/2.0±2.6 in the post-indwelled, 0.4±0.4/3.4±2.1/2.3±2.6 in 5th, and 0.5±0.5/1.7±2.2/1.9±3.1 in 30th fraction’s CT scan with a similar data distribution. There were 6 patients with 5-mm-over displacement in AP and/or CC directions. Conclusion: Reproducibility of positioning errors due to temporarily indwelling catheter was observed. Especially in case of patients with unusually large shifts by indwelling catheter at the planning process, treatment planning should be performed by using the pre-indwelled CT images with transferred contour of the urethra identified by post-indwelled CT images.

60 APPLIED LIFE SCIENCES↗