Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distance Based Approximation”

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 325 records · Page 18

Magnetic Pumping as a Source of Particle Heating and Power-Law Distributions in the Solar Wind

Based on the rate of expansion of the solar wind, the plasma should cool rapidly as a function of distance to the Sun. Observations show this is not the case. In this work, a magnetic pumping model is developed as a possible explanation for the heating and the generation of power-law distribution functions observed in the solar wind plasma. Most previous studies in this area focus on the role that the dissipation of turbulent energy on microscopic kinetic scales plays in the overall heating of the plasma. However, with magnetic pumping, particles are energized by the largest-scale turbulent fluctuations, thus bypassing the energy cascade. In contrast to other models, we include the pressure anisotropy term, providing a channel for the large-scale fluctuations to heat the plasma directly. A complete set of coupled differential equations describing the evolution, and energization, of the distribution function are derived, as well as an approximate closed-form solution. Numerical simulations using the VPIC kinetic code are applied to verify the model's analytical predictions. The results of the model for realistic solar wind scenario are computed, where thermal streaming of particles are important for generating a phase shift between the magnetic perturbations and the pressure anisotropy. In turn, averaged over a pump cycle, the phase shift permits mechanical work to be converted directly to heat in the plasma. Here, the results of this scenario show that magnetic pumping may account for a significant portion of the solar wind energization.

79 ASTRONOMY AND ASTROPHYSICS↗

Combined Optical and Radio-Frequency Measurements of a Lightning Megaflash by the FORTE Satellite

The optical and VHF instrumentation on the Fast On-Orbit Recording of Transient Events (FORTE) satellite is used to document the combined phenomenology evolution of a lightning “megaflash” – mesoscale lightning that propagates laterally over exceptional distances. We identify a FORTE flash whose maximum extent was 82 km and inferred length over multiple distinct branches exceeded 100 km. This flash lasted 1.2 s and produced 250 optical and 591 RF events. We find that the channel development mapped by FORTE’s pixelated lightning imager (LLS) occurred at a typical speed of 2.6x10 5 m s -1 and was accompanied by sustained periods of VHF emission that could individually exceed 100 ms in duration. The impulsive IC events generated by the flash indicate that this development occurred at altitudes between 3 and 8 km. Four +CG strokes were identified in the VHF waveform data that are responsible for two of the three highly-radiant LLS groups (two of the +CGs were not as optically bright as the others). These strokes occurred at different locations throughout the flash footprint with the most distant strokes separated by approximately 50 km. These space-based observations match previous observations of megaflashes as well as ground-based measurements of negative leader development during “spider” lightning, suggesting that FORTE is sensing the same phenomena.

Peterson, Michael Jay↗

Human-in-the-loop Sensing and Control for Commercial Building Energy Efficiency and Occupant Comfort

Most of the existing heating, ventilation and air conditioning (HVAC) systems in commercial buildings operate in a conservative manner by assuming maximum occupancy in each room during pre-specified periods of the week, leading to significant energy being wasted as rooms are over-conditioned compared to the actual requirements of the occupants. Though critical, our understanding of occupancy patterns and thermal comfort needs of the occupants in commercial buildings is lacking and it is well known that both of these quantities are stochastic and time-varying, thus requiring sensing solutions to estimate them. This project had the goal of designing, implementing and evaluating a hardware and software solution to ameliorate this challenge. In particular, a depth camera (one whose pixels reveal distance from the camera as opposed to color values) placed on doorways is used to detect entrance and exit events from thermal zones in the building, and thereby estimate their occupancy levels. This information is then fed to a novel control algorithm that can, through interactions with the HVAC system, learn how to provide control inputs that maximize comfort and minimize energy waste. The resulting system represents a significant improvement over existing controllers for commercial HVAC systems and allowed us to improve our understanding of the design of future human-in-the-loop control solutions. For this solution to be feasible, the project had target metrics for its performance and cost. In particular, entrance and exit events for occupants moving about the building would need to be detected with an accuracy higher than 97%; and the resulting control inputs derived from this information would need to lead to approximately 10% energy savings compared to a schedule-based controller. Furthermore, regarding the final hardware design, the project had a target bill of materials (BOM) cost for the sensing solution of less than US$200 per unit while using less than 25W of power on average. All of these target metrics were met or exceeded by our final proposed solution. We performed evaluations by deploying the system in over 20 rooms of different types across 6 commercial buildings in Pittsburgh, PA over the course of three years, and performing targeted controlled experiments to test its performance along the different metrics. The human-in-the-loop control solutions (both hardware and software) developed through this project are expected to lead to significant improvements in the comfort and energy efficiency of HVAC systems used in commercial buildings. The insights we developed through the project pave the way to HVAC systems that can condition interior spaces according to their real-time utilization and the thermal comfort needs of the occupants, thereby reducing energy use. They also open up a new learning-based way of configuring HVAC controllers without having to manually fine-tune them for each building. These innovations can significantly increase the adoption of novel control solutions by the industry and thereby save resources and reduce costs of operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

