Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLANETARY BASE”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

TOI-1634 b: An Ultra-short-period Keystone Planet Sitting inside the M-dwarf Radius Valley

Studies of close-in planets orbiting M dwarfs have suggested that the M-dwarf radius valley may be well explained by distinct formation timescales between enveloped terrestrials and rocky planets that form at late times in a gas-depleted environment. This scenario is at odds with the picture that close-in rocky planets form with a primordial gaseous envelope that is subsequently stripped away by some thermally driven mass-loss process. These two physical scenarios make unique predictions of the rocky/enveloped transition’s dependence on orbital separation such that studying the compositions of planets within the M-dwarf radius valley may be able to establish the dominant physics. Here, we present the discovery of one such keystone planet: the ultra-short-period planet TOI-1634 b (P = 0.989 days, F=121F{sub ⊕}, r{sub p}=1.790{sub -0.081}{sup +0.080} R {sub ⊕}) orbiting a nearby M2 dwarf (K {sub s} = 8.7, R {sub s} = 0.450 R {sub ⊙}, M {sub s} = 0.502 M {sub ⊙}) and whose size and orbital period sit within the M-dwarf radius valley. We confirm the TESS-discovered planet candidate using extensive ground-based follow-up campaigns, including a set of 32 precise radial velocity measurements from HARPS-N. We measure a planetary mass of 4.91{sub -0.70}{sup +0.68} M {sub ⊕}, which makes TOI-1634 b inconsistent with an Earth-like composition at 5.9σ and thus requires either an extended gaseous envelope, a large volatile-rich layer, or a rocky composition that is not dominated by iron and silicates to explain its mass and radius. The discovery that the bulk composition of TOI-1634 b is inconsistent with that of Earth supports the gas-depleted formation mechanism to explain the emergence of the radius valley around M dwarfs with M{sub s}≲0.5 M {sub ⊙}.

36 MATERIALS SCIENCE↗
DOE OSTI2023

MHz free electron laser x-ray diffraction and modeling of pulsed laser heated diamond anvil cell

A new diamond anvil cell experimental approach has been implemented at the European x-ray Free Electron Laser, combining pulsed laser heating with MHz x-ray diffraction. In this report we use this setup to determine liquidus temperatures under extreme conditions, based on the determination of time-resolved crystallization. The focus is on a Fe-Si-O ternary system, relevant for planetary cores. This time-resolved diagnostic is complemented by a finite-element model, reproducing temporal temperature profiles measured experimentally using streaked optical pyrometry. This model calculates the temperature and strain fields by including (i) pressure and temperature dependencies of material properties, and (ii) the heat-induced thermal stress, including feedback effect on material parameter variations. Making our model more realistic, these improvements are critical as they give 7000 K temperature differences compared to previous models. Laser intensities are determined by seeking minimal deviation between measured and modeled temperatures. Combining models and streak optical pyrometry data extends temperature determination below detection limit. The presented approach can be used to infer the liquidus temperature by the appearance of SiO 2 diffraction spots. In addition, temperatures obtained by the model agree with crystallization temperatures reported for Fe–Si alloys. Our model reproduces the planetary relevant experimental conditions, providing temperature, pressure, and volume conditions. Those predictions are then used to determine liquidus temperatures at experimental timescales where chemical migration is limited. This synergy of novel time-resolved experiments and finite-element modeling pushes further the interpretation capabilities in diamond anvil cell experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Introducing MEGASIM: Multitudinous Earth Greek/Trojan Asteroid SIMulation*

Abstract We announce the completion of two new Earth Trojan Asteroid (ETA) simulations. ETAs co-orbit the Sun with Earth and remain loosely bound to the fourth or fifth Lagrange points. To date, two ETAs have been detected (2010 TK 7 and 2020 XL 5 ); however, based on their orbits, these are likely temporary captures rather than members of a primordial population. The motivation for our simulations is multi-faceted: (1) the lack of a detected population despite literature claims of stability over the age of the solar system, (2) ETAs are important for both solar system science and planetary defense, and (3) the proximity of ETAs makes them prime candidates for space exploration missions. These simulations are of extremely high fidelity, large in number and long in simulation time. We will be releasing a series of papers based upon these simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Multiwavelength Look at the GJ 9827 System: No Evidence of Extended Atmospheres in GJ 9827b and d from HST and CARMENES Data

