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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Revealing the collision energy dependence of η/s in RHIC-BES Au+Au collisions using Bayesian statistics

Here, we investigate the collision energy dependence of η/s in a transport + viscous hydrodynamics hybrid model. A Bayesian analysis is performed on RHIC beam energy scan data for Au + Au collisions at √ s NN = 19.6, 39, and 62.4 GeV. The resulting posterior probability distributions for the model parameters show a preference for a larger value of η/s at 19.6 GeV compared to 62.4 GeV, indicating dependence on baryon chemical potential μ B .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Hot QCD phase diagram from holographic Einstein–Maxwell–Dilaton models

In this review, we provide an up-to-date account of quantitative holographic descriptions of the strongly coupled quark-gluon plasma (QGP) produced in heavy-ion collisions, based on the class of gauge-gravity Einstein-Maxwell-Dilaton (EMD) models. Holography is employed to tentatively map the QCD phase diagram at finite temperature onto a dual theory of charged, asymptotically AdS black holes in 5D. With a quantitative focus on the hot QCD phase diagram, the EMD models reviewed are adjusted to describe lattice results for the finite-temperature QCD equation of state, with 2+1 flavors and physical quark masses, at zero chemical potential and vanishing electromagnetic fields. The predictive power of EMD models is tested by quantitatively comparing their predictions for the hot QCD equation of state at nonzero baryon density and the corresponding state-of-the-art lattice QCD results. The shear and bulk viscosities predicted by these EMD models are also compared to the corresponding profiles favored by the latest phenomenological multistage models describing different heavy-ion data. We report preliminary results from a Bayesian analysis which provide systematic evidence that lattice results at finite temperature and zero baryon density strongly constrains the free parameters of EMD models. Remarkably, the set of parameters constrained by lattice results at zero chemical potential produces EMD models in quantitative agreement with lattice QCD results also at finite baryon density. We also review results for equilibrium and transport properties from magnetic EMD models, describing the QGP at finite temperatures and magnetic fields. Lastly, we provide a critical assessment of the main limitations and drawbacks of the holographic models reviewed in the present work, and point out some perspectives we believe are of fundamental importance for future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

Evidence of recombination and positive selection in cetacean papillomaviruses

Papillomaviruses (PVs) are small DNA viruses that have been associated with increased epithelial proliferation. Over one hundred PV types have been identified in humans; however, only three have been identified in bottlenose dolphins (Tursiops truncatus) to date. Using rolling circle amplification and degenerate PCR, we identified four novel PV genomes of bottlenose dolphins. TtPV4, TtPV5 and TtPV6 were identified in genital lesions while TtPV7 was identified in normal genital mucosa. Bayesian analysis of the full-length L1 genes found that TtPV4 and TtPV7 group within the Upsilonpapillomavirus genus while TtPV5 and TtPV6 group with Omikronpapillomavirus. However, analysis of the E1 gene did not distinguish these genera, implying that these genes may not share a common history, consistent with recombination. Recombination analyses identified several probable events. Signals of positive selection were found mostly in the E1 and E2 genes. Recombination and diversifying selection pressures constitute important driving forces of cetacean PV evolution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Magnetostriction of α-RuCl 3 Flakes in the Zigzag Phase

Motivated by the possibility of enhanced magnetic fluctuations in exfoliated α-RuCl 3 flakes, we study magneto-Raman spectra of exfoliated multilayer α-RuCl 3 in out-of-plane magnetic fields of -6 to 6 T at temperatures of 670 mK to 4 K. While the literature currently suggests that bulk α-RuCl 3 is in an antiferromagnetic zigzag phase with R$\bar3$ symmetry at low temperatures, we do not observe R$\bar3$ symmetry in exfoliated α-RuCl 3 at low temperatures. While we saw no magnetic field-driven transitions, the Raman modes exhibit unexpected stochastic shifts in response to the applied magnetic field that are above the uncertainties inferred from Bayesian analysis. Finally, these stochastic shifts are consistent with the emergence of magnetostrictive interactions in exfoliated α-RuCl 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The damping of terahertz acoustic modes in aqueous nanoparticle suspensions

