Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bayesian Analysis”

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 343 records · Page 19

Optimization of artificial viscosity in production codes based on Gaussian Regression surrogate models

To accurately model flows with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity has to be introduced to convert kinetic energy into internal energy, thereby increasing the entropy across shocks. Determining the appropriate strength of the artificial viscosity is an art and strongly depends on the particular problem and experience of the researcher. The objective of this study is to pose the problem of finding the appropriate strength of artificial viscosity as an optimization problem and solve this problem using machine learning (ML) tools, specifically using surrogate models based on Gaussian Process regression and Bayesian analysis. We describe the optimization method and discuss various practical details of its implementation. The shock-containing problems for which we apply this method all have been implemented in the LANL code FLAG. First, we apply ML to find optimal values to isolated shock problems of different strengths. Second, we apply ML to optimize viscosity for a 1D propagating detonation problem based on Zel’dovich-von Neumann-Doring (ZND) detonation theory using a reactive burn model. We compare results for default (currently used values in FLAG) and optimized values of artificial viscosity for these problems demonstrating the potential for significant improvement in the accuracy of computations.

42 ENGINEERING↗

MultiPEM Toolbox: User Manual [Rev. 2]

This document explains use of the Multi-Phenomenology Explosion Monitoring (Multi PEM) Toolbox, a collection of R scripts for estimating the unknown device parameters of a new event with uncertainty quantification. The methodology and application used for illustration in this user manual are fully documented in a Los Alamos National Laboratory technical report hereafter designated “WPA” for reference. Additional details on the application are found in a recent journal article. Two assessment types are available: rapid and complete. Rapid assessments are conducted in two stages, as described in Section 2. In the first stage, calibration data are used to estimate forward and error model parameters (WPA, §5.1) and (if relevant) errors-in-variables yield values for calibration sources (WPA, §3, Equation (3)). In the second stage, new event data are used to estimate the unknown new event device parameters (WPA, §5.2) with uncertainty quantification. Two options for treating the inferred first stage parameters in second stage Bayesian analysis are available: fixing them at their maximum likelihood estimate (default), or multiple imputation. Multiple imputation involves utilizing several posterior samples (imputations) of the first stage parameters as fixed values in the second stage posterior sampling of the new event device parameters. Second stage sampling is conducted across imputations in parallel to improve computational efficiency. This method produces improved uncertainty quantification of the new event device parameters compared with the default treatment of the first stage parameters, at the expense of additional computation. Complete assessments are conducted in a single stage, as described in Section 3. Calibration and (if relevant) new event data are used simultaneously to estimate all forward model, error model, and (if relevant) new event device parameters with uncertainty quantification on the latter. As the name suggests, rapid assessments generally run substantially faster than complete assessments (even with multiple imputation), because the results of first stage analysis can be stored and incorporated into estimating a relatively low-dimensional space of new event device parameters whenever relevant new event data becomes available. On the other hand, complete assessments must be run on the full set of model and device parameters with calibration and new event data every time the latter becomes available.

97 MATHEMATICS AND COMPUTING↗

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↗

3-flavour results with NOvA

The NOvA collaboration has presented the 3-flavour oscillation results using a frequentist approach before. In this contribution, we present the results obtained with a bayesian approach. The conclusions with respect to the frequentist approach are unchanged, that is, the upper octant of $\theta_{23}$ is preferred, $\delta_\text{CP}=\frac{3\pi}{2}$ is not preferred at the $2\sigma$ level with normal ordering and $\delta_\text{CP}=\frac{\pi}{2}$ excluded at $3\sigma$ with inverted ordering. As a novelty of the bayesian analysis, the results obtained when $\theta_{13}$ is left as a free parameter in the fit are also.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the processes generating latitudinal richness gradients: identifying diagnostic patterns and predictions

