Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Inference”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

VERITAS : A density-functional theory-based multiband kinetic model for understanding x-ray spectroscopy of dense plasmas

X-ray spectroscopy has long been a powerful diagnostic tool for hot, dilute plasmas, providing insights into plasma conditions by measuring line shifts and broadenings of atomic transitions. The technique critically depends on the accuracy of atomic physics models used to interpret spectroscopic measurements for inferring plasma properties such as free-electron density and temperature. Over the past decades, the atomic and plasma physics communities have developed robust atomic physics models to account for various processes in hot, dilute classical plasmas. While these models have been successful in that regime, their applicability becomes uncertain when interpreting x-ray spectroscopy experiments of above-solid-density plasmas. Given that finite-temperature density-functional theory (DFT) offers a more accurate description of dense plasma environments, we present the development of a DFT-based multi-band kinetic model, VERITAS, designed to improve the interpretation of x-ray spectroscopic measurements in high-density plasmas produced by laser-driven spherical implosions. This work details the VERITAS model and its application to both time-integrated and time-resolved x-ray spectra from implosion experiments on OMEGA. The advantages and limitations of the VERITAS model will also be discussed, along with potential directions for advancing x-ray spectroscopy of dense and superdense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Differentially Private Generation of Social Networks via Exponential Random Graph Models

Many social networks contain sensitive relational information. One approach to protect the sensitive relational information while offering flexibility for social network research and analysis is to release synthetic social networks at a pre-specified privacy risk level, given the original observed network. In this work, we propose the DP-ERGM procedure that synthesizes networks that satisfy the differential privacy (DP) via the exponential random graph model (EGRM). We apply DP-ERGM to a college student friendship network and compare its original network information preservation in the generated private networks with two other approaches: differentially private DyadWise Randomized Response (DWRR) and Sanitization of the Conditional probability of Edge given Attribute classes (SCEA). The results suggest that DP-EGRM preserves the original information significantly better than DWRR and SCEA in both network statistics and inferences from ERGMs and latent space models. In addition, DP-ERGM satisfies the node DP, a stronger notion of privacy than the edge DP that DWRR and SCEA satisfy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Development of an Austenitization Kinetics Model for 22MnB5 Steel

This paper presents a first-order austenitization kinetics model for 22MnB5 steel, commonly used in hot forming die quenching. Model parameters are derived from constant heating rate dilatometry measurements. Vickers hardness measurements made on coupons that were quenched at intermediate stages of the process were used to verify the model, and the Ac{sub 1} and Ac{sub 3} temperatures inferred from dilatometry are consistent with correlations found in the literature. The austenitization model was extended to consider non-constant heating rates typical of industrial furnaces and again showed reasonable agreement between predictions and measurements. Finally, the model is used to predict latent heat evolution during industrial heating and is shown to be consistent with values inferred from thermocouple measurements of furnace-heated 22MnB5 coupons reported in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hamiltonian parameter inference from resonant inelastic x-ray scattering with active learning

Identifying model Hamiltonians is a vital step toward creating predictive models of materials. Here, in this study, we combine Bayesian optimization with the EDRIXS numerical package to infer Hamiltonian parameters from resonant inelastic x-ray scattering (RIXS) spectra within the single atom approximation. To evaluate the efficacy of our method, we test it on experimental RIXS spectra of NiPS 3 , NiCl 2 , Ca 3 ⁢LiOsO 6 , and Fe 2⁢ O 3 , and demonstrate that it can reproduce results obtained from hand-fitted parameters to a precision similar to expert human analysis while providing a more systematic mapping of parameter space. Our work provides a key first step toward solving the inverse scattering problem to extract effective multi-orbital models from information-dense RIXS measurements, which can be applied to a host of quantum materials. We also propose atomic model parameter sets for two materials, Ca 3⁢ LiOsO 6 and Fe 2⁢ O 3 , that were previously missing from the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Inferring demographic and selective histories from population genomic data using a 2-step approach in species with coding-sparse genomes: an application to human data