In this work, we investigate the possibility of controlling the acoustic damping in a liquid when nanoparticles are suspended in it. To shed light on this topic, we performed Inelastic X-Ray Scattering (IXS) measurements of the terahertz collective dynamics of aqueous suspensions of nanospheres of various materials, size, and relative concentration, either charged or neutral. A Bayesian analysis of measured spectra indicates that the damping of the two acoustic modes of water increases upon nanoparticle immersion. This effect seems particularly pronounced for the longitudinal acoustic mode, which, whenever visible at all, rapidly damps off when increasing the exchanged wavevector. Results also indicate that the observed effect strongly depends on the material the immersed nanoparticles are made of.

36 MATERIALS SCIENCE↗
DOE OSTI2023

BEYONDPLANCK II. CMB mapmaking through Gibbs sampling

We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the BEYONDPLANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A spectrometer for high-precision ion temperature and velocity measurements in low-temperature plasmas

We have developed a low-cost spectrometer with simple optical design that achieves unprecedented precision measurements of ion temperature (±0.01 eV) and velocity (±20 m/s). A Fabry-Pérot étalon provides the simultaneous high resolving power and high throughput needed for the light levels available in singly ionized helium and argon plasmas. Reducing the systematic uncertainty in the absolute wavelength calibration needed for the specified velocity precision motivates a Bayesian analysis method called Nested Sampling to address the nontrivial uncertainty in the diffraction order. An initial emission measurement of a singly charged stationary argon plasma yields a temperature of 0.339 ± 0.007 eV and a velocity of -3 ± 4 m/s with a systematic uncertainty of 20 m/s.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Demonstration of improved laser preheat with a cryogenically cooled magnetized liner inertial fusion platform

We report on progress implementing and testing cryogenically cooled platforms for Magnetized Liner Inertial Fusion (MagLIF) experiments. Two cryogenically cooled experimental platforms were developed: an integrated platform fielded on the Z pulsed power generator that combines magnetization, laser preheat, and pulsed-power-driven fuel compression and a laser-only platform in a separate chamber that enables measurements of the laser preheat energy using shadowgraphy measurements. The laser-only experiments suggest that ~89% ± 10% of the incident energy is coupled to the fuel in cooled targets across the energy range tested, significantly higher than previous warm experiments that achieved at most 67% coupling and in line with simulation predictions. The laser preheat configuration was applied to a cryogenically cooled integrated experiment that used a novel cryostat configuration that cooled the MagLIF liner from both ends. Furthermore, the integrated experiment, z3576, coupled 2.32 ± 0.25 kJ preheat energy to the fuel, the highest to-date, demonstrated excellent temperature control and nominal current delivery, and produced one of the highest pressure stagnations as determined by a Bayesian analysis of the data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2011

SPATIALLY RESOLVING THE HK Tau B EDGE-ON DISK FROM 1.2 TO 4.7 {mu}m: A UNIQUE SCATTERED LIGHT DISK

We present spatially resolved scattered light images of the circumstellar disk around HK Tau B at 3.8 and 4.7 {mu}m taken with the Keck Telescope Laser Guide Star Adaptive Optics (AO) system, and 1.6-2.12 {mu}m images taken with the Very Large Telescope/NACO AO system. Combined with previously published optical Hubble Space Telescope data, we investigate the spatially resolved scattered light properties of this edge-on circumstellar disk and probe for the presence of large grains. The 0.6-3.8 {mu}m scattered light observations reveal strong, and in some cases, unusual, wavelength dependencies in the observed disk morphology. The separation between the two scattered light nebulae, which is directly proportional to the disk-mass-opacity product, decreases by 30% between 0.6 and 3.8 {mu}m. Over the same wavelength range, the FWHM of the disk nebulosity declines by a factor of two, while the flux ratio between the two nebulae increases by a factor of {approx}8. No other disk known to date shows a flux ratio that increases with wavelength. Both the FWHM and nebula flux ratio are affected by the scattering phase function and the observed behavior can most readily be explained by a phase function that becomes more forward throwing with wavelength. The multi-wavelength scattered light observations also confirm the asymmetric nature of the disk and show that the level of asymmetry is a function of wavelength. We use the MCFOST radiative transfer code to model the disk at four wavelengths, corresponding to the I, H, Ks, and L' bandpasses. A single power-law grain size distribution can recreate the observed disk properties simultaneously at all four wavelengths. Bayesian analysis of the dust parameters finds a 99% probability that the maximum grain size is 5.5 {mu}m or larger. We also find that the grain size distribution is steep, with a 99% probability of a power-law index of 4.2 or larger, suggesting that these large grains are a small fraction of the overall dust population. The best-fit dust asymmetry parameter for each individual wavelength shows an unusual behavior, increasing with wavelength from the optical through the near-infrared, peaking at {approx}0.8 between 2.2 and 3.8 {mu}m, then decreasing by a factor of two by {approx}12 {mu}m. Comparing the wavelength dependence of the asymmetry parameter for HK Tau B with those for the interstellar medium (ISM) and dark cloud dust models, we find considerable evolution from an ISM state and argue for the presence of grain growth within the disk. Further, comparing the wavelength dependence of the asymmetry parameter for GG Tau, HV Tau C, and HK Tau B, the three disks that have been spatially resolved in scattered light between 0.8 and 3.8 {mu}m, finds a diverse range of dust properties, indicating differing degrees of grain growth for disks at a similar age.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

