Search NASA⌕ Search

SEARCH · Search NASA

Results for “Source code instrumentation”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Verification and Validation of Open-Source Computational Fluid Dynamics Capabilities in MOOSE

The Multiphysics Object Oriented Simulation Environment (MOOSE) advances predictive modeling in nuclear research by simplifying the application design process and making modeling more efficient and accessible to researchers and the broader community. Because of MOOSE’s critical role in the scientific computation field, verification and validation are important to evaluate the reliability and credibility of MOOSE’s code, simulations, and models. To assess MOOSE’s Navier Stokes module and turbulence models, the following verification and validation University of Manchester ERCOFTAC cases were simulated in MOOSE, and the results were compared: Case 032 Turbulent Flow in a 2D Channel Case 079 Natural Convection in an Enclosed Cavity

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
NASA NTRS2018

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗
DOE OSTI2020

Experimental validation of scintillator detector response to correlated neutrons with MCNP and associated packages

Time, energy, and angle correlated neutrons are emitted in nuclear fission. New, physics-based fission event generators that capture these correlations have been incorporated into MCNP6.2 because correlations can impact nuclear non-proliferation and safeguards systems response. DRiFT is a detector response code that takes MCNP outputs and reproduces laboratory analogous detector quantities. Previous work has compared DRiFT to single detector experiment data and has discussed recent updates to DRiFT specific to correlated fission experiments. In this work, we demonstrate a validation of recent updates to DRiFT for organic scintillator response in correlated fission neutron experiments using a 252Cf source. Furthermore, outputs from MCNP with the FREYA fission event generator were post-processed with DRiFT and compared to experiment data. Simulation and experiment results are in reasonably good agreement for multiple correlated neutron quantities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

A new method for diagnosing effective radiative forcing from aerosol–cloud interactions in climate models

Aerosol–cloud interactions (ACIs) are a leading source of uncertainty in estimates of the historical effective radiative forcing (ERF). One reason for this uncertainty is the difficulty in estimating the ERF from aerosol–cloud interactions (ERFaci) in climate models, which typically requires multiple calls to the radiation code. Most commonly used methods also cannot disentangle the contributions from different processes to ERFaci. Here, we develop a new, computationally efficient method for estimating the shortwave (SW) ERFaci from liquid clouds using histograms of monthly averaged cloud fraction partitioned by cloud droplet effective radius (r e ) and liquid water path (LWP). Multiplying the histograms with SW cloud radiative kernels gives the total SW ERFaci from liquid clouds, which can be decomposed into contributions from the Twomey effect, LWP adjustments, and cloud fraction (CF) adjustments. We test the method with data from five CMIP6-era models, using the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument simulator to generate the histograms. Our method gives similar total SW ERFaci estimates to other established methods in regions of prevalent liquid cloud and indicates that the Twomey effect, LWP adjustments, and CF adjustments have contributed -0.34 ± 0.23, -0.22 ± 0.13, and -0.09 ± 0.11 W m -2 , respectively, to the effective radiative forcing of the climate since 1850 in the ensemble mean (95 % confidence). These results demonstrate that widespread adoption of a MODIS r e –LWP joint histogram diagnostic would allow the SW ERFaci and its components to be quickly and accurately diagnosed from climate model outputs, a crucial step for reducing uncertainty in the historical ERF.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

The AIRES Optical Design

AIRES (Airborne InfraRed Echelle Spectrometer) is the facility spectrometer for SOFIA (Stratospheric Observatory For Infrared Astronomy). AIRES is a long-slit (approximately 160 in) spectrometer designed to cover the 17 to 210-micron range with good sensitivity using three spectroscopic arrays. Initially, only the 30-130 micron, mid-wavelength array will be available. The instrument has a cryogenic K-mirror to perform field rotation and a slit-viewing camera (lambda < 28 microns, FOV = 160 in diameter) to image source morphology and verify telescope pointing. AIRES employs a large echelle grating to achieve a spectral resolving power (lambda/delta lambda) of approximately 1.0 x 10(exp 6)/lambda (sub mu), where lambda (sub mu) is the wavelength in microns. Hyperfine, Inc. has ruled and tested the AIRES' echelle; its wave-front error is 0.028 waves RMS (root mean square) at 10.6 microns. The instrument is housed in a liquid-helium cryostat which is constrained in diameter (approximately 1 m) and length (approximately 2 m) by the observatory. Hence, the length of the echelle (approximately 1.1 m) and the focal length of its collimator (approximately 5.2 m) severely drive the optical design and packaging. The final design uses diamond-turned aluminum optics and has up to 19 reflections inside the cryostat, depending on the optical path. This design was generated, optimized, and toleranced using Code V. The predicted performance is nearly diffraction-limited at 17 microns; the error budget is dominated by design residuals. Light loss due to slit rotation and slit curvature has been minimized. A thorough diffraction analysis with GLAD (G-Level Analysis Drawer) was used to size the mirrors and baffles; the internal light loss is shown to be a strong function of slit width.

