Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “radii”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bayesian inference of fine features of the nuclear equation of state from future neutron star radius measurements to 0.1 km accuracy

To more precisely constrain the equation of state (EOS) of supradense neutron-rich nuclear matter, future high-precision x-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a metamodel EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility K 0 and skewness J 0 of symmetric nuclear matter (SNM) as well as the slope L, curvature K sym , and skewness J sym characterizing the density dependence of nuclear symmetry energy E sym ⁡(ρ), respectively, from mean values of NS radii consistent with existing observations and an expected accuracy Δ⁢R ranging from about 1.0 to 0.1 km. Here, we found that (1) the Δ⁢R has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller Δ⁢R reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density E sym ⁡(ρ), (3) a double-peak feature of the PDF(K sym ) corresponding to the strong K sym – J sym and K sym – L anticorrelations is revealed when Δ⁢R is less than about 0.2 km, and the locations of the two peaks are sensitive to the maximum value of J sym reflecting the stiffness of E sym ⁡(ρ) above about 3 times the saturation density ρ 0 of SNM, and (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2–3)⁢ρ 0 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Model-Independent Determination of Nuclear Weak Form Factors and Implications for Standard Model Precision Tests

We analyze the recoil corrections in superallowed beta decays of T = 1, J P = 0 + nuclei by fixing the mean square charged weak radius model independently using the data of multiple charge radii across the nuclear isotriplet. By comparing to model estimations, we argue that the existing theory uncertainty in the statistical rate function f might have been substantially underestimated. Here, we discuss the implications of our proposed strategy for precision tests of the standard model, including a potential alleviation of the first-row CKM unitarity deficit, and motivate new experiments for charge radii measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Universal parameters of bulk-solvent masks

The bulk solvent is a major component of biomacromolecular crystals that contributes significantly to the observed diffraction intensities. Accurate modelling of the bulk solvent has been recognized as important for many crystallographic calculations. Owing to its simplicity and modelling power, the flat (mask-based) bulk-solvent model is used by most modern crystallographic software packages to account for disordered solvent. In this model, the bulk-solvent contribution is defined by a binary mask and a scale (scattering) function. The mask is calculated on a regular grid using the atomic model coordinates and their chemical types. The grid step and two radii, solvent and shrinkage, are the three parameters that govern the mask calculation. They are highly correlated and their choice is a compromise between the computer time needed to calculate the mask and the accuracy of the mask. It is demonstrated here that this choice can be optimized using a unique value of 0.6 Å for the grid step irrespective of the data resolution, and the radii values adjusted correspondingly. The improved values were tested on a large sample of Protein Data Bank entries derived from X-ray diffraction data and are now used in the computational crystallography toolbox ( CCTBX ) and in Phenix as the default choice.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimising response to an introduction of African swine fever in wild pigs

African swine fever virus (ASFv) is a virulent pathogen that threatens domestic swine industries globally and persists in wild boar populations in some countries. Persistence in wild boar can challenge elimination and prevent disease-free status, making it necessary to address wild swine in proactive response plans. In the United States, invasive wild pigs are abundant and found across a wide range of ecological conditions that could drive different epidemiological dynamics among populations. Information on the size of the control areas required to rapidly eliminate the ASFv in wild pigs and how this area should change with management constraints and local ecology is needed to optimize response planning. We developed a spatially explicit disease transmission model contrasting wild pig movement and contact ecology in two ecosystems in Southeastern United States. We simulated ASFv spread and determined the optimal response area (reported as the radius of a circle) for eliminating ASFv rapidly over a range of detection times (when ASFv was detected relative to the true date of introduction), culling capacities (proportion of wild pigs in the culling zone removed weekly) and wild pig densities. Large radii for response areas (14 km) were needed under most conditions but could be shortened with early detection (≤ 8 weeks) and high culling capacities (≥ 15% weekly). Under most conditions, the ASFv was eliminated in less than 22 weeks using optimal control radii, although ecological conditions with high rates of wild pig movement required higher culling capacities (≥ 10% weekly) for elimination within 1 year. The results highlight the importance of adjusting response plans based on local ecology and show that wild pig movement is a better predictor of the optimal response area than the number of ASFv cases early in the outbreak trajectory. Furthermore, our framework provides a tool for determining optimal control plans in different areas, guiding expectations of response impacts, and planning resources needed for rapid elimination.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Cylindricity Sensitivity Thermal Model of the AGR-5/6/7 Experiment in the Advanced Test Reactor