SU-E-I-79: Effect of Number of Pinholes in Onboard Robotic Multi-Pinhole SPECT System

Purpose: To study the effect of number of pinholes for a novel Single Photon Emission Computed Tomography (SPECT) system for onboard molecular and functional imaging. Methods: Comparison studies were performed using simulation for the 49-pinhole SPECT system and a series of reductions in number of pinholes. Trajectories about the breast of a supine patient were considered. Minimum distances, radii of rotation (RORs), were determined by requirements to fully view the region of interest (ROI) and to avoid collision between the detector and the patient. Reductions in RORs translate into improvements in sensitivity. Starting from the 49-pinhole system, pinholes were removed pod by pod. The furthest two end pods in the Sup-Inf direction were removed first for their higher likelihood of alleviating the collision avoidance criteria. After iterating through different combinations of pinhole pods, and selecting three combinations, the corresponding RORs were used to analytically calculate sensitivities. Results: Based on the Methods procedure, 3 combination of pods removal were identified: 1) Superior peripheral pod 2) Inferior peripheral pod 3) both pods. RORs were reduced at only one multi-pinhole stop. Analytic calculation showed that sensitivities were reduced from 0.032 for the 49-pinhole system to 0.028 for 42-pinhole and to 0.023 for 39-pinhole system respectively. The sensitivity per pinhole detector was approximately the same for all three cases. Conclusion: For the trajectories considered, only minimal improvements in RORs were identified by removing pinhole pods. Consequently, sensitivities decreased in proportion to the number of pinholes. Studies of other anatomical sites are needed to determine if in some cases sensitivity per pinhole can be improved by removing some pinholes. PHS/NIH/NCI grant R21-CA156390-01A1.

60 APPLIED LIFE SCIENCES↗

Brillouin-Mandelstam spectroscopy of standing spin waves in a ferrite waveguide

This article reports results of experimental investigation of the spin wave interference over large distances in the Y 3 Fe 2 (FeO 4 ) 3 waveguide using Brillouin-Mandelstam spectroscopy. Two coherent spin waves are excited by the micro-antennas fabricated at the edges of the waveguide. The amplitudes of the input spin waves are adjusted to provide approximately the same intensity in the central region of the waveguide. The relative phase between the excited spin waves is controlled by the phase shifter. The change of the local intensity distribution in the standing spin wave is monitored using Brillouin-Mandelstam light scattering spectroscopy. Experimental data demonstrate the oscillation of the scattered light intensity depending on the relative phase of the interfering spin waves. The oscillations of the intensity, tunable via the relative phase shift, are observed as far as 7.5 mm away from the spin-wave generating antennas at room temperature. The obtained results are important for developing techniques for remote control of spin currents, with potential applications in spin-based memory and logic devices.

36 MATERIALS SCIENCE↗

