Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scattering”

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 613 records · Page 34

Rotationally inelastic scattering of NO(A{sup 2}Σ{sup +}) + Ar: Differential cross sections and rotational angular momentum polarization

We present the implementation of a new crossed-molecular beam, velocity-map ion-imaging apparatus, optimized for collisions of electronically excited molecules. We have applied this apparatus to rotational energy transfer in NO(A{sup 2}Σ{sup +}, v = 0, N = 0, j = 0.5) + Ar collisions, at an average energy of 525 cm{sup −1}. We report differential cross sections for scattering into NO(A{sup 2}Σ{sup +}, v = 0, N′ = 3, 5, 6, 7, 8, and 9), together with quantum scattering calculations of the differential cross sections and angle dependent rotational alignment. The differential cross sections show dramatic forward scattered peaks, together with oscillatory behavior at larger scattering angles, while the rotational alignment moments are also found to oscillate as a function of scattering angle. In general, the quantum scattering calculations are found to agree well with experiment, reproducing the forward scattering and oscillatory behavior at larger scattering angles. Analysis of the quantum scattering calculations as a function of total rotational angular momentum indicates that the forward scattering peak originates from the attractive minimum in the potential energy surface at the N-end of the NO. Deviations in the quantum scattering predictions from the experimental results, for scattering at angles greater than 10°, are observed to be more significant for scattering to odd final N′. We suggest that this represents inaccuracies in the potential energy surface, and in particular in its representation of the difference between the N- and O-ends of the molecule, as given by the odd-order Legendre moments of the surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

WE-DE-207B-09: Scatter Radiation Measurement From a Digital Breast Tomosynthesis System and Its Impact On Shielding Consideration

Purpose: To accurately measure the scatter radiation from a Hologic digital breast tomosynthesis (DBT) system and to provide updated scatter distribution to guide radiation shielding calculation for DBT rooms. Methods: A high sensitivity GOS-based linear detector was used to measure the angular distribution of scatter radiation from a Hologic Selenia Dimensions DBT system. The linear detector was calibrated for its energy response of typical DBT spectra. Following the NCRP147 approach, the measured scatter intensity was normalized by the primary beam area and primary air kerma at 1m from the scatter phantom center and presented as the scatter fraction. Direct comparison was made against Simpkin’s initial measurement. Key parameters including the phantom size, primary beam area, and kV/anode/target combination were also studied. Results: The measured scatter-to-primary-ratio and scatter fraction data closely matched with previous data from Simpkin. The measured data demonstrated the unique nonisotropic distribution of the scattered radiation around a Hologic DBT system, with two strong peaks around 25° and 160°. The majority scatter radiation (>70%) originated from the imaging detector assembly, instead of the phantom. With a workload from a previous local survey, the scatter air kerma at 1m from the phantom center for wall/door is 0.018mGy/patient, for floor is 0.164mGy/patient, and for ceiling is 0.037mGy/patient. Conclusion: Comparing to Simpkin’s previous data, the scatter air kerma from Holgoic DBT is at least two times higher. The main reasons include the harder primary beam with higher workload, added tomosynthesis acquisition, and strong small angle forward scattering. Due to the highly conservative initial assumptions, the shielding recommendation from NCRP147 is still sufficient for the Hologic DBT system given the workload from a previous local survey. With the data provided from this study, accurate shielding calculation can be performed for Hologic DBT systems with specific workload and barrier distance.

60 APPLIED LIFE SCIENCES↗

Comparison of Generation of Higher-Order Neutron Scattering Cross Sections

The generation of high-order neutron scattering cross sections consistent with high-fidelity simulations remains an area of active research. Popular options include generating cross sections from continuous energy Monte Carlo calculations or generating cross sections from a deterministic neutron transport calculation with high-fidelity tabulated cross sections. Both options present challenges. Monte Carlo simulations can naturally process continuous energy cross section data and allow for the general description of anisotropic neutron scattering given a scattering law. However, Monte Carlo simulations are inherently related to particle weighting and it has been suggested that this may be unacceptable for generating high-order neutron scattering cross sections. Deterministic neutron transport calculations can easily calculate high-order moments of the flux to appropriately calculate higher-order neutron scattering cross sections but are generally limited by discretization of space, energy, and angle. In this work, the trade-offs between generation of high-order neutron scattering cross sections via Monte Carlo and deterministic neutron transport methods are investigated. The methods implemented in the Monte Carlo computer program Serpent and the deterministic fast reactor neutron cross section generator MC2-3 are compared. Cross sections resulting from these methods are used in Rattlesnake, a deterministic neutron transport code developed by Idaho National Laboratory (INL), and results are compared to a reference continuous energy Monte Carlo calculation. Whereas previous work investigating the effects of anisotropic neutron scattering has focused on light water reactor (LWR) simulations, this work focuses on high-order neutron scattering cross sections as they relate to fast reactor simulations. To investigate the consequences of the Serpent and MC2-3 methodologies, a test problem is developed. The test problem is a one-dimensional geometry with fast reactor materials designed to demonstrate deep penetration and exacerbate the effects of high-order neutron scattering. Based on the results of the deep penetration test problem, it is concluded that P3 neutron scattering cross sections are sufficient to describe anisotropic scattering in fast reactor materials. Any neutron scattering of order higher than P3 offers negligible change in the eigenvalue of the test problem. Additionally, it is determined that the methodology as implemented in Serpent is applicable for generating high-order neutron scattering cross sections through at least P3 in fast reactor materials.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