Many processes have been put forward to explain the latitudinal gradient in species richness. Here, we use a simulation model to examine four of the most common hypotheses and identify patterns that might be diagnostic of those four hypotheses. The hypotheses examined include (1) tropical niche conservatism, or the idea that the tropics are more diverse because a tropical clade origin has allowed more time for diversification in the tropics and has resulted in few species adapted to extra-tropical climates. (2) The productivity, or energetic constraints, hypothesis suggests that species richness is limited by the amount of biologically available energy in a region. (3) The tropical stability hypothesis argues that major climatic fluctuations and glacial cycles in extratropical regions have led to greater extinction rates and less opportunity for specialization relative to the tropics. (4) Finally, the speciation rates hypothesis suggests that the latitudinal richness gradient arises from a parallel gradient in rates of speciation. We found that tropical niche conservatism can be distinguished from the other three scenarios by phylogenies which are more balanced than expected, no relationship between mean root distance and richness across regions, and a homogeneous rate of speciation across clades and through time. The energy gradient, speciation gradient, and disturbance gradient scenarios all exhibited phylogenies which were more imbalanced than expected, showed a negative relationship between mean root distance and richness, and diversity-dependence of speciation rate estimates through time. Using Bayesian Analysis of Macroevolutionary Mixtures on the simulated phylogenies, we found that the relationship between speciation rates and latitude could distinguish among these three scenarios. We emphasize the importance of considering multiple hypotheses and focusing on diagnostic predictions instead of predictions that are consistent with more than one hypothesis.

biodiversity↗

DETERMINING AGES OF APOGEE GIANTS WITH KNOWN DISTANCES

We present a sample of 705 local giant stars observed using the New Mexico State University 1 m telescope with the Sloan Digital Sky Survey-III/Apache Point Observatory Galactic Evolution Experiment (APOGEE) spectrograph, for which we estimate stellar ages and the local star formation history (SFH). The high-resolution (R ∼ 22,500), near infrared (1.51–1.7 μm) APOGEE spectra provide measurements of stellar atmospheric parameters (temperature, surface gravity, [M/H], and [α/M]). Due to the smaller uncertainties in surface gravity possible with high-resolution spectra and accurate Hipparcos distance measurements, we are able to calculate the stellar masses to within 30%. For giants, the relatively rapid evolution up the red giant branch allows the age to be constrained by the mass. We examine methods of estimating age using both the mass–age relation directly and a Bayesian isochrone matching of measured parameters, assuming a constant SFH. To improve the SFH prior, we use a hierarchical modeling approach to constrain the parameters of the model SFH using the age probability distribution functions of the data. The results of an α-dependent Gaussian SFH model show a clear age–[α/M] relation at all ages. Using this SFH model as the prior for an empirical Bayesian analysis, we determine ages for individual stars. The resulting age–metallicity relation is flat, with a slight decrease in [M/H] at the oldest ages and a ∼0.5 dex spread in metallicity across most ages. For stars with ages ≲1 Gyr we find a smaller spread, consistent with radial migration having a smaller effect on these young stars than on the older stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Denser sampling of the Rosette Nebula with Faraday rotation measurements: Improved estimates of magnetic fields in H ii regions

We report Faraday rotation measurements of 11 extragalactic radio sources with lines of sight through the Rosette Nebula, a prominent H ii region associated with the star cluster NGC 2244. The goal of these measurements is to better determine the strength and structure of the magnetic field in the nebula. We calculate the rotation measure (RM) through two methods, a least-squares fit to χ(λ{sup 2}) and Rotation Measure Synthesis. In conjunction with our results from Savage et al., we find an excess RM due to the shell of the nebula of +40 to +1200 rad m{sup −2} above a background RM of +147 rad m{sup −2}. We discuss two forms of a simple shell model intended to reproduce the magnitude of the observed RM as a function of distance from the center of the Rosette Nebula. The models represent different physical situations for the magnetic field within the shell of the nebula. The first assumes that there is an increase in the magnetic field strength and plasma density at the outer radius of the H ii region, such as would be produced by a strong magnetohydrodynamic shock wave. The second model assumes that any increase in the RM is due solely to an increase in the density, and the Galactic magnetic field is unaffected in the shell. We employ a Bayesian analysis to compare the two forms of the model. The results of this analysis were inconclusive, although the model without amplification of the interstellar magnetic field is weakly favored.