THE PHYSICAL CONDITIONS IN STARBURSTS DERIVED FROM BAYESIAN FITTING OF MID-INFRARED SPECTRAL ENERGY DISTRIBUTION MODELS: 30 DORADUS AS A TEMPLATE

To understand and interpret the observed spectral energy distributions (SEDs) of starbursts, theoretical or semi-empirical SED models are necessary. Yet, while they are well founded in theory, independent verification and calibration of these models, including the exploration of possible degeneracies between their parameters, are rarely made. As a consequence, a robust fitting method that leads to unique and reproducible results has been lacking. Here we introduce a novel approach based on Bayesian analysis to fit the Spitzer-Infrared Spectrometer spectra of starbursts using the SED models proposed by Groves et al.. We demonstrate its capabilities and verify the agreement between the derived best-fit parameters and actual physical conditions by modeling the nearby, well-studied, giant H II region 30 Doradus in the LMC. The derived physical parameters, such as cluster mass, cluster age, interstellar medium pressure, and covering fraction of photodissociation regions, are representative of the 30 Doradus region. The inclusion of the emission lines in the modeling is crucial to break degeneracies. We investigate the limitations and uncertainties by modeling subregions, which are dominated by single components, within 30 Doradus. A remarkable result for 30 Doradus in particular is a considerable contribution to its mid-infrared spectrum from hot ({approx}300 K) dust. The demonstrated success of our approach will allow us to derive the physical conditions in more distant, spatially unresolved starbursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

YOUNG RADIO PULSARS IN GALACTIC GLOBULAR CLUSTERS

Currently three isolated radio pulsars and one binary radio pulsar with no evidence of any previous recycling are known in 97 surveyed Galactic globular clusters (GCs). As pointed out by Lyne et al., the presence of these pulsars cannot be explained by core-collapse supernovae, as commonly assumed for their counterparts in the Galactic disk. We apply a Bayesian analysis to the results from surveys for radio pulsars in GCs and find the number of potentially observable non-recycled radio pulsars present in all clusters to be <3600. Accounting for beaming and retention considerations, the implied birthrate for any formation scenario for all 97 clusters is <0.25 pulsars century{sup -1} assuming a Maxwellian distribution of velocities with a dispersion of 10 km s{sup -1}. The implied birthrates for higher velocity dispersions are substantially higher than inferred for such pulsars in the Galactic disk. This suggests that the velocity dispersion of young pulsars in GCs is significantly lower than those of disk pulsars. These numbers may be substantial overestimates due to the fact that the currently known sample of young pulsars is observed only in metal-rich clusters. We propose that young pulsars may only be formed in GCs with metallicities with log[Fe/H] > - 0.6. In this case, the potentially observable population of such young pulsars is 447{sup +1420}{sub -399} (the error bars give a 95% confidence interval) and their birthrate is 0.012{sup +0.037}{sub -0.010} pulsars century{sup -1}. The most likely creation scenario to explain these pulsars is the electron capture supernova of an OMgNe white dwarf.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A POSSIBLE BINARY SYSTEM OF A STELLAR REMNANT IN THE HIGH-MAGNIFICATION GRAVITATIONAL MICROLENSING EVENT OGLE-2007-BLG-514