( XMM-Newton XMM-Newton) and Chandra observations of the filled-center supernova remnant G63.7+1.1

Pulsar wind nebulae (PWNe) provide unique opportunities to study relativistic pulsar winds and their interaction with the surrounding medium. G63.7+1.1 is an intriguing supernova remnant studied so far only at radio wavelengths and classified as a PWN candidate based on its nonthermal radio emission. We present the first X-ray study of G63.7+1.1, obtained with XMM-Newton and the Chandra X-ray Observatory, with the goal of detecting and studying the candidate PWN in X-rays and searching for the neutron star candidate. The peak of the X-ray emission lies approximately 1.′2 from the peak of the 21 cm continuum emission and is offset from the center of the X-ray nebula. The diffuse X-ray emission is irregularly shaped (4.′2 × 3.′2), unlike the nearly symmetric 8′-diameter radio nebula. Its X-ray spectrum is fit with a power-law model, yielding a photon index of 1.75{sub −0.3}{sup +0.4}, a column density of (1.6 ± 0.5) × 10{sup 22} cm{sup −2}, and an X-ray luminosity of ∼2 × 10{sup 33} erg s{sup −1} (0.5–10.0 keV) at a kinematic distance of 6 kpc. Several point-like sources were detected within the limit of the X-ray nebula. We present a study of their spectra, identifying 3XMM J194753.4+274357 (CXO J194753.3+274351) as a candidate for the location of the neutron star powering the nebula. Located in a dense environment, we estimate that G63.7+1.1 is old enough (≥8 kyr) for the reverse shock to have interacted with the PWN. This, together with the large size estimated from radio studies (∼14 pc) and the offset between the radio and X-ray peaks, suggests an evolved nebula, with the radio nebula being a “relic” PWN.

79 ASTRONOMY AND ASTROPHYSICS↗

THE VELOCITY ANISOTROPY OF DISTANT MILKY WAY HALO STARS FROM HUBBLE SPACE TELESCOPE PROPER MOTIONS

Based on long baseline (5-7 years) multi-epoch HST/ACS photometry, used previously to measure the proper motion of M31, we present the proper motions (PMs) of 13 main-sequence Milky Way halo stars. The sample lies at an average distance of r {approx_equal} 24 kpc from the Galactic center, with a root-mean-square spread of 6 kpc. At this distance, the median PM accuracy is 5 km s{sup -1}. We devise a maximum likelihood routine to determine the tangential velocity ellipsoid of the stellar halo. The velocity second moments in the directions of the Galactic (l, b) system are {sup 1/2}= 123{sup +29}{sub -23} km s{sup -1}, and {sup 1/2}= 83{sup +24}{sub -16} km s{sup -1}. We combine these results with the known line-of-sight second moment, {sup 1/2}= 105{+-}5 km s{sup -1}, at this (r) to study the velocity anisotropy of the halo. We find approximate isotropy between the radial and tangential velocity distributions, with anisotropy parameter {beta}= 0.0{sup +0.2}{sub -0.4}. Our results suggest that the stellar halo velocity anisotropy out to r {approx} 30 kpc is less radially biased than solar neighborhood measurements. This is opposite to what is expected from violent relaxation, and may indicate the presence of a shell-type structure at r {approx} 24 kpc. With additional multi-epoch HST data, the method presented here has the ability to measure the transverse kinematics of the halo for more stars, and to larger distances. This can yield new improved constraints on the stellar halo formation mechanism, and the mass of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis and modeling of tungsten emission and net erosion in the DIII-D divertor using updated atomic data