Haas, Michael R.↗
DOE OSTI2024

The LSU-Argonne conversion electron spectrometer: A new detector for the X-Array and SATURN decay station

In this study, a new conversion electron detector has been commissioned at the ATLAS/ CARIBU facility at Argonne National Laboratory. The LSU-Argonne Conversion Electron Spectrometer (LACES) is a LN 2 -cooled Si(Li) detector system designed to be incorporated into a decay station that comprises the dedicated HPGe clover array with a box geometry (X -Array) and the Scintillator and Tape Using Radioactive Nuclei (SATURN) device. This integration enables simultaneous measurements of conversion electrons and gamma-rays in decay experiments, yielding novel information on transition multipolarities, electric monopole transitions, and isomeric states that decay mostly via conversion electrons. A measurement of the energy resolution of LACES yielded 2.3-keV FWHM at 975 keV for electrons and 1.3-keV FWHM at 75 keV for X-rays. A detailed study of the absolute detection efficiency (at 5 mm from the source) was performed, where this quantity was determined experimentally in the range of electron transition energies between 25.5 keV and 1047.8 keV and subsequently simulated using the GEANT4 code. Measurement and simulations are found to be in excellent agreement. A precise characterization, for this type of detector system, of the absolute detection efficiency for such a wide energy range is reported for the first time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Calculational Model of the Radiation Calibration Lab (RCL) and the Neutron Measurements performed at the Facility, with a 252 Cf Source

Calibration of detectors measuring gamma, neutron, alpha and beta radiation, and of personal dosimeters used at Lawrence Livermore National Laboratory (LLNL), is performed at the Radiation Calibration Lab (RCL). The calibrations done with a standard calibrated 252 Cf source are a function of the energy spectra, fluence and dose of the source. A detailed calculational model of the RCL has been carried out with the MCNP6 code. The calculations have also reproduced quite well the measurements of the fluence and dose rates for the bare and moderated 252 Cf source carried out by R. Radev using the Rospec Spectrometer and by Clark and Wilson using a Passive Neutron Spectrometer. The importance of having the calculational model of the RCL facility is to calculate the dose rates expected for a given planned measurement, optimizing the choice of detectors, design and shielding materials before the actual measurement are carried out, resulting in a more efficient and time saving measurement. Furthermore, an accurate model of the facility allows to detect possible problems with the measurement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Wayne—A Simulator for HST WFC3 IR Grism Spectroscopy

Wayne is an algorithm that simulates Hubble Space Telescope Wide Field Camera 3 (WFC3) grism spectroscopic frames, including sources of noise and systematics. It can simulate both staring and spatial scan modes, and observations such as the transit and the eclipse of an exoplanet. Unlike many other instrument simulators, the focus of Wayne is on creating frames with realistic systematics in order to test the effectiveness of different data analysis methods in a variety of different scenarios. This approach is critical for method validation and optimizing observing strategies. In this paper we describe the implementation of Wayne for WFC3 in the near-infrared channel with the G102 and G141 grisms. We compare the simulations to real data obtained for the exoplanet HD 209458b, to verify the accuracy of the simulation. The software is now available as open source at https://github.com/ucl-exoplanets/wayne.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Constraining the phase shift of relativistic species in DESI BAOs

