Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compatibility error”

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 217 records · Page 12

Data from: "Warming and provenance limit tree recruitment across and beyond the elevation range of subalpine forest"

This data package contains data used to support conclusions drawn in “Warming and provenance limit tree recruitment across and beyond the elevation range of subalpine forest”, by Kueppers et al. 2017. Data were collected in field sites within the Alpine Treeline Warming Experiment (ATWE), located on Niwot Ridge, on the eastern slope of the Colorado Rocky Mountains, USA. Files containing geospatial data are also included, to provide additional locational context.There are four document formats associated with this archive: three comma-separated values (.csv) files, three Microsoft Excel (.xlsx) files, one .pdf data user’s guide, four keyhole markup language (.kml) files, and a compressed folder containing seven ESRI shapefiles (.shp). The .csv files can be opened using any simple text-editor software, R, or Microsoft Excel. The .xlsx files can only be opened using Microsoft Excel. The .kml file can be opened by Google Earth and Google Maps, and the shapefiles can be opened with any GIS application compatible with the file type, such as ESRI’s ArcGIS, and QGIS.We provide two versions of the seedling data file: “PIEN_PIFLseedlings20150522_20150525rev12222020.csv/.xlsx” (hereafter PIEN_PIFLseedlings2015) and “PIEN_PIFLseedlings20160408rev12222020.csv/.xlsx” (hereafter PIEN_PIFLseedlings2016). PIEN_PIFLseedlings2015 contains the data we used in the paper. PIEN_PIFLseedlings2016 contains an updated version of these data that includes sampling from later years. The main differences between the two files lie in the columns titled “k[YEAR],” which describe the number of seedlings that were killed in a particular year. In PIEN_PIFLseedlings2016, there also is an additional year of data for k2015, and k2014 also has additional data input for the 2014 cohort. Additionally, in years 2010-2014, there are minor differences in the number of seedlings killed -- in as few as 0 plots (in 2011) to as many as 5 plots (in 2014) -- due to errors in data input that were rectified in later years.------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Upslope range shifts by subalpine tree species are a widely anticipated effect of climate change. Climate niche models predict subalpine forests to expand upslope, given more suitable growing conditions for adult trees. However, these models do not take into account climates required for successful seedling recruitment and establishment, an essential element for expansion. Further, localized upper treeline populations are hypothesized to contain favorable traits for colonizing the alpine. To test these expectations and to expand our knowledge of seedling recruitment under climate change, we designed a common garden, climate-warming experiment spread across an elevation gradient at Niwot Ridge in the Colorado Rocky Mountains. We focus on two widespread Western North American species, Engelmann spruce (Picea engelmannii Parry ex. Engelm) and limber pine (Pinus flexilis James), which occur at treeline. While the former is considered a late-seral species more tolerant of shade, limber pine is a shade-intolerant pioneer species able to establish on infertile sites.Every autumn, seeds of the two species were collected from high- (3370 m–3570 m) and low-provenance (2910–3240 m) sources close to the experimental sites and sown in our plots. A subset of plots were heated and another subset watered over the summer months to offset the effects of warming. Across five years, we found that seeds originating from low elevation recruited more strongly for both species, although this provenance difference diminished by the fourth year for Engelmann spruce, likely due to small sample sizes. Despite the recruitment of low-provenance seed, warming treatments decreased recruitment at all elevations. Combining this with the likeliness and availability of lower-quality, high provenance seed moving upslope at the treeline, tree migration into the alpine may be slowed. Overall, our findings suggest that the hardier limber pine is likely to become a more significant species in subalpine forest communities in the future, while the more sensitive Engelmann spruce may experience range contraction.

54 ENVIRONMENTAL SCIENCES↗

MORPHOLOGY OF GALAXY CLUSTERS: A COSMOLOGICAL MODEL-INDEPENDENT TEST OF THE COSMIC DISTANCE-DUALITY RELATION