Tungsten (W) is one of the leading candidate materials for plasma-facing components. However, its main drawback is its high radiative efficiency; if W penetrates the plasma, it can lead to core degradation or even collapse. Since eroded tungsten tends to ionize in the sheath and redeposit promptly, the net erosion flux that escapes prompt redeposition can differ significantly from the gross erosion. This work presents a modeling framework to estimate net erosion and photon emission from W coatings exposed to the lower divertor of DIII-D using the DiMES material exposure probe. The approach couples RustBCA for sputtering yields with a Monte Carlo transport code (LPTMC) that models redeposition and W emission. Computation is carried out with new atomic data, based on R-matrix and Mons calculations, leading to lower ionization probabilities and a twofold increase in net erosion estimates compared to calculations done with OPEN-ADAS atomic data. The model results are benchmarked against experimental measurements, showing quantitative agreement for erosion, although the trends in W emission are reproduced only qualitatively. The model is also used to assess whether W II emission can serve as a direct measurement of the net erosion of W in the lower divertor of DIII-D. Simulations show that this is not valid if the electron pressure is above ~120 Pa or if the toroidal length of the eroded material is smaller than the parallel-to-B distance traveled by impurity ions before steady-state conditions are reached. Finally, simulations suggest that when W is sputtered by carbon ions with high impact energies (≳300 eV) in DIII-D, W net erosion scales with W gross erosion and can be numerically approximated using W I flux alone as input.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrogen Refueling Reference Station Lot Size Analysis for Urban Sites

Hydrogen Fueling Infrastructure Research and Station Technology (H2FIRST) is a project initiated by the DOE in 2015 and executed by Sandia National Laboratories and the National Renewable Energy Laboratory to address R&D barriers to the deployment of hydrogen fueling infrastructure. One key barrier to the deployment of fueling stations is the land area they require (i.e. "footprint"). Space is particularly a constraint in dense urban areas where hydrogen demand is high but space for fueling stations is limited. This work presents current fire code requirements that inform station footprint, then identifies and quantifies opportunities to reduce footprint without altering the safety profile of fueling stations. Opportunities analyzed include potential new methods of hydrogen delivery, as well as alternative placements of station technologies (i.e. rooftop/underground fuel storage). As interest in heavy-duty fueling stations and other markets for hydrogen grows, this study can inform techniques to reduce the footprint of heavy-duty stations as well. This work characterizes generic designs for stations with a capacity of 600 kg/day hydrogen dispensed and 4 dispenser hoses. Three base case designs (delivered gas, delivered liquid, and on-site electrolysis production) have been modified in 5 different ways to study the impacts of recently released fire code changes, colocation with gasoline refueling, alternate delivery assumptions, underground storage of hydrogen, and rooftop storage of hydrogen, resulting in a total of 32 different station designs. The footprints of the base case stations range from 13,000 to 21,000 ft 2 . A significant focus of this study is the NFPA 2 requirements, especially the prescribed setback distances for bulk gaseous or liquid hydrogen storage. While the prescribed distances are large in some cases, these setback distances are found to have a nuanced impact on station lot size; considerations of the delivery truck path, traffic flow, parking, and convenience store location are also important. Station designs that utilize underground and rooftop storage can reduce footprint but may not be practical or economical. For example, burying hydrogen storage tanks underground can reduce footprint, but the cost savings they enable depend on the cost of burial and the cost land. Siting and economic analysis of station lot sizes illustrate the benefit of smaller station footprints in the flexibility and cost savings they can provide. This study can be used as a reference that provides examples of the key design differences that fueling stations can incorporate, the approximate sizes of generic station lots, and considerations that might be unique to particular designs.

08 HYDROGEN↗

Turbulent entrainment in finite-length wind farms