TH-A-18C-04: Ultrafast Cone-Beam CT Scatter Correction with GPU-Based Monte Carlo Simulation

Purpose: Scatter artifacts severely degrade image quality of cone-beam CT (CBCT). We present an ultrafast scatter correction framework by using GPU-based Monte Carlo (MC) simulation and prior patient CT image, aiming at automatically finish the whole process including both scatter correction and reconstructions within 30 seconds. Methods: The method consists of six steps: 1) FDK reconstruction using raw projection data; 2) Rigid Registration of planning CT to the FDK results; 3) MC scatter calculation at sparse view angles using the planning CT; 4) Interpolation of the calculated scatter signals to other angles; 5) Removal of scatter from the raw projections; 6) FDK reconstruction using the scatter-corrected projections. In addition to using GPU to accelerate MC photon simulations, we also use a small number of photons and a down-sampled CT image in simulation to further reduce computation time. A novel denoising algorithm is used to eliminate MC scatter noise caused by low photon numbers. The method is validated on head-and-neck cases with simulated and clinical data. Results: We have studied impacts of photo histories, volume down sampling factors on the accuracy of scatter estimation. The Fourier analysis was conducted to show that scatter images calculated at 31 angles are sufficient to restore those at all angles with <0.1% error. For the simulated case with a resolution of 512×512×100, we simulated 10M photons per angle. The total computation time is 23.77 seconds on a Nvidia GTX Titan GPU. The scatter-induced shading/cupping artifacts are substantially reduced, and the average HU error of a region-of-interest is reduced from 75.9 to 19.0 HU. Similar results were found for a real patient case. Conclusion: A practical ultrafast MC-based CBCT scatter correction scheme is developed. The whole process of scatter correction and reconstruction is accomplished within 30 seconds. This study is supported in part by NIH (1R01CA154747-01), The Core Technology Research in Strategic Emerging Industry, Guangdong, China (2011A081402003)

60 APPLIED LIFE SCIENCES↗

WE-G-204-06: Grid-Line Artifact Minimization for High Resolution Detectors Using Iterative Residual Scatter Correction

Purpose: Anti-scatter grid-line artifacts are more prominent for high-resolution x-ray detectors since the fraction of a pixel blocked by the grid septa is large. Direct logarithmic subtraction of the artifact pattern is limited by residual scattered radiation and we investigate an iterative method for scatter correction. Methods: A stationary Smit-Rοntgen anti-scatter grid was used with a high resolution Dexela 1207 CMOS X-ray detector (75 µm pixel size) to image an artery block (Nuclear Associates, Model 76-705) placed within a uniform head equivalent phantom as the scattering source. The image of the phantom was divided by a flat-field image obtained without scatter but with the grid to eliminate grid-line artifacts. Constant scatter values were subtracted from the phantom image before dividing by the averaged flat-field-with-grid image. The standard deviation of pixel values for a fixed region of the resultant images with different subtracted scatter values provided a measure of the remaining grid-line artifacts. Results: A plot of the standard deviation of image pixel values versus the subtracted scatter value shows that the image structure noise reaches a minimum before going up again as the scatter value is increased. This minimum corresponds to a minimization of the grid-line artifacts as demonstrated in line profile plots obtained through each of the images perpendicular to the grid lines. Artifact-free images of the artery block were obtained with the optimal scatter value obtained by this iterative approach. Conclusion: Residual scatter subtraction can provide improved grid-line artifact elimination when using the flat-field with grid “subtraction” technique. The standard deviation of image pixel values can be used to determine the optimal scatter value to subtract to obtain a minimization of grid line artifacts with high resolution x-ray imaging detectors. This study was supported by NIH Grant R01EB002873 and an equipment grant from Toshiba Medical Systems Corp.