Aiming at comparing different morphological models of galaxy clusters, we use two new methods to make a cosmological model-independent test of the distance-duality (DD) relation. The luminosity distances come from the Union2 compilation of Supernovae Type Ia. The angular diameter distances are given by two cluster models (De Filippis et al. and Bonamente et al.). The advantage of our methods is that they can reduce statistical errors. Concerning the morphological hypotheses for cluster models, it is mainly focused on the comparison between the elliptical {beta}-model and spherical {beta}-model. The spherical {beta}-model is divided into two groups in terms of different reduction methods of angular diameter distances, i.e., the conservative spherical {beta}-model and corrected spherical {beta}-model. Our results show that the DD relation is consistent with the elliptical {beta}-model at 1{sigma} confidence level (CL) for both methods, whereas for almost all spherical {beta}-model parameterizations, the DD relation can only be accommodated at 3{sigma} CL, particularly for the conservative spherical {beta}-model. In order to minimize systematic uncertainties, we also apply the test to the overlap sample, i.e., the same set of clusters modeled by both De Filippis et al. and Bonamente et al. It is found that the DD relation is compatible with the elliptically modeled overlap sample at 1{sigma} CL; however, for most of the parameterizations the DD relation cannot be accommodated even at 3{sigma} CL for any of the two spherical {beta}-models. Therefore, it is reasonable that the marked triaxial ellipsoidal model is a better geometrical hypothesis describing the structure of the galaxy cluster compared with the spherical {beta}-model if the DD relation is valid in cosmological observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Phase curves of small bodies from the SLOAN Moving Objects Catalog

Extensive photometric surveys continue to produce enormous stores of data on small bodies. These data are typically sparsely obtained at arbitrary (or unknown) rotational phases. Therefore, new methods for processing such data need to be developed to make the most of these vast catalogs. We aim to produce a method of recreating the phase curves of small bodies by considering the uncertainties introduced by the nominal errors in the magnitudes and the effect introduced by rotational variations. Here, we use the SLOAN Moving Objects Catalog data as a benchmark to construct phase curves of all small bodies in u', g', r', i', and z' filters. From the phase curves, we obtain the absolute magnitudes and we use them to set up the absolute colors, which are the colors of the asteroids that are not affected by changes in the phase angle. We selected objects with ≥3 observations taken in at least one filter and spanning over a minimum of 5 degrees in the phase angle. We developed a method that combines Monte Carlo simulations and Bayesian inference to estimate the absolute magnitudes using the HG 12 * photometric system. We obtained almost 15 000 phase curves, with about 12 000 of these including all five filters. The absolute magnitudes and absolute colors are compatible with previously published data that support our method. The method we developed is fully automatic and well suited for a run based on large amounts of data. Moreover, it includes the nominal uncertainties in the magnitudes and the whole distribution of possible rotational states of the objects producing what are possibly less precise values, that is, larger uncertainties, but more accurate, namely, closer to the actual value. To our knowledge, this work is the first to include the effect of rotational variations in such a manner.

79 ASTRONOMY AND ASTROPHYSICS↗

‘Universal Injector’ at LERF - Layout and Optics Architecture

Here, we present a next level design of a compact injector within the LERF vault, which would serve both the 22 GeV CEBAF and the positron program, while being compatible with the electron source required to produce positrons for Ce+BAF. The baseline design of a 3-pass recirculator features a main linac configured with three C-75 cryo-modules, five isochronous return arcs and three straight sections, facilitating electron beam acceleration up to 650 MeV. The Universal Injector offers flexibility of extraction 1-pass and 2-pass energy electrons, as needed for positron production, with a final 3-pass extraction required by the 650 MeV injector for 22 GeV CEBAF. This note provides an overview of the baseline optics design of the Universal Injector complex, including the 8 MeV injector merger and the recirculator racetrack, including individual extraction lines for all three passes. A comprehensive suite of beam dynamics studies to validate the design is under way; starting with the orbit correction scheme, followed by start-to-end tracking with lattice misalignments and magnet errors.