GJ 9827 is a bright star hosting a planetary system with three transiting planets. As a multiplanet system with planets that sprawl within the boundaries of the radius gap between terrestrial and gaseous planets, GJ 9827 is an optimal target to study the evolution of the atmospheres of close-in planets with a common evolutionary history and their dependence from stellar irradiation. Here we report on the Hubble Space Telescope (HST) and CARMENES transit observations of GJ 9827 planets b and d. We performed a stellar and ISM characterization from the ultraviolet HST spectra, obtaining fluxes for Lyα and Mg ii of F(Lyα) = (5.42{sub −0.75}{sup +0.96}) × 10{sup −13} erg cm{sup −2} s{sup −1} and F(Mgii) = (5.64 ± 0.24) × 10{sup −14} erg cm{sup −2} s{sup −1}. We also investigated a possible absorption signature in Lyα in the atmosphere of GJ 9827b during a transit event from HST spectra, as well as Hα and He i signature for the atmosphere of GJ 9827b and d from CARMENES spectra. We found no evidence of an extended atmosphere in either of the planets. This result is also supported by our analytical estimations of mass loss based on the measured radiation fields for all three planets of this system, which led to a mass-loss rate of 0.4, 0.3, and 0.1 planetary masses per Gyr for GJ 9827b, c, and d, respectively. These values indicate that the planets could have lost their volatiles quickly in their evolution and probably do not retain an atmosphere at the current stage.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Simulating Atmospheric Processes in Earth System Models and Quantifying Uncertainties With Deep Learning Multi‐Member and Stochastic Parameterizations

Abstract Deep learning is a powerful tool to represent subgrid processes in climate models, but many application cases have so far used idealized settings and deterministic approaches. Here, we develop stochastic parameterizations with calibrated uncertainty quantification to learn subgrid convective and turbulent processes and surface radiative fluxes of a superparameterization embedded in an Earth System Model (ESM). We explore three methods to construct stochastic parameterizations: (a) a single Deep Neural Network (DNN) with Monte Carlo Dropout; (b) a multi‐member parameterization; and (c) a Variational Encoder Decoder with latent space perturbation. We show that the multi‐member parameterization improves the representation of convective processes, especially in the planetary boundary layer, compared to individual DNNs. The respective uncertainty quantification illustrates that methods (b) and (c) are advantageous compared to a dropout‐based DNN parameterization regarding the spread of convective processes. Hybrid simulations with our best‐performing multi‐member parameterizations remained challenging and crash within the first days. Therefore, we develop a pragmatic partial coupling strategy relying on the superparameterization for condensate emulation. Partial coupling reduces the computational efficiency of hybrid Earth‐like simulations but enables model stability over 5 months with our multi‐member parameterizations. However, our hybrid simulations exhibit biases in thermodynamic fields and differences in precipitation patterns. Despite this, the multi‐member parameterizations enable improvements in reproducing tropical extreme precipitation compared to a traditional convection parameterization. Despite these challenges, our results indicate the potential of a new generation of multi‐member machine learning parameterizations leveraging uncertainty quantification to improve the representation of stochasticity of subgrid effects.