In this article, we present an entrainment-based model for predicting the flow and power output of finite-length wind farms. The model is an extension of the three-layer approach of Luzzatto-Fegiz & Caulfield ( Phys. Rev. Fluids , vol. 3, 2018, 093802) for wind farms of infinite length, and assumes dependence of key flow quantities, such as the wind farm bulk velocity, on the streamwise distance from the farm entrance. To assist our analysis and validate the proposed model, we undertake a series of large-eddy simulations with different turbine spacing arrangements and layouts. Comparisons are also made with the top-down model with entrance effects of Meneveau ( J. Turbul. , vol. 13, 2012, N7) and data from the literature. The finite-length entrainment model is shown to be capable of capturing the power drop between contiguous rows of turbines as well as describing the advection and turbulent transport of kinetic energy in both the entrance and fully developed regions. The fully developed regime is approximated only deep in the wind farm, after approximately 15 rows of turbines. Our data suggest that for the cases considered in this study, the empirical coefficients that can be used to describe turbulent entrainment and transfers above the wind farm exhibit little dependence on the farm layout and may be considered constant for modelling purposes. However, the flow field within the wind farm layer can be strongly modulated by the turbine density (spacing) as well as the array layout, and to that extent it can be argued that they are both primary factors determining the wind farm power output.

17 WIND ENERGY↗

Synchrotron Based Microstructural Characterization Method Development for Pu Oxides

Identification of new plutonium oxide signatures is a high priority focus area of the technical nuclear forensics and nonproliferation communities. Comprehensive materials science studies of plutonium oxide are needed to advance the understanding of the microcrystalline properties and their relationship to macro-scale morphologies and signatures. Plutonium oxide particle growth is depenedent on the microcrystalline strucutures, grain size, and orientation. Microcrystalline structures have been found to vary for plutonium oxide based on the oxidation state of plutonium, concentration of reactants, and the physiochemical conditions of precipitation. Empirical observation confirms the relationship between processing conditions and particle morphology. High sensitivity, non destructive analysis of individual microcrystalline grains requires the use of bright sources of X-rays produced by synchrotron radiation, such as those at the Advanced Photon Source. The light source can reveal higher resolution crystallographic mapping, more sensitive chemical analysis, and higher fidelity morphological information than any current methods at SRNL Progress of FY19 was the analysis of the synchrotron data. Areas of focus included XANES oxidation state and EXAFS to measure bond distances with modeling. Most measurements for XANES indicated that the plutonium was in the +4 oxidation state, however, some indicated that there could also be a +5 oxidation state. The EXAFS measurements showed typical Pu – O bonding distance of 2.33 and Pu – Pu bonding distances of 3.81 A; however, some samples also showed that there could be shorter Pu – O distances of approximately 1.81 A, indicative of plutonium in the +5 oxidation state.

36 MATERIALS SCIENCE↗

Phonon gas model for thermal conductivity of dense, strongly interacting liquids

Developing predictive thermal property models for liquids based on microscopic principles has been elusive. The difficulty is that liquids have gas-like and solid-like attributes that are at odds when considering the frameworks of microscopic models: Models for gases are simple due to randomness and low density, whereas models for crystalline solids rely on symmetry and long-range order for easier calculation. The short-range order in liquids does, however, provide structure to neighboring molecules similar to amorphous solids, and there have been recent advances indicating that collective vibrational modes store heat in liquids. Models combining Debye approximations from solid-state physics and Frenkel’s theory of liquids can accurately predict the heat capacity of liquids. Phonon-like dispersions in liquids have also been widely observed in neutron scattering experiments. These developments motivate us to propose a model where high-frequency vibrational modes, which travel at the speed of sound and have a mean free path on the order of the average intermolecular distance, conduct heat in liquids. We use this liquid phonon gas model to calculate the thermal conductivity of liquids with varying intermolecular interaction energies from strongest to weakest—Coulomb, hydrogen-bonding, Keesom, and London dispersion energy. Generally, the model is more accurate as the intermolecular interaction energy and density of liquids increase. The calculated thermal conductivity of Coulombic-bound molten sodium nitrate and hydrogen-bonded water is within 1.46% and 2.98% of the experimentally measured values, respectively, across their entire temperature ranges. Further modal analysis of the velocity and the mean free path of collective vibrations could establish the liquid phonon gas model as an accurate model for weakly interacting liquids as well.

36 MATERIALS SCIENCE↗