In the early Universe, neutrinos decouple quickly from the primordial plasma and propagate without further interactions. The impact of free-streaming neutrinos is to create a temporal shift in the gravitational potential that impacts the acoustic waves known as baryon acoustic oscillations (BAOs), resulting in a non-linear spatial shift in the Fourier-space BAO signal. In this work, we make use of and extend upon an existing methodology to measure the phase shift amplitude $\beta _{\phi }$ and apply it to the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) BAOs with an anisotropic BAO fitting pipeline. We validate the fitting methodology by testing the pipeline with two publicly available fitting codes applied to highly precise cubic box simulations and realistic simulations representative of the DESI DR1 data. We find further study towards the methods used in fitting the BAO signal will be necessary to ensure accurate constraints on $\beta _{\phi }$ in future DESI data releases. Using DESI DR1, we present individual measurements of the anisotropic BAO distortion parameters and the $\beta _{\phi }$ for the different tracers, and additionally a combined fit to $\beta _{\phi }$ resulting in $\beta _{\phi } = 2.7 \pm 1.7$. After including a prior on the distortion parameters from constraints using Planck we find $\beta _{\phi } = 2.7^{+0.60}_{-0.67}$ suggesting $\beta _{\phi } > 0$ at 4.3$\sigma$ significance. This result may hint at a phase shift that is not purely sourced from the standard model expectation for $N_{\rm {eff}}$ or could be a upwards statistical fluctuation in the measured $\beta _{\phi }$; this result relaxes in models with additional freedom beyond Lambda-cold dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Upgrade of the Lyman-alpha diagnostic system on DIII-D for main chamber edge neutral studies

The LLAMA (Lyman Alpha Measurement Apparatus) pinhole camera diagnostic had previously been deployed on DIII-D to measure radial profiles of the Lyman-α (Ly-α) deuterium neutral line brightness across the plasma boundary in the lower chamber to infer neutral deuterium density and ionization rate profiles. This system has recently been upgraded with a new diagnostic head, named ALPACA, that also encloses two pinhole cameras and duplicates the LLAMA views in the upper chamber. Similar to LLAMA, ALPACA provides two times 20 lines of sight, viewing the plasma edge on the inboard and outboard sides with a radial resolution of ~2.5 cm (FWHM) and an effective time resolution of ~1 ms that allows for the investigation of inter-ELM dynamics. The extended Ly-α system provides better coverage to study neutrals in experiments with various plasma shapes utilizing both the upper and lower divertors. Furthermore, post-campaign calibration of the LLAMA diagnostic has successfully been demonstrated for the first time. This was facilitated by various upgrades to the calibration set-up and detailed measurements of the emissivity distribution of the Ly-α calibration source using a pinhole collimator. It was found that the sensitivity of the inboard LLAMA pinhole camera was reduced by a factor of 2.0 ± 0.2 over the course of six months of plasma operation in 2021. In conclusion, the upgraded Ly-α system, equipped with improved absolute calibration, will provide key input for neutral fueling and pedestal particle transport studies and for 2D edge transport code validation on the DIII-D tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

(U) Energy-Dependent Perturbations of Scattering Cross Sections Using the PERT Card in MCNP6.2

The PERT card is used for first- and second-order perturbations of material densities, isotope densities, and isotope cross sections in the MCNP6.2 Monte Carlo code. Using the PERT card for firstand second-order sensitivities is straightforward. Reference 2 only looked at density perturbations with no energy dependence. This report looks at the sensitivity of a response to scattering cross sections as a function of energy and finds that there is a bug in the PERT capability. The SENSMG multigroup neutron sensitivity code was recently given the capability of computing sensitivities to user-specified reactions whose cross sections are available in a user-supplied NJOY output file. SENSMG uses the PARTISN multigroup discrete ordinates code for the transport. The NJOY-reaction capability is similar to that of SUSD3D, which also uses PARTISN for the multigroup transport. The PERT bug was identified by comparing MCNP6.2 outputs with SENSMG outputs. The next section of this report describes the test problem. Section III presents the MCNP6.2 PERT bug and proposes a solution. Section IV is a summary and conclusions. The SENSMG and MCNP input files are listed in the appendix.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Three-phase Four-leg Inverter LabVIEW FPGA Control Code [SWR-16-26]