Bogacz, Alex [Thomas Jefferson National Accelerato↗

Active Cancellation of Servo-Induced Noise on Stabilized Lasers via Feedforward

Many precision laser applications require active frequency stabilization. However, such stabilization loops operate by pushing noise to frequencies outside their bandwidth, leading to large “servo bumps” that can have deleterious effects for certain applications. The prevailing approach to filtering this noise is to pass the laser through a high-finesse optical cavity, which places constraints on the system design. Here, we propose and demonstrate a different approach where a frequency error signal is derived from a beat note between the laser and the light that passes through the reference cavity. The phase noise derived from this beat note is fed forward to an electro-optic modulator after the laser, carefully accounting for relative delay, for real-time frequency correction. With a hertz-line-width laser, we show ≳20 dB noise suppression at the peak of the servo bump (approximately 250 kHz) and a noise-suppression bandwidth of approximately 5 MHz—well beyond the servo bump. By simulating the Rabi dynamics of a two-level atom with our measured data, we demonstrate substantial improvements to the pulse fidelity over a wide range of Rabi frequencies. Here, our approach offers a simple and versatile method for obtaining a clean spectrum of a narrow-line-width laser, as required in many emerging applications of cold atoms, and is readily compatible with commercial systems that may even include wavelength conversion.

47 OTHER INSTRUMENTATION↗

BioSig-Air-Force

1) Configured servers: In coordination with the INSIGHT team, a hardware configuration was selected. Two nodes were purchased, configured, and shipped with compatible OS and database installation. The servers have been stress tested for reliability as they use leading edge technologies. Each node has two CPUs and 12 cores per CPU with maximum onboard memory for high performance. 2) LIM and Experimental module: The original BioSig system was developed for cancer research. Accordingly, the LIM system its corresponding web pages are being modified to facilitate (i) pathogene-donor interactions, (ii) media composition, (iii) chemical and siRNA plate configurations. The LIM system has been redesigned. The revised system allows design of new media and tracking it from lot-to-lot so that variations in the phenotypic responses can be tracked to a specific media and lot number. Similar associations are also possible with other experimental factors (e.g., donor-pathoge, siRNA, and chemical). Furthermore, the design of the experimental variables has also been revised to (i) interact with the newly developed LIM system, (ii) simplify experimental specifications, and (iii) test for potential operator's error during the data entry. Part of the complication has been due to the handshake between multiple teams that provide the small molecule plates and the team that creates assay plates. Our efforts have focused to harmonize these interactions (e.g., various data formats) so that each assay plate can be mapped to its source so that a correct set of experimental variables can be associated with each image. For example, depending upon the source of the chemical plates, they may have different formats. We have developed a canonical representation that registers SMILES code, for each chemical compound, along with its physiochemical properties. The schema for LIM conjunction with customized Web pages. 3) Import of Images and computed descriptors module: In coordination with the INSIGHT team, policies were designed to route images and computed representation into BioSig. This module (i) examines for completion of image analysis, and imports images, computed masks, and descriptors into BioSig. A database API for efficient retrieval of selection of descriptors (among thousands) was designed and implemented. 4) Computed segmentation masks from external software were imported, boundaries computed, and overlaid on images for quality control.

Source record↗

ECO2M: A TOUGH2 Fluid Property Module for Mixtures of Water, NaCl, and CO2, Including Super- and Sub-Critical Conditions, and Phase Change Between Liquid and Gaseous CO2