Abstract The demographic history of a population, and the distribution of fitness effects (DFE) of newly arising mutations in functional genomic regions, are fundamental factors dictating both genetic variation and evolutionary trajectories. Although both demographic and DFE inference has been performed extensively in humans, these approaches have generally either been limited to simple demographic models involving a single population, or, where a complex population history has been inferred, without accounting for the potentially confounding effects of selection at linked sites. Taking advantage of the coding-sparse nature of the genome, we propose a 2-step approach in which coalescent simulations are first used to infer a complex multi-population demographic model, utilizing large non-functional regions that are likely free from the effects of background selection. We then use forward-in-time simulations to perform DFE inference in functional regions, conditional on the complex demography inferred and utilizing expected background selection effects in the estimation procedure. Throughout, recombination and mutation rate maps were used to account for the underlying empirical rate heterogeneity across the human genome. Importantly, within this framework it is possible to utilize and fit multiple aspects of the data, and this inference scheme represents a generalized approach for such large-scale inference in species with coding-sparse genomes.

Soni, Vivak (ORCID:0000000294969562)↗
DOE OSTI2016

Nitrate vulnerability projections from Bayesian inference of multiple groundwater age tracers

Nitrate is a major source of contamination of groundwater in the United States and around the world. We tested the applicability of multiple groundwater age tracers ( 3 H, 3 He, 4 He, 14 C, 13 C, and 85 Kr) in projecting future trends of nitrate concentration in 9 long-screened, public drinking water wells in Turlock, California, where nitrate concentrations are increasing toward the regulatory limit. Very low 85 Kr concentrations and apparent 3 H/ 3 He ages point to a relatively old modern fraction (40–50 years), diluted with pre-modern groundwater, corroborated by the onset and slope of increasing nitrate concentrations. An inverse Gaussian–Dirac model was chosen to represent the age distribution of the sampled groundwater at each well. Model parameters were estimated using a Bayesian inference, resulting in the posterior probability distribution – including the associated uncertainty – of the parameters and projected nitrate concentrations. Three scenarios were considered, including combined historic nitrate and age tracer data, the sole use of nitrate and the sole use of age tracer data. Each scenario was evaluated based on the ability of the model to reproduce the data and the level of reliability of the nitrate projections. The tracer-only scenario closely reproduced tracer concentrations, but not observed trends in the nitrate concentration. Both cases that included nitrate data resulted in good agreement with historical nitrate trends. Use of combined tracers and nitrate data resulted in a narrower range of projections of future nitrate levels. However, use of combined tracer and nitrate resulted in a larger discrepancy between modeled and measured tracers for some of the tracers. In conclusion, despite nitrate trend slopes between 0.56 and 1.73 mg/L/year in 7 of the 9 wells, the probability that concentrations will increase to levels above the MCL by 2040 are over 95% for only two of the wells, and below 15% in the other wells, due to a leveling off of reconstructed historical nitrate loadings to groundwater since about 1990.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Mapping Soil Cation‐Exchange Capacity using Bayesian Modeling and Proximal Sensors at the Field Scale

Core Ideas A Bayesian inference approach (INLA‐SPDE) was used to map topsoil and subsoil CEC. DEM, gamma‐ray spectrometer and EM induction data were combined to map CEC. Posterior marginal distributions of the model parameters and responses were estimated. Gamma‐ray data performed best in the topsoil, followed by DUALEM421S and elevation. Elevation data performed best in the subsoil, followed by gamma‐ray and DUALEM421S. Cation exchange capacity (CEC) affects soil fertility, acidity, and structural resilience. This is particularly the case in sugarcane growing areas of Australia because the soil there is sandy (>60%), strongly acidic (pH < 5.5), and strongly sodic (exchangeable sodium percentage [ESP] > 15%). Unfortunately, obtaining information on CEC at the field extent is time‐consuming and expensive. Here, we used a digital soil mapping approach to add value to limited (40) topsoil (0–0.3 m) and subsoil (0.6–0.9 m) CEC information. We first collected proximally sensed ancillary data from three sources, including a digital elevation model (DEM), γ‐ray (γ‐ray) spectrometer (RS700) and electromagnetic (EM) induction instruments. We then use a Bayesian inference approach (Integrated Nested Laplace Approximation with Stochastic Partial Differential Equation, INLA‐SPDE) implemented in R software to model the CEC and ancillary data. Accuracy (RMSE), bias (ME), and concordance (Lin's) of models were also generated from the different sources of ancillary data, either in combination or alone. We concluded, overall, that the INLA–SPDE approach could provide estimations of the posterior marginal distributions of the model parameters as well as the model responses as reported by other researchers. We also concluded that using the ancillary data sources in combination was most accurate (e.g., RMSE = 0.72) to predict CEC, least biased (e.g., ME = 0.07) and had the highest concordance (e.g., Lin's = 0.69) in both the topsoil and subsoil than using the ancillary data alone. The best ancillary data, when used alone for mapping CEC in the topsoil, was γ‐ray spectrometry, followed by EM data and elevation. For subsoil CEC, it was elevation, followed by γ‐ray spectrometry and then soil electrical conductivity (EC a ) data. The maps of the credibility interval (CI) indicated that better predictions were achieved in the topsoil and indicated where improvements in prediction could be achieved in the subsoil.