The AGR-5/6/7 experiment is currently being irradiated in the Advanced Test Reactor (ATR) at the Idaho National Laboratory and is approximately 70% complete. Several fuel and material irradiation experiments have been planned for the U.S. Department of Energy Advanced Gas Reactor Fuel Development and Qualification Program, which supports the development and qualification of tristructural isotropic (TRISO) coated particle fuel for use in high-temperature gas-cooled reactors. The goals of these experiments are to provide irradiation performance data to support fuel process development, qualify fuel for normal operating conditions, support development of fuel performance models and codes, and provide irradiated fuel and materials for post-irradiation examination and safety testing. Originally planned and named as separate fuel experiments, but subsequently combined into a single test train, AGR-5/6/7 is testing low-enriched uranium oxycarbide TRISO fuel. The AGR-5/6/7 test train has 5 capsules with thermocouples and independent gas control mixtures. Unique to this paper is a sensitivity study concerning the cylindricity of the graphite holders containing the fuel compacts and their eccentricity in relation to the stainless-steel capsule walls. Each capsule has small nubs on the outside used for centering the graphite holder inside the stainless-steel capsule with a small gas gap used to control temperature. Due to machining tolerances of these nubs and vibration wearing the nubs down when the experiment is running in the reactor, the possibility exists that the holder may move around radially. Each capsule is equipped with several thermocouples placed at various radii and depths within each graphite holder. This paper will show the sensitivity of offsetting the graphite holder for various radii in 45-degree increments around the circle with the objective of minimizing the difference between the measured thermocouples and the modeled thermocouple temperatures. Separate gas mixtures of helium/neon are introduced into this gas gap between the holder and capsule wall and changed as necessary to maintain the desired thermocouple temperatures to keep the fuel compacts at constant temperature as the nuclear reactor conditions change. The goal of the sensitivity study is to find a radius and an angle to offset the holder from perfectly centered for each of the five capsules separately. The complex thermal model includes fission heating, gamma heating, radiation heat transfer, and heat transfer via conduction and radiation across the control gaps. Subroutines linked to the thermal model offer an easy method to offset the graphite holder from the capsule walls without remeshing the entire model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Charging and ion ejection dynamics of large helium nanodroplets exposed to intense femtosecond soft xray pulses

Ion ejection from charged helium nanodroplets exposed to intense femtosecond soft X-ray pulses is studied by single-pulse ion time-of-flight (TOF) spectroscopy in coincidence with small-angle X-ray scattering. Scattering images encode the droplet size and absolute photon flux incident on each droplet, while ion TOF spectra are used to determine the maximum ion kinetic energy, Ekin, of He+ j fragments (j = 1-4). Measurements span HeN droplet sizes between N similar to 10(7) and similar to 10(10) (radii R-0 = 78-578 nm), and droplet charges between similar to 9x10(-5) and similar to 3x10(-3) e/atom. Conditions encompass a wide range of ionization and expansion regimes, from departure of all photoelectrons from the droplet, leading to pure Coulomb explosion, to substantial electron trapping by the electrostatic potential of the charged droplet, indicating the onset of hydrodynamic expansion. The unique combination of absolute Xray intensities, droplet sizes, and ion Ekin on an event-by-event basis reveals a detailed picture of the correlations between the ionization conditions and the ejection dynamics of the ionic fragments. The maximum E-kin of He+ is found to be governed by Coulomb repulsion from unscreened cations across all expansion regimes. The impact of ion-atom interactions resulting from the relatively low charge densities is increasingly relevant with less electron trapping. The findings are consistent with the emergence of a charged spherical shell around a quasineutral plasma core as the degree of ionization increases. The results demonstrate a complex relationship between measured ion E-kin and droplet ionization conditions that can only be disentangled through the use of coincident single-pulse TOF and scattering data.