ECO2M is a fluid property module for the TOUGH2 simulator (Version 2.0) that was designed for applications to geologic storage of CO{sub 2} in saline aquifers. It includes a comprehensive description of the thermodynamics and thermophysical properties of H{sub 2}O - NaCl - CO{sub 2} mixtures, that reproduces fluid properties largely within experimental error for temperature, pressure and salinity conditions in the range of 10 C {le} T {le} 110 C, P {le} 600 bar, and salinity from zero up to full halite saturation. The fluid property correlations used in ECO2M are identical to the earlier ECO2N fluid property package, but whereas ECO2N could represent only a single CO{sub 2}-rich phase, ECO2M can describe all possible phase conditions for brine-CO{sub 2} mixtures, including transitions between super- and sub-critical conditions, and phase change between liquid and gaseous CO{sub 2}. This allows for seamless modeling of CO{sub 2} storage and leakage. Flow processes can be modeled isothermally or non-isothermally, and phase conditions represented may include a single (aqueous or CO{sub 2}-rich) phase, as well as two-and three-phase mixtures of aqueous, liquid CO{sub 2} and gaseous CO{sub 2} phases. Fluid phases may appear or disappear in the course of a simulation, and solid salt may precipitate or dissolve. TOUGH2/ECO2M is upwardly compatible with ECO2N and accepts ECO2N-style inputs. This report gives technical specifications of ECO2M and includes instructions for preparing input data. Code applications are illustrated by means of several sample problems, including problems that had been previously solved with TOUGH2/ECO2N.

54 ENVIRONMENTAL SCIENCES↗

A spatially adaptive high-order meshless method for fluid–structure interactions

We present a scheme implementing an a posteriori refinement strategy in the context of a high-order meshless method for problems involving point singularities and fluid–solid interfaces. The generalized moving least squares (GMLS) discretization used in this work has been previously demonstrated to provide high-order compatible discretization of the Stokes and Darcy problems, offering a high-fidelity simulation tool for problems with moving boundaries. The meshless nature of the discretization is particularly attractive for adaptive h -refinement, especially when resolving the near-field aspects of variables and point singularities governing lubrication effects in fluid–structure interactions. We demonstrate that the resulting spatially adaptive GMLS method is able to achieve optimal convergence in the presence of singularities for both the div-grad and Stokes problems. Furthermore, we present a series of simulations for flows of colloid suspensions, in which the refinement strategy efficiently achieved highly accurate solutions, particularly for colloids with complex geometries.

42 ENGINEERING↗

Matlab Based LOCO

The LOCO algorithm has been used by many accelerators around the world. Although the uses for LOCO vary, the most common use has been to find calibration errors and correct the optics functions. The light source community in particular has made extensive use of the LOCO algorithms to tightly control the beta function and coupling. Maintaining high quality beam parameters requires constant attention so a relatively large effort was put into software development for the LOCO application. The LOCO code was originally written in FORTRAN. This code worked fine but it was somewhat awkward to use. For instance, the FORTRAN code itself did not calculate the model response matrix. It required a separate modeling code such as MAD to calculate the model matrix then one manually loads the data into the LOCO code. As the number of people interested in LOCO grew, it required making it easier to use. The decision to port LOCO to Matlab was relatively easy. It's best to use a matrix programming language with good graphics capability; Matlab was also being used for high level machine control; and the accelerator modeling code AT, [5], was already developed for Matlab. Since LOCO requires collecting and processing a relative large amount of data, it is very helpful to have the LOCO code compatible with the high level machine control, [3]. A number of new features were added while porting the code from FORTRAN and new methods continue to evolve, [7][9]. Although Matlab LOCO was written with AT as the underlying tracking code, a mechanism to connect to other modeling codes has been provided.

43 PARTICLE ACCELERATORS↗

SU-E-P-26: Oncospace: A Shared Radiation Oncology Database System Designed for Personalized Medicine, Decision Support, and Research