We report the extremely high-magnification (A > 1000) binary microlensing event OGLE-2007-BLG-514. We obtained good coverage around the double peak structure in the light curve via follow-up observations from different observatories. The binary lens model that includes the effects of parallax (known orbital motion of the Earth) and orbital motion of the lens yields a binary lens mass ratio of q = 0.321 {+-} 0.007 and a projected separation of s = 0.072 {+-} 0.001 in units of the Einstein radius. The parallax parameters allow us to determine the lens distance D{sub L} = 3.11 {+-} 0.39 kpc and total mass M{sub L} = 1.40 {+-} 0.18 M{sub Sun }; this leads to the primary and secondary components having masses of M{sub 1} = 1.06 {+-} 0.13 M{sub Sun} and M{sub 2} = 0.34 {+-} 0.04 M{sub Sun }, respectively. The parallax model indicates that the binary lens system is likely constructed by the main-sequence stars. On the other hand, we used a Bayesian analysis to estimate probability distributions by the model that includes the effects of xallarap (possible orbital motion of the source around a companion) and parallax (q = 0.270 {+-} 0.005, s = 0.083 {+-} 0.001). The primary component of the binary lens is relatively massive, with M{sub 1} = 0.9{sup +4.6}{sub -0.3} M{sub Sun} and it is at a distance of D{sub L} = 2.6{sup +3.8}{sub -0.9} kpc. Given the secure mass ratio measurement, the companion mass is therefore M{sub 2} = 0.2{sup +1.2}{sub -0.1} M{sub Sun }. The xallarap model implies that the primary lens is likely a stellar remnant, such as a white dwarf, a neutron star, or a black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

ATMOSPHERIC RETRIEVAL FOR SUPER-EARTHS: UNIQUELY CONSTRAINING THE ATMOSPHERIC COMPOSITION WITH TRANSMISSION SPECTROSCOPY

We present a retrieval method based on Bayesian analysis to infer the atmospheric compositions and surface or cloud-top pressures from transmission spectra of exoplanets with general compositions. In this study, we identify what can unambiguously be determined about the atmospheres of exoplanets from their transmission spectra by applying the retrieval method to synthetic observations of the super-Earth GJ 1214b. Our approach to inferring constraints on atmospheric parameters is to compute their joint and marginal posterior probability distributions using the Markov Chain Monte Carlo technique in a parallel tempering scheme. A new atmospheric parameterization is introduced that is applicable to general atmospheres in which the main constituent is not known a priori and clouds may be present. Our main finding is that a unique constraint of the mixing ratios of the absorbers and two spectrally inactive gases (such as N{sub 2} and primordial H{sub 2}+ He) is possible if the observations are sufficient to quantify both (1) the broadband transit depths in at least one absorption feature for each absorber and (2) the slope and strength of the molecular Rayleigh scattering signature. A second finding is that the surface pressure or cloud-top pressure can be quantified if a surface or cloud deck is present at low optical depth. A third finding is that the mean molecular mass can be constrained by measuring either the Rayleigh scattering slope or the shapes of the absorption features, thus enabling one to distinguish between cloudy hydrogen-rich atmospheres and high mean molecular mass atmospheres. We conclude, however, that without the signature of molecular Rayleigh scattering-even with robustly detected infrared absorption features (>10{sigma})-there is no reliable way to tell from the transmission spectrum whether the absorber is a main constituent of the atmosphere or just a minor species with a mixing ratio of X{sub abs} < 0.1%. The retrieval method leads us to a conceptual picture of which details in transmission spectra are essential for unique characterizations of well-mixed exoplanet atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

AN INTERFEROMETRIC STUDY OF THE FOMALHAUT INNER DEBRIS DISK. II. KECK NULLER MID-INFRARED OBSERVATIONS

We report on high-contrast mid-infrared observations of Fomalhaut obtained with the Keck Interferometer Nuller (KIN) showing a small resolved excess over the level expected from the stellar photosphere. The measured null excess has a mean value of 0.35% {+-} 0.10% between 8 and 11 {mu}m and increases from 8 to 13 {mu}m. Given the small field of view of the instrument, the source of this marginal excess must be contained within 2 AU of Fomalhaut. This result is reminiscent of previous VLTI K-band ({approx_equal}2{mu}m) observations, which implied the presence of a {approx}0.88% excess, and argued that thermal emission from hot dusty grains located within 6 AU from Fomalhaut was the most plausible explanation. Using a parametric two-dimensional radiative transfer code and a Bayesian analysis, we examine different dust disk structures to reproduce both the near- and mid-infrared data simultaneously. While not a definitive explanation of the hot excess of Fomalhaut, our model suggests that the most likely inner few AU disk geometry consists of a two-component structure, with two different and spatially distinct grain populations. The 2-11 {mu}m data are consistent with an inner hot ring of very small ({approx_equal}10-300 nm) carbon-rich grains concentrating around 0.1 AU. The second dust population-inferred from the KIN data at longer mid-infrared wavelengths-consists of larger grains (size of a few microns to a few tens of microns) located further out in a colder region where regular astronomical silicates could survive, with an inner edge around 0.4 AU-1 AU. From a dynamical point of view, the presence of the inner concentration of submicron-sized grains is surprising, as such grains should be expelled from the inner planetary system by radiation pressure within only a few years. This could either point to some inordinate replenishment rates (e.g., many grazing comets coming from an outer reservoir) or to the existence of some braking mechanism preventing the grains from moving out.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

