BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Adam Burrows

Adam Burrows appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-phastro-ph.IMastro-ph.HEhep-phgr-qcnucl-thastro-ph.COastro-ph.GAphysics.hist-ph

Catalog footprint

What is connected

102works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

102 published item(s)

preprint2022arXiv

The Collapse and Three-Dimensional Explosion of Three-Dimensional, vis à vis One-Dimensional, Massive-star Supernova Progenitor Models

The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically-symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M$_{\odot}$ model, evolved in three-dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code F{\sc{ornax}}. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically-symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multi-dimensional initial models serve as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Physical Effects of Progenitor Rotation: Comparing Two Long-Duration 3D Core-Collapse Supernova Simulations

We analyse and determine the effects of modest progenitor rotation in the context of core-collapse supernovae by comparing two separate long-duration three-dimensional simulations of 9 M$_{\odot}$ progenitors, one rotating with an initial spin period of $\sim$60 seconds and the other non-rotating. We determine that both models explode early, though the rotating model explodes a bit earlier. Despite this difference, the asymptotic explosion energies ($\sim$10$^{50}$ ergs) and residual neutron star baryon masses ($\sim$1.3 M$_{\odot}$) are similar. We find that the proto-neutron star (PNS) core can deleptonize and cool significantly more quickly. Soon into the evolution of the rotating model, we witness more vigorous and extended PNS core convection that early in its evolution envelopes the entire inner sphere, not just a shell. Moreover, we see a corresponding excursion in both the $ν_e$ luminosity and gravitational-wave strain that may be diagnostic of this observed dramatic phenomenon. In addition, after bounce the innermost region of the rotating model seems to execute meridional circulation. The rotationally-induced growth of the convective PNS region may facilitate the growth of core B-fields by the dynamo mechanism by facilitating the achievement of the critical Rossby number condition for substantial growth of a dipole field, obviating the need for rapid rotation rates to create dipole fields of significance. The next step is to explore the progenitor-mass and spin dependencies across the progenitor continuum of the supernova explosion, dynamics, and evolution of PNS convection and its potential role in the generation of magnetar and pulsar magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Core-Collapse Supernova Explosion Theory

Most supernova explosions accompany the death of a massive star. These explosions give birth to neutron stars and black holes and eject solar masses of heavy elements. However, determining the mechanism of explosion has been a half-century journey of great complexity. In this paper, we present our perspective of the status of this theoretical quest and the physics and astrophysics upon which its resolution seems to depend. The delayed neutrino-heating mechanism is emerging as a robust solution, but there remain many issues to address, not the least of which involves the chaos of the dynamics, before victory can unambiguously be declared. It is impossible to review in detail all aspects of this multi-faceted, more-than-half-century-long theoretical quest. Rather, we here map out the major ingredients of explosion and the emerging systematics of the observables with progenitor mass, as we currently see them. Our discussion will of necessity be speculative in parts, and many of the ideas may not survive future scrutiny. Some statements may be viewed as informed predictions concerning the numerous observables that rightly exercise astronomers witnessing and diagnosing the supernova Universe. Importantly, the same explosion in the inside, by the same mechanism, can look very different in photons, depending upon the mass and radius of the star upon explosion. A 10$^{51}$-erg (one "Bethe") explosion of a red supergiant with a massive hydrogen-rich envelope, a diminished hydrogen envelope, no hydrogen envelope, and, perhaps, no hydrogen envelope or helium shell all look very different, yet might have the same core and explosion evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the Origin of Pulsar and Magnetar Magnetic Fields

In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue that these expectations should apply to the convective proto-neutron stars at the centers of core-collapse supernovae. We analyze a suite of three-dimensional simulations of core collapse, featuring a realistic equation of state and full neutrino transport, in this context. All our progenitor models, ranging from 9 solar masses to 25 solar masses, including one with initial rotation, have sufficiently vigorous proto-neutron-star convection to generate dipole fields of order ~10^15 gauss, if the modified Rossby number resides in the critical range. Thus, the magnetar/radio pulsar dichotomy may arise naturally in part from the distribution of core rotation rates in massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Thermal Phase Curves of XO-3b: an Eccentric Hot Jupiter at the Deuterium Burning Limit

We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $μ$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $μ$m and $1610 \pm 70$ ppm at 4.5 $μ$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3's parallax measurement and find a planetary mass $M_p=11.79 \pm 0.98 ~M_{\rm{Jup}}$ and radius $R_p=1.295 \pm 0.066 ~R_{\rm{Jup}}$. We compare our \textit{Spitzer} observations with multiple atmospheric models to constrain the radiative and advective properties of XO-3b. While the decorrelated 4.5 $μ$m observations are pristine, the 3.6 $μ$m phase curve remains polluted with detector systematics due to larger amplitude intrapixel sensitivity variations in this channel. We focus our analysis on the more reliable 4.5 $μ$m phase curve and fit an energy balance model with solid body rotation to estimate the zonal wind speed and the pressure of the bottom of the mixed layer. Our energy balance model fit suggests an eastward equatorial wind speed of $3.13 ^{+0.26} _{-0.83}$ km/s, an atmospheric mixed layer down to $2.40 ^{+0.92} _{-0.16}$ bar, and Bond albedo of $0.106 ^{+0.008} _{-0.106}$. We assume that the wind speed and mixed layer depth are constant throughout the orbit. We compare our observations with a 1D planet-averaged model predictions at apoapse and periapse and 3D general circulation model (GCM) predictions for XO-3b. We also investigate the inflated radius of XO-3b and find that it would require an unusually large amount of internal heating to explain the observed planetary radius.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A systematic study of proto-neutron star convection in three-dimensional core-collapse supernova simulations

This paper presents the first systematic study of proto-neutron star (PNS) convection in three dimensions (3D) based on our latest numerical Fornax models of core-collapse supernova (CCSN). We confirm that PNS convection commonly occurs, and then quantify the basic physical characteristics of the convection. By virtue of the large number of long-term models, the diversity of PNS convective behavior emerges. We find that the vigor of PNS convection is not correlated with CCSN dynamics at large radii, but rather with the mass of PNS $-$ heavier masses are associated with stronger PNS convection. We find that PNS convection boosts the luminosities of $ν_μ$, $ν_τ$, $\barν_μ$, and $\barν_τ$ neutrinos, while the impact on other species is complex due to a competition of factors. Finally, we assess the consequent impact on CCSN dynamics and the potential for PNS convection to generate pulsar magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neutrino oscillations in supernovae: angular moments and fast instabilities

Recent theoretical work indicates that the neutrino radiation in core-collapse supernovae may be susceptible to flavor instabilities that set in far behind the shock, grow extremely rapidly, and have the potential to profoundly affect supernova dynamics and composition. Here we analyze the nonlinear collective oscillations that are prefigured by these instabilities. We demonstrate that a zero-crossing in $n_{ν_e} - n_{\barν_e}$ as a function of propagation angle is not sufficient to generate instability. Our analysis accounts for this fact and allows us to formulate complementary criteria. Using Fornax simulation data, we show that fast collective oscillations qualitatively depend on how forward-peaked the neutrino angular distributions are.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Profiles of James Peebles, Michel Mayor, and Didier Queloz: 2019 Nobel Laureates in Physics

The 2019 Nobel Prize in Physics honors three pioneering scientists for their fundamental contributions to basic cosmic questions - Professor James Peebles (Princeton University), Michel Mayor (University of Geneva), and Didier Queloz (University of Geneva and the University of Cambridge) - "for contributions to our understanding of the evolution of the universe and Earth's place in the cosmos," with one half to James Peebles "for theoretical discoveries in physical cosmology," and the other half jointly to Michel Mayor and Didier Queloz "for the discovery of an exoplanet orbiting a solar-type star." We summarize the historical and scientific backdrop to this year's Physics Nobel.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Statistical Characterization of Hot Jupiter Atmospheres using Spitzer's Secondary Eclipses

We report 78 secondary eclipse depths for a sample of 36 transiting hot Jupiters observed at 3.6- and 4.5 microns using the Spitzer Space Telescope. Our eclipse results for 27 of these planets are new, and include highly irradiated worlds such as KELT-7b, WASP-87b, WASP-76b, and WASP-64b, and important targets for JWST such as WASP-62b. We find that WASP-62b has a slightly eccentric orbit e cos(omega) = 0.00614+/- 0.00064, and we confirm the eccentricity of HAT-P-13b and WASP-14b. The remainder are individually consistent with circular orbits, but we find statistical evidence for eccentricity increasing with orbital period in our range from 1 to 5 days. Our day-side brightness temperatures for the planets yield information on albedo and heat redistribution, following Cowan and Agol (2011). Planets having maximum day side temperatures exceeding ~ 2200K are consistent with zero albedo and distribution of stellar irradiance uniformly over the day-side hemisphere. Our most intriguing result is that we detect a systematic difference between the emergent spectra of these hot Jupiters as compared to blackbodies. The ratio of observed brightness temperatures, Tb(4.5)/Tb(3.6), increases with equilibrium temperature by 100 +/- 24 parts-per-million per Kelvin, over the entire temperature range in our sample (800K to 2500K). No existing model predicts this trend over such a large range of temperature. We suggest that this may be due to a structural difference in the atmospheric temperature profile between the real planetary atmospheres as compared to models.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Prospects for Directly Imaging Young Giant Planets at Optical Wavelengths

In this work, we investigate the properties of young giant planet spectra in the optical and suggest that future space-based direct imaging missions should be considering young planets as a valuable and informative science case. While young planets are dimmer in the optical than in the infrared, they can still be brighter in the optical than a mature planet of similar mass. Therefore, an instrument designed to characterize mature planets should also be suitable for high-precision photometric imaging and spectroscopy of young self-luminous planets in a wavelength range and at a contrast ratio not attainable from the ground. We identify known young self-luminous companions which are feasible targets for WFIRST-CGI and compute spectra for them, including a treatment of scattering and reflected light at optical wavelengths. Using these results, we highlight potentially diagnostic spectral features that will be present in the WFIRST-CGI wavelengths. Expanding to direct imaging missions beyond WFIRST-CGI, we also use evolutionary models across a grid of masses and planet-star separations as inputs to compute spectra of hypothetical objects, exploring when reflected light may contribute to a degree comparable to that of thermal emission from the residual heat of formation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Overarching Framework of Core-Collapse Supernova Explosions as Revealed by 3D Fornax Simulations

We have conducted nineteen state-of-the-art 3D core-collapse supernova simulations spanning a broad range of progenitor masses. This is the largest collection of sophisticated 3D supernova simulations ever performed. We have found that while the majority of these models explode, not all do, and that even models in the middle of the available progenitor mass range may be less explodable. This does not mean that those models for which we did not witness explosion would not explode in Nature, but that they are less prone to explosion than others. One consequence is that the "compactness" measure is not a metric for explodability. We find that lower-mass massive star progenitors likely experience lower-energy explosions, while the higher-mass massive stars likely experience higher-energy explosions. Moreover, most 3D explosions have a dominant dipole morphology, have a pinched, wasp-waist structure, and experience simultaneous accretion and explosion. We reproduce the general range of residual neutron-star masses inferred for the galactic neutron-star population. The most massive progenitor models, however, in particular vis à vis explosion energy, need to be continued for longer physical times to asymptote to their final states. We find that while the majority of the inner ejecta have Y$_e = 0.5$, there is a substantial proton-rich tail. This result has important implications for the nucleosynthetic yields as a function of progenitor. Finally, we find that the non-exploding models eventually evolve into compact inner configurations that experience a quasi-periodic spiral SASI mode. We otherwise see little evidence of the SASI in the exploding models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

3.6 and 4.5 $μ$m ${\it Spitzer}$ Phase Curves of the Highly-Irradiated Hot Jupiters WASP-19b and HAT-P-7b

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 $μ$m obtained using the Spitzer Space Telescope. For WASP-19b, we measure secondary eclipse depths of $0.485\%\pm 0.024\%$ and $0.584\%\pm 0.029\%$ at 3.6 and 4.5 $μ$m, which are consistent with a single blackbody with effective temperature $2372 \pm 60$ K. The measured 3.6 and 4.5 $μ$m secondary eclipse depths for HAT-P-7b are $0.156\%\pm 0.009\%$ and $0.190\%\pm 0.006\%$, which are well-described by a single blackbody with effective temperature $2667\pm 57$ K. Comparing the phase curves to the predictions of one-dimensional and three-dimensional atmospheric models, we find that WASP-19b's dayside emission is consistent with a model atmosphere with no dayside thermal inversion and moderately efficient day-night circulation. We also detect an eastward-shifted hotspot, suggesting the presence of a superrotating equatorial jet. In contrast, HAT-P-7b's dayside emission suggests a dayside thermal inversion and relatively inefficient day-night circulation; no hotspot shift is detected. For both planets, these same models do not agree with the measured nightside emission. The discrepancies in the model-data comparisons for WASP-19b might be explained by high-altitude silicate clouds on the nightside and/or high atmospheric metallicity, while the very low 3.6 $μ$m nightside planetary brightness for HAT-P-7b may be indicative of an enhanced global C/O ratio. We compute Bond albedos of 0 ($<0.08$ at $1σ$) and $0.38\pm 0.06$ for WASP-19b and HAT-P-7b, respectively. In the context of other planets with thermal phase curve measurements, we show that WASP-19b and HAT-P-7b fit the general trend of decreasing day-night heat recirculation with increasing irradiation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detecting the Supernova Breakout Burst in Terrestrial Neutrino Detectors

We calculate the distance-dependent performance of a few representative terrestrial neutrino detectors in detecting and measuring the properties of the $ν_e$ breakout burst light curve in a Galactic core-collapse supernova. The breakout burst is a signature phenomenon of core collapse and offers a probe into the stellar core through collapse and bounce. We examine cases of no neutrino oscillations and oscillations due to normal and inverted neutrino-mass hierarchies. For the normal hierarchy, other neutrino flavors emitted by the supernova overwhelm the $ν_e$ signal, making a detection of the breakout burst difficult. For the inverted hierarchy (IH), some detectors at some distances should be able to see the $ν_e$ breakout burst peak and measure its properties. For the IH, the maximum luminosity of the breakout burst can be measured at 10 kpc to accuracies of $\sim$30% for Hyper-Kamiokande (Hyper-K) and $\sim$60% for the Deep Underground Neutrino Experiment (DUNE). Super-Kamiokande (Super-K) and Jiangmen Underground Neutrino Observatory (JUNO) lack the mass needed to make an accurate measurement. For the IH, the time of the maximum luminosity of the breakout burst can be measured in Hyper-K to an accuracy of $\sim$3 ms at 7 kpc, in DUNE $\sim$2 ms at 4 kpc, and JUNO and Super-K can measure the time of maximum luminosity to an accuracy of $\sim$2 ms at 1 kpc. Detector backgrounds in IceCube render a measurement of the $ν_e$ breakout burst unlikely. For the inverted hierarchy, a measurement of the maximum luminosity of the breakout burst could be used to differentiate between nuclear equations of state.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical Constraints on the Core Mass of Hot Jupiter HAT-P-13b