79 ASTRONOMY AND ASTROPHYSICS↗

The next generation virgo cluster survey (NGVS). XIII. The luminosity and mass function of galaxies in the core of the virgo cluster and the contribution from disrupted satellites

We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg{sup 2} (0.3 Mpc{sup 2}) area in the core of the Virgo Cluster, based on u{sup ∗} griz data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample—which consists of 352 objects brighter than M {sub g} = −9.13 mag, the 50% completeness limit of the survey—reaches 2.2 mag deeper than the widely used Virgo Cluster Catalog and at least 1.2 mag deeper than any sample previously used to measure the luminosity function in Virgo. Using a Bayesian analysis, we find a best-fit faint-end slope of α = −1.33 ± 0.02 for the g-band luminosity function; consistent results are found for the stellar mass function and the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultracompact dwarfs (UCDs)—previously compiled by the NGVS in this region—to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, M {sub g} = −9.6 mag) increases dramatically, up to α = −1.60 ± 0.06 when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo owing to the disruption of satellites, both nucleated and non-nucleated. We estimate that ∼150 objects with M {sub g} ≲ −9.6 mag and that are currently classified as globular clusters might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the Virgo core might once have belonged to such satellites; these same disrupted satellites might have contributed ∼40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, α = −1.21 ± 0.05, which is only 1.7σ shallower than measured in the core of the Virgo Cluster.

79 ASTRONOMY AND ASTROPHYSICS↗

OGLE-2012-BLG-0724Lb: a Saturn-mass planet around an M dwarf

We report the discovery of a planet by the microlensing method, OGLE-2012-BLG-0724Lb. Although the duration of the planetary signal for this event was one of the shortest seen for a planetary event, the anomaly was well covered thanks to high-cadence observations taken by the survey groups OGLE and MOA. By analyzing the light curve, this planetary system is found to have a mass ratio q=(1.58±0.15)×10{sup −3}. By conducting a Bayesian analysis, we estimate that the host star is an M dwarf with a mass of M{sub L}=0.29{sub −0.16}{sup +0.33} M{sub ⊙} located at D{sub L}=6.7{sub −1.2}{sup +1.1} kpc away from the Earth and the companion’s mass is m{sub P}=0.47{sub −0.26}{sup +0.54} M{sub Jup}. The projected planet–host separation is a{sub ⊥}=1.6{sub −0.3}{sup +0.4} AU. Because the lens–source relative proper motion is relatively high, future high-resolution images would detect the lens host star and determine the lens properties uniquely. This system is likely a Saturn-mass exoplanet around an M dwarf, and such systems are commonly detected by gravitational microlensing. This adds another example of a possible pileup of sub-Jupiters (0.2<m{sub P}/M{sub Jup}<1) in contrast to a lack of Jupiters (∼1--2 M{sub Jup}) around M dwarfs, supporting the prediction by core accretion models that Jupiter-mass or more massive planets are unlikely to form around M dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian mass estimates of the Milky Way: the dark and light sides of parameter assumptions

We present mass and mass profile estimates for the Milky Way (MW) Galaxy using the Bayesian analysis developed by Eadie et al. and using globular clusters (GCs) as tracers of the Galactic potential. The dark matter and GCs are assumed to follow different spatial distributions; we assume power-law model profiles and use the model distribution functions described in Evans et al. and Deason et al. We explore the relationships between assumptions about model parameters and how these assumptions affect mass profile estimates. We also explore how using subsamples of the GC population beyond certain radii affect mass estimates. After exploring the posterior distributions of different parameter assumption scenarios, we conclude that a conservative estimate of the Galaxy’s mass within 125 kpc is 5.22×10{sup 11} M{sub ⊙}, with a 50% probability region of (4.79,5.63)×10{sup 11} M{sub ⊙}. Extrapolating out to the virial radius, we obtain a virial mass for the MW of 6.82×10{sup 11} M{sub ⊙} with 50% credible region of (6.06,7.53)×10{sup 11} M{sub ⊙}(r{sub vir}=185{sub −7}{sup +7} kpc). If we consider only the GCs beyond 10 kpc, then the virial mass is 9.02 (5.69,10.86)×10{sup 11} M{sub ⊙} (r{sub vir}=198{sub −24}{sup +19} kpc). We also arrive at an estimate of the velocity anisotropy parameter β of the GC population, which is β=0.28 with a 50% credible region (0.21, 0.35). Interestingly, the mass estimates are sensitive to both the dark matter halo potential and visible matter tracer parameters, but are not very sensitive to the anisotropy parameter.