AND YET IT MOVES: THE DANGERS OF ARTIFICIALLY FIXING THE MILKY WAY CENTER OF MASS IN THE PRESENCE OF A MASSIVE LARGE MAGELLANIC CLOUD

Motivated by recent studies suggesting that the Large Magellanic Cloud (LMC) could be significantly more massive than previously thought, we explore whether the approximation of an inertial Galactocentric reference frame is still valid in the presence of such a massive LMC. We find that previous estimates of the LMC’s orbital period and apocentric distance derived assuming a fixed Milky Way (MW) are significantly shortened for models where the MW is allowed to move freely in response to the gravitational pull of the LMC. Holding other parameters fixed, the fraction of models favoring first infall is reduced. Due to this interaction, the MW center of mass within the inner 50 kpc can be significantly displaced in phase-space in a very short period of time that ranges from 0.3 to 0.5 Gyr by as much as 30 kpc and 75 km s{sup −1}. Furthermore, we show that the gravitational pull of the LMC and response of the MW are likely to significantly affect the orbit and phase space distribution of tidal debris from the Sagittarius dwarf galaxy (Sgr). Such effects are larger than previous estimates based on the torque of the LMC alone. As a result, Sgr deposits debris in regions of the sky that are not aligned with the present-day Sgr orbital plane. In addition, we find that properly accounting for the movement of the MW around its common center of mass with the LMC significantly modifies the angular distance between apocenters and tilts its orbital pole, alleviating tensions between previous models and observations. While these models are preliminary in nature, they highlight the central importance of accounting for the mutual gravitational interaction between the MW and LMC when modeling the kinematics of objects in the MW and Local Group.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-J-45: Design and Study of An In-House Respiratory Gating Phantom Platform for Gated Radiotherapy

Purpose: The main purpose of this work was to develop an in-house low cost respiratory motion phantom platform for testing the accuracy of the gated radiotherapy system and analyze the dosimetric difference during gated radiotherapy. Methods: An in-house respiratory motion platform(RMP) was designed and constructed for testing the targeting accuracy of respiratory tracking system. The RMP consist of acrylic Chest Wall Platform, 2 DC motors, 4 IR sensors, speed controller circuit, 2 LED and 2 moving rods inside the RMP. The velocity of the movement can be varied from 0 to 30 cycles per minute. The platform mounted to a base using precision linear bearings. The base and platform are made of clear, 15mm thick polycarbonate plastic and the linear ball bearings are oriented to restrict the platform to a movement of approximately 50mm up and down with very little friction. Results: The targeting accuracy of the respiratory tracking system was evaluated using phantom with and without respiratory movement with varied amplitude. We have found the 5% dose difference to the PTV during the movement in comparison with stable PTV. The RMP can perform sinusoidal motion in 1D with fixed peak to peak motion of 5 to 50mm and cycle interval from 2 to 6 seconds. The RMP was designed to be able to simulate the gross anatomical anterior posterior motion attributable to respiration-induced motion of the thoracic region. Conclusion: The unique RMP simulates breathing providing the means to create a comprehensive program for commissioning, training, quality assurance and dose verification of gated radiotherapy treatments. Create the anterior/posterior movement of a target over a 5 to 50 mm distance to replicate tumor movement. The targeting error of the respiratory tracking system is less than 1.0 mm which shows suitable for clinical treatment with highly performance.

60 APPLIED LIFE SCIENCES↗

Recombination smooths the time-signal disrupted by latency in within-host HIV phylogenies