Helium↗
DOE OSTI2021

Neutron Stars and the Nuclear Matter Equation of State

Neutron stars provide a window into the properties of dense nuclear matter. Several recent observational and theoretical developments provide powerful constraints on their structure and internal composition. Among these are the first observed binary neutron star merger, GW170817, whose gravitational radiation was accompanied by electromagnetic radiation from a short γ-ray burst and an optical afterglow believed to be due to the radioactive decay of newly minted heavy r-process nuclei. These observations give important constraints on the radii of typical neutron stars and on the upper limit to the neutron star maximum mass and complement recent pulsar observations that established a lower limit. Pulse-profile observations by the Neutron Star Interior Composition Explorer (NICER) X-ray telescope provide an independent, consistent measure of the neutron star radius. Theoretical many-body studies of neutron matter reinforce these estimates of neutron star radii. Studies using parameterized dense matter equations of state (EOSs) reveal several EOS-independent relations connecting global neutron star properties.

Physics↗
DOE OSTI2021

Characterization of Nanoscale Pores in Tight Gas Sandstones Using Complex Techniques: A Case Study of a Linxing Tight Gas Sandstone Reservoir

Pore structures with rich nanopores and permeability in tight gas reservoirs are poorly understood up to date. Advanced techniques are needed to be employed to accurately characterize pore structures, especially tiny pores which include micron and nanopores. In this study, various experimental techniques such as scanning electron microscopy (SEM), nuclear magnetic resonance (NMR) T 2 , nitrogen adsorption method, and NMR cryoporometry (NMRC) are combined to interrogate the complex pore systems of the tight gas reservoir in the Linxing formation, Ordos Basin, China. Results show that tight gas sandstones are primarily comprised of residual interparticle and clay-dominated pores. Clay and quartz are two dominate minerals while pyrite occupies a nontrivial amount as well. The permeability of tight gas sandstones is very low, exhibiting an extremely poor positive correlation with porosity. While pore types and relative pore contents are more influential factors on the permeability, accurate characterization of pore size distribution is critical for the permeability of tight gas sandstones. Therefore, complementary characterization methods are carried out, indicating that neither small pores with radii < 100 nm (around peak 1 in NMR T 2 distribution) nor large pores with radii > 5 μ m (around peak 3 in NMR T 2 distribution) control the permeability by analyzing the connectivity of the pores in various size ranges, but rather pores averaging approximately 350 ± X nm (around peak 2 in NMR T 2 distribution) have sufficient connectivity to host and transmit hydrocarbons. The pore size of tight gas sandstones is dominated by the clay-rich mineral assemblage. The study shows that the NMRC technique can be a very promising method, especially when referred to as a promising “roadmap” on how to interrogate tight formations such as the tight gas sands or even shale especially for the nanopore characterization.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Cloud microphysical structure analysis based on high-resolution insitu measurements

This study investigates the microphysical parameters and shapes of droplet size distributions (DSDs) along three aircraft traverses of developing convective clouds during Cloud Aerosol Interactions and Precipitation Enhancement Experiment (CAIPEEX) Phase III 2015 at a sampling frequency of 25 Hz. The droplet number concentration (Nc, cm -3 ), and liquid water content (LWC, g m -3 ) present steep gradients within zones of a few tens of meters near the cloud edges and relatively gentle gradients in the strong updraft zones. Sometimes, the horizontal LWC distribution resembles a trapezoidal shape with steep LWC and Nc gradients near the cloud edges. The LWC maximums (LWC max ) are lower than the adiabatic LWCs, but the high adiabatic fractions in the cloud core indicate low dilution. High LWC/LWC max , largest droplets, and narrow and similarly shaped DSDs are found in the regions of high updrafts. Zones of low LWC/LWC max are found close to the cloud edges, where DSDs are highly diverse, containing both large and small droplets. Finally, we analyze the mixing diagrams. Significant in-phase turbulent fluctuations in LWC and N c were found. The effective radii change slightly across cloud updraft zones but decrease at the zone of low LWC/LWC max ratio close to cloud edges. The spectra of LWC and N c obey Kolmogorov -5/3 turbulence law. The radii of the correlation of LWC and N c in updraft zones are of several tens of meters. Filaments up to 120–175 m in size are also noticed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Climatology of Marine Boundary Layer Cloud and Drizzle Properties derived from Ground-based Observations over the Azores