In the area of power electronics control, Field Programmable Gate Arrays (FPGAs) have the capability to outperform their Digital Signal Processor (DSP) counterparts due to the FPGA’s ability to implement true parallel processing and therefore facilitate higher switching frequencies, higher control bandwidth, and/or enhanced functionality. National Instruments (NI) has developed two platforms, Compact RIO (cRIO) and Single Board RIO (sbRIO), which combine a real-time processor with an FPGA. The FPGA can be programmed with a subset of the well-known LabVIEW graphical programming language. The use of cRIO and sbRIO for power electronics control has developed over the last few years to include control of three-phase inverters. Most three-phase inverter topologies include three switching legs. The addition of a fourth-leg to natively generate the neutral connection allows the inverter to serve single-phase loads in a microgrid or stand-alone power system and to balance the three-phase voltages in the presence of significant load imbalance. However, the control of a four-leg inverter is much more complex. In particular, instead of standard two-dimensional space vector modulation (SVM), the inverter requires three-dimensional space vector modulation (3D-SVM). The candidate software implements complete control algorithms in LabVIEW FPGA for a three-phase four-leg inverter. The software includes feedback control loops, three-dimensional space vector modulation gate-drive algorithms, advanced alarm handling capabilities, contactor control, power measurements, and debugging and tuning tools. The feedback control loops allow inverter operation in AC voltage control, AC current control, or DC bus voltage control modes based on external mode selection by a user or supervisory controller. The software includes the ability to synchronize its AC output to the grid or other voltage-source before connection. The software also includes provisions to allow inverter operation in parallel with other voltage regulating devices on the AC or DC buses. This flexibility allows the Inverter to operate as a stand-alone voltage source, connected to the grid, or in parallel with other controllable voltage sources as part of a microgrid or remote power system. In addition, as the inverter is expected to operate under severe unbalanced conditions, the software includes algorithms to accurately compute real and reactive power for each phase based on definitions provided in the IEEE Standard 1459: IEEE Standard Definitions for the Measurement of Electric Power Quantities Under Sinusoidal, Nonsinusoidal, Balanced, or Unbalanced Conditions. Finally, the software includes code to output analog signals for debugging and for tuning of control loops. The software fits on the Xilinx Virtex V LX110 FPGA embedded in the NI cRIO-9118 FPGA chassis, and with a 40 MHz base clock, supports a modulation update rate of 40 MHz, user-settable switching frequencies and synchronized control loop update rates of tens of kHz, and reference waveform generation, including Phase Lock Loop (PLL), update rate of 100 kHz.

Martin, Gregory↗
DOE OSTI2016

The Earth System Prediction Suite: Toward a Coordinated U.S. Modeling Capability

The Earth System Prediction Suite (ESPS) is a collection of flagship U.S. weather and climate models and model components that are being instrumented to conform to interoperability conventions, documented to follow metadata standards, and made available either under open-source terms or to credentialed users. Furthermore, the ESPS represents a culmination of efforts to create a common Earth system model architecture, and the advent of increasingly coordinated model development activities in the United States. ESPS component interfaces are based on the Earth System Modeling Framework (ESMF), community-developed software for building and coupling models, and the National Unified Operational Prediction Capability (NUOPC) Layer, a set of ESMF-based component templates and interoperability conventions. Our shared infrastructure simplifies the process of model coupling by guaranteeing that components conform to a set of technical and semantic behaviors. The ESPS encourages distributed, multiagency development of coupled modeling systems; controlled experimentation and testing; and exploration of novel model configurations, such as those motivated by research involving managed and interactive ensembles. ESPS codes include the Navy Global Environmental Model (NAVGEM), the Hybrid Coordinate Ocean Model (HYCOM), and the Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS); the NOAA Environmental Modeling System (NEMS) and the Modular Ocean Model (MOM); the Community Earth System Model (CESM); and the NASA ModelE climate model and the Goddard Earth Observing System Model, version 5 (GEOS-5), atmospheric general circulation model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

CLC: A cross-level program characterization method

Characterization of program execution plays a key role in performance improvement. There are numerous transformations applied to each step that a program takes on its lowering from source code to a compiler intermediate representation to machine language to microarchitecture-specific execution. The unpredictable benefit of each transformation step could lead a notionally superior algorithm to exhibit inferior performance once actually run, and it can be hard to discern which step in the transformation path contradicted the code developer’s assumptions. Conventional approaches to program-execution characterization consider the behavior after only a single one of those steps, which limits the information that can be provided to the user. To help address the issue of myopic views of program execution, this paper presents a novel cross-level characterization approach for understanding the behavior of program execution at different levels in the process of writing, compiling, and running a program. We show that this approach provides a richer view of the sources of performance gains and losses and helps identify program execution in a more accurate manner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Resonator Stabilization Architecture to Suppress Switching Transient Crosstalk in I-CDM