Purpose: Advancement in Radiation Oncology (RO) practice develops through evidence based medicine and clinical trial. Knowledge usable for treatment planning, decision support and research is contained in our clinical data, stored in an Oncospace database. This data store and the tools for populating and analyzing it are compatible with standard RO practice and are shared with collaborating institutions. The question is - what protocol for system development and data sharing within an Oncospace Consortium? We focus our example on the technology and data meaning necessary to share across the Consortium. Methods: Oncospace consists of a database schema, planning and outcome data import and web based analysis tools.1) Database: The Consortium implements a federated data store; each member collects and maintains its own data within an Oncospace schema. For privacy, PHI is contained within a single table, accessible to the database owner.2) Import: Spatial dose data from treatment plans (Pinnacle or DICOM) is imported via Oncolink. Treatment outcomes are imported from an OIS (MOSAIQ).3) Analysis: JHU has built a number of webpages to answer analysis questions. Oncospace data can also be analyzed via MATLAB or SAS queries.These materials are available to Consortium members, who contribute enhancements and improvements. Results: 1) The Oncospace Consortium now consists of RO centers at JHU, UVA, UW and the University of Toronto. These members have successfully installed and populated Oncospace databases with over 1000 patients collectively.2) Members contributing code and getting updates via SVN repository. Errors are reported and tracked via Redmine. Teleconferences include strategizing design and code reviews.3) Successfully remotely queried federated databases to combine multiple institutions’ DVH data for dose-toxicity analysis (see below – data combined from JHU and UW Oncospace). Conclusion: RO data sharing can and has been effected according to the Oncospace Consortium model: http://oncospace.radonc.jhmi.edu/ . John Wong - SRA from Elekta; Todd McNutt - SRA from Elekta; Michael Bowers - funded by Elekta.

60 APPLIED LIFE SCIENCES↗

Development and implementation of high-throughput proteomic and metabolomics assays by using advanced chromatographic and mass spectrometric systems (CRADA Final Report)

The mission of this CRADA with Agilent was to couple powerful MS platforms (QQQ, IM-QTOFMS) with Agilent’s novel Ultra-High-Performance Liquid Chromatography (UHPLC) fast metabolomic workflows and perform ABF Machine Learning (ML) to generated datasets. Agilent transferred UHPLC methods to PNNL and LBNL and methods were implemented and demonstrated in both labs, achieving total acquisition times of < 10 min. Metabolites analyzed using Agilent’s shared methods included metabolites from central carbon metabolism, common across hosts, and metabolites unique to engineered strains. Standards were acquired in an UHPLC-Drift Tube Ion Mobility Mass Spectrometer (DTIMS) system for the first time within the context of ABF and methods were optimized based on Agilent’s protocols. Samples from ABF hosts Pseudomonas putida, Aspergillus pseudoterreus, Aspergillus niger and Rhodosporidium toruloides were analyzed using the UHPLC-DTIMS platform for a total of 276 runs. A data analysis workflow compatible with the Experimental Data Depot (EDD) and completely shareable was developed for the acquired UHPLC-DTIMS data. Samples were analyzed using a Data Independent Acquisition Approach (DIA), which for most of the standards provided more transitions therefore increasing detection confidence. Using the data acquired by PNNL, LBNL, and Agilent’s specifications from previous ML projects, SNL applied an ensemble ML strategy to pick the best performing model for automated LC-method selection. Finally, with the contribution of the participant labs and Agilent, SNL developed an Automated Method Selection (AMS) software tool to predict the best liquid chromatography method for analysis of any new molecules of interest. Samples with novel pathways and new metabolite targets of interest are generated at a high pace in the ABF. Overall, the project advanced rapid metabolomics by combining liquid chromatography, ion mobility spectrometry, and data-independent mass spectrometry with machine learning. This multidimensional approach uses retention time, collision cross-section, precursor mass, and fragment-ion information to distinguish chemically similar metabolites that can be difficult to resolve using conventional liquid- or gas-chromatography methods. The resulting workflow also provided automated metabolite-identification error estimates, addressing a recognized need for statistical confidence measures in metabolomics.