60 APPLIED LIFE SCIENCES↗

TH-AB-209-10: Breast Cancer Identification Through X-Ray Coherent Scatter Spectral Imaging

Purpose: We have previously described the development and testing of a coherent-scatter spectral imaging system for identification of cancer. Our prior evaluations were performed using either tissue surrogate phantoms or formalin-fixed tissue obtained from pathology. Here we present the first results from a scatter imaging study using fresh breast tumor tissues obtained through surgical excision. Methods: A coherent-scatter imaging system was built using a clinical X-ray tube, photon counting detectors, and custom-designed coded-apertures. System performance was characterized using calibration phantoms of biological materials. Fresh breast tumors were obtained from patients undergoing mastectomy and lumpectomy surgeries for breast cancer. Each specimen was vacuum-sealed, scanned using the scatter imaging system, and then sent to pathology for histological workup. Scatter images were generated separately for each tissue specimen and analyzed to identify voxels containing malignant tissue. The images were compared against histological analysis (H&E + pathologist identification of tumors) to assess the match between scatter-based and histological diagnosis. Results: In all specimens scanned, the scatter images showed the location of cancerous regions within the specimen. The detection and classification was performed through automated spectral matching without the need for manual intervention. The scatter spectra corresponding to cancer tissue were found to be in agreement with those reported in literature. Inter-patient variability was found to be within limits reported in literature. The scatter images showed agreement with pathologist-identified regions of cancer. Spatial resolution for this configuration of the scanner was determined to be 2–3 mm, and the total scan time for each specimen was under 15 minutes. Conclusion: This work demonstrates the utility of coherent scatter imaging in identifying cancer based on the scatter properties of the tissue. It presents the first results from coherent scatter imaging of fresh (unfixed) breast tissue using our coded-aperture scatter imaging approach for cancer identification.

60 APPLIED LIFE SCIENCES↗

TH-CD-207A-08: Simulated Real-Time Image Guidance for Lung SBRT Patients Using Scatter Imaging

Purpose: To develop a comprehensive Monte Carlo-based model for the acquisition of scatter images of patient anatomy in real-time, during lung SBRT treatment. Methods: During SBRT treatment, images of patient anatomy can be acquired from scattered radiation. To rigorously examine the utility of scatter images for image guidance, a model is developed using MCNP code to simulate scatter images of phantoms and lung cancer patients. The model is validated by comparing experimental and simulated images of phantoms of different complexity. The differentiation between tissue types is investigated by imaging objects of known compositions (water, lung, and bone equivalent). A lung tumor phantom, simulating materials and geometry encountered during lung SBRT treatments, is used to investigate image noise properties for various quantities of delivered radiation (monitor units(MU)). Patient scatter images are simulated using the validated simulation model. 4DCT patient data is converted to an MCNP input geometry accounting for different tissue composition and densities. Lung tumor phantom images acquired with decreasing imaging time (decreasing MU) are used to model the expected noise amplitude in patient scatter images, producing realistic simulated patient scatter images with varying temporal resolution. Results: Image intensity in simulated and experimental scatter images of tissue equivalent objects (water, lung, bone) match within the uncertainty (∼3%). Lung tumor phantom images agree as well. Specifically, tumor-to-lung contrast matches within the uncertainty. The addition of random noise approximating quantum noise in experimental images to simulated patient images shows that scatter images of lung tumors can provide images in as fast as 0.5 seconds with CNR∼2.7. Conclusions: A scatter imaging simulation model is developed and validated using experimental phantom scatter images. Following validation, lung cancer patient scatter images are simulated. These simulated patient images demonstrate the clinical utility of scatter imaging for real-time tumor tracking during lung SBRT.

60 APPLIED LIFE SCIENCES↗

Optimization-based scatter estimation using primary modulation for computed tomography