In this study, more than four years of ground-based observations and retrievals were collected and analyzed to investigate the seasonal and diurnal variations of single-layered MBL (with three subsets: non-drizzling, virga and rain) cloud and drizzle properties, as well as their vertical and horizontal variations. The annual mean drizzle frequency was ~55%, with ~70% in winter and ~45% in summer. The cloud-top (-base) height for rain clouds was the highest (lowest), resulted in the deepest cloud layer (0.8 km), which is four (two) times of non-drizzling (virga) clouds. The retrieved cloud-droplet effective radii (rc) were the largest (smallest) for rain (non-drizzling) clouds, and the nighttime values were greater than the daytime values. Drizzle number concentration (Nd) and liquid water content (LWCd) were three and one order(s) lower than their cloud counterparts. The rc and LWCc increased from the cloud base to ???? = ~0.75 by condensational growth, while drizzle median radii (rd) increased from the cloud top downwards the cloud base by collision-coalescence. The adiabaticity values monotonically increased from the cloud top to the cloud base with maxima of ~0.7 (0.3) for non-drizzling (rain) clouds. The drizzling process decreases the adiabaticity by 0.25 to 0.4, and the cloud-top entrainment mixing impacts as deep as upper 40% of the cloud layers. Cloud and drizzle homogeneities decreased with increased horizontal sampling lengths. Cloud homogeneity increases with increasing cloud fraction. These results can serve as baselines for studying MBL cloud-to-rain conversion and growth processes over the Azores.

Wu, Peng↗
DOE OSTI2022

Analysis of Defects in Metal Additive Manufacturing with Augmented Data Generation

Laser powder bed fusion (LPBF) is a method of additive manufacturing (AM) that selectively melts and fuses together microscopic metallic powder. LPBF offers the benefit of producing custom structures out of high strength metals that can be difficult to fabricate with conventional methods. The challenge of LPBF is that 3D printed structures often have internal pores due to process flaws. Pulsed thermal tomography (PTT) is a method for reconstructing the depth profile of materials, allowing the visualization internal voids in solids. In prior work, we developed a convolutional neural network (CNN) which, having been trained on simulated 2D PTT images of subsurface elliptical defects, was able to classify the semi-major radii, semi-minor radii, and angular orientation of the best-fit ellipses in previously unseen PTT images. The unseen PTT images contained subsurface irregular defect shapes imported from scanning electron microscopy (SEM) images of metallic LPBF-printed specimens. Training the CNN on irregular defect shapes instead of on elliptical shapes would make the resulting classifications more descriptive of actual defect shapes. However, this requires a much higher volume of SEM images of material defects, which are difficult to obtain because of random occurrence of defects in LPBF. To address this challenge, we developed a generative adversarial network (GAN) to augment the existing dataset of SEM defect images. The GAN model is demonstrated to create novel yet realistic defect shapes that can be used as input for simulated PTT images to train CNN.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
DOE OSTI2026

Bimetallic Coinage Metal Complexes of Tetradentate ortho Ter arylene Ligands