Li, Nan↗
DOE OSTI2025

Demonstrate new plasticity models for doped UO 2 that capture dislocation mechanisms

In light water reactors, fuel vendors are investigating the use of dopants to modify the properties of UO 2 pellets, with the goal of improving pellet-cladding mechanical interactions during operation. Dopants are expected to ‘soften’ the pellets; that is, the doped pellets have higher plastic deformation than conventional UO 2 . This leads to a reduction in the severity of mechanical pellet-cladding interactions, helping to reduce the hoop strain on the cladding. By minimizing the strain exerted by the pellet on the cladding, it is anticipated that cladding performance under accident conditions can be enhanced (i.e., lowering the risk of burst during a LOCA). Dopants such as chromium (Cr) promote grain growth during pellet fabrication, leading to larger grains; therefore, understanding the link between chemistry, microstructure and mechanical deformation (enhanced creep rates) behavior of UO 2 is critical to helping operators further substantiate the benefits of doping UO 2 . Historically, the nuclear energy industry has relied on empirical models to make assessments of performance. Compared to empirical models, mechanistic physics-based models provide benefits, such as, fewer data points for validation and better extrapolation where experimental data is scarce or non-existent. In this report, Bayesian inference techniques have been applied to a previously developed lower length-scale-informed diffusional creep model. The objective is to i) infer lower-length-scale parameter distributions from available experiment and then ii) determine the uncertainties in the measurable quantity (in this case creep rates) after propagating the inferred lower length scale parameter uncertainties. The approach requires many evaluations of the model, which becomes computationally insurmountable; therefore, a neural-network model is trained to data obtained by sampling the full model over the most important parameters. This neural-network is then used in the Bayesian inference approach to determine probability distributions in the parameter values that represent the uncertainty in the model given what is known from the experiments (posterior). A significant reduction compared to conservative initial (prior) uncertainties is achieved through inference against the experimental data, demonstrating the efficacy of this approach. Furthermore, by accounting for uncertainties in the experimental conditions and sample non-stoichiometry, it is possible to resolve apparent discrepancies in experimental measurements within a self-consistent grain boundary (Coble) creep model that is sensitive to chemistry. This work has been written up and submitted to Nuclear Technology for a special issue on accelerated fuel qualification (AFQ). This uncertainty quantification (UQ) work not only improves the diffusional model, while accounting for uncertainty, but also establishes a framework which can readily be applied to the mechanistic models of dislocation deformation developed in this study. The most likely values from the Bayesian analysis are incorporated into our UO 2 diffusional creep model and a lower length scale-informed irradiation UO 2 creep mechanistic model to generate a dataset. This dataset has been provided to our INL collaborators for training an artificial neural network surrogate model, which will be implemented in the BISON fuel performance code to assess how the results differ from those currently obtained using a fully empirical model and that of using the nominal (uncalibrated) atomic scale parameters in our mechanistic model. Plastic deformation (creep and glide) in UO 2 is a complex phenomenon, governed by multiple underlying processes such as local defect concentrations, applied stresses, and microstructural characteristics. Consequently, there is a need for a meso-scale model with polycrystalline resolution capable of extrapolating to large grain sizes applicable to doped UO 2 , where data is limited and the model can help bridge the knowledge gap. By integrating atomistic data into the polycrystal LApx code, it becomes possible to predict dislocation climb and glide plasticity that simple analytical models cannot accurately represent. The application of atomic-scale data within LApx demonstrated the importance of climb and glide mechanisms in reproducing high-stress UO 2 behavior. Behaviors such as this are crucial to capture and implement in BISON, as parts of the fuel pellet can reach temperatures where glide can occur before pellet cracking. This model which captures dislocation based mechanisms for UO 2 is then used to stand up the doped model accounting for larger grain sizes. It was found that larger grain sizes can lead to enhanced deformation rates in the glide regime, and therefore can help with the pellet cladding mechanical interaction. Therefore if the fuel pellet reaches conditions (stress/temperature) where glide is active, the enhanced creep rates for larger grains in the glide regime (doped UO 2 ) can help with pellet cladding mechanical interactions. Plastic deformation in UO 2 involves multiple mechanisms, including diffusional creep, dislocation climb, and glide. This milestone contains two parts: (1) UQ of a pre-existing lower length scale informed mechanistic diffusional creep model, and (2) development of a new LApx based model for dislocation-mediated creep mechanisms in UO 2 , with application to large-grain doped UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Non-local thermodynamic equilibrium 1.5D modeling of red giant stars