79 ASTRONOMY AND ASTROPHYSICS↗

The molecular baryon cycle of M82

Baryons cycle into galaxies from the intergalactic medium and are converted into stars; a fraction of the baryons are ejected out of galaxies by stellar feedback. Here we present new high-resolution (3.″9; 68 pc) {sup 12}CO(2–1) and {sup 12}CO(3–2) images that probe these three stages of the baryon cycle in the nearby starburst M82. We combine these new observations with previous {sup 12}CO(1–0) and [Fe ii] images to study the physical conditions within the molecular gas. Using a Bayesian analysis and the radiative transfer code RADEX, we model temperatures and densities of molecular hydrogen, as well as column densities of CO. Besides the disk, we concentrate on two regions within the galaxy: an expanding super-bubble and the base of a molecular streamer. Shock diagnostics, kinematics, and optical extinction suggest that the streamer is an inflowing filament, with a mass inflow rate of molecular gas of 3.5 M{sub ⊙} yr{sup −1}. We measure the mass outflow rate of molecular gas of the expanding super-bubble to be 17 M{sub ⊙} yr{sup −1}, five times higher than the inferred inflow rate and 1.3 times the star formation rate of the galaxy. The high mass outflow rate and large star formation rate will deplete the galaxy of molecular gas within eight million years, unless there are additional sources of molecular gas.

79 ASTRONOMY AND ASTROPHYSICS↗

MOA-2012-BLG-505Lb: A Super-Earth-mass Planet That Probably Resides in the Galactic Bulge

We report the discovery of a super-Earth-mass planet in the microlensing event MOA-2012-BLG-505. This event has the second shortest event timescale of t {sub E} = 10 ± 1 days where the observed data show evidence of a planetary companion. Our 15 minute high cadence survey observation schedule revealed the short subtle planetary signature. The system shows the well known close/wide degeneracy. The planet/host-star mass ratio is q = 2.1 × 10{sup −4} and the projected separation normalized by the Einstein radius is s = 1.1 or 0.9 for the wide and close solutions, respectively. We estimate the physical parameters of the system by using a Bayesian analysis and find that the lens consists of a super-Earth with a mass of 6.7{sub −3.6}{sup +10.7} M{sub ⊕} orbiting around a brown dwarf or late-M-dwarf host with a mass of 0.10{sub −0.05}{sup +0.16} M{sub ⊙} with a projected star–planet separation of 0.9{sub −0.2}{sup +0.3} au. The system is at a distance of 7.2 ± 1.1 kpc, i.e., it is likely to be in the Galactic bulge. The small angular Einstein radius (θ {sub E} = 0.12 ± 0.02 mas) and short event timescale are typical for a low-mass lens in the Galactic bulge. Such low-mass planetary systems in the Bulge are rare because the detection efficiency of planets in short microlensing events is relatively low. This discovery may suggest that such low-mass planetary systems are abundant in the Bulge and currently on-going high cadence survey programs will detect more such events and may reveal an abundance of such planetary systems.

79 ASTRONOMY AND ASTROPHYSICS↗

OGLE-2013-BLG-0132Lb and OGLE-2013-BLG-1721Lb: Two Saturn-mass Planets Discovered around M-dwarfs