The ever-increasing sizes of transition-edge sensor (TES) microcalorimeter arrays motivates improved multiplexed readout with large multiplexing factors, low power dissipation, and low levels of crosstalk. Current-summed code division multiplexing (I-CDM) has been proposed as an alternative to flux-summed code division multiplexing (Φ-CDM) because of its lower power dissipation and greater robustness against the failure of individual readout elements. Simulating I-CDM arrays, we find that unswitched circuit components provide a mechanism for crosstalk, the magnitude of which is determined by their inductance. To mitigate this source of crosstalk, we propose a technique called resonator-stabilized I-CDM (RI-CDM), which the simulations predict will reduce crosstalk by an order of magnitude. RI-CDM reads out dc-biased TESs on an amplitude-modulated carrier wave.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Radiation shielding analysis of the small-angle X-ray scattering flight tube end station of APS upgrade project

The CSSI beamline of the APSU will have a 20.5 m long, 2.75 m diameter SS end station designed to accept pink beam. In this work, the radiation shielding analysis of this end station is analyzed using STAC8 and FLUKA codes. Effect of mirror properties such as reflectivity, coating, inclination and roughness as well as variation in the steel composition are also studied. 8 mm thick SS is found to be adequate if the source is defined by a pinhole and the white beam is reflected with two mirrors inclined at 3.0 mradian or more. The dose rates from monochromatic beams are found to be below the desired levels when the XOP calculated bandwidth (BW) are used while the use of a flat 0.1% BW requires the mirror inclinations to be above 2.5 mradian or more. STAC8 results are consistently higher than the FLUKA results by about 1.5–2.0 times.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Assessment of Model Confidence of a Laser Source Model in xRAGE Using Omega Direct-Drive Implosion Experiments

Los Alamos National Laboratory is interested in developing high-energy-density physics validation capabilities for its multiphysics code xRAGE. xRAGE was recently updated with the laser package Mazinisin to improve predictability. We assess the current implementation and coupling of the laser package via validation of laser-driven, direct-drive spherical capsule experiments from the Omega laser facility. Here, the ASME V&V 20-2009 standard is used to determine the model confidence of xRAGE, and considerations for high-energy-density physics are identified. With current modeling capabilities in xRAGE, the model confidence is overwhelmed by significant systematic errors from the experiment or model. Validation evidence suggests cross-beam energy transfer as a dominant source of the systematic error.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

MOA-2009-BLG-319Lb: A Sub-Saturn Planet inside the Predicted Mass Desert

We present an adaptive optics (AO) analysis of images from the Keck II telescope NIRC2 instrument of the planetary microlensing event MOA-2009-BLG-319. The ∼10 yr baseline between the event and the Keck observations allows the planetary host star to be detected at a separation of 66.5 ± 1.7 mas from the source star, consistent with the light-curve model prediction. The combination of the host star brightness and light-curve parameters yields host star and planet masses of M {sub host} = 0.524 ± 0.048M {sub ⊙} and m {sub p} = 67.3 ± 6.2M {sub ⊕} at a distance of D {sub L} = 7.1 ± 0.7 kpc. The star−planet projected separation is 2.03 ± 0.21 au. The planet-to-star mass ratio of this system, q = (3.857 ± 0.029) × 10{sup −4}, places it in the predicted “planet desert” at 10{sup −4} < q < 4 × 10{sup −4} according to the runaway gas accretion scenario of the core accretion theory. Seven of the 30 planets in the Suzuki et al. sample fall in this mass ratio range, and this is the third with a measured host mass. All three of these host stars have masses of 0.5 ≤ M {sub host}/M {sub ⊙} ≤ 0.7, which implies that this predicted mass ratio gap is filled with planets that have host stars within a factor of two of 1M {sub ⊙}. This suggests that runaway gas accretion does not play a major role in determining giant planet masses for stars somewhat less massive than the Sun. Our analysis has been accomplished with a modified DAOPHOT code that has been designed to measure the brightness and positions of closely blended stars. This will aid in the development of the primary method that the Nancy Grace Roman Space Telescope mission will use to determine the masses of microlens planets and their hosts.

47 OTHER INSTRUMENTATION↗