Purpose: Scatter reduces the image quality in computed tomography (CT), but scatter correction remains a challenge. A previously proposed primary modulation method simultaneously obtains the primary and scatter in a single scan. However, separating the scatter and primary in primary modulation is challenging because it is an underdetermined problem. In this study, an optimization-based scatter estimation (OSE) algorithm is proposed to estimate and correct scatter. Methods: In the concept of primary modulation, the primary is modulated, but the scatter remains smooth by inserting a modulator between the x-ray source and the object. In the proposed algorithm, an objective function is designed for separating the scatter and primary. Prior knowledge is incorporated in the optimization-based framework to improve the accuracy of the estimation: (1) the primary is always positive; (2) the primary is locally smooth and the scatter is smooth; (3) the location of penumbra can be determined; and (4) the scatter-contaminated data provide knowledge about which part is smooth. Results: The simulation study shows that the edge-preserving weighting in OSE improves the estimation accuracy near the object boundary. Simulation study also demonstrates that OSE outperforms the two existing primary modulation algorithms for most regions of interest in terms of the CT number accuracy and noise. The proposed method was tested on a clinical cone beam CT, demonstrating that OSE corrects the scatter even when the modulator is not accurately registered. Conclusions: The proposed OSE algorithm improves the robustness and accuracy in scatter estimation and correction. This method is promising for scatter correction of various kinds of x-ray imaging modalities, such as x-ray radiography, cone beam CT, and the fourth-generation CT.

60 APPLIED LIFE SCIENCES↗

On Spectral Invariance of Single Scattering Albedo for Weakly Absorbing Wavelengths

The single scattering albedo omega (sub 0 lambda) in atmospheric radiative transfer is the ratio of the scattering coefficient to the total extinction coefficient. For cloud water droplets both the scattering and absorption coefficients, thus the single scattering albedo, are functions of wavelength A and droplet size r. In this presentation we will show that for water droplets at weakly absorbing wavelengths, the ratio omega (sub 0 lambda)(r). The slope and intercept of the linear function are wavelength independent and sum to unity. This relationship allows for a representation of any single scattering albedo omega (sub 0 lambda) via one known spectrum omega (sub 0 lambda)(r(sub o)). We will provide a simple physical explanation of the discovered relationship. In addition to water droplets, similar linear relationships were found for the single scattering albedo of non-spherical ice crystals. The single scattering albedo $\omega _ {0\lambda }$ in atmospheric radiative transfer is the ratio of the scattering coefficient to the total extinction coefficient. For cloud water droplets both the scattering and absorption coefficients, and thus the single scattering albedo, are functions of wavelength $\lambda $ and droplet size $r$. We show that for water droplets at weakly absorbing wavelengths, the ratio $\omega _ {0\lambda } (r)$/$\omega _ {0\lambda } (r_{0})$ of two single scattering albedo spectra for two different droplet sizes is a linear function of $\omega _{0\lambda }(r)$. The slope and intercept of the linear function are wavelength independent and sum to unity. This relationship allows for a representation of any single scattering albedo $\omega_{0\lambda }(r)$ via one known spectrum $\omega_{0\lambda }(r_{0})$. We provide a simple physical explanation of the discovered relationship. Similar linear relationships characterize the single scattering albedo of non-spherical ice crystals.

Marshak, Alexander↗

Scattering and Diffraction of Electromagnetic Radiation: An Effective Probe to Material Structure

Scattered electromagnetic waves from material bodies of different forms contain, in an intricate way, precise information on the intrinsic, geometrical and physical properties of the objects. Scattering theories, ever deepening, aim to provide dependable interpretation and prediction to the complicated interaction of electromagnetic radiation with matter. There are well-established multiple-scattering formulations based on classical electromagnetic theories. An example is the Generalized Multi-particle Mie-solution (GMM), which has recently been extended to a special version ̶ the GMM-PA approach, applicable to finite periodic arrays consisting of a huge number (e.g., >>106) of identical scattering centers [1]. The framework of the GMM-PA is nearly complete. When the size of the constituent unit scatterers becomes considerably small in comparison with incident wavelength, an appropriate array of such small element volumes may well be a satisfactory representation of a material entity having an arbitrary structure. X-ray diffraction is a powerful characterization tool used in a variety of scientific and technical fields, including material science. A diffraction pattern is nothing more than the spatial distribution of scattered intensity, determined by the distribution of scattering matter by way of its Fourier transform [1]. Since all linear dimensions entered into Maxwell's equations are normalized by wavelength, an analogy exists between optical and X-ray diffraction patterns. A large set of optical diffraction patterns experimentally obtained can be found in the literature [e.g., 2,3]. Theoretical results from the GMM-PA have been scrutinized using a large collection of publically accessible, experimentally obtained Fraunhofer diffraction patterns. As far as characteristic structures of the patterns are concerned, theoretical and experimental results are in uniform agreement; no exception has been found so far. Closely connected with the spatial distribution of scattered intensities are cross sections, such as for extinction, scattering, absorption, and radiation pressure, as a critical type of key quantity addressed in most theoretical and experimental studies of radiative scattering. Cross sections predicted from different scattering theories are supposed to be in general agreement. For objects of irregular shape, the GMM-PA solutions can be compared with the highly flexible Discrete Dipole Approximation (DDA) [4,5] when dividing a target to no more than ~106 unit cells. Also, there are different ways to calculate the cross sections in the GMM-PA, providing an additional means to examine the accuracy of the numerical solutions and to unveil potential issues concerning the theoretical formulations and numerical aspects. To solve multiple scattering by an assembly of material volumes through classical theories such as the GMM-PA, the radiative properties of the component scatterers, the complex refractive index in particular, must be provided as input parameters. When using a PA to characterize a material body, this involves the use of an adequate theoretical tool, an effective medium theory, to connect Maxwell's phenomenogical theory with the atomistic theory of matter. In the atomic theory, one regards matter as composed of interacting particles (atoms and molecules) embedded in the vacuum [6]. However, the radiative properties of atomic-scaled particles are known to be substantially different from bulk materials. Intensive research efforts in the fields of cluster science and nanoscience attempt to bridge the gap between bulk and atom and to understand the transition from classical to quantum physics. The GMM-PA calculations, which place virtually no restriction on the component-particle size, might help to gain certain insight into the transition.