HAT-P-13b is a Jupiter-mass transiting exoplanet that has settled onto a stable, short-period, and mildly eccentric orbit as a consequence of the action of tidal dissipation and perturbations from a second, highly eccentric, outer companion. Due to the special orbital configuration of the HAT-P-13 system, the magnitude of HAT-P-13b's eccentricity ($e_b$) is in part dictated by its Love number ($k_{2_b}$), which is in turn a proxy for the degree of central mass concentration in its interior. Thus, the measurement of $e_b$ constrains $k_{2_b}$ and allows us to place otherwise elusive constraints on the mass of HAT-P-13b's core ($M_{\rm{core,b}}$). In this study we derive new constraints on the value of $e_b$ by observing two secondary eclipses of HAT-P-13b with the Infrared Array Camera on board the $\textit{Spitzer Space Telescope}$. We fit the measured secondary eclipse times simultaneously with radial velocity measurements and find that $e_b = 0.00700 \pm 0.00100$. We then use octupole-order secular perturbation theory to find the corresponding $k_{2_b} = 0.31^{+0.08}_{-0.05}$. Applying structural evolution models, we then find, with 68\% confidence, that $M_{\rm{core,b}}$ is less than 25 Earth masses ($M_{\oplus}$). The most likely value of $M_{\rm{core,b}} = 11 M_{\oplus}$, which is similar to the core mass theoretically required for runaway gas accretion. This is the tightest constraint to date on the core mass of a hot Jupiter. Additionally, we find that the measured secondary eclipse depths, which are in the 3.6 $μ$m and 4.5 $μ$m bands, best match atmospheric model predictions with a dayside temperature inversion and relatively efficient day-night circulation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

HAT-P-65b and HAT-P-66b: Two Transiting Inflated Hot Jupiters and Observational Evidence for the Re-Inflation of Close-In Giant Planets

We present the discovery of the transiting exoplanets HAT-P-65b and HAT-P-66b, with orbital periods of 2.6055 d and 2.9721 d, masses of $0.527 \pm 0.083$ M$_{J}$ and $0.783 \pm 0.057$ M$_{J}$ and inflated radii of $1.89 \pm 0.13$ R$_{J}$ and $1.59^{+0.16}_{-0.10}$ R$_{J}$, respectively. They orbit moderately bright ($V=13.145 \pm 0.029$, and $V=12.993 \pm 0.052$) stars of mass $1.212 \pm 0.050$ M$_{\odot}$ and $1.255^{+0.107}_{-0.054}$ M$_{\odot}$. The stars are at the main sequence turnoff. While it is well known that the radii of close-in giant planets are correlated with their equilibrium temperatures, whether or not the radii of planets increase in time as their hosts evolve and become more luminous is an open question. Looking at the broader sample of well-characterized close-in transiting giant planets, we find that there is a statistically significant correlation between planetary radii and the fractional ages of their host stars, with a false alarm probability of only 0.0041%. We find that the correlation between the radii of planets and the fractional ages of their hosts is fully explained by the known correlation between planetary radii and their present day equilibrium temperatures, however if the zero-age main sequence equilibrium temperature is used in place of the present day equilibrium temperature then a correlation with age must also be included to explain the planetary radii. This suggests that, after contracting during the pre-main-sequence, close-in giant planets are re-inflated over time due to the increasing level of irradiation received from their host stars. Prior theoretical work indicates that such a dynamic response to irradiation requires a significant fraction of the incident energy to be deposited deep within the planetary interiors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Transiting Exoplanet Studies and Community Targets for JWST's Early Release Science Program