Spectra for two-dimensional (2D) stars in the 1.5D approximation are created from synthetic spectra of one-dimensional (1D) non-local thermodynamic equilibrium (NLTE) spherical model atmospheres produced by the PHOENIX code. The 1.5D stars have the spatially averaged Rayleigh-Jeans flux of a K3-4 III star while varying the temperature difference between the two 1D component models (ΔT {sub 1.5D}) and the relative surface area covered. Synthetic observable quantities from the 1.5D stars are fitted with quantities from NLTE and local thermodynamic equilibrium (LTE) 1D models to assess the errors in inferred T {sub eff} values from assuming horizontal homogeneity and LTE. Five different quantities are fit to determine the T {sub eff} of the 1.5D stars: UBVRI photometric colors, absolute surface flux spectral energy distributions (SEDs), relative SEDs, continuum normalized spectra, and TiO band profiles. In all cases except the TiO band profiles, the inferred T {sub eff} value increases with increasing ΔT {sub 1.5D}. In all cases, the inferred T {sub eff} value from fitting 1D LTE quantities is higher than from fitting 1D NLTE quantities and is approximately constant as a function of ΔT {sub 1.5D} within each case. The difference between LTE and NLTE for the TiO bands is caused indirectly by the NLTE temperature structure of the upper atmosphere, as the bands are computed in LTE. We conclude that the difference between T {sub eff} values derived from NLTE and LTE modeling is relatively insensitive to the degree of the horizontal inhomogeneity of the star being modeled and largely depends on the observable quantity being fit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Stomatal response to humidity and CO 2 implicated in recent decline in US evaporation

Abstract Evapotranspiration, defined as the total flux of water from the land surface to the atmosphere, is a major component of the hydrologic cycle and surface energy balance. Although evapotranspiration is expected to intensify with increasing temperatures, long‐term, regional trends in evapotranspiration remain uncertain due to spatially and temporally limited direct measurements. In this study, we utilize an emergent relation between the land surface and atmospheric boundary layer to infer daily evapotranspiration from historical meteorological data collected at 236 weather stations across the United States. Our results suggest a statistically significant ( α = 0.05) decrease in evapotranspiration of approximately 6% from 1961 to 2014, with a significant ( α = 0.05) sharp decline of 13% from 1998 to 2014. We attribute the decrease in evapotranspiration mostly to declines in surface conductance, but also to offsetting changes in longwave radiation, wind speed, and incoming solar radiation. Using an established stomatal conductance model, we explain the changes in inferred surface conductance as a response to increases in carbon dioxide and, more recently, to an abrupt decrease in atmospheric humidity.

Rigden, Angela J.↗
DOE OSTI2022

Differentially Private Synthesis and Sharing of Network Data Via Bayesian Exponential Random Graph Models

Abstract Network data often contain sensitive relational information. One approach to protecting sensitive information while offering flexibility for network analysis is to share synthesized networks based on the information in originally observed networks. We employ differential privacy (DP) and exponential random graph models (ERGMs) and propose the DP-ERGM method to synthesize network data. We apply DP-ERGM to two real-world networks. We then compare the utility of synthesized networks generated by DP-ERGM, the DyadWise Randomized Response (DWRR) approach, and the Synthesis through Conditional distribution of Edge given nodal Attribute (SCEA) approach. In general, the results suggest that DP-ERGM preserves the original information significantly better than two other approaches in network structural statistics and inference for ERGMs and latent space models. Furthermore, DP-ERGM satisfies node DP through modeling the global network structure with ERGM, a stronger notion of privacy than the edge DP under which DWRR and SCEA operate.