Xu, Yu-Lin↗

Implementation of A New Microwave Scattering Database and A Forward Model for Active Microwave Sensors in CRTM

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. CRTM is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk scattering lookup tables in order to perform all-sky RT calculations. However, the current CRTM lookup tables for microwave frequencies were generated based on the Mie theory by assuming spherical frozen particles. The scattering lookup tables generated using the DDA technique has shown to largely improve the RT scattering calculations in the MW region. This presentation targets (i) the implementation and validation of a DDA database that was originally developed for the ARTS RT model into CRTM, and (ii) developing the CRTM active sensor module that takes advantage of the backscattering coefficients computed using the DDA method. The DDA database only provides single scattering properties of different habits, while CRTM requires bulk scattering properties. The CRTM cloud coefficients were previously generated based on the effective radius for representing the size of the particles. However, effective radius is neither measurable nor provided by the NWP models, thus need to be estimated from other geophysical variables such as water content. Therefore, in addition to calculating the CRTM bulk scattering properties from the DDA single scattering database, the CRTM was also largely modified to use cloud water content (kg.m-3), instead of effective radius, for performing the interpolation over size/mass of the particles. CRTM already requires water content as input, thus no extra variables are required for performing scattering calculations using the new ARTS DDA database. The CRTM scattering modules search for effective radius in cloud coefficient files and will use the cloud water content if the effective radius dimension is not found in the cloud coefficient files. Figure 1 shows the CRTM simulated brightness temperatures computed using different cloud coefficients versus ATMS observed values over Hurricane Irma on September 7, 2017 at 18:00 UTC. We used all the cloud water content values included in ERA5 with default CRTM/DDA habits for water, rain, snow, ice, hail, and graupel. ERA5 does not provide separate water content values for ice, hail, and graupel, thus the ice water content values were divided between ice, hail, and graupel clouds similar to what was explained in the previous section. In channels with a frequency lower than 90 GHz, emission from water and rain clouds can compensate for cloud scattering so that cloud contaminated Tbs are larger than corresponding clear sky Tbs. The DDA simulations for channels 1-7 largely perform better than the Mie simulations. The DDA simulations show a mix of small negative and positive simulated minus observed values, while the Mie results show large negative biases. The weighting functions for some of the ATMS temperature sounding channels (channels 9-15) peak mostly above the clouds, therefore the measured Tbs become less sensitive to clouds so that the results of both Mie and DDA become very similar. The Mie lookup tables generate excessive scattering for channel 16, but not enough scattering for the water vapor channels. In the specific case of Hurricane Maria, the DDA lookup tables do not generate enough scattering for channel 16, but the DDA results are much more consistent with observations for water vapor channels than for channel 16. It should be noted that the results may vary if we use other habits to represent snow, hail, and graupel in the DDA simulations. Although these results clearly show the advantage of the DDA database over the Mie dataset, different error sources such as error in the observations, displacement of clouds in the ERA5 reanalysis, and also lack of convective clouds or in general errors in the input atmospheric and cloud profiles contribute to the differences between the simulated and observed values. Aside from the improvements in the simulations, a major advantage of the new dataset is a large number of habits that can be used to tune the data assimilation systems to perform well in different weather conditions.

Isaac Moradi↗

A Theory of Exoplanet Transits with Light Scattering