Petzold, Christopher [Lawrence Berkeley National L↗

A Closed-Loop Proportional-Integral (PI) Control Software for Fully Mechanically Controlled Automated Electron Microscopic Tomography

A closed-loop proportional-integral (PI) control software is provided for fully mechanically controlled automated electron microscopic tomography. The software is developed based on Gatan DigitalMicrograph, and is compatible with Zeiss LIBRA 120 transmission electron microscope. However, it can be expanded to other TEM instrument with modification. The software consists of a graphical user interface, a digital PI controller, an image analyzing unit, and other drive units (i.e.: image acquire unit and goniometer drive unit). During a tomography data collection process, the image analyzing unit analyzes both the accumulated shift and defocus value of the latest acquired image, and provides the results to the digital PI controller. The digital PI control compares the results with the preset values and determines the optimum adjustments of the goniometer. The goniometer drive unit adjusts the spatial position of the specimen according to the instructions given by the digital PI controller for the next tilt angle and image acquisition. The goniometer drive unit achieves high precision positioning by using a backlash elimination method. The major benefits of the software are: 1) the goniometer drive unit keeps pre-aligned/optimized beam conditions unchanged and achieves position tracking solely through mechanical control; 2) the image analyzing unit relies on only historical data and therefore does not require additional images/exposures; 3) the PI controller enables the system to dynamically track the imaging target with extremely low system error.

Ren, Gang↗

Quantifying the Diversification of Hepatitis C Virus (HCV) during Primary Infection: Estimates of the In Vivo Mutation Rate

Hepatitis C virus (HCV) is present in the host with multiple variants generated by its error prone RNA-dependent RNA polymerase. Little is known about the initial viral diversification and the viral life cycle processes that influence diversity. We studied the diversification of HCV during acute infection in 17 plasma donors, with frequent sampling early in infection. To analyze these data, we developed a new stochastic model of the HCV life cycle. We found that the accumulation of mutations is surprisingly slow: at 30 days, the viral population on average is still 46% identical to its transmitted viral genome. Fitting the model to the sequence data, we estimate the median in vivo viral mutation rate is 2.561025 mutations per nucleotide per genome replication (range 1.6–6.261025 ), about 5-fold lower than previous estimates. To confirm these results we analyzed the frequency of stop codons (N = 10) among all possible non-sense mutation targets (M = 898,335), and found a mutation rate of 2.8–3.261025 , consistent with the estimate from the dynamical model. The slow accumulation of mutations is consistent with slow turnover of infected cells and replication complexes within infected cells. This slow turnover is also inferred from the viral load kinetics. Our estimated mutation rate, which is similar to that of other RNA viruses (e.g., HIV and influenza), is also compatible with the accumulation of substitutions seen in HCV at the population level. Our model identifies the relevant processes (long-lived cells and slow turnover of replication complexes) and parameters involved in determining the rate of HCV diversification.

60 APPLIED LIFE SCIENCES↗

TH-E-17A-04: Geometric Validation of K-Space Self-Gated 4D-MRI Vs. 4D-CT Using A Respiratory Motion Phantom

Purpose: 4D-CT is often limited by motion artifacts, low temporal resolution, and poor phase-based target definition. We recently developed a novel k-space self-gated 4D-MRI technique with high spatial and temporal resolution. The goal here is to geometrically validate 4D-MRI using a MRI-CT compatible respiratory motion phantom and comparison to 4D-CT. Methods: 4D-MRI was acquired using 3T spoiled gradient echo-based 3D projection sequences. Respiratory phases were resolved using self-gated k-space lines as the motion surrogate. Images were reconstructed into 10 temporal bins with 1.56×1.56×1.56mm3. A MRI-CT compatible phantom was designed with a 23mm diameter ball target filled with highconcentration gadolinium(Gd) gel embedded in a 35×40×63mm3 plastic box stabilized with low-concentration Gd gel. The whole phantom was driven by an air pump. Human respiratory motion was mimicked using the controller from a commercial dynamic phantom (RSD). Four breathing settings (rates/depths: 10s/20mm, 6s/15mm, 4s/10mm, 3s/7mm) were scanned with 4D-MRI and 4D-CT (slice thickness 1.25mm). Motion ground-truth was obtained from input signals and real-time video recordings. Reconstructed images were imported into Eclipse(Varian) for target contouring. Volumes and target positions were compared with ground-truth. Initial human study was investigated on a liver patient. Results: 4D-MRI and 4D-CT scans for the different breathing cycles were reconstructed with 10 phases. Target volume in each phase was measured for both 4D-CT and 4D-MRI. Volume percentage difference for the 6.37ml target ranged from 6.67±5.33 to 11.63±5.57 for 4D-CT and from 1.47±0.52 to 2.12±1.60 for 4D-MRI. The Mann-Whitney U-test shows the 4D-MRI is significantly superior to 4D-CT (p=0.021) for phase-based target definition. Centroid motion error ranges were 1.35–1.25mm (4D-CT), and 0.31–0.12mm (4D-MRI). Conclusion: The k-space self-gated 4D-MRI we recently developed can accurately determine phase-based target volume while avoiding typical motion artifacts found in 4D-CT, and is being further studied for use in GI targeting and motion management. This work supported in part by grant 1R03CA173273-01.