graph synthesis↗
DOE OSTI2025

Deep Learning for Subsurface Flow: A Comparative Study of U‐Net, Fourier Neural Operators, and Transformers in Underground Hydrogen Storage

Subsurface flow research is essential for the sustainable management of natural resources and the environment. Deep learning (DL) has significantly advanced this field by developing efficient and accurate surrogate models to replace computationally expensive physics‐based simulations. These surrogate models are commonly used to predict the spatiotemporal evolution of state variables, such as gas saturation and reservoir pressure, in heterogeneous geological formations. Despite the various DL models applied to this task, there is a lack of studies systematically comparing their performance. This absence of comparative analysis leads to somewhat arbitrary DL model selection in subsurface flow research, resulting in suboptimal performance and potentially inaccurate predictions. To bridge this gap, we conduct a systematic comparison study of three popular DL architectures—U‐Net, Fourier Neural Operators (FNO), and Segmentation Transformer (SETR)—in surrogate modeling of underground hydrogen storage (UHS). We focus on UHS due to its promise of enhancing clean energy resilience and its cyclic operational conditions that represent common scenarios in various subsurface applications. We evaluate the models based on accuracy, training cost, and inference speed. The comparison shows that U‐Net achieves the highest accuracy, followed by SETR and FNO. Despite its lower accuracy, FNO has the highest inference speed. SETR offers competitive accuracy with the least training memory usage, demonstrating the potential of transformers in learning subsurface flow. Our results provide guidance for selecting DL models for surrogate modeling in a wide range of subsurface flow problems.

42 ENGINEERING↗
DOE OSTI2025

Rapid wavefield forecasting for earthquake early warning via deep sequence to sequence learning

We propose a deep learning model, WaveCastNet, to forecast high-dimensional wavefields. WaveCastNet integrates a convolutional long expressive memory architecture into a sequence-to-sequence forecasting framework, enabling it to model long-term dependencies and multiscale patterns in both space and time. By sharing weights across spatial and temporal dimensions, WaveCastNet requires significantly fewer parameters than more resource-intensive models such as transformers, resulting in faster inference times. Crucially, WaveCastNet also generalizes better than transformers to rare and critical seismic scenarios, such as high-magnitude earthquakes. Here, we show the ability of the model to predict the intensity and timing of destructive ground motions in real time, using simulated data from the San Francisco Bay Area. Furthermore, we demonstrate its zero-shot capabilities by evaluating WaveCastNet on real earthquake data. Our approach does not require estimating earthquake magnitudes and epicenters, steps that are prone to error in conventional methods, nor does it rely on empirical ground-motion models, which often fail to capture strongly heterogeneous wave propagation effects.

Geophysics↗
DOE OSTI2025

PRIME: Protein Representation Inference for Mutation Evaluation

Protein language machine learning models built upon existing ESM-2 model developed by Evolutionary Scale (evolutionaryscale.ai) and an in-house protein language model based on the BERT model developed by Google. The code also includes model training scripts and saved checkpoints from our own training using publicly available SARS-CoV-2 protein sequences.

Gibson, Kaetlyn [Los Alamos National Lab]↗
DOE OSTI2015

Variable protostellar accretion with episodic bursts

We present the latest development of the disk gravitational instability and fragmentation model, originally introduced by us to explain episodic accretion bursts in the early stages of star formation. Using our numerical hydrodynamics model with improved disk thermal balance and star-disk interaction, we computed the evolution of protostellar disks formed from the gravitational collapse of prestellar cores. In agreement with our previous studies, we find that cores of higher initial mass and angular momentum produce disks that are more favorable to gravitational instability and fragmentation, while a higher background irradiation and magnetic fields moderate the disk tendency to fragment. The protostellar accretion in our models is time-variable, thanks to the nonlinear interaction between different spiral modes in the gravitationally unstable disk, and can undergo episodic bursts when fragments migrate onto the star owing to the gravitational interaction with other fragments or spiral arms. Most bursts occur in the partly embedded Class I phase, with a smaller fraction taking place in the deeply embedded Class 0 phase and a few possible bursts in the optically visible Class II phase. The average burst duration and mean luminosity are found to be in good agreement with those inferred from observations of FUors. The model predicts the existence of two types of bursts: the isolated ones, showing well-defined luminosity peaks separated with prolonged periods (∼10{sup 4} yr) of quiescent accretion, and clustered ones, demonstrating several bursts occurring one after another during just a few hundred years. Finally, we estimate that 40%–70% of the star-forming cores can display bursts after forming a star-disk system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Pantheon+ Analysis: Forward Modeling the Dust and Intrinsic Color Distributions of Type Ia Supernovae, and Quantifying Their Impact on Cosmological Inferences