Exoplanet transit spectroscopy enables the characterization of distant worlds, and will yield key results for NASA's James Webb Space Telescope . However, transit spectra models are often simplified, omitting potentially important processes like refraction and multiple scattering. While the former process has seen recent development, the effects of light multiple scattering on exoplanet transit spectra have received little attention. Here, we develop a detailed theory of exoplanet transit spectroscopy that extends to the full refracting and multiple scattering case. We explore the importance of scattering for planet-wide cloud layers, where the relevant parameters are the slant scattering optical depth, the scattering asymmetry parameter, and the angular size of the host star. The latter determines the size of the “target” for a photon that is back-mapped from an observer. We provide results that straightforwardly indicate the potential importance of multiple scattering for transit spectra. When the orbital distance is smaller than 10–20 times the stellar radius, multiple scattering effects for aerosols with asymmetry parameters larger than 0.8–0.9 can become significant. We provide examples of the impacts of cloud/haze multiple scattering on transit spectra of a hot Jupiter-like exoplanet. For cases with a forward and conservatively scattering cloud/haze, differences due to multiple scattering effects can exceed 200 ppm, but shrink to zero at wavelength ranges corresponding to strong gas absorption or when the slant optical depth of the cloud exceeds several tens. We conclude with a discussion of types of aerosols for which multiple scattering in transit spectra may be important.

79 ASTRONOMY AND ASTROPHYSICS↗

Distributions of off-diagonal scattering matrix elements: Exact results

Scattering is a ubiquitous phenomenon which is observed in a variety of physical systems which span a wide range of length scales. The scattering matrix is the key quantity which provides a complete description of the scattering process. The universal features of scattering in chaotic systems is most generally modeled by the Heidelberg approach which introduces stochasticity to the scattering matrix at the level of the Hamiltonian describing the scattering center. The statistics of the scattering matrix is obtained by averaging over the ensemble of random Hamiltonians of appropriate symmetry. We derive exact results for the distributions of the real and imaginary parts of the off-diagonal scattering matrix elements applicable to orthogonally-invariant and unitarily-invariant Hamiltonians, thereby solving a long standing problem. -- Highlights: •Scattering problem in complex or chaotic systems. •Heidelberg approach to model the chaotic nature of the scattering center. •A novel route to the nonlinear sigma model based on the characteristic function. •Exact results for the distributions of off-diagonal scattering-matrix elements. •Universal aspects of the scattering-matrix fluctuations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A global study of hygroscopicity-driven light-scattering enhancement in the context of other in situ aerosol optical properties

The scattering and backscattering enhancement factors (ƒ(RH) and ƒ b (RH)) describe how aerosol particle light scattering and backscattering, respectively, change with relative humidity (RH). They are important parameters in estimating direct aerosol radiative forcing (DARF). In this study we use the dataset presented in Burgos et al. (2019) that compiles ƒ(RH) and ƒ b (RH) measurements at three wavelengths (i.e., 450, 550 and 700 nm) performed with tandem nephelometer systems at multiple sites around the world. We present an overview of ƒ(RH) and ƒ b (RH) based on both long-term and campaign observations from 23 sites representing a range of aerosol types. The scattering enhancement shows a strong variability from site to site, with no clear pattern with respect to the total scattering coefficient. In general, higher ƒ(RH) is observed at Arctic and marine sites, while lower values are found at urban and desert sites, although a consistent pattern as a function of site type is not observed. The backscattering enhancement ƒ b (RH) is consistently lower than ƒ(RH) at all sites, with the difference between ƒ(RH) and ƒ b (RH) increasing for aerosol with higher ƒ(RH). This is consistent with Mie theory, which predicts higher enhancement of the light scattering in the forward than in the backward direction as the particle takes up water. Our results show that the scattering enhancement is higher for PM 1 than PM 10 at most sites, which is also supported by theory due to the change in scattering efficiency with the size parameter that relates particle size and the wavelength of incident light. At marine-influenced sites this difference is enhanced when coarse particles (likely sea salt) predominate. For most sites, ƒ(RH) is observed to increase with increasing wavelength, except at sites with a known dust influence where the spectral dependence of ƒ(RH) is found to be low or even exhibit the opposite pattern. The impact of RH on aerosol properties used to calculate radiative forcing (e.g., single-scattering albedo, ω 0 , and backscattered fraction, b) is evaluated. The single-scattering albedo generally increases with RH, while b decreases. The net effect of aerosol hygroscopicity on radiative forcing efficiency (RFE) is an increase in the absolute forcing effect (negative sign) by a factor of up to 4 at RH = 90 % compared to dry conditions (RH < 40 %). Because of the scarcity of scattering enhancement measurements, an attempt was made to use other more commonly available aerosol parameters (i.e., ω 0 and scattering Ångström exponent, α sp ) to parameterize ƒ(RH). The majority of sites (75 %) showed a consistent trend with ω 0 (higher ƒ(RH = 85 %) for higher ω 0 ), while no clear pattern was observed between ƒ(RH = 85 %) and αsp. This suggests that aerosol ω 0 is more promising than α sp as a surrogate for the scattering enhancement factor, although neither parameter is ideal. Nonetheless, the qualitative relationship observed between ω 0 and ƒ(RH) could serve as a constraint on global model simulations.