60 APPLIED LIFE SCIENCES↗

SPACE VELOCITIES OF SOUTHERN GLOBULAR CLUSTERS. VI. NINE CLUSTERS IN THE INNER MILKY WAY

We have measured the absolute proper motions of nine low-latitude, inner Galaxy globular clusters, namely, NGC 6273 (M 19), NGC 6284, NGC 6287, NGC 6293, NGC 6333 (M 9), NGC 6342, NGC 6356, NGC 6388, and NGC 6441. These are the first determinations ever made for these clusters. The proper motions are on the International Celestial Reference System via Hipparcos. The proper-motion errors range between 0.4 and 0.9 mas yr{sup -1} and are dominated by the number of measurable cluster members in these regions which are very crowded by the bulge/bar and the thick disk. This sample contains five metal-poor ([Fe/H] < -1.0) and four metal-rich ([Fe/H] {>=} -1.0) clusters; seven clusters are located within {approx}4 kpc from the Galactic center, while the remaining two, namely NGC 6356 and NGC 6284, are in the background of the bulge at {approx}7.5 kpc from the Galactic center. By combining proper motions with radial velocities and distances from the literature we derive three-dimensional velocities. In a number of cases, distance uncertainties make the kinematical classification ambiguous. For the metal-poor group of clusters, we find that three clusters, namely NGC 6273, NGC 6287, and NGC 6293 are members of a kinematically hot system, the inner halo. As for the remaining two metal-poor clusters, NGC 6284 is located at {approx}7.5 kpc from the Galactic center and kinematically belongs to the thick disk, while NGC 6333, located in the inner {approx}2 kpc, has an uncertain membership (between halo and thick disk) due to the distance uncertainty. Within the metal-rich group of clusters, NGC 6356 and NGC 6342 have velocities compatible with membership in the thick disk; however, velocity uncertainties do not allow us to rule out their belonging to a hotter system. NGC 6342 is within {approx}2 kpc from the Galactic center, and thus it may belong to the bulge. NGC 6356 is at {approx}7.5 kpc from the Galactic center and thus its metallicity, kinematics, and location argue together in favor of thick-disk membership. The remaining two metal-rich clusters, NGC 6388 and NGC 6441, have velocities incompatible with membership in the thick disk or the bar of the Milky Way. They can be thought of as members of a kinematically hot system in the inner Galaxy. Curiously, both clusters have similar velocity components. Together with their similar Galactic location and peculiar but similar stellar-population characteristics, these two clusters may share a common origin. Their velocities are also very low indicating that the two clusters are now at apocenter, i.e., they will not leave the inner {approx}4 kpc of the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

THE HELIUM CONTENT OF GLOBULAR CLUSTERS: NGC 6121 (M4)