We present the discovery of two planetary systems consisting of a Saturn-mass planet orbiting an M-dwarf, which were detected in faint microlensing events OGLE-2013-BLG-0132 and OGLE-2013-BLG-1721. The planetary anomalies were covered with high cadence by Optical Gravitational Lensing Experiment (OGLE) and Microlensing Observations in Astrophysics (MOA) photometric surveys. The light curve modeling indicates that the planet-to-host mass ratios are (5.15±0.28)×10{sup -4} and (13.18±0.72)×10{sup -4}, respectively. Both events were too short and too faint to measure a reliable parallax signal and hence the lens mass. We therefore used a Bayesian analysis to estimate the masses of both planets: 0.29{sub -0.13}{sup +0.16} M{sub Jup} (OGLE-2013-BLG-0132Lb) and 0.64{sub -0.31}{sup +0.35} M{sub Jup} (OGLE-2013-BLG-1721Lb). Thanks to a high relative proper motion, OGLE-2013-BLG-0132 is a promising candidate for the high-resolution imaging follow-up. Both planets belong to an increasing sample of sub-Jupiter-mass planets orbiting M-dwarfs beyond the snow line.

79 ASTRONOMY AND ASTROPHYSICS↗

Two Jupiter-mass Planets Discovered by the KMTNet Survey in 2017

We report two microlensing events, KMT-2017-BLG-1038 and KMT-2017-BLG-1146, that are caused by planetary systems. These events were discovered by Korea Microlensing Telescope Network survey observations from the 2017 bulge season. The discovered systems consist of a planet and host star with mass ratios of 5.3{sub −0.4}{sup +0.2}×10{sup −3} and 2.0{sub −0.1}{sup +0.6}×10{sup −3}, respectively. Based on a Bayesian analysis assuming a Galactic model without stellar remnant hosts, we find that the planet KMT-2017-BLG-1038Lb is a super-Jupiter-mass planet (M{sub p} =2.04{sub −1.15}{sup +2.02} M{sub J}) orbiting a mid-M dwarf host (M{sub h}=0.37{sub −0.20}{sup +0.36} M{sub ⊙}) that is located at 6.01{sub −1.72}{sup +1.27} kpc toward the Galactic bulge. The other planet, KMT-2017-BLG-1146Lb, is a sub-Jupiter-mass planet (M{sub p}=0.71{sub −0.42}{sup +0.80} M{sub J}) orbiting a mid-M dwarf host (M{sub h}=0.33{sub −0.20}{sup +0.36} M{sub ⊙}) at a distance of 6.50{sub −2.00}{sup +1.38} kpc toward the Galactic bulge. Both are potentially gaseous planets that are beyond their hosts’ snow lines. These typical microlensing planets will be routinely discovered by second-generation microlensing surveys, rapidly increasing the number of detections.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian Characterization of Main-sequence Binaries in the Old Open Cluster NGC 188

The binary fractions of open and globular clusters yield powerful constraints on their dynamical state and evolutionary history. We apply publicly available Bayesian analysis tools to a UBV RIJHK {sub S} photometric catalog of the open cluster NGC 188 to detect and characterize photometric binaries along the cluster main sequence. This technique has the advantage of self-consistently handling photometric errors, missing data in various bandpasses, and star-by-star prior constraints on cluster membership. Simulations are used to verify uncertainties and quantify selection biases in our analysis, illustrating that among binaries with mass ratios >0.5, we recover the binary fraction to better than 7% in the mean, with no significant dependence on binary fraction and a mild dependence on assumed mass-ratio distribution. Using our photometric catalog, we recover the majority (65% ± 11%) of spectroscopically identified main-sequence binaries, including eight of the nine with spectroscopically measured mass ratios. Accounting for incompleteness and systematics, we derive a mass-ratio distribution that rises toward lower mass ratios (within our q > 0.5 analysis domain). We observe a raw binary fraction for solar-type main-sequence stars with mass ratios q > 0.5 of 42% ± 4%, independent of the assumed mass-ratio distribution to within its uncertainties, consistent with literature values for old open clusters but significantly higher than the field solar-type binary fraction. We confirm that the binaries identified by our method are more concentrated than single stars, in agreement with previous studies, and we demonstrate that the binary nature of those candidates that remain unidentified spectroscopically is strongly supported by photometry from Gaia DR2.

79 ASTRONOMY AND ASTROPHYSICS↗

KMT-2016-BLG-1836Lb: A Super-Jovian Planet from a High-cadence Microlensing Field