The synthesis, structural and photophysical properties of five bimetallic coinage metal complexes of three tetradentate charge-transfer (CT) ligands are presented. The ligand structures incorporate donor/acceptor moieties electronically decoupled by ortho connection to a bridging biphenylene group. The M M distances observed in the crystal structures of the bimetallic complexes are shorter than the sum of the Van der Waals radii of the metal ions, suggesting stabilizing interactions between the d10 metal centers. Additionally, short distances between the bridging biphenylene and metal atoms suggest arylene metal interactions. Bimetallic Cu, Ag, and Au complexes of the ligand incorporating a carbazole donor and pyridyl acceptor (N2R22-) show absorption and emission energies that trend with the ionization potential of each metal, indicating metal character in the excited state. These properties contrast with our previously reported mononuclear Zn complex of N2R22-The structural and photophysical properties of five bimetallic d10 coinage metal complexes were presented. The M M distances observed in the crystal structures of the bimetallic complexes are shorter than the sum of the metals’ Van der Waals radii, suggesting attractive interactions between the d10 metal centers. Additionally, short distances between the bridging biphenylene and metal atoms suggest arylene metal interactions. The absorption and emission spectra of the bimetallic complexes redshift based on the oxidation potential of the metal ion. This suggests involvement of metal based orbitals in the excited state transitions of the bimetallic complexes, which was corroborated by molecular modeling. High photoluminescent quantum yields were observed in solution and in rigid media (up to PL = 0.75), and the excited state lifetimes in solution were found to exceed several microseconds.