In the context of the multiple stellar population scenario in globular clusters, helium (He) has been proposed as a key element to interpret the observed multiple main sequences, subgiant branches, and red giant branches, as well as the complex horizontal branch (HB) morphology. In particular, second-generation stars belonging to the bluer part of the HB are thought to be more He-rich ({Delta}Y = 0.03 or more) but also more Na-rich/O-poor than those located in the redder part that should have Y equal to the cosmological value. Up to now this hypothesis was only partially confirmed in NGC 6752, where stars of the redder zero-age HB showed an He content of Y = 0.25 {+-} 0.01, fully compatible with the primordial He content of the universe, and were all Na-poor/O-rich. Here we study hot blue horizontal branch (BHB) stars in the GC NGC 6121 (M4) to measure their He plus O/Na content. Our goal is to complete the partial results obtained for NGC 6752, focusing our attention on targets located on the bluer part of the HB of M4. We observed six BHB stars using the VLT2/UVES spectroscopic facility. Spectra of signal-to-noise ratio {approx} 150 were obtained and the very weak He line at 5875 A measured for all our targets. We compared this feature with synthetic spectra to obtain He abundances. In addition O, Na, and Fe abundances were estimated. Stars turned out to be all Na-rich and O-poor and to have a homogeneous He content with a mean value of Y = 0.29 {+-} 0.01(random) {+-} 0.01(systematic), which is enhanced by {Delta}Y {approx} 0.04 with respect to the most recent measurements of the primordial He content of the universe (Y {approx} 0.24/0.25). The high He content of blue HB stars in M4 is also confirmed by the fact that they are brighter than red HB stars (RHB). Theoretical models suggest the BHB stars are He-enhanced by {Delta}(Y) = 0.02/0.03 with respect to the RHB stars. The whole sample of stars has a metallicity of [Fe/H] = -1.06 {+-} 0.02 (internal error), in agreement with other studies available in the literature. This is a rare direct measurement of the (primordial) He abundance for stars belonging to the Na-rich/O-poor population of GC stars in a temperature regime where the He content is not altered by sedimentation or extreme mixing as suggested for the hottest, late helium flash HB stars. Our results support theoretical predictions that the Na-rich/O-poor population is also more He-rich than the Na-poor/O-rich generation and that a leading contender for the second parameter is the He abundance.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-180: Fano Cavity Test of Proton Transport in Monte Carlo Codes Running On GPU and Xeon Phi

Purpose: In proton dose calculation, clinically compatible speeds are now achieved with Monte Carlo codes (MC) that combine 1) adequate simplifications in the physics of transport and 2) the use of hardware architectures enabling massive parallel computing (like GPUs). However, the uncertainties related to the transport algorithms used in these codes must be kept minimal. Such algorithms can be checked with the so-called “Fano cavity test”. We implemented the test in two codes that run on specific hardware: gPMC on an nVidia GPU and MCsquare on an Intel Xeon Phi (60 cores). Methods: gPMC and MCsquare are designed for transporting protons in CT geometries. Both codes use the method of fictitious interaction to sample the step-length for each transport step. The considered geometry is a water cavity (2×2×0.2 cm{sup 3}, 0.001 g/cm{sup 3}) in a 10×10×50 cm{sup 3} water phantom (1 g/cm{sup 3}). CPE in the cavity is established by generating protons over the phantom volume with a uniform momentum (energy E) and a uniform intensity per unit mass I. Assuming no nuclear reactions and no generation of other secondaries, the computed cavity dose should equal IE, according to Fano's theorem. Both codes were tested for initial proton energies of 50, 100, and 200 MeV. Results: For all energies, gPMC and MCsquare are within 0.3 and 0.2 % of the theoretical value IE, respectively (0.1% standard deviation). Single-precision computations (instead of double) increased the error by about 0.1% in MCsquare. Conclusion: Despite the simplifications in the physics of transport, both gPMC and MCsquare successfully pass the Fano test. This ensures optimal accuracy of the codes for clinical applications within the uncertainties on the underlying physical models. It also opens the path to other applications of these codes, like the simulation of ion chamber response.

60 APPLIED LIFE SCIENCES↗