We report the discovery of a super-Jovian planet in the microlensing event KMT-2016-BLG-1836, which was found by the Korea Microlensing Telescope Network (KMTNet) high-cadence observations (Γ~4 hr -1 ). The planet–host mass ratio q ~ 0.004. A Bayesian analysis indicates that the planetary system is composed of a super-Jovian M{ planet =2.2 -1.1 +1.9 M J planet orbiting an M or K dwarf, M{ host =0.49 -0.25 0.38 M ⊙ , at a distance of D L =7.1 -2.4} +0.8 kpc. The projected planet–host separation is 3.5 -0.9 +1.1 au, implying that the planet is located beyond the snow line of the host star. Future high-resolution images can potentially strongly constrain the lens brightness and thus the mass and distance of the planetary system. Without considering detailed detection efficiency, selection, or publication biases, we find a potential mass-ratio desert at -3.7 ≲ log q ≲ -3.0 for the 31 published KMTNet planets.

79 ASTRONOMY AND ASTROPHYSICS↗

OGLE-2018-BLG-1269Lb: A Jovian Planet with a Bright I = 16 Host

We report the discovery of a planet in the microlensing event OGLE-2018-BLG-1269 with a planet–host mass ratio q ∼ 6 × 10{sup −4}, i.e., 0.6 times smaller than the Jupiter/Sun mass ratio. Combined with the Gaia parallax and proper motion, a strong one-dimensional constraint on the microlens parallax vector allows us to significantly reduce the uncertainties of lens physical parameters. A Bayesian analysis that ignores any information about light from the host yields that the planet is a cold giant (M{sub 2}=0.69{sub −0.22}{sup +0.44} M{sub J}) orbiting a Sun-like star (M{sub 1}=1.13{sub −0.35}{sup +0.72} M{sub ⊙}) at a distance of D{sub L}=2.56{sub −0.62}{sup +0.92} kpc. The projected planet–host separation is a{sub ⊥}=4.61{sub −1.17}{sup +1.70} au. Using Gaia astrometry, we show that the blended light lies ≲12 mas from the host and therefore must be either the host star or a stellar companion to the host. An isochrone analysis favors the former possibility at >99.6%. The host is therefore a subgiant. For host metallicities in the range of 0.0⩽[Fe/H]⩽+0.3, the host and planet masses are then in the range of 1.16⩽M{sub 1}/M{sub ⊙}⩽1.38 and 0.74⩽M{sub 2}/M{sub J}⩽0.89, respectively. Low host metallicities are excluded. The brightness and proximity of the lens make the event a strong candidate for spectroscopic follow-up both to test the microlensing solution and to further characterize the system.

79 ASTRONOMY AND ASTROPHYSICS↗

New Determination of Fundamental Properties of Palomar 5 Using Deep DESI Imaging Data

The legacy imaging surveys for the Dark Energy Spectroscopic Instrument project provides multiple-color photometric data, which are about 2 mag deeper than those from the SDSS. In this study, we redetermine the fundamental properties for an old halo globular cluster of Palomar 5 based on these new imaging data, including structure parameters, stellar population parameters, and luminosity and mass functions. These characteristics, together with its tidal tails, are key for dynamical studies of the cluster and constraining the mass model of the Milky Way. By fitting the King model to the radial surface density profile of Palomar 5, we derive the core radius of [[Formula]], tidal radius of [[Formula]], and concentration parameter of c = 0.78 ± 0.04. We apply a Bayesian analysis method to derive the stellar population properties and get an age of 11.508 ± 0.027 Gyr, metallicity of [Fe/H] = −1.798 ± 0.014, reddening of E(B − V) = 0.0552 ± 0.0005, and distance modulus of (m−M){sub 0} = 16.835±0.006. The main-sequence luminosity and mass functions for both the cluster center and tidal tails are investigated. The luminosity and mass functions at different distances from the cluster center suggest that there is obvious spatial mass segregation. Many faint low-mass stars have been evaporated at the cluster center, and the tidal tails are enhanced by low-mass stars. Both the concentration and relaxation times suggest that Palomar 5 is a totally relaxed system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