Shariaty, Darius Allen [University of Southern Cal↗
DOE OSTI2023

Density functional modeling of the binding energies between aluminosilicate oligomers and different metal cations

Interactions between negatively charged aluminosilicate species and positively charged metal cations are critical to many important engineering processes and applications, including sustainable cements and aluminosilicate glasses. In an effort to probe these interactions, here we have calculated the pair-wise interaction energies (i.e., binding energies) between aluminosilicate dimer/trimer and 17 different metal cations M n+ (M n+ = Li + , Na + , K + , Cu + , Cu 2+ , Co 2+ , Zn 2+ , Ni 2+ , Mg 2+ , Ca 2+ , Ti 2+ , Fe 2+ , Fe 3+ , Co 3+ , Cr 3+ , Ti 4+ and Cr 6+ ) using a density functional theory (DFT) approach. Analysis of the DFT-optimized structural representations for the clusters (dimer/trimer + M n+ ) shows that their structural attributes (e.g., interatomic distances) are generally consistent with literature observations on aluminosilicate glasses. The DFT-derived binding energies are seen to vary considerably depending on the type of cations (i.e., charge and ionic radii) and aluminosilicate species (i.e., dimer or trimer). A survey of the literature reveals that the difference in the calculated binding energies between different M n+ can be used to explain many literature observations associated with the impact of metal cations on materials properties (e.g., glass corrosion, mineral dissolution, and ionic transport). Analysis of all the DFT-derived binding energies reveals that the correlation between these energy values and the ionic potential and field strength of the metal cations are well captured by 2nd order polynomial functions ( R 2 values of 0.99–1.00 are achieved for regressions). Given that the ionic potential and field strength of a given metal cation can be readily estimated using well-tabulated ionic radii available in the literature, these simple polynomial functions would enable rapid estimation of the binding energies of a much wider range of cations with the aluminosilicate dimer/trimer, providing guidance on the design and optimization of sustainable cements and aluminosilicate glasses and their associated applications. Finally, the limitations associated with using these simple model systems to model complex interactions are also discussed.

Gong, Kai↗
DOE OSTI2023

Bayes goes fast: Uncertainty quantification for a covariant energy density functional emulated by the reduced basis method

A covariant energy density functional is calibrated using a principled Bayesian statistical framework informed by experimental binding energies and charge radii of several magic and semi-magic nuclei. The Bayesian sampling required for the calibration is enabled by the emulation of the high-fidelity model through the implementation of a reduced basis method (RBM)—a set of dimensionality reduction techniques that can speed up demanding calculations involving partial differential equations by several orders of magnitude. The RBM emulator we build—using only 100 evaluations of the high-fidelity model—is able to accurately reproduce the model calculations in tens of milliseconds on a personal computer, an increase in speed of nearly a factor of 3,300 when compared to the original solver. Besides the analysis of the posterior distribution of parameters, we present model calculations for masses and radii with properly estimated uncertainties. We also analyze the model correlation between the slope of the symmetry energy L and the neutron skin of 48 Ca and 208 Pb. The straightforward implementation and outstanding performance of the RBM makes it an ideal tool for assisting the nuclear theory community in providing reliable estimates with properly quantified uncertainties of physical observables. Such uncertainty quantification tools will become essential given the expected abundance of data from the recently inaugurated and future experimental and observational facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Visual Orbits of Spectroscopic Binaries with the CHARA Array. III. HD 8374 and HD 24546

We present the visual orbits of two long-period spectroscopic binary stars, HD 8374 and HD 24546, using interferometric observations acquired with the CHARA Array and the Palomar Testbed Interferometer. We also obtained new radial velocities from echelle spectra using the APO 3.5 m and Fairborn 2.0 m telescopes. By combining the visual and spectroscopic observations, we solve for the full, three-dimensional orbits and determine the stellar masses and distances to within 3% uncertainty. We then estimate the effective temperature and radius of each component star through Doppler tomography and spectral energy distribution analyses, in order to compare the observed stellar parameters to the predictions of stellar evolution models. For HD 8374, we find masses of M {sub 1} = 1.636 ± 0.050M {sub ⊙} and M {sub 2} = 1.587 ± 0.049M {sub ⊙}, radii of R {sub 1} = 1.84 ± 0.05R {sub ⊙} and R {sub 2} = 1.66 ± 0.12R {sub ⊙}, temperatures of T{sub eff1}=7280±110 K and T{sub eff2}=7280±120 K, and an estimated age of 1.0 Gyr. For HD 24546, we find masses of M {sub 1} = 1.434 ± 0.014M {sub ⊙} and M {sub 2} = 1.409 ± 0.014M {sub ⊙}, radii of R {sub 1} = 1.67 ± 0.06R {sub ⊙} and R {sub 2} = 1.60 ± 0.10R {sub ⊙}, temperatures of T{sub eff1}=6790±120 K and T{sub eff2}=6770±90 K, and an estimated age of 1.4 Gyr. HD 24546 is therefore too old to be a member of the Hyades cluster, despite its physical proximity to the group.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Gaia–Kepler Stellar Properties Catalog. II. Planet Radius Demographics as a Function of Stellar Mass and Age

Studies of exoplanet demographics require large samples and precise constraints on exoplanet host stars. Using the homogeneous Kepler stellar properties derived using the Gaia Data Release 2 by Berger et al., we recompute Kepler planet radii and incident fluxes and investigate their distributions with stellar mass and age. We measure the stellar mass dependence of the planet radius valley to be dlogR{sub p}/dlogM{sub ⋆} = 0.26{sub −0.16}{sup +0.21}, consistent with the slope predicted by a planet mass dependence on stellar mass (0.24–0.35) and core-powered mass loss (0.33). We also find the first evidence of a stellar age dependence of the planet populations straddling the radius valley. Specifically, we determine that the fraction of super-Earths (1–1.8 R{sub ⊕}) to sub-Neptunes (1.8–3.5 R{sub ⊕}) increases from 0.61 ± 0.09 at young ages (<1 Gyr) to 1.00 ± 0.10 at old ages (>1 Gyr), consistent with the prediction by core-powered mass loss that the mechanism shaping the radius valley operates over Gyr timescales. Additionally, we find a tentative decrease in the radii of relatively cool (F {sub p} < 150 F{sub ⊕}) sub-Neptunes over Gyr timescales, which suggests that these planets may possess H/He envelopes instead of higher mean molecular weight atmospheres. We confirm the existence of planets within the hot sub-Neptunian “desert” (2.2 R {sub ⊕} < R {sub p} < 3.8 R{sub ⊕}, F {sub p} > 650 F{sub ⊕}) and show that these planets are preferentially orbiting more evolved stars compared to other planets at similar incident fluxes. In addition, we identify candidates for cool (F {sub p} < 20 F{sub ⊕}) inflated Jupiters, present a revised list of habitable zone candidates, and find that the ages of single and multiple transiting planet systems are statistically indistinguishable.

79 ASTRONOMY AND ASTROPHYSICS↗