Behrens, Gunnar [Deutsches Zentrum für Luft‐ und R↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Field validation of dynamic mechanical torque measurements using fiber-optic strain sensors for geared wind turbines

Abstract Accurate knowledge of the mechanical loads of wind turbine gearboxes has become essential in modern, highly loaded gearbox designs, as maintaining or even improving gearbox reliability with increasing torque density demands is proving to be challenging. Unfortunately, the traditional method of measuring dynamic mechanical torque using strain gauges placed on the outer surface of a rotating shaft and transmitting the resulting signal is unsuitable for serial deployment due to technical and economic constraints. An alternative method based on fiber-optic strain sensors placed on the stationary outer surface of the gearbox ring gear has been proposed. Like shaft torsion, the radial deformation of the ring gear is proportionate to the rotor torque. Placing the sensors on a stationary component is a cost-effective alternative for serial implementation because the need for complex and expensive data transfer via wireless transmission or a slip ring is eliminated. In this paper, we present the results of an extensive field experiment conducted to evaluate the torque measurement accuracy of this novel sensing solution installed on the gearbox of a Gamesa G97 2-MW wind turbine at the National Renewable Energy Laboratory’s Flatirons Campus. Torque measurements derived from fiber-optic strain sensors placed on the ring gear of the planetary stage are compared to conventional torque measurements from strain gauges placed on the main shaft. Two different torque estimation data processing methods were evaluated, with the method based on operational deflection shapes providing the most accurate results with an average normalized root mean square error below 0.7% for a load revolution distribution analysis. The effect of operating conditions on the torque estimate was also investigated, and the third planet-passing operational deflection shape was found to be the least sensitive to nontorque load-related effects. The fiber-optic strain sensors’ successful operation during the complete test campaign has demonstrated a robust and accurate solution for fleet-wide enhanced gearbox remaining useful life estimation.

17 WIND ENERGY↗
DOE OSTI2025

Field Validation of Dynamic Mechanical Torque Measurements for Geared Wind Turbines

Accurate knowledge of the mechanical loads of wind turbine gearboxes has become essential in modern, highly loaded gearbox designs, as maintaining or even improving gearbox reliability with increasing torque density demands is proving to be challenging. Unfortunately, the traditional method of measuring dynamic mechanical torque using strain gauges placed on the outer surface of a rotating shaft and transmitting the resulting signal is unsuitable for serial deployment due to technical and economic constraints. An alternative method based on fiber-optic strain sensors placed on the stationary outer surface of the gearbox ring gear has been proposed. Like shaft torsion, the radial deformation of the ring gear is proportionate to the rotor torque. Placing the sensors on a stationary component is a cost-effective alternative for serial implementation because the need for complex and expensive data transfer via wireless transmission or a slip ring is eliminated. In this paper, we present the results of an extensive field experiment conducted to evaluate the torque measurement accuracy of this novel sensing solution installed on the gearbox of a Gamesa G97 2-MW wind turbine at the National Renewable Energy Laboratory's Flatirons Campus. Torque measurements derived from fiber-optic strain sensors placed on the ring gear of the planetary stage are compared to conventional torque measurements from strain gauges placed on the main shaft. Two different torque estimation data processing methods were evaluated, with the method based on operational deflection shapes providing the most accurate results with an average normalized root mean square error below 0.7% for a load revolution distribution analysis. The effect of operating conditions on the torque estimate was also investigated, and the third planet-passing operational deflection shape was found to be the least sensitive to nontorque load-related effects. The fiber-optic strain sensors' successful operation during the complete test campaign has demonstrated a robust and accurate solution for fleet-wide enhanced gearbox remaining useful life estimation.

17 WIND ENERGY↗
DOE OSTI2022

The Influence of Dust-Smoke Mixtures on Boundary Layer Processes and Nocturnal Warming in the Sahel

Using a single column model with ground-based, aircraft, and satellite data sets we assess the combined role of smoke and dust aerosols, land degradation/aridization (LDA), and their impact on the planetary boundary layer (PBL) in influencing near-surface air temperature over the Sahel. Our study is unique because it assesses the combined role of smoke and dust aerosols on PBL evolution and near-surface air temperatures during both day and nighttime. More importantly, using a theoretical framework, we provide a careful explanation of the geophysical processes responsible for the changes in PBL and near-surface air temperature. Our results indicate that during northern hemisphere winter months, dust, and smoke over Sahel radiatively combine to impact the PBL. We show that aerosol mixtures dominated by dust modify PBL height in a manner that minimizes/maximizes surface layer cooling/warming at times when daytime maximum/nocturnal minimum temperatures occur. Furthermore, we find that increasing smoke contribution to total column aerosol optical extinction counteracts nighttime warming through daytime cooling. When smoke constitutes half or more of to the total column aerosol optical extinction, the ratio of longwave to shortwave radiative forcing is less than 10%, and nighttime cooling ensues. Minimum temperature is most sensitive to changes in mid-visible aerosol optical depth (AOD) values <1 and doubling of dust AOD within this range during the 1950–1980 Sahelian LDA event is estimated to have a nocturnal warming potential of 0.6°C.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Novel Spatio-Temporal Regime Tracking Method for Impact Simulations

In this proposal, we present a novel method of tracking the rheological regimes activated during impact cratering events that will allow researchers to gain new insights into cratering mechanics. Rheology describes the stress-strain response of rocks to different conditions. Planetary impact cratering events often occur on too large of a scale to be feasibly captured in controlled experiments. Instead, these dynamic events are primarily studied using multi-physics codes equipped with complex material models that enable calculations of the impact event at scale. However, determining which physical processes are required for the problem of interest is challenging. Because the dominant rheological regimes change with space and time during crater formation, it is difficult to link numerical simulations with observable features of craters at the end of the event. The basis of this work is the implementation of numerical flags that track the activation of each rheological regime throughout impact simulations. We demonstrate this with the ‘Rock Model’ implemented in the CTH shock-physics code. We will use this rheology tracking method to ’zoom in’ on a specific region within an event and track the conditions the rock experiences over time. This work will develop community benchmarks to validate and distribute our implemented rheological models. We will focus on improving the melt models used by the planetary impact modeling community by developing an EOS-aware rheological transition from solid to melt. Through analysis of cell and tracer-particle based tracking data, we will study the effects of different rheologies on modeled outcomes, particularly on the volume and distributions of impacts melts. We will also use this method to link observable features with the rheological mechanisms responsible for them. The deliverables (peer-reviewed papers) from the proposed work are (i) tests of the implemented rheologic processes and demonstrations of the tracking flags; (ii) the first calculations of the spatio-temporal evolution of the dominant rheologies during impact cratering events; and (iii) application to delivery of impactor iron during basin-scale impacts.

58 GEOSCIENCES↗
DOE OSTI2022

High-Resolution Wind Resource Data Set of the Greater Puerto Rico Region

In February 2022, the U.S. Department of Energy and six national laboratories launched the Puerto Rico Grid Resilience and Transitions to 100% Renewable Energy Study (PR100). PR100 aims to provide a comprehensive analysis of possible pathways for Puerto Rico's energy future, with a goal of 100% renewable energy by 2050. As a part of the renewable energy potential assessment in this project, we developed 20 years (2001-2020) of data using a numerical weather prediction (NWP) model for onshore and offshore wind resource assessment for the Puerto Rico. The research steps in developing the long-term wind resource data sets based on the NWP model were: 1. Model wind resource based on the Weather Research and Forecasting (WRF) model. 2. Develop WRF model configurations for Puerto Rico. 3. Test WRF with 11 different physics parameterizations for planetary boundary layer (PBL). 4. Assess WRF output from the different PBL schemes against observations. 5. Select a final model configuration which can produce the modeled wind speed with sufficient accuracy. 6. Produce 20 years wind resource data sets for Puerto Rico region. In the first stage of our framework for developing wind resource data, we developed a WRF model configuration using two nested domains (9 km and 3 km) to cover Puerto Rico and U.S. Virgin Islands and downscale the ERA5 reanalysis data (0.25 degrees x 0.25 degrees; hourly interval) to a 3-km domain. For the second stage, we implemented one-year simulations focused on using 11 different PBL physics parameterizations to find a combination of WRF physics parameterizations that could provide accurately modeled wind speed for Puerto Rico. We also analyzed the sensitivity of the modeled wind speed to PBL schemes for onshore and offshore locations. The WRF output resulting from the 11 WRF experiments using different PBL parameterizations were evaluated against observations obtained from the National Data Buoy Center (NDBC) as well as at hub height for a location for which measurements were available. A final model setup selected through the validation with observational data was used to produce 20 years of data with 3-km spatial and 5-minute temporal resolution. The WRF model output was post-processed to include wind profiles and basic atmospheric variables in a format that can be easily used for downstream modeling. The 20 years of wind resource data will be made available through NREL and support the estimation of wind energy development costs for the PR100 study.

17 WIND ENERGY↗
DOE OSTI2021

TOI-532b: The Habitable-zone Planet Finder confirms a Large Super Neptune in the Neptune Desert orbiting a metal-rich M-dwarf host

We confirm the planetary nature of TOI-532b, using a combination of precise near-infrared radial velocities with the Habitable-zone Planet Finder, Transiting Exoplanet Survey Satellite (TESS) light curves, ground-based photometric follow up, and high-contrast imaging. TOI-532 is a faint (J ~ 11.5) metal-rich M dwarf with T{sub eff} = 3957 ± 69 K and [Fe/H] = 0.38 ± 0.04; it hosts a transiting gaseous planet with a period of ~2.3 days. Joint fitting of the radial velocities with the TESS and ground-based transits reveal a planet with radius of 5.82 ± 0.19 R {sub ⊕}, and a mass of 61.5{sub -9.3}{sup +9.7} M {sub ⊕}. TOI-532b is the largest and most massive super Neptune detected around an M dwarf with both mass and radius measurements, and it bridges the gap between the Neptune-sized planets and the heavier Jovian planets known to orbit M dwarfs. It also follows the previously noted trend between gas giants and host-star metallicity for M-dwarf planets. In addition, it is situated at the edge of the Neptune desert in the Radius–Insolation plane, helping place constraints on the mechanisms responsible for sculpting this region of planetary parameter space.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring the Meteorological Impacts of Surface and Rooftop Heat Mitigation Strategies Over a Tropical City

Different heat mitigation technologies have been developed to improve the thermal environment in cities. However, the regional impacts of such technologies, especially in the context of a tropical city, remain unclear. The deployment of heat mitigation technologies at city-scale can change the radiation balance, advective flow, and energy balance between urban areas and the overlying atmosphere. We used the mesoscale Weather Research and Forecasting model coupled with a physically based single-layer urban canopy model to assess the impacts of five different heat mitigation technologies on surface energy balance, standard surface meteorological fields, and planetary boundary layer (PBL) dynamics for premonsoon typical hot summer days over a tropical coastal city in the month of April in 2018, 2019, and 2020. Results indicate that the regional impacts of cool materials (CMs), super-cool broadband radiative coolers, green roofs (GRs), vegetation fraction change, and a combination of CMs and GRs (i.e., “Cool city (CC)”) on the lower atmosphere are different at diurnal scale. Results showed that super-cool materials have the maximum potential of ambient temperature reduction of 1.6°C during peak hour (14:00 LT) compared to other technologies in the study. During the daytime hours, the PBL height was considerably lower than the reference scenario with no implementation of strategies by 700 m for super-cool materials and 500 m for both CMs and CC cases; however, the green roofing system underwent nominal changes over the urban area. During the nighttime hours, the PBL height increased by CMs and the CC strategies compared to the reference scenario, but minimal changes were evident for super-cool materials. The changes of temperature on the vertical profile of the heat mitigation implemented city reveal a stable PBL over the urban domain and a reduction of the vertical mixing associated with a pollution dome. This would lead to crossover phenomena above the PBL due to the decrease in vertical wind speed. Therefore, assessing the coupled regional impact of urban heat mitigation over the lower atmosphere at city-scale is urgent for sustainable urban planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Constraining the Aerosol Effects on Deep Convective Clouds by Considering the Coupling Between Clouds and the Planetary Boundary Layer

Several mechanisms have been proposed for the aerosol invigoration effect. Although their principles are well established, their actual magnitudes and roles in cloud development remain uncertain and debatable. This uncertainty partly stems from observational-based studies, in which it has been challenging to separate the co-variability between aerosols and meteorology. Addressing this problem requires large data samples. To this end, this study employs the Atmospheric Radiation Measurement data set expanding to 16 years (some up to 17 years, compared to 10 years in previous work) in the U.S. Southern Great Plains. It also conducts a more careful and rigorous analysis to isolate the influences of convective available potential energy (CAPE) and synoptic patterns to address a previously raised concern. We incorporated a new key process affecting aerosol-cloud interaction: cloud-surface coupling. The state/degree of the coupling relationship turns out to play an important role in the invigoration effect. Our analysis reinforces earlier findings of a robust positive relationship between cloud thickness and aerosol loading across CAPE percentiles—but only under cloud-surface coupled conditions. The increase in cloud thickness with aerosol loading is most pronounced in coupled clouds with high CAPE and bases below 1 km. Coupled clouds with bases below 1 km thicken between 1 and 4 km, depending on the CAPE percentile. Decoupled clouds show no such systematic changes. Synoptic patterns also lead to different strengths of the invigoration effect. Clean and polluted air masses are predominantly associated with northerly and southerly winds, respectively, with a stronger invigoration effect in cleaner air masses.

Geosciences↗
DOE OSTI2026

Investigating ultrafine aerosol turbulent fluxes during atmospheric new particle formation events

New particle formation (NPF) is an important source of atmospheric aerosols, yet quantitatively predicting its occurrence remains challenging, partly because ground-based measurements cannot capture NPF accurately if they occur in the upper atmosphere. While observations have confirmed the presence of new particles near the top of the planetary boundary layer (PBL), their origins and the direction of their vertical transport have remained ambiguous. Here, we propose and validate an analytical framework using airborne eddy-covariance (EC) and continuous wavelet transform (CWT) to directly quantify the vertical flux of newly formed particles and determine its dominant vertical transportation direction. Analyzing data from a dedicated airborne campaign over the Southern Great Plains, we observed a persistent and strong downward particle flux during NPF events, with a mean value of 133.8 cm?3 m s?¹ downward flux in the entrainment zone, whereas fluxes on non-NPF days were negligible. Our spectral analysis further confirms that these directional fluxes can be reliably captured using standard 1 Hz aerosol instrumentation. These findings suggest that new particle formation, driven by the entrainment of air from the overlying residual/stable layer during PBL growth, is a significant and potentially understudied source of boundary layer aerosols. The framework presented here provides a methodology to correctly attribute NPF events to specific altitudes, thereby improving the mechanistic understanding required for accurate atmospheric models.

Zhang, Ruoyu↗
DOE OSTI2020

Orbital Refinement and Stellar Properties for the HD 9446, HD 43691, and HD 179079 Planetary Systems

The Transit Ephemeris Refinement and Monitoring Survey is a project that aims to detect transits of intermediate-long period planets by refining orbital parameters of the known radial velocity planets using additional data from ground-based telescopes, calculating a revised transit ephemeris for the planet, then monitoring the planet host star during the predicted transit window. Here we present the results from three systems that had high probabilities of transiting planets: HD 9446 b and c, HD 43691 b, and HD 179079 b. We provide new radial velocity (RV) measurements that are then used to improve the orbital solution for the known planets. We search the RV data for indications of additional planets in orbit and find that HD 9446 shows a strong linear trend of 4.8σ. Using the newly refined planet orbital solutions, which include a new best-fit solution for the orbital period of HD 9446 c, and an improved transit ephemerides, we found no evidence of transiting planets in the photometry for each system. Transits of HD 9446 b can be ruled out completely and transits HD 9446 c and HD 43691 b can be ruled out for impact parameters up to b = 0.5778 and b = 0.898, respectively, due to gaps in the photometry. A transit of HD 179079 b cannot be ruled out, however, due to the relatively small size of this planet compared to the large star and thus low signal to noise. We determine properties of the three host stars through spectroscopic analysis and find through photometric analysis that HD 9446 exhibits periodic variability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Fast, Two-dimensional Gaussian Process Method Based on Celerite: Applications to Transiting Exoplanet Discovery and Characterization

Gaussian processes (GPs) are commonly used as a model of stochastic variability in astrophysical time series. In particular, GPs are frequently employed to account for correlated stellar variability in planetary transit light curves. The efficient application of GPs to light curves containing thousands to tens of thousands of data points has been made possible by recent advances in GP methods, including the celerite method. Here we present an extension of the celerite method to two input dimensions where, typically, the second dimension is small. This method scales linearly with the total number of data points when the noise in each large dimension is proportional to the same celerite kernel and only the amplitude of the correlated noise varies in the second dimension. We demonstrate the application of this method to the problem of measuring precise transit parameters from multiwavelength light curves and show that it has the potential to improve transit parameters measurements by orders of magnitude. Applications of this method include transit spectroscopy and exomoon detection, as well a broader set of astronomical problems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Noble gas incorporation into silicate glasses: implications for planetary volatile storage

Incorporation of small molecules in silicate melts may provide an important mechanism for storing noble gases in the deep Earth, yet the means by which chemically inert noble gases enter and are retained in silica-based materials is not understood. High pressure, room temperature sound velocity measurements on silica and natural basalt glasses in different pressure-transmitting media reveal that neon enters the structure of silicate glasses and enhances their elastic strengths, whereas an ethanol-methanol mixture does not. Combined with literature data, we found the incorporation of small molecules into silica and basalt glasses is controlled by the void size distribution of the glass and size of the molecules. Pressure primarily reduces the size of noble gases, thereby increasing their solubilities in silicate melts and glasses.

58 GEOSCIENCES↗