REVISITING THE SCALING RELATIONS OF BLACK HOLE MASSES AND HOST GALAXY PROPERTIES

New kinematic data and modeling efforts in the past few years have substantially expanded and revised dynamical measurements of black hole masses (M {sub .}) at the centers of nearby galaxies. Here we compile an updated sample of 72 black holes and their host galaxies, and present revised scaling relations between M {sub .} and stellar velocity dispersion ({sigma}), V-band luminosity (L), and bulge stellar mass (M {sub bulge}), for different galaxy subsamples. Our best-fitting power-law relations for the full galaxy sample are log{sub 10}(M {sub .}) = 8.32 + 5.64log{sub 10}({sigma}/200 km s{sup -1}), log{sub 10}(M {sub .}) = 9.23 + 1.11log{sub 10}(L/10{sup 11} L {sub Sun }), and log{sub 10}(M {sub .}) = 8.46 + 1.05log{sub 10}(M {sub bulge}/10{sup 11} M {sub Sun }). A log-quadratic fit to the M {sub .}-{sigma} relation with an additional term of {beta}{sub 2} [log{sub 10}({sigma}/200 km s{sup -1})]{sup 2} gives {beta}{sub 2} = 1.68 {+-} 1.82 and does not decrease the intrinsic scatter in M {sub .}. Including 92 additional upper limits on M {sub .} does not change the slope of the M {sub .}-{sigma} relation. When the early- and late-type galaxies are fit separately, we obtain similar slopes of 5.20 and 5.06 for the M {sub .}-{sigma} relation but significantly different intercepts-M {sub .} in early-type galaxies are about two times higher than in late types at a given sigma. Within early-type galaxies, our fits to M {sub .}({sigma}) give M {sub .} that is about two times higher in galaxies with central core profiles than those with central power-law profiles. Our M {sub .}-L and M {sub .}-M {sub bulge} relations for early-type galaxies are similar to those from earlier compilations, and core and power-law galaxies yield similar L- and M {sub bulge}-based predictions for M {sub .}. When the conventional quadrature method is used to determine the intrinsic scatter in M {sub .}, our data set shows weak evidence for increased scatter at M {sub bulge} < 10{sup 11} M {sub Sun} or L{sub V} < 10{sup 10.3} L {sub Sun }, while the scatter stays constant for 10{sup 11} < M {sub bulge} < 10{sup 12.3} M {sub Sun} and 10{sup 10.3} < L{sub V} < 10{sup 11.5} L {sub Sun }. A Bayesian analysis indicates that a larger sample of M {sub .} measurements would be needed to detect any statistically significant trend in the scatter with galaxy properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

A HOT URANUS ORBITING THE SUPER METAL-RICH STAR HD 77338 AND THE METALLICITY-MASS CONNECTION

We announce the discovery of a low-mass planet orbiting the super metal-rich K0V star HD 77338 as part of our ongoing Calan-Hertfordshire Extrasolar Planet Search. The best-fit planet solution has an orbital period of 5.7361 {+-} 0.0015 days and with a radial velocity semi-amplitude of only 5.96 {+-} 1.74 ms{sup -1}, we find a minimum mass of 15.9{sup +4.7}{sub -5.3} M{sub Circled-Plus }. The best-fit eccentricity from this solution is 0.09{sup +0.25}{sub -0.09}, and we find agreement for this data set using a Bayesian analysis and a periodogram analysis. We measure a metallicity for the star of +0.35 {+-} 0.06 dex, whereas another recent work finds +0.47 {+-} 0.05 dex. Thus HD 77338b is one of the most metal-rich planet-host stars known and the most metal-rich star hosting a sub-Neptune-mass planet. We searched for a transit signature of HD 77338b but none was detected. We also highlight an emerging trend where metallicity and mass seem to correlate at very low masses, a discovery that would be in agreement with the core accretion model of planet formation. The trend appears to show that for Neptune-mass planets and below, higher masses are preferred when the host star is more metal-rich. Also a lower boundary is apparent in the super metal-rich regime where there are no very low mass planets yet discovered in comparison to the sub-solar metallicity regime. A Monte Carlo analysis shows that this low-mass planet desert is statistically significant with the current sample of 36 planets at the {approx}4.5{sigma} level. In addition, results from Kepler strengthen the claim for this paucity of the lowest-mass planets in super metal-rich systems. Finally, this discovery adds to the growing population of low-mass planets around low-mass and metal-rich stars and shows that very low mass planets can now be discovered with a relatively small number of data points using stable instrumentation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