The James Webb Space Telescope will revolutionize transiting exoplanet atmospheric science due to its capability for continuous, long-duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful future transiting exoplanet characterization programs. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed "community targets") that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations.(Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

3.6 and 4.5 $μ$m Phase Curves of the Highly-Irradiated Eccentric Hot Jupiter WASP-14b

We present full-orbit phase curve observations of the eccentric ($e\sim 0.08$) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 $μ$m bands using the \textit{Spitzer Space Telescope}. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of $0.1882\%\pm 0.0048\%$ and $0.2247\%\pm 0.0086\%$ at 3.6 and 4.5 $μ$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2σ$ upper limit on the nightside flux at 3.6 $μ$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $μ$m, the minimum planet flux is $30\%\pm 5\%$ of the maximum flux, corresponding to a brightness temperature of $1380\pm 65$ K. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b's measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models tend to over-predict the nightside emission at 3.6 $μ$m, while under-predicting the nightside emission at 4.5 $μ$m. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b ($7.8 M_{\mathrm{Jup}}$) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolving the HD 100546 Protoplanetary System with the Gemini Planet Imager: Evidence for Multiple Forming, Accreting Planets

We report Gemini Planet Imager H band high-contrast imaging/integral field spectroscopy and polarimetry of the HD 100546, a 10 $Myr$-old early-type star recently confirmed to host a thermal infrared bright (super)jovian protoplanet at wide separation, HD 100546 b. We resolve the inner disk cavity in polarized light, recover the thermal-infrared (IR) bright arm, and identify one additional spiral arm. We easily recover HD 100546 b and show that much of its emission originates an unresolved, point source. HD 100546 b likely has extremely red infrared colors compared to field brown dwarfs, qualitatively similar to young cloudy superjovian planets, however, these colors may instead indicate that HD 100546 b is still accreting material from a circumplanetary disk. Additionally, we identify a second point source-like peak at $r_{proj}$ $\sim$ 14 AU, located just interior to or at inner disk wall consistent with being a 10--20 $M_{J}$ candidate second protoplanet-- "HD 100546 c" -- and lying within a weakly polarized region of the disk but along an extension of the thermal IR bright spiral arm. Alternatively, it is equally plausible that this feature is a weakly polarized but locally bright region of the inner disk wall. Astrometric monitoring of this feature over the next 2 years and emission line measurements could confirm its status as a protoplanet, rotating disk hot spot that is possibly a signpost of a protoplanet, or a stationary emission source from within the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spitzer Secondary Eclipse Observations of Five Cool Gas Giant Planets and Empirical Trends in Cool Planet Emission Spectra

In this work we present Spitzer 3.6 and 4.5 micron secondary eclipse observations of five new cool (<1200 K) transiting gas giant planets: HAT-P-19b, WASP-6b, WASP-10b, WASP-39b, and WASP-67b. We compare our measured eclipse depths to the predictions of a suite of atmosphere models and to eclipse depths for planets with previously published observations in order to constrain the temperature- and mass-dependent properties of gas giant planet atmospheres. We find that the dayside emission spectra of planets less massive than Jupiter require models with efficient circulation of energy to the night side and/or increased albedos, while those with masses greater than that of Jupiter are consistently best-matched by models with inefficient circulation and low albedos. At these relatively low temperatures we expect the atmospheric methane to CO ratio to vary as a function of metallicity, and we therefore use our observations of these planets to constrain their atmospheric metallicities. We find that the most massive planets have dayside emission spectra that are best-matched by solar metallicity atmosphere models, but we are not able to place strong constraints on metallicities of the smaller planets in our sample. Interestingly, we find that the ratio of the 3.6 and 4.5 micron brightness temperatures for these cool transiting planets is independent of planet temperature, and instead exhibits a tentative correlation with planet mass. If this trend can be confirmed, it would suggest that the shape of these planets' emission spectra depends primarily on their masses, consistent with the hypothesis that lower-mass planets are more likely to have metal-rich atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spitzer Secondary Eclipses of the Dense, Modestly-irradiated, Giant Exoplanet HAT-P-20b Using Pixel-Level Decorrelation

HAT-P-20b is a giant exoplanet that orbits a metal-rich star. The planet itself has a high total density, suggesting that it may also have a high metallicity in its atmosphere. We analyze two eclipses of the planet in each of the 3.6- and 4.5 micron bands of Warm Spitzer. These data exhibit intra-pixel detector sensitivity fluctuations that were resistant to traditional decorrelation methods. We have developed a simple, powerful, and radically different method to correct the intra-pixel effect for Warm Spitzer data, which we call pixel-level decorrelation (PLD). PLD corrects the intra-pixel effect very effectively, but without explicitly using - or even measuring - the fluctuations in the apparent position of the stellar image. We illustrate and validate PLD using synthetic and real data, and comparing the results to previous analyses. PLD can significantly reduce or eliminate red noise in Spitzer secondary eclipse photometry, even for eclipses that have proven to be intractable using other methods. Our successful PLD analysis of four HAT-P-20b eclipses shows a best-fit blackbody temperature of 1134 +/-29K, indicating inefficient longitudinal transfer of heat, but lacking evidence for strong molecular absorption. We find sufficient evidence for variability in the 4.5 micron band that the eclipses should be monitored at that wavelength by Spitzer, and this planet should be a high priority for JWST spectroscopy. All four eclipses occur about 35 minutes after orbital phase 0.5, indicating a slightly eccentric orbit. A joint fit of the eclipse and transit times with extant RV data yields e(cos{omega}) = 0.01352 (+0.00054, -0.00057), and establishes the small eccentricity of the orbit to high statistical confidence. Given the existence of a bound stellar companion, HAT-P-20b is another excellent candidate for orbital evolution via Kozai migration or other three-body mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Direct Detectability of Giant Exoplanets in the Optical

Motivated by the possibility that a coronagraph will be put on WFIRST/AFTA, we explore the direct detectability of extrasolar giant planets (EGPs) in the optical. We quantify a planet's detectability by the fraction of its orbit for which it is in an observable configuration ($f_\mathrm{obs}$). Using a suite of Monte Carlo experiments, we study the dependence of $f_\mathrm{obs}$ upon the inner working angle (IWA) and minimum achievable contrast ($C_\mathrm{min}$) of the direct-imaging observatory; the planet's phase function, geometric albedo, single-scattering albedo, radius, and distance from Earth; and the semi-major axis distribution of EGPs. We calculate phase functions for a given geometric or single-scattering albedo, assuming various scattering mechanisms. We find that the Lambertian phase function can predict significantly larger $f_\mathrm{obs}$'s with respect to the more realistic Rayleigh phase function. For observations made with WFIRST/AFTA's baseline capabilities ($C_\mathrm{min}\sim10^{-9}$, $\mathrm{IWA}\sim0.2''$), Jupiter-like planets orbiting stars within 10, 30, and 50 parsecs of Earth have volume-averaged observability fractions of ${\sim}$12%, 3%, and 0.5%, respectively. At 10 parsecs, such observations yield $f_\mathrm{obs}>1\%$ for low- to modest-eccentricity planets with semi-major axes in the range ${\sim}2 - 10$ AU. If $C_\mathrm{min}=10^{-10}$, this range extends to ${\sim}35$ AU. We find that, in all but the most optimistic configurations, the probability for detection in a blind search is low (${<}\,5\%$). However, with orbital parameter constraints from long-term radial-velocity campaigns and Gaia astrometry, the tools we develop in this work can be used to determine both the most promising systems to target and when to observe them.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Two-Dimensional Core-Collapse Supernova Models with Multi-Dimensional Transport

We present new two-dimensional (2D) axisymmetric neutrino radiation/hydrodynamic models of core-collapse supernova (CCSN) cores. We use the CASTRO code, which incorporates truly multi-dimensional, multi-group, flux-limited diffusion (MGFLD) neutrino transport, including all relevant $\mathcal{O}(v/c)$ terms. Our main motivation for carrying out this study is to compare with recent 2D models produced by other groups who have obtained explosions for some progenitor stars and with recent 2D VULCAN results that did not incorporate $\mathcal{O}(v/c)$ terms. We follow the evolution of 12, 15, 20, and 25 solar-mass progenitors to approximately 600 milliseconds after bounce and do not obtain an explosion in any of these models. Though the reason for the qualitative disagreement among the groups engaged in CCSN modeling remains unclear, we speculate that the simplifying ``ray-by-ray' approach employed by all other groups may be compromising their results. We show that ``ray-by-ray' calculations greatly exaggerate the angular and temporal variations of the neutrino fluxes, which we argue are better captured by our multi-dimensional MGFLD approach. On the other hand, our 2D models also make approximations, making it difficult to draw definitive conclusions concerning the root of the differences between groups. We discuss some of the diagnostics often employed in the analyses of CCSN simulations and highlight the intimate relationship between the various explosion conditions that have been proposed. Finally, we explore the ingredients that may be missing in current calculations that may be important in reproducing the properties of the average CCSNe, should the delayed neutrino-heating mechanism be the correct mechanism of explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A First-Look Atmospheric Modeling Study of the Young Directly-Imaged Planet-Mass Companion, ROXs 42Bb

We present and analyze $J$$K_{s}$$L^\prime$ photometry and our previously published $H$-band photometry and $K$-band spectroscopy for ROXs 42Bb, an object Currie et al. (2014) first reported as a young directly imaged planet-mass companion. ROXs 42Bb exhibits IR colors redder than field L dwarfs but consistent with other planet-mass companions. From the H$_{2}$O-2 spectral index, we estimate a spectral type of L0 $\pm$ 1; weak detections/non-detections of the $CO$ bandheads, $Na I$, and $Ca I$ support evidence for a young, low surface gravity object primarily derived from the $H_{2}$(K) index. ROXs 42Bb's photometry/K-band spectrum are inconsistent with limiting cases of dust-free atmospheres ($COND$) and marginally inconsistent with the AMES/DUSTY models and the BT-SETTL models. However, ROXS 42Bb data are simultaneously fit by atmosphere models incorporating several micron-sized dust grains entrained in thick clouds, although further modifications are needed to better reproduce the $K$-band spectral shape. ROXs 42Bb's best-estimated temperature is $T_{eff}$ $\sim$ 1950--2000 $K$, near the low end of the empirically-derived range in Currie et al. (2014). For an age of $\sim$ 1--3 $Myr$ and considering the lifetime of the protostar phase, ROXs 42Bb's luminosity of log($L$/$L_{\odot}$) $\sim$ -3.07 $\pm$ 0.07 implies a mass of 9$^{+3}_{-3}$ $M_{J}$, making it one of the lightest planetary mass objects yet imaged.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Precise Water Abundance Measurement for the Hot Jupiter WASP-43b

The water abundance in a planetary atmosphere provides a key constraint on the planet's primordial origins because water ice is expected to play an important role in the core accretion model of planet formation. However, the water content of the Solar System giant planets is not well known because water is sequestered in clouds deep in their atmospheres. By contrast, short-period exoplanets have such high temperatures that their atmospheres have water in the gas phase, making it possible to measure the water abundance for these objects. We present a precise determination of the water abundance in the atmosphere of the 2 $M_\mathrm{Jup}$ short-period exoplanet WASP-43b based on thermal emission and transmission spectroscopy measurements obtained with the Hubble Space Telescope. We find the water content is consistent with the value expected in a solar composition gas at planetary temperatures (0.4-3.5x solar at 1 $σ$ confidence). The metallicity of WASP-43b's atmosphere suggested by this result extends the trend observed in the Solar System of lower metal enrichment for higher planet masses.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Atmospheric Characterization of the Hot Jupiter Kepler-13Ab

Kepler-13Ab (= KOI-13.01) is a unique transiting hot Jupiter. It is one of very few known short-period planets orbiting a hot A-type star, making it one of the hottest planets currently known. The availability of Kepler data allows us to measure the planet's occultation (secondary eclipse) and phase curve in the optical, which we combine with occultations observed by warm Spitzer at 4.5 mic and 3.6 mic and a ground-based occultation observation in the Ks band (2.1 mic). We derive a day-side hemisphere temperature of 2,750 +- 160 K as the effective temperature of a black body showing the same occultation depths. Comparing the occultation depths with one-dimensional planetary atmosphere models suggests the presence of an atmospheric temperature inversion. Our analysis shows evidence for a relatively high geometric albedo, Ag= 0.33 +0.04 -0.06. While measured with a simplistic method, a high Ag is supported also by the fact that the one-dimensional atmosphere models underestimate the occultation depth in the optical. We use stellar spectra to determine the dilution, in the four wide bands where occultation was measured, due to the visual stellar binary companion 1.15 +- 0.05" away. The revised stellar parameters measured using these spectra are combined with other measurements leading to revised planetary mass and radius estimates of Mp = 4.94 - 8.09 Mjup and Rp = 1.406 +- 0.038 Rjup. Finally, we measure a Kepler mid-occultation time that is 34.0 +- 6.9 s earlier than expected based on the mid-transit time and the delay due to light travel time, and discuss possible scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Characterization of the atmosphere of the hot Jupiter HAT-P-32Ab and the M-dwarf companion HAT-P-32B

We report secondary eclipse photometry of the hot Jupiter HAT-P-32Ab, taken with Hale/WIRC in H and Ks bands and with Spitzer/IRAC at 3.6 and 4.5 micron. We carried out adaptive optics imaging of the planet host star HAT-P-32A and its companion HAT-P-32B in the near-IR and the visible. We clearly resolve the two stars from each other and find a separation of 2.923" +/- 0. 004" and a position angle 110.64 deg +/- 0.12 deg. We measure the flux ratios of the binary in g' r' i' z' and H & Ks bands, and determine Teff = 3565 +/- 82 K for the companion star, corresponding to an M1.5 dwarf. We use PHOENIX stellar atmosphere models to correct the dilution of the secondary eclipse depths of the hot Jupiter due to the presence of the M1.5 companion. We also improve the secondary eclipse photometry by accounting for the non-classical, flux-dependent nonlinearity of the WIRC IR detector in the H band. We measure planet-to-star flux ratios of 0.090 +/- 0.033%, 0.178 +/- 0.057%, 0.364 +/- 0.016%, and 0.438 +/- 0.020% in the H, Ks, 3.6 and 4.5 micron bands, respectively. We compare these with planetary atmospheric models, and find they prefer an atmosphere with a temperature inversion and inefficient heat redistribution. However, we also find that the data are equally well-described by a blackbody model for the planet with Tp = 2042 +/- 50 K. Finally, we measure a secondary eclipse timing offset of 0.3 +/- 1.3 min from the predicted mid-eclipse time, which constrains e = 0.0072 +0.0700/-0.0064 when combined with RV data and is more consistent with a circular orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

CHARIS Science: Performance Simulations for the Subaru Telescope's Third-Generation of Exoplanet Imaging Instrumentation

We describe the expected scientific capabilities of CHARIS, a high-contrast integral-field spectrograph (IFS) currently under construction for the Subaru telescope. CHARIS is part of a new generation of instruments, enabled by extreme adaptive optics (AO) systems (including SCExAO at Subaru), that promise greatly improved contrasts at small angular separation thanks to their ability to use spectral information to distinguish planets from quasistatic speckles in the stellar point-spread function (PSF). CHARIS is similar in concept to GPI and SPHERE, on Gemini South and the Very Large Telescope, respectively, but will be unique in its ability to simultaneously cover the entire near-infrared $J$, $H$, and $K$ bands with a low-resolution mode. This extraordinarily broad wavelength coverage will enable spectral differential imaging down to angular separations of a few $λ/D$, corresponding to $\sim$$0.\!\!''1$. SCExAO will also offer contrast approaching $10^{-5}$ at similar separations, $\sim$$0.\!\!''1$--$0.\!\!''2$. The discovery yield of a CHARIS survey will depend on the exoplanet distribution function at around 10 AU. If the distribution of planets discovered by radial velocity surveys extends unchanged to $\sim$20 AU, observations of $\sim$200 mostly young, nearby stars targeted by existing high-contrast instruments might find $\sim$1--3 planets. Carefully optimizing the target sample could improve this yield by a factor of a few, while an upturn in frequency at a few AU could also increase the number of detections. CHARIS, with a higher spectral resolution mode of $R \sim 75$, will also be among the best instruments to characterize planets and brown dwarfs like HR 8799 cde and $κ$ And b.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraints on the Atmospheric Circulation and Variability of the Eccentric Hot Jupiter XO-3b

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5~$μ$m band taken with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We measure individual eclipse depths and center of eclipse times for a total of twelve secondary eclipses. We fit these data simultaneously with two transits observed in the same band in order to obtain a global best-fit secondary eclipse depth of $0.1580\pm 0.0036\%$ and a center of eclipse phase of $0.67004\pm 0.00013 $. We assess the relative magnitude of variations in the dayside brightness of the planet by measuring the size of the residuals during ingress and egress from fitting the combined eclipse light curve with a uniform disk model and place an upper limit of 0.05$\%$. The new secondary eclipse observations extend the total baseline from one and a half years to nearly three years, allowing us to place an upper limit on the periastron precession rate of $2.9\times 10^{-3}$ degrees/day the tightest constraint to date on the periastron precession rate of a hot Jupiter. We use the new transit observations to calculate improved estimates for the system properties, including an updated orbital ephemeris. We also use the large number of secondary eclipses to obtain the most stringent limits to date on the orbit-to-orbit variability of an eccentric hot Jupiter and demonstrate the consistency of multiple-epoch Spitzer observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deep Thermal Infrared Imaging of HR 8799 bcde: New Atmospheric Constraints and Limits on a Fifth Planet

We present new $L^\prime$ (3.8 $μm$) and Br-$α$ (4.05 $μm$) data and reprocessed archival $L^\prime$ data for the young, planet-hosting star HR 8799 obtained with Keck/NIRC2, VLT/NaCo and Subaru/IRCS. We detect all four HR 8799 planets in each dataset at a moderate to high signal-to-noise (SNR $\gtrsim$ 6-15). We fail to identify a fifth planet, "HR 8799 f", at $r$ $<$ 15 $AU$ at a 5-$σ$ confidence level: one suggestive, marginally significant residual at 0.2" is most likely a PSF artifact. Assuming companion ages of 30 $Myr$ and the Baraffe (Spiegel \& Burrows) planet cooling models, we rule out an HR 8799 f with mass of 5 $M_{J}$ (7 $M_{J}$), 7 $M_{J}$ (10 $M_{J}$), and 12 $M_{J}$ (13 $M_{J}$) at $r_{proj}$ $\sim$ 12 $AU$, 9 $AU$, and 5 $AU$, respectively. All four HR 8799 planets have red early T dwarf-like $L^\prime$ - [4.05] colors, suggesting that their SEDs peak in between the $L^\prime$ and $M^\prime$ broadband filters. We find no statistically significant difference in HR 8799 cde's colors. Atmosphere models assuming thick, patchy clouds appear to better match HR 8799 bcde's photometry than models assuming a uniform cloud layer. While non-equilibrium carbon chemistry is required to explain HR 8799 bc's photometry/spectra, evidence for it from HR 8799 de's photometry is weaker. Future, deep IR spectroscopy/spectrophotometry with the Gemini Planet Imager, SCExAO/CHARIS, and other facilities may clarify whether the planets are chemically similar or heterogeneous.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Expectations for the Hard X-ray Continuum and Gamma-ray Line Fluxes from the Type Ia supernova SN 2014J in M82

The hard X-ray continuum and gamma-ray lines from a Type Ia supernova dominate its integrated photon emissions and can provide unique diagnostics of the mass of the ejecta, the $^{56}$Ni yield and spatial distribution, its kinetic energy and expansion speed, and the mechanism of explosion. Such signatures and their time behavior "X-ray" the bulk debris field in direct fashion, and do not depend upon the oftimes problematic and elaborate UV, optical, and near-infrared spectroscopy and radiative transfer that have informed the study of these events for decades. However, to date no hard photons have ever been detected from a Type Ia supernova in explosion. With the advent of the supernova SN 2014J in M82, at a distance of $\sim$3.5 Mpc, this situation may soon change. Both NuSTAR and INTEGRAL have the potential to detect SN 2014J, and, if spectra and light curves can be measured, would usefully constrain the various explosion models published during the last $\sim$thirty years. In support of these observational campaigns, we provide predictions for the hard X-ray continuum and gamma-line emissions for fifteen Type Ia explosion models gleaned from the literature. The model set, containing as it does deflagration, delayed detonation, merger detonation, pulsational delayed detonation, and sub-Chandrasekhar helium detonation models, collectively spans a wide range of properties and, hence, signatures. We provide a brief discussion of various diagnostics (with examples), but importantly make the spectral and line results available electronically to aid in the interpretation of the anticipated data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gemini Planet Imager Spectroscopy of the HR 8799 planets c and d

During the first-light run of the Gemini Planet Imager (GPI) we obtained K-band spectra of exoplanets HR 8799 c and d. Analysis of the spectra indicates that planet d may be warmer than planet c. Comparisons to recent patchy cloud models and previously obtained observations over multiple wavelengths confirm that thick clouds combined with horizontal variation in the cloud cover generally reproduce the planets' spectral energy distributions. When combined with the 3 to 4 um photometric data points, the observations provide strong constraints on the atmospheric methane content for both planets. The data also provide further evidence that future modeling efforts must include cloud opacity, possibly including cloud holes, disequilibrium chemistry, and super-solar metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Highlights in the Study of Exoplanet Atmospheres

Exoplanets are now being discovered in profusion. However, to understand their character requires spectral models and data. These elements of remote sensing can yield temperatures, compositions, and even weather patterns, but only if significant improvements in both the parameter retrieval process and measurements are achieved. Despite heroic efforts to garner constraining data on exoplanet atmospheres and dynamics, reliable interpretation has oftimes lagged ambition. I summarize the most productive, and at times novel, methods employed to probe exoplanet atmospheres, highlight some of the most interesting results obtained, and suggest various broad theoretical topics in which further work could pay significant dividends.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recovery of the Candidate Protoplanet HD 100546 b with Gemini/NICI and Detection of Additional (Planet-Induced?) Disk Structure at Small Separations

We report the first independent, second-epoch (re-)detection of a directly-imaged protoplanet candidate. Using $L^\prime$ high-contrast imaging of HD 100546 taken with the Near-Infrared Coronagraph and Imager (NICI) on Gemini South, we recover `HD 100546 b' with a position and brightness consistent with the original VLT/NaCo detection from Quanz et al, although data obtained after 2013 will be required to decisively demonstrate common proper motion. HD 100546 b may be spatially resolved, up to $\approx$ 12-13 AU in diameter, and is embedded in a finger of thermal IR bright, polarized emission extending inwards to at least 0.3". Standard hot-start models imply a mass of $\approx$ 15 $M_{J}$. But if HD 100546 b is newly formed or made visible by a circumplanetary disk, both of which are plausible, its mass is significantly lower (e.g. 1--7 $M_{J}$). Additionally, we discover a thermal IR-bright disk feature, possibly a spiral density wave, at roughly the same angular separation as HD 100546 b but 90 degrees away. Our interpretation of this feature as a spiral arm is not decisive, but modeling analyses using spiral density wave theory implies a wave launching point exterior to $\approx$ 0.45" embedded within the visible disk structure: plausibly evidence for a second, hitherto unseen wide-separation planet. With one confirmed protoplanet candidate and evidence for 1--2 others, HD 100546 is an important evolutionary precursor to intermediate-mass stars with multiple super-jovian planets at moderate/wide separations like HR 8799.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scientific Return of Coronagraphic Exoplanet Imaging and Spectroscopy Using WFIRST

In this study, we explore and review the scientific potential for exoplanet characterization by a high-contrast optical coronagraph on WFIRST/AFTA. We suggest that the heterogeneity in albedo spectra and planet/star flux ratios as a function of orbital distance, planet mass, and composition expected for the giant exoplanets at separations from their primaries accessible to WFIRST will provide unique constraints on giant planet atmospheres, evolution, aerosol and cloud properties, and general theory. Such exoplanets are not merely extrapolations of Jupiter and Saturn, but are likely to occupy a rich continuum of varied behaviors. Each in themselves and jointly, optical spectra, photometry, and polarization measurements of a diverse population of giant exoplanets in the solar neighborhood has the potential to reveal a multitude of fundamental features of their gas-giant chemistry, atmospheres, and formation. Such a campaign will enrich our understanding of this class of planets beyond what is possible with even a detailed exploration of the giants in our own solar system, and will compliment ongoing studies of exoplanets in the infrared and on close-in orbits inaccessible to coronagraphy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The 4.5 $μ$m full-orbit phase curve of the hot Jupiter HD 209458b

The hot Jupiter HD 209458b is particularly amenable to detailed study as it is among the brightest transiting exoplanet systems currently known (V-mag = 7.65; K-mag = 6.308) and has a large planet-to-star contrast ratio. HD 209458b is predicted to be in synchronous rotation about its host star with a hot spot that is shifted eastward of the substellar point by superrotating equatorial winds. Here we present the first full-orbit observations of HD 209458b, in which its 4.5 $μ$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $μ$m measurement of HD 209458b's secondary eclipse emission downward by $\sim$35% to $0.1391%^{+0.0072%}_{-0.0069%}$, changing our interpretation of the properties of its dayside atmosphere. We find that the hot spot on the planet's dayside is shifted eastward of the substellar point by $40.9^{\circ}\pm{6.0^{\circ}}$, in agreement with circulation models predicting equatorial superrotation. HD 209458b's dayside (T$_{bright}$ = 1499 $\pm$ 15 K) and nightside (T$_{bright}$ = 972 $\pm$ 44 K) emission indicates a day-to-night brightness temperature contrast smaller than that observed for more highly irradiated exoplanets, suggesting that the day-to-night temperature contrast may be partially a function of the incident stellar radiation. The observed phase curve shape deviates modestly from global circulation model predictions potentially due to disequilibrium chemistry or deficiencies in the current hot CH$_{4}$ line lists used in these models. Observations of the phase curve at additional wavelengths are needed in order to determine the possible presence and spatial extent of a dayside temperature inversion, as well as to improve our overall understanding of this planet's atmospheric circulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Astronomical Reach of Fundamental Physics

Using basic physical arguments, we derive by dimensional and physical analysis the characteristic masses and sizes of important objects in the Universe in terms of just a few fundamental constants. This exercise illustrates the unifying power of physics and the profound connections between the small and the large in the Cosmos we inhabit. We focus on the minimum and maximum masses of normal stars, the corresponding quantities for neutron stars, the maximum mass of a rocky planet, the maximum mass of a white dwarf, and the mass of a typical galaxy. To zeroth order, we show that all these masses can be expressed in terms of either the Planck mass or the Chandrasekar mass, in combination with various dimensionless quantities. With these examples we expose the deep interrelationships imposed by Nature between disparate realms of the Universe and the amazing consequences of the unifying character of physical law.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Emergent 1.1-1.7 Micron Spectrum of the Exoplanet CoRoT-2b as Measured Using the Hubble Space Telescope

We have used Hubble/WFC3 and the G141 grism to measure the secondary eclipse of the transiting very hot Jupiter CoRoT-2b in the 1.1-1.7$μ$m spectral region. We find an eclipse depth averaged over this band equal to $395^{+69}_{-45}$ parts per million, equivalent to a blackbody temperature of $1788\pm18$K. We study and characterize several WFC3 instrumental effects, especially the "hook" phenomenon described by Deming et al. (2013). We use data from several transiting exoplanet systems to find a quantitative relation between the amplitude of the hook and the exposure level of a given pixel. Although the uncertainties in this relation are too large to allow us to develop an empirical correction for our data, our study provides a useful guide for optimizing exposure levels in future WFC3 observations. We derive the planet's spectrum using a differential method. The planet-to-star contrast increases to longer wavelength within the WFC3 bandpass, but without water absorption or emission to a $3σ$ limit of 85 ppm. The slope of the WFC3 spectrum is significantly less than the slope of the best-fit blackbody. We compare all existing eclipse data for this planet to a blackbody spectrum, and to spectra from both solar abundance and carbon-rich (C/O=1) models. A blackbody spectrum is an acceptable fit to the full dataset. Extra continuous opacity due to clouds or haze, and flattened temperature profiles, are strong candidates to produce quasi-blackbody spectra, and to account for the amplitude of the optical eclipses. Our results show ambiguous evidence for a temperature inversion in this planet.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Gemini Planet Imager: First Light

The Gemini Planet Imager (GPI) is a dedicated facility for directly imaging and spectroscopically characterizing extrasolar planets. It combines a very high-order adaptive optics system, a diffraction-suppressing coronagraph, and an integral field spectrograph with low spectral resolution but high spatial resolution. Every aspect of GPI has been tuned for maximum sensitivity to faint planets near bright stars. During first light observations, we achieved an estimated H band Strehl ratio of 0.89 and a 5-sigma contrast of $10^6$ at 0.75 arcseconds and $10^5$ at 0.35 arcseconds. Observations of Beta Pictoris clearly detect the planet, Beta Pictoris b, in a single 60-second exposure with minimal post-processing. Beta Pictoris b is observed at a separation of $434 \pm 6$ milli-arcseconds and position angle $211.8 \pm 0.5$ deg. Fitting the Keplerian orbit of Beta Pic b using the new position together with previous astrometry gives a factor of three improvement in most parameters over previous solutions. The planet orbits at a semi-major axis of $9.0^{+0.8}_{-0.4}$ AU near the 3:2 resonance with the previously-known 6 AU asteroidal belt and is aligned with the inner warped disk. The observations give a 4% posterior probability of a transit of the planet in late 2017.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Thermal structure of an exoplanet atmosphere from phase-resolved emission spectroscopy

Exoplanets that orbit close to their host stars are much more highly irradiated than their Solar System counterparts. Understanding the thermal structures and appearances of these planets requires investigating how their atmospheres respond to such extreme stellar forcing. We present spectroscopic thermal emission measurements as a function of orbital phase ("phase-curve observations") for the highly-irradiated exoplanet WASP-43b spanning three full planet rotations using the Hubble Space Telescope. With these data, we construct a map of the planet's atmospheric thermal structure, from which we find large day-night temperature variations at all measured altitudes and a monotonically decreasing temperature with pressure at all longitudes. We also derive a Bond albedo of 0.18 +0.07,-0.12 and an altitude dependence in the hot-spot offset relative to the substellar point.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Combined VLT and Gemini Study of the Atmosphere of the Directly-Imaged Planet, beta Pictoris b

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $β$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $β$ Pic b's near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects -- i.e. $κ$ And b and 1RXJ 1609B -- match $β$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $μm$ and 5 $μm$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $μm$) dust grains fail to reproduce the $β$ Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde but also with small (a few microns) modal particle sizes yield fits consistent with the data within uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ($<a>$ = 4 $μm$) we derive atmosphere parameters of log(g) = 3.8 $\pm$ 0.2 and $T_{eff}$ = 1575--1650 $K$, an inferred mass of 7$^{+4}_{-3}$ $M_{J}$, and a luminosity of log(L/L$_{\odot}$) $\sim$ -3.80 $\pm$ 0.02. The best-estimated planet radius, $\approx$ 1.65 $\pm$ 0.06 $R_{J}$, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges with constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if $β$ Pic b is younger than $\approx$ 7 Myr, consistent with a late formation, well after its host star's birth $\sim$ 12$^{+8}_{-4}$ Myr ago.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exoplanet Transit Spectroscopy Using WFC3: WASP-12 b, WASP-17 b, and WASP-19 b

We report analysis of transit spectroscopy of the extrasolar planets WASP-12 b, WASP-17 b, and WASP-19 b using the Wide Field Camera 3 on the HST. We analyze the data for a single transit for each planet using a strategy similar in certain aspects to the techniques used by Berta et al. (2012), but we extend their methodology to allow us to correct for channel- or wavelength-dependent instrumental effects by utilizing the band-integrated time series and measurements of the drift of the spectrum on the detector over time. We achieve almost photon-limited results for individual spectral bins, but the uncertainties in the transit depth for the the band-integrated data are exacerbated by the uneven sampling of the light curve imposed by the orbital phasing of HST's observations. Our final transit spectra for all three objects are consistent with the presence of a broad absorption feature at 1.4 microns potentially due to water. However, the amplitude of the absorption is less than that expected based on previous observations with Spitzer, possibly due to hazes absorbing in the NIR or non-solar compositions. The degeneracy of models with different compositions and temperature structures combined with the low amplitude of any features in the data preclude our ability to place unambiguous constraints on the atmospheric composition without additional observations with WFC3 to improve the S/N and/or a comprehensive multi-wavelength analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

HST Spectral Mapping of L/T Transition Brown Dwarfs Reveals Cloud Thickness Variations

Most directly imaged giant exoplanets are fainter than brown dwarfs with similar spectra. To explain their relative underluminosity unusually cloudy atmospheres have been proposed. However, with multiple parameters varying between any two objects, it remained difficult to observationally test this idea. We present a new method, sensitive time-resolved Hubble Space Telescope near-infrared spectroscopy, to study two rotating L/T transition brown dwarfs (2M2139 and SIMP0136). The observations provide spatially and spectrally resolved mapping of the cloud decks of the brown dwarfs. The data allow the study of cloud structure variations while other parameters are unchanged. We find that both brown dwarfs display variations of identical nature: J- and H-band brightness variations with minimal color and spectral changes. Our light curve models show that even the simplest surface brightness distributions require at least three elliptical spots. We show that for each source the spectral changes can be reproduced with a linear combination of only two different spectra, i.e. the entire surface is covered by two distinct types of regions. Modeling the color changes and spectral variations together reveal patchy cloud covers consisting of a spatially heterogenous mix of low-brightness, low-temperature thick clouds and brighter, thin and warm clouds. We show that the same thick cloud patches seen in our varying brown dwarf targets, if extended to the entire photosphere, predict near-infrared colors/magnitudes matching the range occupied by the directly imaged exoplanets that are cooler and less luminous than brown dwarfs with similar spectral types. This supports the models in which thick clouds are responsible for the near infrared properties of these underluminous exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Infrared Transmission Spectroscopy of the Exoplanets HD209458b and XO-1b Using the Wide Field Camera-3 on the Hubble Space Telescope

Exoplanetary transmission spectroscopy in the near-infrared using Hubble/NICMOS is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with Hubble/WFC3 provides an opportunity to resolve this ambiguity. We here report WFC3 spectroscopy of the giant planets HD209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6-percent (XO-1) and 26-percent (HD209458b) of the photon-limit at a spectral resolving power of 70, and are better than 0.01-percent per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 microns. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. (2002). Model atmospheres having uniformly-distributed extra opacity of 0.012 cm^2 per gram account approximately for both our water measurement and the sodium absorption in this planet. Our results for HD209458b support the picture advocated by Pont et al. (2013) in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and/or dust. However, the extra opacity needed for HD209458b is grayer than for HD189733b, with a weaker Rayleigh component.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Orbital Phase Variations of the Eccentric Giant Planet HAT-P-2b

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 μm bands of the Spitzer Space Telescope. The 3.6 and 4.5 μm data sets span an entire orbital period of HAT-P-2b, making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 μm that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39+/-0.28, 5.84+/-0.39, and 4.68+/-0.37 hours after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138%+/-0.0089%, 0.1162%+/-0.0080%, and 0.1888%+/-0.0072% in the 3.6, 4.5, and 8.0 μm bands respectively. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity and argument of periapse of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Secondary Eclipse Photometry of the Exoplanet WASP-5b with Warm Spitzer

We present photometry of the extrasolar planet WASP-5b in the 3.6 and 4.5 micron bands taken with the Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By examining the depth of the planet's secondary eclipse at these two wavelengths, we can place joint constraints on the planet's atmospheric pressure-temperature profile and chemistry. We measure secondary eclipse depths of 0.197% +/- 0.028% and 0.227% +/- 0.025% in the 3.6 micron and 4.5 micron bands, respectively. Our observations are best matched by models showing a hot dayside and, depending on our choice of model, a weak thermal inversion or no inversion at all. We measure a mean offset from the predicted center of eclipse of 0.078 +/- 0.032 hours, translating to ecos(omega) = 0.0031 +/- 0.0013 and consistent with a circular orbit. We see no evidence for any eclipse timing variations comparable to those reported in a previous transit study.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spectra as Windows into Exoplanet Atmospheres

Understanding a planet's atmosphere is a necessary condition for understanding not only the planet itself, but also its formation, structure, evolution, and habitability, This puts a premium on obtaining spectra, and developing credible interpretative tools with which to retrieve vital planetary information. However, for exoplanets these twin goals are far from being realized. In this paper, I provide a personal perspective on exoplanet theory and remote sensing via photometry and low-resolution spectroscopy. Though not a review in any sense, this paper highlights the limitations in our knowledge of compositions, thermal profiles, and the effects of stellar irradiation, focussing on, but not restricted to, transiting giant planets. I suggest that the true function of the recent past of exoplanet atmospheric research has been not to constrain planet properties for all time, but to train a new generation of scientists that, by rapid trial and error, is fast establishing a solid future foundation for a robust science of exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spectral and Photometric Diagnostics of Giant Planet Formation Scenarios

Gas-giant planets that form via core accretion might have very different characteristics from those that form via disk-instability. Disk-instability objects are typically thought to have higher entropies, larger radii, and (generally) higher effective temperatures than core-accretion objects. We provide a large set of models exploring the observational consequences of high-entropy (hot) and low-entropy (cold) initial conditions, in the hope that this will ultimately help to distinguish between different physical mechanisms of planet formation. However, the exact entropies and radii of newly-formed planets due to these two modes of formation cannot, at present, be precisely predicted. We introduce a broad range of "Warm Start" gas-giant planet models. Between the hottest and the coldest models that we consider, differences in radii, temperatures, luminosities, and spectra persist for only a few million to a few tens of millions of years for planets that are a few times Jupiter's mass or less. For planets that are ~five times Jupiter's mass or more, significant differences between hottest-start and coldest-start models persist for on the order of 100 Myrs. We find that out of the standard infrared bands (J, H, K, L', M, N) the K and H bands are the most diagnostic of the initial conditions. A hottest-start model can be from ~4.5 magnitudes brighter (at Jupiter's mass) to ~9 magnitudes brighter (at ten times Jupiter's mass) than a coldest-start model in the first few million years. In more massive objects, these large differences in luminosity and spectrum persist for much longer than in less massive objects. We consider the influence of atmospheric conditions on spectra, and find that the presence or absence of clouds, and the metallicity of an atmosphere, can affect an object's apparent brightness in different bands by up to several magnitudes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Dominance of Neutrino-Driven Convection in Core-Collapse Supernovae

Multi-dimensional instabilities have become an important ingredient in core-collapse supernova (CCSN) theory. Therefore, it is necessary to understand the driving mechanism of the dominant instability. We compare our parameterized three-dimensional CCSN simulations with other buoyancy-driven simulations and propose scaling relations for neutrino-driven convection. Through these comparisons, we infer that buoyancy-driven convection dominates post-shock turbulence in our simulations. In support of this inference, we present four major results. First, the convective fluxes and kinetic energies in the neutrino-heated region are consistent with expectations of buoyancy-driven convection. Second, the convective flux is positive where buoyancy actively drives convection, and the radial and tangential components of the kinetic energy are in rough equipartition (i.e. K_r ~ K_θ + K_ϕ). Both results are natural consequences of buoyancy-driven convection, and are commonly observed in simulations of convection. Third, buoyant driving is balanced by turbulent dissipation. Fourth, the convective luminosity and turbulent dissipation scale with the driving neutrino power. In all, these four results suggest that in neutrino-driven explosions, the multi-dimensional motions are consistent with neutrino-driven convection.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermal Processes Governing Hot-Jupiter Radii

There have been many proposed explanations for the larger-than-expected radii of some transiting hot Jupiters, including either stellar or orbital energy deposition deep in the atmosphere or deep in the interior. In this paper, we explore the important influences on hot-Jupiter radius evolution of (i) additional heat sources in the high atmosphere, the deep atmosphere, and deep in the convective interior; (ii) consistent cooling of the deep interior through the planetary dayside, nightside, and poles; (iii) the degree of heat redistribution to the nightside; and (iv) the presence of an upper atmosphere absorber inferred to produce anomalously hot upper atmospheres and inversions in some close-in giant planets. In particular, we compare the radius expansion effects of atmospheric and deep-interior heating at the same power levels and derive the power required to achieve a given radius increase when night-side cooling is incorporated. We find that models that include consistent day/night cooling are more similar to isotropically irradiated models when there is more heat redistributed from the dayside to the nightside. In addition, we consider the efficacy of ohmic heating in the atmosphere and/or convective interior in inflating hot Jupiters. Among our conclusions are that (i) the most highly irradiated planets cannot stably have uB > (10 km/s Gauss) over a large fraction of their daysides, where u is the zonal wind speed and B is the dipolar magnetic field strength in the atmosphere, and (ii) that ohmic heating cannot in and of itself lead to a runaway in planet radius.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transport Equations for Oscillating Neutrinos

We derive a suite of generalized Boltzmann equations, based on the density-matrix formalism, that incorporates the physics of neutrino oscillations for two- and three-flavor oscillations, matter refraction, and self-refraction. The resulting equations are straightforward extensions of the classical transport equations that nevertheless contain the full physics of quantum oscillation phenomena. In this way, our broadened formalism provides a bridge between the familiar neutrino transport algorithms employed by supernova modelers and the more quantum-heavy approaches frequently employed to illuminate the various neutrino oscillation effects. We also provide the corresponding angular-moment versions of this generalized equation set. Our goal is to make it easier for astrophysicists to address oscillation phenomena in a language with which they are familiar. The equations we derive are simple and practical, and are intended to facilitate progress concerning oscillation phenomena in the context of core-collapse supernova theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Warm Spitzer and Palomar Near-IR Secondary Eclipse Photometry of Two Hot Jupiters: WASP-48b and HAT-P-23b

We report secondary eclipse photometry of two hot Jupiters, WASP-48b and HAT-P-23b, at 3.6 and 4.5 um taken with the InfraRed Array Camera aboard the Spitzer Space Telescope during the warm Spitzer mission and in the H and Ks bands with the Wide Field IR Camera at the Palomar 200-inch Hale Telescope. WASP-48b and HAT-P-23b are Jupiter-mass and twice Jupiter-mass objects orbiting an old, slightly evolved F star and an early G dwarf star, respectively. In the H, Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.047% +/- 0.016%, 0.109% +/- 0.027%, 0.176% +/- 0.013%, and 0.214% +/- 0.020% for WASP-48b. In the Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.234% +/- 0.046%, 0.248% +/- 0.019%, and 0.309% +/- 0.026% for HAT-P-23b. For WASP-48b and HAT-P-23b, respectively, we measure delays of 2.6 +/- 3.9 minutes and 4.0 +/- 2.4 minutes relative to the predicted times of secondary eclipse for circular orbits, placing 2-sigma upper limits on |e cos(omega)| of 0.0053 and 0.0080, both of which are consistent with circular orbits. The dayside emission spectra of these planets are well-described by blackbodies with effective temperatures of 2158 +/- 100 K (WASP-48b) and 2154 +/- 90 K (HAT-P-23b), corresponding to moderate recirculation in the zero albedo case. Our measured eclipse depths are also consistent with one-dimensional radiative transfer models featuring varying degrees of recirculation and weak thermal inversions or no inversions at all. We discuss how the absence of strong temperature inversions on these planets may be related to the activity levels and metallicities of their host stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Warm Spitzer Photometry of Three Hot Jupiters: HAT-P-3b, HAT-P-4b and HAT-P-12b

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5 micron bands. HAT-P-3b and HAT-P-4b are Jupiter-mass, objects orbiting an early K and an early G dwarf stars, respectively. For HAT-P-3b we find eclipse depths of 0.112%+0.015%-0.030% (3.6 micron) and 0.094%+0.016%-0.009% (4.5 micron). The HAT-P-4b values are 0.142%+0.014%-0.016% (3.6 micron) and 0.122%+0.012%-0.014% (4.5micron). The two planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT-P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, along with hot Neptune GJ 436b and hot Saturn WASP-29b. We are able to place 3$σ$ upper limits on the secondary eclipse depth of HAT-P-12b in both wavelengths: < 0.042% (3.6 micron) and <0.085% (4.5 micron). We discuss these results in the context of the {\it Spitzer} secondary eclipse measurements of GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. We place 3 sigma upper limits on the |e cos(omega)| quantity (where e is eccentricity and omega is the argument of periapsis) for HAT-P-3b (<0.0081) and HAT-P-4b (<0.0042), based on the secondary eclipse timings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Monte Carlo Method for Time-Dependent Neutrino Radiation Transport

Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck & Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Investigation into the Character of Pre-Explosion Core-Collapse Supernova Shock Motion

We investigate the structure of the stalled supernova shock in both 2D and 3D and explore the differences in the effects of neutrino heating and the standing accretion shock instability (SASI). We find that early on the amplitude of the dipolar mode of the shock is factors of 2 to 3 smaller in 3D than in 2D. However, later in both 3D and 2D the monopole and dipole modes start to grow until explosion. Whereas in 2D the (l,m) = (1,0) mode changes sign quasi-periodically, producing the "up-and-down" motion always seen in modern 2D simulations, in 3D this almost never happens. Rather, in 3D when the dipolar mode starts to grow, it grows in magnitude and wanders stochastically in direction until settling before explosion to a particular patch of solid angle. In 2D we find that the amplitude of the dipolar shock deformation separates into two classes. For the first, identified with the SASI and for a wide range of "low" neutrino luminosities, this amplitude remains small and roughly constant. For the other, identified with higher luminosities and neutrino-driven convection, the dipolar amplitude grows sharply. Importantly, it is only for this higher luminosity class that we see neutrino-driven explosions within ~1 second of bounce. Moreover, for the "low" luminosity runs, the power spectra of these dipolar oscillations peak in the 30-50 Hz range associated with advection timescales, while for the high-luminosity runs the power spectra at lower frequencies are significantly more prominent. We associate this enhanced power at lower frequencies with slower convective effects and the secular growth of the dipolar shock amplitude. On the basis of our study, we hypothesize that neutrino-driven buoyant convection should almost always dominate the SASI when the supernova explosion is neutrino-driven.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dimensional Dependence of the Hydrodynamics of Core-Collapse Supernovae

The multidimensional character of the hydrodynamics in core-collapse supernova (CCSN) cores is a key facilitator of explosions. Unfortunately, much of this work has necessarily been performed assuming axisymmetry and it remains unclear whether or not this compromises those results. In this work, we present analyses of simplified two- and three-dimensional CCSN models with the goal of comparing the multidimensional hydrodynamics in setups that differ only in dimension. Not surprisingly, we find many differences between 2D and 3D models. While some differences are subtle and perhaps not crucial to understanding the explosion mechanism, others are quite dramatic and make interpreting 2D CCSN models problematic. In particular, we find that imposing axisymmetry artificially produces excess power at the largest spatial scales, power that has been deemed critical in the success of previous explosion models and has been attributed solely to the standing accretion shock instability. Nevertheless, our 3D models, which have an order of magnitude less power on large scales compared to 2D models, explode earlier. Since we see explosions earlier in 3D than in 2D, the vigorous sloshing associated with the large scale power in 2D models is either not critical in any dimension or the explosion mechanism operates differently in 2D and 3D. Possibly related to the earlier explosions in 3D, we find that about 25% of the accreted material spends more time in the gain region in 3D than in 2D, being exposed to more integrated heating and reaching higher peak entropies, an effect we associate with the differing characters of turbulence in 2D and 3D. Finally, we discuss a simple model for the runaway growth of buoyant bubbles that is able to quantitatively account for the growth of the shock radius and predicts a critical luminosity relation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Direct Imaging Confirmation and Characterization of a Dust-Enshrouded Candidate Exoplanet Orbiting Fomalhaut

We present Subaru/IRCS J band data for Fomalhaut and a (re)reduction of archival 2004--2006 HST/ACS data first presented by Kalas et al. (2008). We confirm the existence of a candidate exoplanet, Fomalhaut b, in both the 2004 and 2006 F606W data sets at a high signal-to-noise. Additionally, we confirm the detection at F814W and present a new detection in F435W. Fomalhaut b's space motion may be consistent with it being in an apsidally-aligned, non debris ring-crossing orbit, although new astrometry is required for firmer conclusions. We cannot confirm that Fomalhaut b exhibits 0.7-0.8 mag variability cited as evidence for planet accretion or a semi-transient dust cloud. The new, combined optical SED and IR upper limits confirm that emission identifying Fomalhaut b originates from starlight scattered by small dust, but this dust is most likely associated with a massive body. The Subaru and IRAC/4.5 micron upper limits imply M < 2 Mj, still consistent with the range of Fomalhaut b masses needed to sculpt the disk. Fomalhaut b is very plausibly "a planet identified from direct imaging" even if current images of it do not, strictly speaking, show thermal emission from a directly imaged planet.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

First Light LBT AO Images of HR 8799 bcde at 1.65 and 3.3 Microns: New Discrepancies between Young Planets and Old Brown Dwarfs

As the only directly imaged multiple planet system, HR 8799 provides a unique opportunity to study the physical properties of several planets in parallel. In this paper, we image all four of the HR 8799 planets at H-band and 3.3 microns with the new LBT adaptive optics system, PISCES, and LBTI/LMIRCam. Our images offer an unprecedented view of the system, allowing us to obtain H and 3.3$ micron photometry of the innermost planet (for the first time) and put strong upper-limits on the presence of a hypothetical fifth companion. We find that all four planets are unexpectedly bright at 3.3 microns compared to the equilibrium chemistry models used for field brown dwarfs, which predict that planets should be faint at 3.3 microns due to CH4 opacity. We attempt to model the planets with thick-cloudy, non-equilibrium chemistry atmospheres, but find that removing CH4 to fit the 3.3 micron photometry increases the predicted L' (3.8 microns) flux enough that it is inconsistent with observations. In an effort to fit the SED of the HR 8799 planets, we construct mixtures of cloudy atmospheres, which are intended to represent planets covered by clouds of varying opacity. In this scenario, regions with low opacity look hot and bright, while regions with high opacity look faint, similar to the patchy cloud structures on Jupiter and L/T transition brown-dwarfs. Our mixed cloud models reproduce all of the available data, but self-consistent models are still necessary to demonstrate their viability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

NIMBUS: The Near-Infrared Multi-Band Ultraprecise Spectroimager for SOFIA

We present a new and innovative near-infrared multi-band ultraprecise spectroimager (NIMBUS) for SOFIA. This design is capable of characterizing a large sample of extrasolar planet atmospheres by measuring elemental and molecular abundances during primary transit and occultation. This wide-field spectroimager would also provide new insights into Trans-Neptunian Objects (TNO), Solar System occultations, brown dwarf atmospheres, carbon chemistry in globular clusters, chemical gradients in nearby galaxies, and galaxy photometric redshifts. NIMBUS would be the premier ultraprecise spectroimager by taking advantage of the SOFIA observatory and state of the art infrared technologies. This optical design splits the beam into eight separate spectral bandpasses, centered around key molecular bands from 1 to 4 microns. Each spectral channel has a wide field of view for simultaneous observations of a reference star that can decorrelate time-variable atmospheric and optical assembly effects, allowing the instrument to achieve ultraprecise calibration for imaging and photometry for a wide variety of astrophysical sources. NIMBUS produces the same data products as a low-resolution integral field spectrograph over a large spectral bandpass, but this design obviates many of the problems that preclude high-precision measurements with traditional slit and integral field spectrographs. This instrument concept is currently not funded for development.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Perspectives on Core-Collapse Supernova Theory

Core-collapse theory brings together many facets of high-energy and nuclear astrophysics and the numerical arts to present theorists with one of the most important, yet frustrating, astronomical questions: "What is the mechanism of core-collapse supernova explosions?" A review of all the physics and the fifty-year history involved would soon bury the reader in minutiae that could easily obscure the essential elements of the phenomenon, as we understand it today. Moreover, much remains to be discovered and explained, and a complicated review of an unresolved subject in flux could grow stale fast. Therefore, in this paper I describe what I think are various important facts and perspectives that may have escaped the attention of those interested in this puzzle. Furthermore, I attempt to describe the modern theory's physical underpinnings and briefly summarize the current state of play. In the process, I identify a few myths (as I see them) that have crept into modern discourse. However, there is much more to do and humility in the face of this age-old challenge is clearly the most prudent stance as we seek its eventual resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pulsation Frequencies and Modes of Giant Exoplanets

We calculate the eigenfrequencies and eigenfunctions of the acoustic oscillations of giant exoplanets and explore the dependence of the characteristic frequency and the eigenfrequencies on several parameters: the planet mass, the planet radius, the core mass, and the heavy element mass fraction in the envelope. We provide the eigenvalues for degree $l$ up to 8 and radial order n up to 12. For the selected values of l and n, we find that the pulsation eigenfrequencies depend strongly on the planet mass and radius, especially at high frequency. We quantify this dependence through the calculation of the characteristic frequency which gives us an estimate of the scale of the eigenvalue spectrum at high frequency. For the mass range 0.5 < M_P < 15 M_J, and fixing the planet radius to the Jovian value, we find that the characteristic frequency is ~164.0 * (M_P/M_J)^(0.48) microHz, where M_P is the planet mass and M_J is Jupiter's mass. For the radius range from 0.9 to 2.0 R_J, and fixing the planet's mass to the Jovian value, we find that the characteristic frequency is ~164.0 * (R_P/R_J)^(-2.09) microHz, where R_P is the planet radius and R_J is Jupiter's radius. We explore the influence of the presence of a dense core on the pulsation frequencies and on the characteristic frequency of giant exoplanets. We find that the presence of heavy elements in the envelope affects the eigenvalue distribution in ways similar to the presence of a dense core. Additionally, we apply our formalism to Jupiter and Saturn and find results consistent with both the observationnal data of Gaulme et al. (2011) and previous theoretical work.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scientific Design of a High Contrast Integral Field Spectrograph for the Subaru Telescope

Ground-based telescopes equipped with adaptive-optics (AO) systems and specialized science cameras are now capable of directly detecting extrasolar planets. We present the expected scientific capabilities of CHARIS, the Coronagraphic High Angular Resolution Imaging Spectrograph, which is being built for the Subaru 8.2 m telescope of the National Astronomical Observatory of Japan. CHARIS will be implemented behind the new extreme adaptive optics system at Subaru, SCExAO, and the existing 188-actuator system AO188. CHARIS will offer three observing modes over near-infrared wavelengths from 0.9 to 2.4 microns (the y-, J-, H-, and K-bands), including a low-spectral-resolution mode covering this entire wavelength range and a high-resolution mode within a single band. With these capabilities, CHARIS will offer exceptional sensitivity for discovering giant exoplanets, and will enable detailed characterization of their atmospheres. CHARIS, the only planned high-contrast integral field spectrograph on an 8m-class telescope in the Northern Hemisphere, will complement the similar instruments such as Project 1640 at Palomar, and GPI and SPHERE in Chile.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Gemini NICI Planet-Finding Campaign: Discovery of a Multiple System Orbiting the Young A Star HD 1160

We report the discovery by the Gemini NICI Planet-Finding Campaign of two low-mass companions to the young A0V star HD 1160 at projected separations of 81 +/- 5 AU (HD 1160 B) and 533 +/- 25 AU (HD 1160 C). VLT images of the system taken over a decade for the purpose of using HD 1160 A as a photometric calibrator confirm that both companions are physically associated. By comparing the system to members of young moving groups and open clusters with well-established ages, we estimate an age of 50 (+50,-40) Myr for HD 1160 ABC. While the UVW motion of the system does not match any known moving group, the small magnitude of the space velocity is consistent with youth. Near-IR spectroscopy shows HD 1160 C to be an M3.5 +/- 0.5 star with an estimated mass of 0.22 (+0.03,-0.04) M_Sun, while NIR photometry of HD 1160 B suggests a brown dwarf with a mass of 33 (+12,-9) M_Jup. The very small mass ratio (0.014) between the A and B components of the system is rare for A star binaries, and would represent a planetary-mass companion were HD 1160 A to be slightly less massive than the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Vertical Atmospheric Structure in a Variable Brown Dwarf: Pressure-dependent Phase Shifts in Simultaneous Hubble Space Telescope-Spitzer Light Curves

Heterogeneous clouds or temperature perturbations in rotating brown dwarfs produce variability in the observed flux. We report time-resolved simultaneous observations of the variable T6.5 brown dwarf 2MASSJ22282889-431026 over the wavelength ranges 1.1-1.7 microns and broadband 4.5 microns. Spectroscopic observations were taken with Wide Field Camera 3 on board the Hubble Space Telescope and photometry with the Spitzer Space Telescope. The object shows sinusoidal infrared variability with a period of 1.4 hr at most wavelengths with peak-to-peak amplitudes between 1.45% and 5.3% of the mean flux. While the light curve shapes are similar at all wavelengths, their phases differ from wavelength to wavelength with a maximum difference of more than half of a rotational period. We compare the spectra with atmospheric models of different cloud prescriptions, from which we determine the pressure levels probed at different wavelengths. We find that the phase lag increases with decreasing pressure level, or higher altitude. We discuss a number of plausible scenarios that could cause this trend of light curve phase with probed pressure level. These observations are the first to probe heterogeneity in an ultracool atmosphere in both horizontal and vertical directions, and thus are an ideal test case for realistic three dimensional simulations of the atmospheric structure with clouds in brown dwarfs and extrasolar planets.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A 5 Micron Image of beta Pictoris b at a Sub-Jupiter Projected Separation: Evidence for a Misalignment Between the Planet and the Inner, Warped Disk

We present and analyze a new M' detection of the young exoplanet beta Pictoris b from 2008 VLT/NaCo data at a separation of ~ 4 AU and a high signal-to-noise rereduction of L' data taken in December 2009. Based on our orbital analysis, the planet's orbit is viewed almost perfectly edge-on (i ~ 89 degrees) and has a Saturn-like semimajor axis of 9.50 (+3.93, -1.7) AU. Intriguingly, the planet's orbit is aligned with the major axis of the outer disk (Omega ~ 31 degrees) but probably misaligned with the warp/inclined disk at 80 AU often cited as a signpost for the planet's existence. Our results motivate new studies to clarify how $β$ Pic b sculpts debris disk structures and whether a second planet is required to explain the warp/inclined disk.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analytic models for albedos, phase curves and polarization of reflected light from exoplanets

New observational facilities are becoming increasingly capable of observing reflected light from transiting and directly imaged extrasolar planets. In this study, we provide an analytic framework to interpret observed phase curves, geometric albedos, and polarization of giant planet atmospheres. We compute the observables for non-conservative Rayleigh scattering in homogeneous semi-infinite atmospheres using both scalar and vector formalisms. In addition, we compute phase curves and albedos for Lambertian, isotropic, and anisotropic scattering phase functions. We provide analytic expressions for geometric albedos and spherical albedos as a function of the scattering albedo for Rayleigh scattering in semi-infinite atmospheres. Given an observed geometric albedo our prescriptions can be used to estimate the underlying scattering albedo of the atmosphere, which in turn is indicative of the scattering and absorptive properties of the atmosphere. We also study the dependence of polarization in Rayleigh scattering atmospheres on the orbital parameters of the planet-star system, particularly on the orbital inclination. We show how the orbital inclination of non-transiting exoplanets can be constrained from their observed polarization parameters. We consolidate the formalism, solution techniques, and results from analytic models available in the literature, often scattered in various sources, and present a systematic procedure to compute albedos, phase curves, and polarization of reflected light.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Model atmospheres for massive gas giants with thick clouds: Application to the HR 8799 planets and predictions for future detections

We have generated an extensive new suite of massive giant planet atmosphere models and used it to obtain fits to photometric data for the planets HR 8799b, c, and d. We consider a wide range of cloudy and cloud-free models. The cloudy models incorporate different geometrical and optical thicknesses, modal particle sizes, and metallicities. For each planet and set of cloud parameters, we explore grids in gravity and effective temperature, with which we determine constraints on the planet's mass and age. Our new models yield statistically significant fits to the data, and conclusively confirm that the HR 8799 planets have much thicker clouds than those required to explain data for typical L and T dwarfs. Both models with 1) physically thick forsterite clouds and a 60-micron modal particle size and 2) clouds made of 1 micron-sized pure iron droplets and 1% supersaturation fit the data. Current data are insufficient to accurately constrain the microscopic cloud properties, such as composition and particle size. The range of best-estimated masses for HR 8799b, HR 8799c, and HR 8799d conservatively span 2-12 M_J, 6-13 M_J, and 3-11 M_J, respectively and imply coeval ages between ~10 and ~150 Myr, consistent with previously reported stellar age. The best-fit temperatures and gravities are slightly lower than values obtained by Currie et al. (2011) using even thicker cloud models. Finally, we use these models to predict the near-to-mid IR colors of soon-to-be imaged planets. Our models predict that planet-mass objects follow a locus in some near-to-mid IR color-magnitude diagrams that is clearly separable from the standard L/T dwarf locus for field brown dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Dependence of Brown Dwarf Radii on Atmospheric Metallicity and Clouds: Theory and Comparison with Observations

Employing realistic and consistent atmosphere boundary conditions, we have generated evolutionary models for brown dwarfs and very-low-mass stars (VLMs) for different metallicities ([Fe/H]), with and without clouds. We find that the spread in radius at a given mass and age can be as large as $\sim$10% to $\sim$25%, with higher-metallicity, higher-cloud-thickness atmospheres resulting quite naturally in larger radii. For each 0.1 dex increase in [Fe/H], radii increase by $\sim$1% to $\sim$2.5%, depending upon age and mass. We also find that, while for smaller masses and older ages brown dwarf radii decrease with increasing helium fraction ($Y$) (as expected), for more massive brown dwarfs and a wide range of ages they increase with helium fraction. The increase in radius in going from $Y=0.25$ to $Y=0.28$ can be as large as $\sim$0.025 \rj\ ($\sim$2.5%). Furthermore, we find that for VLMs an increase in atmospheric metallicity from 0.0 to 0.5 dex, increases radii by $\sim$4%, and from -0.5 to 0.5 dex by $\sim$10%. Therefore, we suggest that opacity due to higher metallicity might naturally account for the apparent radius anomalies in some eclipsing VLM systems. Ten to twenty-five percent variations in radius exceed errors stemming from uncertainities in the equation of state alone. This serves to emphasize that transit and eclipse measurements of brown dwarf radii constrain numerous effects collectively, importantly including the atmosphere and condensate cloud models, and not just the equation of state. At all times, one is testing a multi-parameter theory, and not a universal radius$-$mass relation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Impact of Circumplantary Jets on Transit Spectra and Timing Offsets for Hot-Jupiters

We present theoretical wavelength-dependent transit light curves for the giant planet HD209458b based on a number of state of the art 3D radiative hydrodynamical models. By varying the kinematic viscosity in the model we calculate observable signatures associated with the emergence of a super-rotating circumplanetary jet that strengthens with decreased viscosity. We obtain excellent agreement between our mid-transit transit spectra and existing data from Hubble and Spitzer, finding the best fit for intermediate values of viscosity. We further exploit dynamically driven differences between eastern and western hemispheres to extract the spectral signal imparted by a circumplanetary jet. We predict that: (i) the transit depth should decrease as the jet becomes stronger; (ii) the measured transit times should show timing offsets of up to 6 seconds at wavelengths with higher opacity, which increases with jet strength; (iii) wavelength-dependent differences between ingress and egress spectra increase with jet strength; (iv) the color-dependent transit shape should vary more strongly with wavelength for stronger jets. These techniques and trends should be valid for other hot Jupiters as well. Observations of transit timing offsets may be accessible with current instrumentation, though the other predictions may require the capabilities JWST and other future missions. Hydrodynamical models utilized solve the 3D Navier-Stokes equations together with decoupled thermal and radiative energy equations and wavelength dependent stellar heating.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Warm Spitzer Photometry of XO-4b, HAT-P-6b and HAT-P-8b

We have analyzed Warm Spitzer/IRAC observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6μm and 4.5μm for each target. XO-4b exhibits a stronger eclipse depth at 4.5μm than at 3.6μm, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6μm, and is best-described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversions and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. (2010) and Cowan & Agol (2011), respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for e cos(ω) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Atmosphere and Spectral Models of the Kepler-Field Planets HAT-P-7b and TrES-2

We develop atmosphere models of two of the three Kepler-field planets that were known prior to the start of the Kepler mission (HAT-P-7b and TrES-2). We find that published Kepler and Spitzer data for HAT-P-7b appear to require an extremely hot upper atmosphere on the dayside, with a strong thermal inversion and little day-night redistribution. The Spitzer data for TrES-2 suggest a mild thermal inversion with moderate day-night redistribution. We examine the effect of nonequilibrium chemistry on TrES-2 model atmospheres and find that methane levels must be adjusted by extreme amounts in order to cause even mild changes in atmospheric structure and emergent spectra. Our best-fit models to the Spitzer data for TrES-2 lead us to predict a low secondary eclipse planet-star flux ratio (~2 x 10^-5) in the Kepler bandpass, which is consistent with what very recent observations have found. Finally, we consider how the Kepler-band optical flux from a hot exoplanet depends on the strength of a possible extra optical absorber in the upper atmosphere. We find that the optical flux is not monotonic in optical opacity, and the non-monotonicity is greater for brighter, hotter stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Day-side z'-band emission and eccentricity of Wasp-12b

We report the detection of the eclipse of the very-hot Jupiter WASP-12b via z'-band time-series photometry obtained with the 3.5-meter ARC telescope at Apache Point Observatory. We measure a decrease in flux of 0.082+/-0.015% during the passage of the planet behind the star. That planetary flux is equally well reproduced by atmospheric models with and without extra absorbers, and blackbody models with f > 0.585+/-0.080. It is therefore necessary to measure the planet at other wavelengths to further constrain its atmospheric properties. The eclipse appears centered at phase = 0.5100 (+0.0072,-0.0061), consistent with an orbital eccentricity of |e cos w| = 0.016 (+0.011,-0.009) (see note at end of Section 4). If the orbit of the planet is indeed eccentric, the large radius of WASP-12b can be explained by tidal heating.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Models of Neptune-Mass Exoplanets: Emergent Fluxes and Albedos

There are now many known exoplanets with Msin(i) within a factor of two of Neptune's, including the transiting planets GJ436b and HAT-P-11b. Planets in this mass-range are different from their more massive cousins in several ways that are relevant to their radiative properties and thermal structures. By analogy with Neptune and Uranus, they are likely to have metal abundances that are an order of magnitude or more greater than those of larger, more massive planets. This increases their opacity, decreases Rayleigh scattering, and changes their equation of state. Furthermore, their smaller radii mean that fluxes from these planets are roughly an order of magnitude lower than those of otherwise identical gas giant planets. Here, we compute a range of plausible radiative equilibrium models of GJ436b and HAT-P-11b. In addition, we explore the dependence of generic Neptune-mass planets on a range of physical properties, including their distance from their host stars, their metallicity, the spectral type of their stars, the redistribution of heat in their atmospheres, and the possible presence of additional optical opacity in their upper atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Photometric and Spectral Signatures of 3D Models of Transiting Giant Exoplanets

Using a 3D GCM, we create dynamical model atmospheres of a representative transiting giant exoplanet, HD 209458b. We post-process these atmospheres with an opacity code to obtain transit radius spectra during the primary transit. Using a spectral atmosphere code, we integrate over the face of the planet seen by an observer at various orbital phases and calculate light curves as a function of wavelength and for different photometric bands. The products of this study are generic predictions for the phase variations of a zero-eccentricity giant planet's transit spectrum and of its light curves. We find that for these models the temporal variations in all quantities and the ingress/egress contrasts in the transit radii are small ($< 1.0$\%). Moreover, we determine that the day/night contrasts and phase shifts of the brightness peaks relative to the ephemeris are functions of photometric band. The $J$, $H$, and $K$ bands are shifted most, while the IRAC bands are shifted least. Therefore, we verify that the magnitude of the downwind shift in the planetary ``hot spot" due to equatorial winds is strongly wavelength-dependent. The phase and wavelength dependence of light curves, and the associated day/night contrasts, can be used to constrain the circulation regime of irradiated giant planets and to probe different pressure levels of a hot Jupiter atmosphere. We posit that though our calculations focus on models of HD 209458b similar calculations for other transiting hot Jupiters in low-eccentricity orbits should yield transit spectra and light curves of a similar character.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Results From Core-Collapse Simulations with Multi-Dimensional, Multi-Angle Neutrino Transport

We present new results from the only 2D multi-group, multi-angle calculations of core-collapse supernova evolution. The first set of results from these calculations was published in Ott et al. (2008). We have followed a nonrotating and a rapidly rotating 20 solar mass model for ~400 ms after bounce. We show that the radiation fields vary much less with angle than the matter quantities in the region of net neutrino heating. This obtains because most neutrinos are emitted from inner radiative regions and because the specific intensity is an integral over sources from many angles at depth. The latter effect can only be captured by multi-angle transport. We then compute the phase relationship between dipolar oscillations in the shock radius and in matter and radiation quantities throughout the postshock region. We demonstrate a connection between variations in neutrino flux and the hydrodynamical shock oscillations, and use a variant of the Rayleigh test to estimate the detectability of these neutrino fluctuations in IceCube and Super-K. Neglecting flavor oscillations, fluctuations in our nonrotating model would be detectable to ~10 kpc in IceCube, and a detailed power spectrum could be measured out to ~5 kpc. These distances are considerably lower in our rapidly rotating model or with significant flavor oscillations. Finally, we measure the impact of rapid rotation on detectable neutrino signals. Our rapidly rotating model has strong, species-dependent asymmetries in both its peak neutrino flux and its light curves. The peak flux and decline rate show pole-equator ratios of up to ~3 and ~2, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

We present photometry of the giant extrasolar planet WASP-4b at 3.6 and 4.5 micron taken with the Infrared Array Camera on board the Spitzer Space Telescope as part of Spitzer's extended warm mission. We find secondary eclipse depths of 0.319+/-0.031% and 0.343+/-0.027% for the 3.6 and 4.5 micron bands, respectively and show model emission spectra and pressure-temperature profiles for the planetary atmosphere. These eclipse depths are well fit by model emission spectra with water and other molecules in absorption, similar to those used for TrES-3 and HD 189733b. Depending on our choice of model, these results indicate that this planet has either a weak dayside temperature inversion or no inversion at all. The absence of a strong thermal inversion on this highly irradiated planet is contrary to the idea that highly irradiated planets are expected to have inversions, perhaps due the presence of an unknown absorber in the upper atmosphere. This result might be explained by the modestly enhanced activity level of WASP-4b's G7V host star, which could increase the amount of UV flux received by the planet, therefore reducing the abundance of the unknown stratospheric absorber in the planetary atmosphere as suggested in Knutson et al. (2010). We also find no evidence for an offset in the timing of the secondary eclipse and place a 2 sigma upper limit on |ecos(omega)| of 0.0024, which constrains the range of tidal heating models that could explain this planet's inflated radius.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Deuterium-Burning Mass Limit for Brown Dwarfs and Giant Planets

There is no universally acknowledged criterion to distinguish brown dwarfs from planets. Numerous studies have used or suggested a definition based on an object's mass, taking the ~13-Jupiter mass (M_J) limit for the ignition of deuterium. Here, we investigate various deuterium-burning masses for a range of models. We find that, while 13 M_J is generally a reasonable rule of thumb, the deuterium fusion mass depends on the helium abundance, the initial deuterium abundance, the metallicity of the model, and on what fraction of an object's initial deuterium abundance must combust in order for the object to qualify as having burned deuterium. Even though, for most proto-brown dwarf conditions, 50% of the initial deuterium will burn if the object's mass is ~(13.0 +/- 0.8)M_J, the full range of possibilities is significantly broader. For models ranging from zero-metallicity to more than three times solar metallicity, the deuterium burning mass ranges from ~11.0 M_J (for 3-times solar metallicity, 10% of initial deuterium burned) to ~16.3 M_J (for zero metallicity, 90% of initial deuterium burned).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Gemini NICI Planet-Finding Campaign

Our team is carrying out a multi-year observing program to directly image and characterize young extrasolar planets using the Near-Infrared Coronagraphic Imager (NICI) on the Gemini-South 8.1-meter telescope. NICI is the first instrument on a large telescope designed from the outset for high-contrast imaging, comprising a high-performance curvature adaptive optics system with a simultaneous dual-channel coronagraphic imager. Combined with state-of-the-art observing methods and data processing, NICI typically achieves ~2 magnitudes better contrast compared to previous ground-based or space-based programs, at separations inside of ~2 arcsec. In preparation for the Campaign, we carried out efforts to identify previously unrecognized young stars, to rigorously construct our observing strategy, and to optimize the combination of angular and spectral differential imaging. The Planet-Finding Campaign is in its second year, with first-epoch imaging of 174 stars already obtained out of a total sample of 300 stars. We describe the Campaign's goals, design, implementation, performance, and preliminary results. The NICI Campaign represents the largest and most sensitive imaging survey to date for massive (~1 Mjup) planets around other stars. Upon completion, the Campaign will establish the best measurements to date on the properties of young gas-giant planets at ~5-10 AU separations. Finally, Campaign discoveries will be well-suited to long-term orbital monitoring and detailed spectrophotometric followup with next-generation planet-finding instruments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Gemini NICI Planet-Finding Campaign: Discovery of a Close Substellar Companion to the Young Debris Disk Star PZ Tel

We report the discovery of a tight substellar companion to the young solar analog PZ Tel, a member of the Beta Pictoris moving group observed with high contrast adaptive optics imaging as part of the Gemini NICI Planet-Finding Campaign. The companion was detected at a projected separation of 16.4 +/- 1.0 AU (0.33 +/- 0.01") in April 2009. Second-epoch observations in May 2010 demonstrate that the companion is physically associated and shows significant orbital motion. Monte Carlo modeling constrains the orbit of PZ Tel B to eccentricities > 0.6. The near-IR colors of PZ Tel B indicate a spectral type of M7+/-2 and thus this object will be a new benchmark companion for studies of ultracool, low-gravity photospheres. Adopting an age of 12 +8 -4 Myr for the system, we estimate a mass of 36 +/- 6 Mjup based on the Lyon/DUSTY evolutionary models. PZ Tel B is one of few young substellar companions directly imaged at orbital separations similar to those of giant planets in our own solar system. Additionally, the primary star PZ Tel A shows a 70 um emission excess, evidence for a significant quantity of circumstellar dust that has not been disrupted by the orbital motion of the companion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermal Emission and Tidal Heating of the Heavy and Eccentric Planet XO-3b

We determined the flux ratios of the heavy and eccentric planet XO-3b to its parent star in the four IRAC bands of the Spitzer Space Telescope: 0.101% +- 0.004% at 3.6 micron; 0.143% +- 0.006% at 4.5 micron; 0.134% +- 0.049% at 5.8 micron and 0.150% +- 0.036% at 8.0 micron. The flux ratios are within [-2.2,0.3, -0.8, -1.7]-sigma of the model of XO-3b with a thermally inverted stratosphere in the 3.6 micron, 4.5 micron, 5.8 micron and 8.0 micron channels, respectively. XO-3b has a high illumination from its parent star (Fp ~(1.9 - 4.2) x 10^9 ergs cm^-2 s^-1) and is thus expected to have a thermal inversion, which we indeed observe. When combined with existing data for other planets, the correlation between the presence of an atmospheric temperature inversion and the substellar flux is insufficient to explain why some high insolation planets like TrES-3 do not have stratospheric inversions and some low insolation planets like XO-1b do have inversions. Secondary factors such as sulfur chemistry, atmospheric metallicity, amounts of macroscopic mixing in the stratosphere or even dynamical weather effects likely play a role. Using the secondary eclipse timing centroids we determined the orbital eccentricity of XO-3b as e = 0.277 +- 0.009. The model radius-age trajectories for XO-3b imply that at least some amount of tidal-heating is required to inflate the radius of XO-3b, and the tidal heating parameter of the planet is constrained to Qp < 10^6 .

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tidal Heating Models for the Radii of the Inflated Transiting Giant Planets WASP-4b, WASP-6b, WASP-12b, and TrES-4

In order to explain the inflated radii of some transiting extrasolar giant planets, we investigate a tidal heating scenario for the inflated planets WASP-4b, WASP-6b, WASP-12b, WASP-15b, and TrES-4. To do so, we assume that they retain a nonzero eccentricity, possibly by dint of continuing interaction with a third body. We calculate the amount of extra heating in the envelope that is then required to fit the radius of each planet, and we explore how this additional power depends on the planetary atmospheric opacity and on the mass of a heavy-element central core. There is a degeneracy between the core mass $M_{\rm core}$ and the heating $\dot{E}_{\rm heating}$. Therefore, in the case of tidal heating, there is for each planet a range of the couple $\{M_{\rm core},e^2/Q'_p\}$ that can lead to the same radius, where $Q'_p$ is the tidal dissipation factor and $e$ is the eccentricity. With this in mind, we also investigate the case of the non-inflated planet HAT-P-12b, which can admit solutions combining a heavy-element core and tidal heating. A substantial improvement of the measured eccentricities of such planetary systems could simplify this degeneracy by linking the two unknown parameters $\{M_{\rm core},Q'_p\}$. Further independent constraints on either of these parameters would, through our calculations, constrain the other.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Warm Spitzer Photometry of the Transiting Exoplanets CoRoT-1 and CoRoT-2 at Secondary Eclipse

We measure secondary eclipses of the hot giant exoplanets CoRoT-1 at 3.6 and 4.5 microns, and CoRoT-2 at 3.6 microns, both using Warm Spitzer. We find that the Warm Spitzer mission is working very well for exoplanet science. For consistency of our analysis we also re-analyze archival cryogenic Spitzer data for secondary eclipses of CoRoT-2 at 4.5 and 8 microns. We compare the total data for both planets, including optical eclipse measurements by the CoRoT mission, and ground-based eclipse measurements at 2 microns, to existing models. Both planets exhibit stronger eclipses at 4.5 than at 3.6 microns, which is often indicative of an atmospheric temperature inversion. The spectrum of CoRoT-1 is best reproduced by a 2460K blackbody, due either to a high altitude layer that strongly absorbs stellar irradiance, or an isothermal region in the planetary atmosphere. The spectrum of CoRoT-2 is unusual because the 8 micron contrast is anomalously low. Non-inverted atmospheres could potentially produce the CoRoT-2 spectrum if the planet exhibits line emission from CO at 4.5 microns, caused by tidal-induced mass loss. However, the viability of that hypothesis is questionable because the emitting region cannot be more than about 30-percent larger than the planet's transit radius, based on the ingress and egress times at eclipse. An alternative possibility to account for the spectrum of CoRoT-2 is an additional opacity source that acts strongly at wavelengths less than 5 microns, heating the upper atmosphere while allowing the deeper atmosphere seen at 8 microns to remain cooler. We obtain a similar result as Gillon et al. for the phase of the secondary eclipse of CoRoT-2, implying an eccentric orbit with e*cos(omega)=-0.0030 +/- 0.0004.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Model for Gravitational Wave Emission from Neutrino-Driven Core-Collapse Supernovae

Abridged: Using a suite of progenitor models, neutrino luminosities, and two- dimensional (2D) simulations, we investigate the matter gravitational-wave (GW) emission from postbounce phases of neutrino-driven core-collapse supernovae (CCSNe). The relevant phases are prompt and steady-state convection, the standing accretion shock instability (SASI), and asymmetric explosions. For the stages before explosion, we propose a model for the source of GW emission. Downdrafts of the postshock-convection/SASI region strike the protoneutron star "surface" with large speeds and are decelerated by buoyancy forces. We find that the GW amplitude is set by the magnitude of deceleration and, by extension, the downdraft's speed and the vigor of postshock-convective/SASI motions. However, the characteristic frequencies, which evolve from ~100 Hz after bounce to ~300-400 Hz, are primarily independent of these speeds, but are set by the deceleration timescale, which is in turn set by the buoyancy frequency at the lower boundary of postshock convection. Consequently, the GW characteristic frequencies are dependent upon a combination of core structure attributes, specifically the dense-matter equation of state (EOS) and details that determine the gradients at the boundary, including the accretion-rate history, the EOS at subnuclear densities, and neutrino transport. During explosion, the high frequency signal wanes and is replaced by a strong low frequency, ~10s of Hz, signal that reveals the general morphology of the explosion (i.e. prolate, oblate, or spherical). However, current and near-future GW detectors are sensitive to GW power at frequencies >50 Hz. Therefore, the signature of explosion will be the abrupt reduction of detectable GW emission.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Detection of Thermal Emission of XO-2b: Evidence for a Weak Temperature Inversion

We estimate flux ratios of the extrasolar planet XO-2b to its host star XO-2 at 3.6, 4.5, 5.8 and 8.0 micron with IRAC on the Spitzer Space Telescope to be 0.00081 +- 0.00017, 0.00098 +- 0.00020, 0.00167 +- 0.00036 and 0.00133 +- 0.00049, respectively. The fluxes provide tentative evidence for a weak temperature inversion in the upper atmosphere, the precise nature of which would need to be confirmed by longer wavelength observations. XO-2b substellar flux of 0.76 x 10^9 ergs cm^-2 s^-1 lies in the predicted transition region between atmospheres with and without upper atmospheric temperature inversion.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Explorations into the Viability of Coupled Radius-Orbit Evolutionary Models for Inflated Planets

The radii of some transiting extrasolar giant planets are larger than would be expected by the standard theory. We address this puzzle with the model of coupled radius-orbit tidal evolution developed by \citet{Ibgui_and_Burrows_2009}. The planetary radius is evolved self-consistently with orbital parameters, under the influence of tidal torques and tidal dissipation in the interior of the planet. A general feature of this model, which we have previously demonstrated in the generic case, is that a possible transient inflation of the planetary radius can temporarily interrupt its standard monotonic shrinking and can lead to the inflated radii that we observe. In particular, a bloated planet with even a circular orbit may still be inflated due to an earlier episode of tidal heating. We have modified our model to include an orbital period dependence of the tidal dissipation factor in the star, $Q'_{\ast} \propto P^γ$, $-1 \leqslant γ\leqslant 1$. With this model, we search, for a tidally heated planet, orbital and radius evolutionary tracks that fall within the observational limits of the radius, the semimajor axis, and the eccentricity of the planet in its current estimated age range. We find that, for some inflated planets (WASP-6b and WASP-15b), there are such tracks; for another (TrES-4), there are none; and for still others (WASP-4b and WASP-12b), there are such tracks, but our model might imply that we are observing the planets at a special time. Finally, we stress that there is a two to three order-of-magnitude timescale uncertainty of the inspiraling phase of the planet into its host star, arising from uncertainties in the tidal dissipation factor in the star $Q'_{\ast}$.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Ks-band detection of thermal emission and color constraints to CoRoT-1b: A low-albedo planet with inefficient atmospheric energy redistribution and a temperature inversion

We report the detection in Ks-band of the secondary eclipse of the hot Jupiter CoRoT-1b, from time series photometry with the ARC 3.5-m telescope at Apache Point Observatory. The eclipse shows a depth of 0.336+/-0.042 percent and is centered at phase 0.5022 (+0.0023,-0.0027), consistent with a zero eccentricity orbit ecosω = 0.0035 (+0.0036,-0.0042). We perform the first optical to near-infrared multi-band photometric analysis of an exoplanet's atmosphere and constrain the reflected and thermal emissions by combining our result with the recent 0.6, 0.71, and 2.09 micron secondary eclipse detections by Snellen et al. (2009), Gillon et al. (2009), and Alonso et al. (2009a). Comparing the multi-wavelength detections to state-of-the-art radiative-convective chemical-equilibrium atmosphere models, we find the near-infrared fluxes difficult to reproduce. The closest blackbody-based and physical models provide the following atmosphere parameters: a temperature T = 2454 (+84,-170) K, a very low Bond albedo A_B = 0.000 (+0.087,-0.000), and an energy redistribution parameter P_n = 0.1, indicating a small but nonzero amount of heat transfer from the day- to night-side. The best physical model suggests a thermal inversion layer with an extra optical absorber of opacity kappa_e =0.05cm^2g^-1, placed near the 0.1-bar atmospheric pressure level. This inversion layer is located ten times deeper in the atmosphere than the absorbers used in models to fit mid-infrared Spitzer detections of other irradiated hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Search for very low-mass brown dwarfs and free-floating planetary-mass objects in Taurus

The number of low-mass brown dwarfs and even free floating planetary mass objects in young nearby star-forming regions and associations is continuously increasing, offering the possibility to study the low-mass end of the IMF in greater detail. In this paper, we present six new candidates for (very) low-mass objects in the Taurus star-forming region one of which was recently discovered in parallel by Luhman et al. (2009). The underlying data we use is part of a new database from a deep near-infrared survey at the Calar Alto observatory. The survey is more than four magnitudes deeper than the 2MASS survey and covers currently ~1.5 square degree. Complementary optical photometry from SDSS were available for roughly 1.0 square degree. After selection of the candidates using different color indices, additional photometry from Spitzer/IRAC was included in the analysis. In greater detail we focus on two very faint objects for which we obtained J-band spectra. Based on comparison with reference spectra we derive a spectral type of L2+/-0.5 for one object, making it the object with the latest spectral type in Taurus known today. From models we find the effective temperature to be 2080+/-140 K and the mass 5-15 Jupiter masses. For the second source the J-band spectrum does not provide a definite proof of the young, low-mass nature of the object as the expected steep water vapor absorption at 1.33 micron is not present in the data. We discuss the probability that this object might be a background giant or carbon star. If it were a young Taurus member, however, a comparison to theoretical models suggests that it lies close to or even below the deuterium burning limit (<13 Jupiter masses) as well. A first proper motion analysis for both objects shows that they are good candidates for being Taurus members.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Testing the models: NIR imaging and spectroscopy of the benchmark T-dwarf binary Eps Indi B

The relative roles of metallicity and surface gravity on the near-infrared spectra of late-T brown dwarfs are not yet fully understood, and evolutionary models still need to be calibrated in order to provide accurate estimates of brown dwarf physical parameters from measured spectra. The T-type brown dwarfs Eps Indi Ba and Bb forming the tightly bound binary Eps Indi B, which orbits the K4V star Eps Indi A, are nowadays the only such benchmark T dwarfs for which all important physical parameters such as metallicity, age and mass are (or soon will be) known. We present spatially resolved VLT/NACO images and low resolution spectra of Eps Indi B in the J, H and K near-infrared bands. The spectral types of Eps Indi Ba and Bb are determined by direct comparison of the flux-calibrated JHK spectra with T dwarf standard template spectra and also by NIR spectral indices. Eps Indi Bb is confirmed as a T6 while the spectral type of Eps Indi Ba is T1.5 so somewhat later than the previously reported T1. Constrained values for surface gravity and effective temperature are derived by comparison with model spectra. The evolutionary models predict masses around about 53 M_J for Eps Indi Ba and about 34 M_J for Eps Indi Bb, slightly higher than previously reported values. The suppressed J-band and enhanced K-band flux of Eps Indi Ba indicates that a noticeable cloud layer is still present in a T1.5 dwarf while no clouds are needed to model the spectrum of Eps Indi Bb.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Broadband Infrared Emission Spectrum of the Exoplanet TrES-3

We use the Spitzer Space Telescope to estimate the dayside thermal emission of the exoplanet TrES-3 integrated in the 3.6, 4.5, 5.8, and 8.0 micron bandpasses of the Infrared Array Camera (IRAC) instrument. We observe two secondary eclipses and find relative eclipse depths of 0.00346 +/- 0.00035, 0.00372 +/- 0.00054, 0.00449 +/- 0.00097, and 0.00475 +/- 0.00046, respectively in the 4 IRAC bandpasses. We combine our results with the earlier K band measurement of De Mooij et al. (2009), and compare them with models of the planetary emission. We find that the planet does not require the presence of an inversion layer in the high atmosphere. This is the first very strongly irradiated planet that does not have a temperature inversion, which indicates that stellar or planetary characteristics other than temperature have an important impact on temperature inversion. De Mooij & Snellen (2009) also detected a possible slight offset in the timing of the secondary eclipse in K band. However, based on our 4 Spitzer channels, we place a 3sigma upper limit of |ecos(w)| < 0.0056 where e is the planets orbital eccentricity and w is the longitude of the periastron. This result strongly indicates that the orbit is circular, as expected from tidal circularization theory.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

2D Multi-Angle, Multi-Group Neutrino Radiation-Hydrodynamic Simulations of Postbounce Supernova Cores

We perform axisymmetric (2D) multi-angle, multi-group neutrino radiation-hydrodynamic calculations of the postbounce phase of core-collapse supernovae using a genuinely 2D discrete-ordinate (S_n) method. We follow the long-term postbounce evolution of the cores of one nonrotating and one rapidly-rotating 20-solar-mass stellar model for ~400 milliseconds from 160 ms to ~550 ms after bounce. We present a multi-D analysis of the multi-angle neutrino radiation fields and compare in detail with counterpart simulations carried out in the 2D multi-group flux-limited diffusion (MGFLD) approximation to neutrino transport. We find that 2D multi-angle transport is superior in capturing the global and local radiation-field variations associated with rotation-induced and SASI-induced aspherical hydrodynamic configurations. In the rotating model, multi-angle transport predicts much larger asymptotic neutrino flux asymmetries with pole to equator ratios of up to ~2.5, while MGFLD tends to sphericize the radiation fields already in the optically semi-transparent postshock regions. Along the poles, the multi-angle calculation predicts a dramatic enhancement of the neutrino heating by up to a factor of 3, which alters the postbounce evolution and results in greater polar shock radii and an earlier onset of the initially rotationally weakened SASI. In the nonrotating model, differences between multi-angle and MGFLD calculations remain small at early times when the postshock region does not depart significantly from spherical symmetry. At later times, however, the growing SASI leads to large-scale asymmetries and the multi-angle calculation predicts up to 30% higher average integral neutrino energy deposition rates than MGFLD.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Detection of a Temperature Inversion in the Broadband Infrared Emission Spectrum of TrES-4

We estimate the strength of the bandpass-integrated thermal emission from the extrasolar planet TrES-4 at 3.6, 4.5, 5.8, and 8.0 micron using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We find relative eclipse depths of 0.137 +/- 0.011%, 0.148 +/- 0.016%, 0.261 +/- 0.059%, and 0.318 +/- 0.044% in these four bandpasses, respectively. We also place a 2 sigma upper limit of 0.37% on the depth of the secondary eclipse in the 16 micron IRS peak-up array. These eclipse depths reveal that TrES-4 has an emission spectrum similar to that of HD 209458b, which requires the presence of water emission bands created by an thermal inversion layer high in the atmosphere in order to explain the observed features. TrES-4 receives more radiation from its star than HD 209458b and has a correspondingly higher effective temperature, therefore the presence of a temperature inversion in this planet's atmosphere lends support to the idea that inversions might be correlated with the irradiance received by the planet. We find no evidence for any offset in the timing of the secondary eclipse, and place a 3 sigma upper limit of |ecos(omega)|<0.0058 where e is the planet's orbital eccentricity and omega is the argument of pericenter. From this we conclude that tidal heating from ongoing orbital circulatization is unlikely to be the explanation for TrES-4's inflated radius.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Spitzer Planet Limits around the Pulsating White Dwarf GD66

We present infrared observations in search of a planet around the white dwarf, GD66. Time-series photometry of GD66 shows a variation in the arrival time of stellar pulsations consistent with the presence of a planet with mass > 2.4Mj. Any such planet is too close to the star to be resolved, but the planet's light can be directly detected as an excess flux at 4.5um. We observed GD66 with the two shorter wavelength channels of IRAC on Spitzer but did not find strong evidence of a companion, placing an upper limit of 5--7Mj on the mass of the companion, assuming an age of 1.2--1.7Gyr.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Maximum Entropy for Gravitational Wave Data Analysis: Inferring the Physical Parameters of Core-Collapse Supernovae

The gravitational wave signal arising from the collapsing iron core of a Type II supernova progenitor star carries with it the imprint of the progenitor's mass, rotation rate, degree of differential rotation, and the bounce depth. Here, we show how to infer the gravitational radiation waveform of a core collapse event from noisy observations in a network of two or more LIGO-like gravitational wave detectors and, from the recovered signal, constrain these source properties. Using these techniques, predictions from recent core collapse modeling efforts, and the LIGO performance during its S4 science run, we also show that gravitational wave observations by LIGO might have been sufficient to provide reasonable estimates of the progenitor mass, angular momentum and differential angular momentum, and depth of the core at bounce, for a rotating core collapse event at a distance of a few kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Simulations of Magnetically-Driven Supernova and Hypernova Explosions in the Context of Rapid Rotation

We present here the first 2D rotating, multi-group, radiation magnetohydrodynamics (RMHD) simulations of supernova core collapse, bounce, and explosion. In the context of rapid rotation, we focus on the dynamical effects of magnetic stresses and the creation and propagation of MHD jets. We find that a quasi-steady state can be quickly established after bounce, during which a well-collimated MHD jet is maintained by continuous pumping of power from the differentially rotating core. If the initial spin period of the progenitor core is $\sles$ 2 seconds, the free energy reservoir in the secularly evolving protoneutron star is adequate to power a supernova explosion, and may be enough for a hypernova. The jets are well collimated by the infalling material and magnetic hoop stresses, and maintain a small opening angle. We see evidence of sausage instabilities in the emerging jet stream. Neutrino heating is sub-dominant in the rapidly rotating models we explore, but can contribute 10$-$25% to the final explosion energy. Our simulations suggest that even in the case of modest or slow rotation, a supernova explosion might be followed by a secondary, weak MHD jet explosion, which, because of its weakness may to date have gone unnoticed in supernova debris. Furthermore, we suggest that the generation of a non-relativistic MHD precursor jet during the early protoneutron star/supernova phase is implicit in both the collapsar and "millisecond magnetar" models of GRBs. The multi-D, multi-group, rapidly rotating RMHD simulations we describe here are a start along the path towards more realistic simulations of the possible role of magnetic fields in some of Nature's most dramatic events.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The 3.6-8.0 Micron Broadband Emission Spectrum of HD 209458b: Evidence for an Atmospheric Temperature Inversion

We estimate the strength of the bandpass-integrated thermal emission from the extrasolar planet HD 209458b at 3.6, 4.5, 5.8, and 8.0 microns using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We observe a single secondary eclipse simultaneously in all four bandpasses and find relative eclipse depths of 0.00094 +/- 0.00009, 0.00213 +/- 0.00015, 0.00301 +/- 0.00043, and 0.00240 +/- 0.00026, respectively. These eclipse depths reveal that the shape of the inferred emission spectrum for the planet differs significantly from the predictions of standard atmosphere models; instead the most plausible explanation would require the presence of an inversion layer high in the atmosphere leading to significant water emission in the 4.5 and 5.8 micron bandpasses. This is the first clear indication of such a temperature inversion in the atmosphere of a hot Jupiter, as previous observations of other planets appeared to be in reasonably good agreement with the predictions of models without such an inversion layer.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A New Algorithm for 2-D Transport for Astrophysical Simulations: I. General Formulation and Tests for the 1-D Spherical Case

We derive new equations using the mixed-frame approach for one- and two-dimensional (axisymmetric) time-dependent radiation transport and the associated couplings with matter. Our formulation is multi-group and multi-angle and includes anisotropic scattering, frequency(energy)-dependent scattering and absorption, complete velocity dependence to order v/c, rotation, and energy redistribution due to inelastic scattering. Hence, the "2D" realization is actually "6 1/2"-dimensional. The effects of radiation viscosity are automatically incorporated. Moreover, we develop Accelerated-Lambda-Iteration, Krylov subspace (GMRES), Discontinuous-Finite-Element, and Feautrier numerical methods for solving the equations and present the results of one-dimensional numerical tests of the new formalism. The virtues of the mixed-frame approach include simple velocity dependence with no velocity derivatives, straight characteristics, simple physical interpretation, and clear generalization to higher dimensions. Our treatment can be used for both photon and neutrino transport, but we focus on neutrino transport and applications to core-collapse supernova theory in the discussions and examples.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Telescope to Observe Planetary Systems (TOPS): a high throughput 1.2-m visible telescope with a small inner working angle

The Telescope to Observe Planetary Systems (TOPS) is a proposed space mission to image in the visible (0.4-0.9 micron) planetary systems of nearby stars simultaneously in 16 spectral bands (resolution R~20). For the ~10 most favorable stars, it will have the sensitivity to discover 2 R_E rocky planets within habitable zones and characterize their surfaces or atmospheres through spectrophotometry. Many more massive planets and debris discs will be imaged and characterized for the first time. With a 1.2m visible telescope, the proposed mission achieves its power by exploiting the most efficient and robust coronagraphic and wavefront control techniques. The Phase-Induced Amplitude Apodization (PIAA) coronagraph used by TOPS allows planet detection at 2 lambda/d with nearly 100% throughput and preserves the telescope angular resolution. An efficient focal plane wavefront sensing scheme accurately measures wavefront aberrations which are fed back to the telescope active primary mirror. Fine wavefront control is also performed independently in each of 4 spectral channels, resulting in a system that is robust to wavefront chromaticity.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Spin Periods and Rotational Profiles of Neutron Stars at Birth

We present results from an extensive set of one- and two-dimensional radiation-hydrodynamic simulations of the supernova core collapse, bounce, and postbounce phases, and focus on the protoneutron star (PNS) spin periods and rotational profiles as a function of initial iron core angular velocity, degree of differential rotation, and progenitor mass. For the models considered, we find a roughly linear mapping between initial iron core rotation rate and PNS spin. The results indicate that the magnitude of the precollapse iron core angular velocities is the single most important factor in determining the PNS spin. Differences in progenitor mass and degree of differential rotation lead only to small variations in the PNS rotational period and profile. Based on our calculated PNS spins, at ~ 200-300 milliseconds after bounce, and assuming angular momentum conservation, we estimate final neutron star rotation periods. We find periods of one millisecond and shorter for initial central iron core periods of below ~ 10 s. This is appreciably shorter than what previous studies have predicted and is in disagreement with current observational data from pulsar astronomy. After considering possible spindown mechanisms that could lead to longer periods we conclude that there is no mechanism that can robustly spin down a neutron star from ~ 1 ms periods to the "injection" periods of tens to hundreds of milliseconds observed for young pulsars. Our results indicate that, given current knowledge of the limitations of neutron star spindown mechanisms, precollapse iron cores must rotate with periods around 50-100 seconds to form neutron stars with periods generically near those inferred for the radio pulsar population.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Line Intensities and Molecular Opacities of the FeH $F^4Δ_i-X^4Δ_i$ Transition

We calculate new line lists and opacities for the $F^4Δ_i-X^4Δ_i$ transition of FeH. The 0-0 band of this transition is responsible for the Wing-Ford band seen in M-type stars, sunspots and brown dwarfs. The new Einstein A values for each line are based on a high level ab initio calculation of the electronic transition dipole moment. The necessary rotational line strength factors (Hönl-London factors) are derived for both the Hund's case (a) and (b) coupling limits. A new set of spectroscopic constants were derived from the existing FeH term values for v=0, 1 and 2 levels of the $X$ and $F$ states. Using these constants extrapolated term values were generated for v=3 and 4 and for $J$ values up to 50.5. The line lists (including Einstein A values) for the 25 vibrational bands with v$\leq$4 were generated using a merged list of experimental and extrapolated term values. The FeH line lists were use to compute the molecular opacities for a range of temperatures and pressures encountered in L and M dwarf atmospheres. Good agreement was found between the computed and observed spectral energy distribution of the L5 dwarf 2MASS-1507.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The Spectra of T Dwarfs. II. Red Optical Data

We present 6300--10100 Å spectra for a sample of 13 T dwarfs observed using LRIS mounted on the Keck I 10m Telescope. A variety of features are identified and analyzed, including pressure-broadened K I and Na I doublets; narrow Cs I and Rb I lines; weak CaH, CrH, and FeH bands; strong H$_2$O absorption; and a possible weak CH$_4$ band. H$α$ emission is detected in three of the T dwarfs, strong in the previously reported active T dwarf 2MASS 1237+6526 and weak in SDSS 1254$-$0122 and 2MASS 1047+2124. None of the T dwarfs exhibit Li I absorption. Guided by the evolution of optical spectral features with near-infrared spectral type, we derive a parallel optical classification scheme, focusing on spectral types T5 to T8, anchored to select spectral standards. We find general agreement between optical and near-infrared types for nearly all of the T dwarfs so far observed, including two earlier-type T dwarfs, within our classification uncertainties ($\sim$1 subtype). These results suggest that competing gravity and temperature effects compensate for each other over the 0.6--2.5 $\micron$ spectral region. We identify one possible means of disentangling these effects by comparing the strength of the K I red wing to the 9250 Å H$_2$O band. One of our objects, 2MASS 0937+2931, exhibits a peculiar spectrum, with a substantial red slope and relatively strong FeH absorption, both consequences of a metal-deficient atmosphere. Based on its near-infrared properties and substantial space motion, this object may be a thick disk or halo brown dwarf.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The Physics of Protoneutron Star Winds: Implications for r-Process Nucleosynthesis

We solve the general-relativistic steady-state eigenvalue problem of neutrino-driven protoneutron star winds, which immediately follow core-collapse supernova explosions. We provide velocity, density, temperature, and composition profiles and explore the systematics and structures generic to such a wind for a variety of protoneutron star characteristics. Furthermore, we derive the entropy, dynamical timescale, and neutron-to-seed ratio in the general relativistic framework essential in assessing this site as a candidate for $r$-process nucleosynthesis. Generally, we find that for a given mass outflow rate ($\dot{M}$), the dynamical timescale of the wind is significantly shorter than previously thought. We argue against the existence or viability of a high entropy ($\gtrsim300$ per k$_{B}$ per baryon), long dynamical timescale $r$-process epoch. In support of this conclusion, we model the protoneutron star cooling phase, calculate nucleosynthetic yields in our steady-state profiles, and estimate the integrated mass loss. We find that transonic winds enter a high entropy phase only with very low $\dot{M}$ ($\lesssim1\times10^{-9}$ M$_\odot$ s$^{-1}$) and extremely long dynamical timescale ($τ_ρ\gtrsim0.5$ seconds). Our results support the possible existence of an early $r$-process epoch at modest entropy ($\sim150$) and very short dynamical timescale, consistent in our calculations with a very massive or very compact protoneutron star that contracts rapidly after the preceding supernova. We explore possible modifications to our models, which might yield significant $r$-process nucleosynthesis generically. Finally, we speculate on the effect of fallback and shocks on both the wind physics and nucleosynthesis.

0 citations0 reviewsNo code linkedOpen access