Within-host HIV evolution involves several features that may disrupt standard phylogenetic reconstruction. One important feature is re-activation of latently integrated provirus, which has the potential to disrupt the temporal signal, leading to variation in the branch lengths and apparent evolutionary rates in a tree. Yet, real within-host HIV phylogenies tend to show clear, ladder-like trees structured by the time of sampling. Another important feature is recombination, which violates the fundamental assumption that evolutionary history can be represented by a single bifurcating tree. Thus, recombination complicates the within-host HIV dynamic by mixing genomes and creating evolutionary loop structures that cannot be represented in bifurcating trees. In this paper, we develop a coalescent-based simulator of within-host HIV evolution that includes latency, recombination, and effective population size dynamics that allows us to study the relationship between the true, complex genealogy of within-host HIV evolution, encoded as an Ancestral Recombination Graph (ARG), and the observed phylogenetic tree. To compare our ARG results to the familiar phylogeny format, we calculate the expected bifurcating tree after decomposing the ARG into all unique site trees, their combined distance matrix, and the overall corresponding bifurcating tree. While latency and recombination separately disrupt the phylogenetic signal, remarkably, we find that recombination recovers the temporal signal of within-host HIV evolution caused by latency by mixing fragments of old, latent genomes into the contemporary population. In effect, recombination averages over extant heterogeneity, whether it stems from mixed time-signals or population bottlenecks. Further, we establish that the signals of latency and recombination can be observed in phylogenetic trees despite being an incorrect representation of the true evolutionary history. Using an Approximate Bayesian Computation method, we develop a set of statistical probes to tune our simulation model to nine longitudinally-sampled within-host HIV phylogenies. Because ARGs are exceedingly difficult to infer from real HIV data, our simulation system allows investigating effects of latency, recombination, and population size bottlenecks by matching decomposed ARGs to real data as observed in standard phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗

LIFETIME AND SPECTRAL EVOLUTION OF A MAGMA OCEAN WITH A STEAM ATMOSPHERE: ITS DETECTABILITY BY FUTURE DIRECT IMAGING

We present the thermal evolution and emergent spectra of solidifying terrestrial planets along with the formation of steam atmospheres. The lifetime of a magma ocean and its spectra through a steam atmosphere depends on the orbital distance of the planet from the host star. For a Type I planet, which is formed beyond a certain critical distance from the host star, the thermal emission declines on a timescale shorter than approximately 10{sup 6} years. Therefore, young stars should be targets when searching for molten planets in this orbital region. In contrast, a Type II planet, which is formed inside the critical distance, will emit significant thermal radiation from near-infrared atmospheric windows during the entire lifetime of the magma ocean. The K{sub s} and L bands will be favorable for future direct imaging because the planet-to-star contrasts of these bands are higher than approximately 10{sup −7}–10{sup −8}. Our model predicts that, in the Type II orbital region, molten planets would be present over the main sequence of the G-type host star if the initial bulk content of water exceeds approximately 1 wt%. In visible atmospheric windows, the contrasts of the thermal emission drop below 10{sup −10} in less than 10{sup 5} years, whereas those of the reflected light remain 10{sup −10} for both types of planets. Since the contrast level is comparable to those of reflected light from Earth-sized planets in the habitable zone, the visible reflected light from molten planets also provides a promising target for direct imaging with future ground- and space-based telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗

A novel approach to assess the treatment response using Gaussian random field in PET