54 ENVIRONMENTAL SCIENCES↗

On-the-fly generation of differential resonance scattering probability distribution functions for Monte Carlo codes

Current Monte Carlo codes use one of three models to model neutron scattering in the epithermal energy range: (1) the asymptotic scattering model, (2) the free gas scattering model, or (3) the S({alpha},{beta}) model, depending on the neutron energy and the specific Monte Carlo code. The free gas scattering model assumes the scattering cross section is constant over the neutron energy range, which is usually a good approximation for light nuclei, but not for heavy nuclei where the scattering cross section may have several resonances in the epithermal region. Several researchers in the field have shown that using the free gas scattering model in the vicinity of the resonances in the lower epithermal range can under-predict resonance absorption due to the up-scattering phenomenon. Existing methods all involve performing the collision analysis in the center-of-mass frame, followed by a conversion back to the laboratory frame. In this paper, we will present a new sampling methodology that (1) accounts for the energy-dependent scattering cross sections in the collision analysis and (2) acts in the laboratory frame, avoiding the conversion to the center-of-mass frame. The energy dependence of the scattering cross section was modeled with even-ordered polynomials to approximate the scattering cross section in Blackshaw's equations for the moments of the differential scattering PDFs. These moments were used to sample the outgoing neutron speed and angle in the laboratory frame on-the-fly during the random walk of the neutron. Results for criticality studies on fuel pin and fuel assembly calculations using these methods showed very close comparison to results using the reference Doppler-broadened rejection correction (DBRC) scheme. (authors)

97 MATHEMATICS AND COMPUTING↗

Constituent Components of Out-of-Field Scatter Dose for 18-MV Intensity Modulated Radiation Therapy Versus 3-Dimensional Conformal Radiation Therapy: A Comparison With 6-MV and Implications for Carcinogenesis

Purpose: To characterize and compare the components of out-of-field dose for 18-MV intensity modulated radiation therapy (IMRT) versus 3-dimensional conformal radiation therapy (3D-CRT) and their 6-MV counterparts and consider implications for second cancer induction. Methods and Materials: Comparable plans for each technique/energy were delivered to a water phantom with a sloping wall; under full scatter conditions; with field edge abutting but outside the bath to prevent internal/phantom scatter; and with shielding below the linear accelerator head to attenuate head leakage. Neutron measurements were obtained from published studies. Results: Eighteen-megavolt IMRT produces 1.7 times more out-of-field scatter than 18-MV 3D-CRT. In absolute terms, however, differences are just approximately 0.1% of central axis dose. Eighteen-megavolt IMRT reduces internal/patient scatter by 13%, but collimator scatter (C) is 2.6 times greater than 18-MV 3D-CRT. Head leakage (L) is minimal. Increased out-of-field photon scatter from 18-MV IMRT carries out-of-field second cancer risks of approximately 0.2% over and above the 0.4% from 18-MV 3D-CRT. Greater photoneutron dose from 18-MV IMRT may result in further maximal, absolute increased risk to peripheral tissue of approximately 1.2% over 18-MV 3D-CRT. Out-of-field photon scatter remains comparable for the same modality irrespective of beam energy. Machine scatter (C+L) from 18 versus 6 MV is 1.2 times higher for IMRT and 1.8 times for 3D-CRT. It is 4 times higher for 6-MV IMRT versus 3D-CRT. Reduction in internal scatter with 18 MV versus 6 MV is 27% for 3D-CRT and 29% for IMRT. Compared with 6-MV 3D-CRT, 18-MV IMRT increases out-of-field second cancer risk by 0.2% from photons and adds 0.28-2.2% from neutrons. Conclusions: Out-of-field photon dose seems to be independent of beam energy for both techniques. Eighteen-megavolt IMRT increases out-of-field scatter 1.7-fold over 3D-CRT because of greater collimator scatter despite reducing internal/patient scatter. Out-of-field carcinogenic risk is thus increased (but improved in-field dose conformity may offset this). Potentially increased carcinogenic risk should be weighed against any benefit 18-MV IMRT may provide.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Demonstration of a novel technique to measure two-photon exchange effects in elastic e±p scattering