Abstract Recent studies have shown that the observed color distributions of Type Ia supernovae (SNe Ia) can be well described by a combination of a dust distribution and an intrinsic color distribution. Using the Pantheon+ sample of 1701 SN Ia, we apply a new forward-modeling fitting method (Dust2Dust) to measure the parent dust and color distributions, including their dependence on host-galaxy mass. At each fit step, the SN Ia selection efficiency is determined from a large simulated sample that is reweighted to reflect the proposed distributions. We use five separate metrics to describe the goodness of fit: distribution of fitted light-curve color c , cosmological residual trends with c , cosmological residual scatter with c , fitted color–luminosity relationship β SALT2 , and intrinsic scatter σ int . We present the results and the uncertainty in 12-dimensional space. Furthermore, we measure that the uncertainty on this modeling propagates to an upper threshold uncertainty in the equation of state of dark energy w of 0.014(1) for the Pantheon+ cosmology analysis and contributes negligible uncertainty to the Hubble constant H 0 . The Dust2Dust code is made publicly available at https://github.com/djbrout/dustdriver .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Emulator-based Bayesian calibration of a subglacial drainage model

Subglacial drainage models, often motivated by the relationship between hydrology and ice flow, sensitively depend on numerous unconstrained parameters. We explore using borehole water-pressure time series to calibrate the uncertain parameters of a popular subglacial drainage model, taking a Bayesian perspective to quantify the uncertainty in parameter estimates and in the calibrated model predictions. To reduce the computation time associated with Markov Chain Monte Carlo sampling, we construct a fast Gaussian process emulator to stand in for the subglacial drainage model. We first carry out a calibration experiment using synthetic observations consisting of model simulations with hidden parameter values as a demonstration of the method. Using real borehole water pressures measured in western Greenland, we find meaningful constraints on four of the eight model parameters and a factor-of-three reduction in uncertainty of the calibrated model predictions. These experiments illustrate Gaussian process-based Bayesian inference as a useful tool for calibration and uncertainty quantification of complex glaciological models using field data. However, significant differences between the calibrated model and the borehole data suggest that structural limitations of the model, rather than poorly constrained parameters or computational cost, remain the most important constraint on subglacial drainage modelling.

58 GEOSCIENCES↗
DOE OSTI2020

Revisiting the Architecture of the KOI-89 System

While high stellar obliquities observed in exoplanetary systems may be attributed to processes that tilt the planetary orbits, it is also possible that they reflect misalignments between protoplanetary disks and stellar spins. This latter hypothesis predicts the presence of coplanar multiplanetary systems misaligned with their central stars. Here we reevaluate the evidence of such an architecture that has been claimed for the KOI-89 system. An early-type star, KOI-89 has one validated transiting planet, KOI-89.01/Kepler-462b (period 84.7 days, radius 3.0 R {sub ⊕}), and one transiting planet candidate, KOI-89.02 (period 207.6 days, radius 4.0 R {sub ⊕}), where the latter exhibits transit timing variations (TTVs). A previous modeling of the stellar gravity-darkening effect in the transit light curves inferred a high stellar obliquity of ≈70°. We perform photodynamical modeling of the Kepler transit light curves and use the resulting constraints on the orbital configuration and transit times to update the gravity-darkened transit model. As a result, we find no firm evidence for the gravity-darkening effect in the transit shapes and conclude that stellar obliquity is not constrained by the data. Given the evidence for low orbital eccentricities from the dynamical analysis, the system architecture can thus be consistent with many other multitransiting systems with flat, near-circular orbits aligned with the stellar spin. We find that the TTVs imparted on its neighbor imply that KOI-89.01 has a mass ≳20 M {sub ⊕}. This would render it one of the densest known sub-Neptunes, mostly composed of a solid core. Lower masses are possible if the TTVs are instead due to an unseen third planet.

79 ASTRONOMY AND ASTROPHYSICS↗