STARS DO NOT EAT THEIR YOUNG MIGRATING PLANETS: EMPIRICAL CONSTRAINTS ON PLANET MIGRATION HALTING MECHANISMS

The discovery of ''hot Jupiters'' very close to their parent stars confirmed that Jovian planets migrate inward via several potential mechanisms. We present empirical constraints on planet migration halting mechanisms. We compute model density functions of close-in exoplanets in the orbital semi-major axis-stellar mass plane to represent planet migration that is halted via several mechanisms, including the interior 1:2 resonance with the magnetospheric disk truncation radius, the interior 1:2 resonance with the dust sublimation radius, and several scenarios for tidal halting. The models differ in the predicted power-law dependence of the exoplanet orbital semi-major axis as a function of stellar mass, and thus we also include a power-law model with the exponent as a free parameter. We use a Bayesian analysis to assess the model success in reproducing empirical distributions of confirmed exoplanets and Kepler candidates that orbit interior to 0.1 AU. Our results confirm a correlation of the halting distance with stellar mass. Tidal halting provides the best fit to the empirical distribution of confirmed Jovian exoplanets at a statistically robust level, consistent with the Kozai mechanism and the spin-orbit misalignment of a substantial fraction of hot Jupiters. We can rule out migration halting at the interior 1:2 resonances with the magnetospheric disk truncation radius and the interior 1:2 resonance with the dust disk sublimation radius, a uniform random distribution, and a distribution with no dependence on stellar mass. Note that our results do not rule out Type-II migration, but rather eliminate the role of a circumstellar disk in stopping exoplanet migration. For Kepler candidates, which have a more restricted range in stellar mass compared to confirmed planets, we are unable to discern between the tidal dissipation and magnetospheric disk truncation braking mechanisms at a statistically significant level. The power-law model favors exponents in the range of 0.38-0.9. This is larger than that predicted for tidal halting (0.23-0.33), which suggests that additional physics may be missing in the tidal halting theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

CANDIDATE PLANETS IN THE HABITABLE ZONES OF KEPLER STARS

A key goal of the Kepler mission is the discovery of Earth-size transiting planets in ''habitable zones'' where stellar irradiance maintains a temperate climate on an Earth-like planet. Robust estimates of planet radius and irradiance require accurate stellar parameters, but most Kepler systems are faint, making spectroscopy difficult and prioritization of targets desirable. The parameters of 2035 host stars were estimated by Bayesian analysis and the probabilities p{sub HZ} that 2738 candidate or confirmed planets orbit in the habitable zone were calculated. Dartmouth Stellar Evolution Program models were compared to photometry from the Kepler Input Catalog, priors for stellar mass, age, metallicity and distance, and planet transit duration. The analysis yielded probability density functions for calculating confidence intervals of planet radius and stellar irradiance, as well as p{sub HZ}. Sixty-two planets have p{sub HZ} > 0.5 and a most probable stellar irradiance within habitable zone limits. Fourteen of these have radii less than twice the Earth; the objects most resembling Earth in terms of radius and irradiance are KOIs 2626.01 and 3010.01, which orbit late K/M-type dwarf stars. The fraction of Kepler dwarf stars with Earth-size planets in the habitable zone ({eta}{sub Circled-Plus }) is 0.46, with a 95% confidence interval of 0.31-0.64. Parallaxes from the Gaia mission will reduce uncertainties by more than a factor of five and permit definitive assignments of transiting planets to the habitable zones of Kepler stars.

79 ASTRONOMY AND ASTROPHYSICS↗