The discrepancy between proton electromagnetic form factors extracted using unpolarized and polarized scattering data is believed to be a consequence of two-photon exchange (TPE) effects. However, the calculations of TPE corrections have significant model dependence, and there is limited direct experimental evidence for such corrections. The TPE contributions depend on the sign of the lepton charge in e±p scattering, but the luminosities of secondary positron beams limited past measurement at large scattering angles, where the TPE effects are believe to be most significant. We present the results of a new experimental technique for making direct e±p comparisons, which has the potential to make precise measurements over a broad range in Q 2 and scattering angles. We use the Jefferson Laboratory electron beam and the Hall B photon tagger to generate a clean but untagged photon beam. The photon beam impinges on a converter foil to generate a mixed beam of electrons, positrons, and photons. A chicane is used to separate and recombine the electron and positron beams while the photon beam is stopped by a photon blocker. This provides a combined electron and positron beam, with energies from 0.5 to 3.2 GeV, which impinges on a liquid hydrogen target. The large acceptance CLAS detector is used to identify and reconstruct elastic scattering events, determining both the initial lepton energy and the sign of the scattered lepton. The data were collected in two days with a primary electron beam energy of only 3.3 GeV, limiting the data from this run to smaller values of Q 2 and scattering angle. Nonetheless, this measurement yields a data sample for e±p with statistics comparable to those of the best previous measurements. We have shown that we can cleanly identify elastic scattering events and correct for the difference in acceptance for electron and positron scattering. Because we ran with only one polarity for the chicane, we are unable to study the difference between the incoming electron and positron beams. This systematic effect leads to the largest uncertainty in the final ratio of positron to electron scattering: R=1.027±0.005±0.05 for < Q 2 >=0.206 GeV 2 and 0.830 ≤ ε ≤ 0.943. We have demonstrated that the tertiary e ± beam generated using this technique provides the opportunity for dramatically improved comparisons of e±p scattering, covering a significant range in both Q 2 and scattering angle. Combining data with different chicane polarities will allow for detailed studies of the difference between the incoming e + and e - beams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Characterization of scatter in digital mammography from use of Monte Carlo simulations and comparison to physical measurements

Purpose: Monte Carlo simulations were performed with the goal of verifying previously published physical measurements characterizing scatter as a function of apparent thickness. A secondary goal was to provide a way of determining what effect tissue glandularity might have on the scatter characteristics of breast tissue. The overall reason for characterizing mammography scatter in this research is the application of these data to an image processing-based scatter-correction program. Methods: MCNPX was used to simulate scatter from an infinitesimal pencil beam using typical mammography geometries and techniques. The spreading of the pencil beam was characterized by two parameters: mean radial extent (MRE) and scatter fraction (SF). The SF and MRE were found as functions of target, filter, tube potential, phantom thickness, and the presence or absence of a grid. The SF was determined by separating scatter and primary by the angle of incidence on the detector, then finding the ratio of the measured scatter to the total number of detected events. The accuracy of the MRE was determined by placing ring-shaped tallies around the impulse and fitting those data to the point-spread function (PSF) equation using the value for MRE derived from the physical measurements. The goodness-of-fit was determined for each data set as a means of assessing the accuracy of the physical MRE data. The effect of breast glandularity on the SF, MRE, and apparent tissue thickness was also considered for a limited number of techniques. Results: The agreement between the physical measurements and the results of the Monte Carlo simulations was assessed. With a grid, the SFs ranged from 0.065 to 0.089, with absolute differences between the measured and simulated SFs averaging 0.02. Without a grid, the range was 0.28–0.51, with absolute differences averaging −0.01. The goodness-of-fit values comparing the Monte Carlo data to the PSF from the physical measurements ranged from 0.96 to 1.00 with a grid and 0.65 to 0.86 without a grid. Analysis of the data suggested that the nongrid data could be better described by a biexponential function than the single exponential used here. The simulations assessing the effect of breast composition on SF and MRE showed only a slight impact on these quantities. When compared to a mix of 50% glandular/50% adipose tissue, the impact of substituting adipose or glandular breast compositions on the apparent thickness of the tissue was about 5%. Conclusions: The findings show agreement between the physical measurements published previously and the Monte Carlo simulations presented here; the resulting data can therefore be used more confidently for an application such as image processing-based scatter correction. The findings also suggest that breast composition does not have a major impact on the scatter characteristics of breast tissue. Application of the scatter data to the development of a scatter-correction software program can be simplified by ignoring the variations in density among breast tissues.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