Purpose: The assessment of early therapeutic response to anticancer therapy is vital for treatment planning and patient management in clinic. With the development of personal treatment plan, the early treatment response, especially before any anatomically apparent changes after treatment, becomes urgent need in clinic. Positron emission tomography (PET) imaging serves an important role in clinical oncology for tumor detection, staging, and therapy response assessment. Many studies on therapy response involve interpretation of differences between two PET images, usually in terms of standardized uptake values (SUVs). However, the quantitative accuracy of this measurement is limited. This work proposes a statistically robust approach for therapy response assessment based on Gaussian random field (GRF) to provide a statistically more meaningful scale to evaluate therapy effects. Methods: The authors propose a new criterion for therapeutic assessment by incorporating image noise into traditional SUV method. An analytical method based on the approximate expressions of the Fisher information matrix was applied to model the variance of individual pixels in reconstructed images. A zero mean unit variance GRF under the null hypothesis (no response to therapy) was obtained by normalizing each pixel of the post-therapy image with the mean and standard deviation of the pretherapy image. The performance of the proposed method was evaluated by Monte Carlo simulation, where XCAT phantoms (128{sup 2} pixels) with lesions of various diameters (2–6 mm), multiple tumor-to-background contrasts (3–10), and different changes in intensity (6.25%–30%) were used. The receiver operating characteristic curves and the corresponding areas under the curve were computed for both the proposed method and the traditional methods whose figure of merit is the percentage change of SUVs. The formula for the false positive rate (FPR) estimation was developed for the proposed therapy response assessment utilizing local average method based on random field. The accuracy of the estimation was validated in terms of Euler distance and correlation coefficient. Results: It is shown that the performance of therapy response assessment is significantly improved by the introduction of variance with a higher area under the curve (97.3%) than SUVmean (91.4%) and SUVmax (82.0%). In addition, the FPR estimation serves as a good prediction for the specificity of the proposed method, consistent with simulation outcome with ∼1 correlation coefficient. Conclusions: In this work, the authors developed a method to evaluate therapy response from PET images, which were modeled as Gaussian random field. The digital phantom simulations demonstrated that the proposed method achieved a large reduction in statistical variability through incorporating knowledge of the variance of the original Gaussian random field. The proposed method has the potential to enable prediction of early treatment response and shows promise for application to clinical practice. In future work, the authors will report on the robustness of the estimation theory for application to clinical practice of therapy response evaluation, which pertains to binary discrimination tasks at a fixed location in the image such as detection of small and weak lesion.

61 RADIATION PROTECTION AND DOSIMETRY↗

Development of a novel high quantum efficiency MV x-ray detector for image-guided radiotherapy: A feasibility study

Purpose: To develop a new scintillating fiber-based electronic portal imaging device (EPID) with a high quantum efficiency (QE) while preserving an adequate spatial resolution. Methods: Two prototypes were built: one with a single pixel readout and the other with an active matrix flat-panel imager (AMFPI) for readout. The energy conversion layer of both prototypes was made of scintillating fiber layers interleaved with corrugated lead sheets to form a honeycomb pattern. The scintillating fibers have a diameter of 1 mm and the distance between the centers of neighboring fibers on the same layer is 1.35 mm. The layers have 1.22 mm spacing between them. The energy conversion layer has a thickness of 2 cm. The modulation transfer function (MTF), antiscatter properties and sensitivity of the detector with a single pixel readout were measured using a 6-MV beam on a LINAC machine. In addition, a Monte Carlo simulation was conducted to calculate the zero-frequency detective quantum efficiency (DQE(0)) of the proposed detector with an active matrix flat-panel imager for readout. Results: The DQE(0) of the proposed detector can be 11.5%, which is about an order of magnitude higher than that of current EPIDs. The frequency of 50% modulation (f 50 ) of the measured MTF is 0.2 mm –1 at 6 MV, which is comparable to that of video-based EPIDs. The scatter to primary ratio (SPR) measured with the detector at 10 cm air gap and 20 x 20 cm 2 field size is approximately 30% lower than that of ionization chamber–based detectors with a comparable QE. The detector noise which includes the x-ray quantum noise and absorption noise is much larger than the electronic noise per pixel of the flat-panel imager at a dose of less than two LINAC pulses. Thus, the proposed detector is quantum noise limited down to very low doses (~a couple of radiation pulses of the LINAC). A proof-of-concept image has been obtained using a 6-MV beam. Conclusions: This work indicates that by using scintillating fibers and lead layers it is possible to increase the thickness of the detecting materials, and therefore the QE or the DQE(0) of the detector, while maintaining an adequate spatial resolution for MV x-ray imaging. Due to the use of lead as the spacing material, the new detector also has antiscatter property, which will help improve the signal-to-noise ratio of the images. Further investigation to optimize the design of the detector and achieve a better combination of DQE and spatial resolution is warranted.

62 RADIOLOGY AND NUCLEAR MEDICINE↗