Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star (PNS) convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral Standing Accretion Shock Instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate ($\gtrsim 50 \%$), but that the time-integrated signal has a relatively modest asymmetry ($\lesssim 20 \%$). Both features are associated with the lepton-number emission self-sustained asymmetry (LESA) and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavors of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of $\sim$20\% if the distance to the CCSN is $\lesssim 6$ kpc.

JWST will provide moderate resolution transit spectra with continuous wavelength coverage from the optical to the mid-infrared for the first time. In this paper, we illustrate how different aerosol species, size-distributions, and spatial distributions encode information in JWST transit spectra. We use the transit spectral modeling code METIS, along with Mie theory and several flexible treatments of aerosol size and spatial distributions to perform parameter sensitivity studies, calculate transit contribution functions, compute Jacobians, and retrieve parameters with Markov Chain Monte Carlo. The broader wavelength coverage of JWST will likely encompass enough non-gray aerosol behavior to recover information about the species and size-distribution of particles, especially if distinct resonance features arising from the aerosols are present. Within the JWST wavelength range, the optical and mid-infrared typically provide information about 0.1-1 $\mu$m sized aerosols, while the near-infrared to mid-infrared wavelengths usually provide information about gaseous absorption, even if aerosols are present. Strong gaseous absorption features in the infrared often remain visible, even when clouds and hazes are flattening the optical and NIR portion of the spectrum that is currently observable. For some combinations of aerosol properties, temperature, and surface gravity, one can make a precise measure of metallicity despite the presence of aerosols, but more often the retrieved metallicity of a cloudy or hazy atmosphere has significantly lower precision than for a clear atmosphere with otherwise similar properties. Future efforts to securely link aerosol properties to atmospheric metallicity and temperature in a physically motivated manner will ultimately enable a robust physical understanding of the processes at play in cloudy, hazy exoplanet atmospheres.

Based on the Kompaneets approximation, we develop a robust methodology to calculate spectral redistribution via inelastic neutrino-nucleon scattering in the context of core-collapse supernova simulations. The resulting equations conserve lepton number to machine precision and scale linearly, not quadratically, with number of energy groups. The formalism also provides an elegant means to derive the rate of energy transfer to matter which, as it must, automatically goes to zero when the neutrino radiation field is in thermal equilibrium. Furthermore, we derive the next-higher-order in {\epsilon}/mc2 correction to the neutrino Kompaneets equation. Unlike other Kompaneets schema, ours also generalizes to the case of anisotropic angular distributions, while retaining the conservative form that is a hallmark of the classical Kompaneets equation. Our formalism enables immediate incorporation into supernova codes that follow the spectral angular moments of the neutrino radiation fields.

Many of the exoplanets for which we can obtain the highest SNR transit spectra are tidally locked. The atmospheres on tidally-locked planets likely exhibit large differences between the day and night side of the planet, the poles, and the morning versus evening terminators. In this paper, we illustrate how the combined effects of aerosols and day-night temperature gradients shape transit spectra of tidally-locked exoplanets when full 3D structures are taken into account and evaluate the implications for retrievals of atmospheric properties. To do this, we have developed a new code, METIS, which can compute transit spectra for an arbitrary longitude-latitude-altitude grid of temperatures and pressures. Using METIS, we pair flexible treatments of clouds and hazes with simple parameterized day-night temperature gradients to compute transit spectra and perform retrieval experiments across a wide array of possible exoplanet atmospheric properties. Our key findings are that: (1) the presence of aerosols typically increases the effects of day-night temperature gradients on transit spectra; (2) ignoring day-night temperature gradients when attempting to perform Bayesian parameter estimation will still return biased results when aerosols are present, as has already been shown for clear atmospheres in the literature; (3) when a day-night temperature gradient is present and accounted for in the retrieval model, some transit spectra can provide sufficient information to constrain temperatures and the width of the transition from day to night. The presence of clouds and hazes can actually tighten such constraints, but also weaken constraints on metallicity. This paper represents a step towards the larger goal of developing models and theory of adequate complexity to match the superior quality data that will soon be available.

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

Background: Dataset skills are used in STEM fields from healthcare work to astronomy research. Few fields explicitly teach students the skills to analyze datasets, and yet the increasing push for authentic science implies these skills should be taught. Purpose: The overarching motivation is to understand learning of dataset skills within an astronomy context. Specifically, when participants work with a 200-entry Google Sheets dataset of astronomical data about quasars, what are they learning, how are they learning it, and who is doing the learning? Sample: The authors studied a matched set of participants (n=87) consisting of 54 university undergraduate students (34 male, 18 female), and 33 science educators (16 male, 17 female). Design and methods: Participants explored a three-phase dataset activity and were given an eight-question multiple-choice pre/post-test covering skills of analyzing datasets and astronomy content, with questions spanning Bloom's Taxonomy. Pre/post-test scores were compared and a t-test performed for subsamples by population. Results: Participants exhibited learning of both dataset skills and astronomy content, indicating that dataset skills can be learned through this astronomy activity. Participants exhibited gains in both recall and synthesis questions, indicating learning is non-sequential. Female undergraduate students exhibited lower levels of learning than other populations. Conclusions: Implications of the study include a stronger dataset focus in post-secondary STEM education and among science educators, and the need for further investigation into how instructors can ameliorate the challenges faced by female undergraduate students.

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_{\nu_e} - n_{\bar{\nu}_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.

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 $\nu_{\mu}$, $\nu_{\tau}$, $\bar{\nu}_{\mu}$, and $\bar{\nu}_{\tau}$ 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.

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.

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.

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.

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 \`a 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.

Using our new state-of-the-art core-collapse supernova (CCSN) code Fornax, we explore the dependence upon spatial resolution of the outcome and character of three-dimensional (3D) supernova simulations. For the same 19-M$_{\odot}$ progenitor star, energy and radial binning, neutrino microphysics, and nuclear equation of state, changing only the number of angular bins in the $\theta$ and $\phi$ directions, we witness that our lowest resolution 3D simulation does not explode. However, when jumping progressively up in resolution by factors of two in each angular direction on our spherical-polar grid, models then explode, and explode slightly more vigorously with increasing resolution. This suggests that there can be a qualitative dependence of the outcome of 3D CCSN simulations upon spatial resolution. The critical aspect of higher spatial resolution is the adequate capturing of the physics of neutrino-driven turbulence, in particular its Reynolds stress. The greater numerical viscosity of lower-resolution simulations results in greater drag on the turbulent eddies that embody turbulent stress, and, hence, in a diminution of their vigor. Turbulent stress not only pushes the temporarily stalled shock further out, but bootstraps a concomitant increase in the deposited neutrino power. Both effects together lie at the core of the resolution dependence we observe.

We provide the time series and angular distributions of the neutrino and gravitational-wave emissions of eleven state-of-the-art three-dimensional non-rotating core-collapse supernova models and explore correlations between these signatures and the real-time dynamics of the shock and the proto-neutron-star core. The neutrino emissions are roughly isotropic on average, with instantaneous excursions about the mean inferred luminosity of as much as $\pm$20%. The deviation from isotropy is least for the "$\nu_{\mu}$"-type neutrinos and the lowest-mass progenitors. Instantaneous temporal luminosity variations along a given direction for exploding models average $\sim$2$-$4%, but can be as high as $\sim$10%. For non-exploding models, they can achieve $\sim$25%. The temporal variations in the neutrino emissions correlate with the temporal and angular variations in the mass accretion rate. We witness the LESA phenomenon in all our models and find that the vector direction of the LESA dipole and that of the inner Y$_\mathrm{e}$ distribution are highly correlated. For our entire set of 3D models, we find strong connections between the cumulative neutrino energy losses, the radius of the proto-neutron star, and the $f$-mode frequency of the gravitational wave emissions. When physically normalized, the progenitor-to-progenitor variation in any of these quantities is no more than $\sim$10%. Moreover, the reduced $f$-mode frequency is independent of time after bounce to better than $\sim$10%. Therefore, simultaneous measurement of gravitational waves and neutrinos from a given supernova event can be used synergistically to extract real physical quantities of the supernova core.

LOX is a lunar-orbiting astrophysics mission that will probe the cosmos at MeV energies. It is guided by open questions regarding thermonuclear, or Type-Ia, supernovae (SNeIa) and will characterize these inherently radioactive objects by enabling a systematic survey of SNeIa at gamma-ray energies for the first time. Astronomical investigations from lunar orbit afford new opportunities to advance our understanding of the cosmos. The foundation of LOX is an observational approach well suited to the all-sky monitoring demands of supernova investigations and time-domain astronomy. Its inherently wide field-of-view and continuous all-sky monitoring provides an innovative way of addressing decadal survey questions at MeV energies (0.1-10 MeV). The LOX approach achieves high sensitivity with a simple, high-heritage instrument design that eliminates the need for complex, position-sensitive detectors, kinematic event reconstruction, masks, or other insensitive detector mass, while also mitigating technology development, implementation complexity, and their associated costs. LOX can be realized within existing programs, like Explorer.

We optimize the third-generation gravitational-wave detector to maximize the range to detect core-collapse supernovae. Based on three-dimensional simulations for core-collapse and the corresponding gravitational-wave waveform emitted, the corresponding detection range for these waveforms is limited to within our galaxy even in the era of third-generation detectors. The corresponding event rate is two per century. We find from the waveforms that to detect core-collapse supernovae with an event rate of one per year, the gravitational-wave detectors need a strain sensitivity of 3$\times10^{-27}~$Hz$^{-1/2}$ in a frequency range from 100~Hz to 1500~Hz. We also explore detector configurations technologically beyond the scope of third-generation detectors. We find with these improvements, the event rate for gravitational-wave observations from CCSN is still low, but is improved to one in twenty years.

The 2017 detection of the inspiral and merger of two neutron stars in gravitational waves and gamma rays was accompanied by a quickly-reddening transient. Such a transient was predicted to occur following a rapid neutron capture (r-process) nucleosynthesis event, which synthesizes neutron-rich, radioactive nuclei and can take place in both dynamical ejecta and in the wind driven off the accretion torus formed after a neutron star merger. We present the first three-dimensional general relativistic, full transport neutrino radiation magnetohydrodynamics (GRRMHD) simulations of the black hole-accretion disk-wind system produced by the GW170817 merger. We show that the small but non-negligible optical depths lead to neutrino transport globally coupling the disk electron fraction, which we capture by solving the transport equation with a Monte Carlo method. The resulting absorption drives up the electron fraction in a structured, continuous outflow, with electron fraction as high as $Y_e\sim 0.4$ in the extreme polar region. We show via nuclear reaction network and radiative transfer calculations that nucleosynthesis in the disk wind will produce a blue kilonova.

We study the gravitational wave signal from eight new 3D core-collapse supernova simulations. We show that the signal is dominated by $f$- and $g$-mode oscillations of the protoneutron star and its frequency evolution encodes the contraction rate of the latter, which, in turn, is known to depend on the star's mass, on the equation of state, and on transport properties in warm nuclear matter. A lower-frequency component of the signal, associated with the standing accretion shock instability, is found in only one of our models. Finally, we show that the energy radiated in gravitational waves is proportional to the amount of turbulent energy accreted by the protoneutron star.

WFIRST-CGI is a NASA technology demonstration mission that is charged with demonstrating key technologies for future exo-Earth imaging missions in space. In the process, it will obtain images and low-resolution spectra of a handful to a dozen extrasolar planets and possibly protoplanetary disks. Its unprecedented contrast levels in the optical will provide astronomers with their first direct look at mature, Jupiter sized planets at moderate separations. This paper addresses the question: what science can be done with such data? An analytic noise model, which is informed by the ongoing engineering developments, is used to compute maximum achievable signal-to-noise ratios and scientifically viable integration times for hypothetical star planet systems, as well as to investigate the constraining power of various combinations of WFIRST-CGI photometric and spectral observations. This work introduces two simple models for planetary geometric albedos, which are inspired largely by the solar system's gas giants. The first planet model is a hybrid Jupiter-Neptune model, which separately treats the short and long wavelengths where chromophores and methane dominate absorption, respectively. The second planet model fixes cloud and haze properties in CoolTLusty to match Jupiter's albedo spectrum, it then perturbs only the metallicity. MCMC retrievals performed on simulated observations are used to assess the precision with which planet model parameters can be measured subject to different exposure times and observing cases. Fit results for both models' parameterizations of geometric albedo spectra demonstrate that a rough indication of the metallicity or methane content should be possible for some WFIRST-CGI targets. We conclude that real observations will likely be able to differentiate between extreme cases using these models, but will lack the precision necessary to uncover subtle trends.

Observations with next-generation ground-based detectors further enhanced with multi-messenger (electromagnetic and neutrino) detections will allow us to probe new extreme astrophysics. Target sources included: core-collapse supernovae, continuous emission from isolated or accreting neutron stars, and bursts from magnetars and other pulsars.

Using the new state-of-the-art core-collapse supernova (CCSN) code F{\sc{ornax}}, we have simulated the three-dimensional dynamical evolution of the cores of 9-, 10-, 11-, 12-, and 13-M$_{\odot}$ stars from the onset of collapse. Stars from 8-M$_{\odot}$ to 13-M$_{\odot}$ constitute roughly 50% of all massive stars, so the explosive potential for this mass range is important to the overall theory of CCSNe. We find that the 9-, 10-, 11-, and 12-M$_{\odot}$ models explode in 3D easily, but that the 13-M$_{\odot}$ model does not. From these findings, and the fact that slightly more massive progenitors seem to explode \citep{vartanyan2019}, we suggest that there is a gap in explodability near 12-M$_{\odot}$ to 14-M$_{\odot}$ for non-rotating progenitor stars. Factors conducive to explosion are turbulence behind the stalled shock, energy transfer due to neutrino-matter absorption and neutrino-matter scattering, many-body corrections to the neutrino-nucleon scattering rate, and the presence of a sharp silicon-oxygen interface in the progenitor. Our 3D exploding models frequently have a dipolar structure, with the two asymmetrical exploding lobes separated by a pinched waist where matter temporarily continues to accrete. This process maintains the driving neutrino luminosty, while partially shunting matter out of the way of the expanding lobes, thereby modestly facilitating explosion. The morphology of all 3D explosions is characterized by multiple bubble structures with a range of low-order harmonic modes. Though much remains to be done in CCSN theory, these and other results in the literature suggest that, at least for these lower-mass progenitors, supernova theory is converging on a credible solution.

The lowest-mass stars, brown dwarfs and giant exoplanets span a minimum in the mass-radius relationship that probes the fundamental physics of extreme states of matter, magnetism, and fusion. This White Paper outlines scientific opportunities and the necessary resources for modeling and measuring the mass-radius relationship in this regime.

This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes. Instead of focusing exclusively on tests of neutrino radiation hydrodynamics relevant to the core-collapse supernova problem for which Fornax is primarily intended, we present here classical and rigorous demonstrations of code performance relevant to a broad range of multi-dimensional hydrodynamic and multi-group radiation hydrodynamic problems. Our code solves the comoving-frame radiation moment equations using the M1 closure, utilizes conservative high-order reconstruction, employs semi-explicit matter and radiation transport via a high-order time stepping scheme, and is suitable for application to a wide range of astrophysical problems. To this end, we first describe the philosophy, algorithms, and methodologies of Fornax and then perform numerous stringent code tests, that collectively and vigorously exercise the code, demonstrate the excellent numerical fidelity with which it captures the many physical effects of radiation hydrodynamics, and show excellent strong scaling well above 100k MPI tasks.

The histories of core-collapse supernova theory and of neutrino physics have paralleled one another for more than seventy years. Almost every development in neutrino physics necessitated modifications in supernova models. What has emerged is a complex and rich dynamical scenario for stellar death that is being progressively better tested by increasingly sophisiticated computer simulations. Though there is still much to learn about the agency and details of supernova explosions, whatever final theory emerges will have the neutrino at its core. I summarize in this brief contribution some of the salient developments in neutrino physics as they related to supernova theory, while avoiding any attempt to review the hundreds of pivotal papers that have pushed supernova theory forward. My goal has been merely to highlight the debt of supernova astrophysics to neutrino physics.

We observed eclipses of the transiting brown dwarf CWW 89Ab at 3.6um and 4.5um using Spitzer/IRAC. The CWW 89 binary system is a member of the $3.0\pm0.25$ Gyr-old open cluster Ruprecht 147, and is composed of a Sun-like primary and an early M-dwarf secondary separated by a projected distance of 25 AU. CWW 89Ab has a radius of $0.937\pm0.042$ RJ and a mass of $36.5\pm0.1$ MJ, and is on a 5.3 day orbit about CWW 89A with a non-zero eccentricity of $e=0.19$ (Curtis et al. 2016). We strongly detect the eclipses of CWW 89Ab in both Spitzer channels as $\delta_{3.6}=1147\pm213$ ppm and $\delta_{4.5}=1097\pm225$ ppm after correcting for the dilution from CWW 89B. After accounting for the irradiation that CWW 89Ab receives from its host star, these measurements imply that the brown dwarf has an internal luminosity of $\log(\mathrm{L_{bol}/\mathrm{L}_\odot})=-4.19\pm0.14$. This is 16 times, or $9.3\,\sigma$, higher than model predictions given the known mass, radius, and age of CWW 89Ab. As we discuss, this over-luminosity is not explainable by an inaccurate age determination, additional stellar heating, nor tidal heating. Instead, we suggest that the anomalous luminosity of CWW 89Ab is caused by a dayside temperature inversion -- though a significant error in the evolutionary models is also a possibility. Importantly, a temperature inversion would require a super-stellar C/O ratio in CWW 89Ab's atmosphere. If this is indeed the case, it implies that CWW 89Ab is a 36.5 MJ object that formed via core accretion processes. Finally, we use our measurement of CWW 89Ab's orbital eccentricity, improved via these observations, to constrain the tidal quality factors of the brown dwarf and the host star CWW 89A to be $Q_{BD}>10^{4.15}$ and $Q_*>10^{9}$, respectively.

In this paper, we present the results of our three-dimensional, multi-group, multi-neutrino-species radiation/hydrodynamic simulation using the state-of-the-art code F{\sc{ornax}} of the terminal dynamics of the core of a non-rotating 16-M$_{\odot}$ stellar progenitor. The calculation incorporates redistribution by inelastic scattering, a correction for the effect of many-body interactions on the neutrino-nucleon scattering rates, approximate general relativity (including the effects of gravitational redshifts), velocity-dependent frequency advection, and an implementation of initial perturbations in the progenitor core. The model explodes within $\sim$100 milliseconds of bounce (near when the silicon-oxygen interface is accreted through the temporarily-stalled shock) and by the end of the simulation (here, $\sim$677 milliseconds after bounce) is accumulating explosion energy at a rate of $\sim$2.5$\times$10$^{50}$ ergs s$^{-1}$. The supernova explosion resembles an asymmetrical multi-plume structure, with one hemisphere predominating. The gravitational mass of the residual proto-neutron star at $\sim$677 milliseconds is $\sim$1.42 M$_{\odot}$. Even at the end of the simulation, explosion in most of the solid angle is accompanied by some accretion in an annular fraction at the wasp-like waist of the debris field. The ejecta electron fraction (Y$_e$) is distributed from $\sim$0.48 to $\sim$0.56, with most of the ejecta mass proton-rich. This may have implications for supernova nucleosynthesis, and could have a bearing on the p- and $\nu$p-processes and on the site of the first peak of the r-process. The ejecta spatial distributions of both Y$_e$ and mass density are predominantly in wide-angle plumes and large-scale structures, but are nevertheless quite patchy.

We present observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Field Camera 3 in the near infrared. We engaged groups across the transiting exoplanet community to participate in the analysis of the data and present results from each. There is general agreement amongst all results that we find an H2O absorption feature with 3.5-4.0 sigma significance. However, the feature is muted in comparison to a clear atmosphere at solar composition. Although the detection of the water feature is robust, the reasons for the muting of this feature are ambiguous due to a degeneracy between clouds and composition. The data does not yield robust detections of any molecular species other than H2O. The group was motivated to perform an additional set of retrieval exercises to investigate an apparent bump in the spectrum at ~ 1.55 um. We explore possible disequilibrium chemistry and find this feature is consistent with super-solar HCN abundance but it is questionable if the required mixing ratio of HCN is chemically and physically plausible. The ultimate goal of this study is to vet WASP-63b as a potential community target to best demonstrate the capabilities and systematics of JWST instruments for transiting exoplanet science. In the case of WASP-63b, the presence of a detectable water feature indicates that WASP-63b remains a plausible target for ERS observations.

For a suite of fourteen core-collapse models during the dynamical first second after bounce, we calculate the detailed neutrino "light" curves expected in the underground neutrino observatories Super-Kamiokande, DUNE, JUNO, and IceCube. These results are given as a function of neutrino-oscillation modality (normal or inverted hierarchy) and progenitor mass (specifically, post-bounce accretion history), and illuminate the differences between the light curves for 1D (spherical) models that don't explode with the corresponding 2D (axisymmetric) models that do. We are able to identify clear signatures of explosion (or non-explosion), the post-bounce accretion phase, and the accretion of the silicon/oxygen interface. In addition, we are able to estimate the supernova detection ranges for various physical diagnostics and the distances out to which various temporal features embedded in the light curves might be discerned. We find that the progenitor mass density profile and supernova dynamics during the dynamical explosion stage should be identifiable for a supernova throughout most of the galaxy in all the facilities studied and that detection by any one of them, but in particular more than one in concert, will speak volumes about the internal dynamics of supernovae.

We present a comparison between several simulation codes designed to study the core-collapse supernova mechanism. We pay close attention to controlling the initial conditions and input physics in order to ensure a meaningful and informative comparison. Our goal is three-fold. First, we aim to demonstrate the current level of agreement between various groups studying the core-collapse supernova central engine. Second, we desire to form a strong basis for future simulation codes and methods to compare to. Lastly, we want this work to be a stepping stone for future work exploring more complex simulations of core-collapse supernovae, i.e., simulations in multiple dimensions and simulations with modern neutrino and nuclear physics. We compare the early (first ~500ms after core bounce) spherically-symmetric evolution of a 20 solar mass progenitor star from six different core-collapse supernovae codes: 3DnSNe-IDSA, AGILE-BOLTZTRAN, FLASH, F{\sc{ornax}}, GR1D, and PROMETHEUS-VERTEX. Given the diversity of neutrino transport and hydrodynamic methods employed, we find excellent agreement in many critical quantities, including the shock radius evolution and the amount of neutrino heating. Our results provide an excellent starting point from which to extend this comparison to higher dimensions and compare the development of hydrodynamic instabilities that are crucial to the supernova explosion mechanism, such as turbulence and convection.

We study gravitational waves (GWs) from a set of two-dimensional multi-group neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Our goal is to systematize the current knowledge about the post-bounce CCSN GW signal and recognize the templatable features that could be used by the ground-based laser interferometers. We demonstrate that starting from ~400ms after core bounce the dominant GW signal represents the fundamental quadrupole (l=2) oscillation mode (f-mode) of the proto-neutron star (PNS), which can be accurately reproduced by a linear perturbation analysis of the angle-averaged PNS profile. Before that, in the time interval between ~200 and ~400ms after bounce, the dominant mode has two radial nodes and represents a g-mode. We associate the high-frequency noise in the GW spectrograms above the main signal with p-modes, while below the dominant frequency there is a region with very little power. The collection of models presented here summarizes the dependence of the CCSN GW signal on the progenitor mass, equation of state, many-body corrections to the neutrino opacity, and rotation. Weak dependence of the dominant GW frequency on the progenitor mass motivates us to provide a simple fit for it as a function of time, which can be used as a prior when looking for CCSN candidates in the LIGO data.

We explore KLIP forward modeling spectral extraction on Gemini Planet Imager coronagraphic data of HR 8799, using PyKLIP and show algorithm stability with varying KLIP parameters. We report new and re-reduced spectrophotometry of HR 8799 c, d, and e in H & K bands. We discuss a strategy for choosing optimal KLIP PSF subtraction parameters by injecting simulated sources and recovering them over a range of parameters. The K1/K2 spectra for HR 8799 c and d are similar to previously published results from the same dataset. We also present a K band spectrum of HR 8799 e for the first time and show that our H-band spectra agree well with previously published spectra from the VLT/SPHERE instrument. We show that HR 8799 c and d show significant differences in their H & K spectra, but do not find any conclusive differences between d and e or c and e, likely due to large error bars in the recovered spectrum of e. Compared to M, L, and T-type field brown dwarfs, all three planets are most consistent with mid and late L spectral types. All objects are consistent with low gravity but a lack of standard spectra for low gravity limit the ability to fit the best spectral type. We discuss how dedicated modeling efforts can better fit HR 8799 planets' near-IR flux and discuss how differences between the properties of these planets can be further explored.

We present results of 2D axisymmetric core-collapse supernova simulations, employing the FORNAX code, of nine progenitor models spanning 12 to 25 M$_{\odot}$ and evolved over a 20,000-km grid. We find that four of the nine models explode with inelastic scattering off electrons and neutrons as well as the many-body correction to neutrino-nucleon scattering opacities. We show that these four models feature sharp Si-O interfaces in their density profiles, and that the corresponding dip in density reduces the accretion rate around the stalled shock and prompts explosion. The non-exploding models lack such a steep feature, suggesting that Si-O interface is one key to explosion. Furthermore, we show that all of the non-exploding models can be nudged to explosion with modest changes to macrophysical inputs, including moderate rotation and perturbations to infall velocities, as well as to microphysical inputs, including changes to neutrino-nucleon interaction rates, suggesting that all the models are perhaps close to criticality. Exploding models have energies of few $\times$10$^{50}$ ergs at the end of our simulation, and are rising, suggesting the need to continue these simulations over larger grids and for longer times to reproduce the energies seen in Nature. We find that the morphology of the explosion contributes to the explosion energy, with more isotropic ejecta producing larger explosion energies. However, we do not find evidence for the Lepton-number Emission Self-Sustained Asymmetry. Finally, we look at PNS properties and explore the role of dimension in our simulations. We find that convection in the proto-neutron star (PNS) produces larger PNS radii as well as greater "$\nu_\mu$" luminosities in 2D compared to 1D.

We explore with self-consistent 2D F{\sc{ornax}} simulations the dependence of the outcome of collapse on many-body corrections to neutrino-nucleon cross sections, the nucleon-nucleon bremsstrahlung rate, electron capture on heavy nuclei, pre-collapse seed perturbations, and inelastic neutrino-electron and neutrino-nucleon scattering. Importantly, proximity to criticality amplifies the role of even small changes in the neutrino-matter couplings, and such changes can together add to produce outsized effects. When close to the critical condition the cumulative result of a few small effects (including seeds) that individually have only modest consequence can convert an anemic into a robust explosion, or even a dud into a blast. Such sensitivity is not seen in one dimension and may explain the apparent heterogeneity in the outcomes of detailed simulations performed internationally. A natural conclusion is that the different groups collectively are closer to a realistic understanding of the mechanism of core-collapse supernovae than might have seemed apparent.

We present new 3.6 and 4.5 $\mu m$ Spitzer phase curves for the highly irradiated hot Jupiter WASP-33b and the unusually dense Saturn-mass planet HD 149026b. As part of this analysis, we develop a new variant of pixel level decorrelation that is effective at removing intrapixel sensitivity variations for long observations (>10 hours) where the position of the star can vary by a significant fraction of a pixel. Using this algorithm, we measure eclipse depths, phase amplitudes, and phase offsets for both planets at 3.6 $\mu m$ and 4.5 $\mu m$. We use a simple toy model to show that WASP-33b's phase offset, albedo, and heat recirculation efficiency are largely similar to those of other hot Jupiters despite its very high irradiation. On the other hand, our fits for HD 149026b prefer a very high albedo and an unusually high recirculation efficiency. We also compare our results to predictions from general circulation models, and find that while neither are a good match to the data, the discrepancies for HD 149026b are especially large. We speculate that this may be related to its high bulk metallicity, which could lead to enhanced atmospheric opacities and the formation of reflective cloud layers in localized regions of the atmosphere. We then place these two planets in a broader context by exploring relationships between the temperatures, albedos, heat transport efficiencies, and phase offsets of all planets with published thermal phase curves. We find a striking relationship between phase offset and irradiation temperature--the former drops with increasing temperature until around 3400 K, and rises thereafter. Although some aspects of this trend are mirrored in the circulation models, there are notable differences that provide important clues for future modeling efforts.

We present new 1D (spherical) and 2D (axisymmetric) simulations of electron-capture (EC) and low-mass iron-core-collapse supernovae (SN). We consider six progenitor models: the ECSN progenitor from Nomoto (1984, 1987); two ECSN-like low-mass low-metallicity iron core progenitors from Heger (private communication); and the 9-, 10-, and 11-$M_\odot$ (zero-age main sequence) progenitors from Sukhbold et al. (2016). We confirm that the ECSN and ESCN-like progenitors explode easily even in 1D with explosion energies of up to a 0.15 Bethes ($1 {\rm B} \equiv 10^{51}\ {\rm erg}$), and are a viable mechanism for the production of very low-mass neutron stars. However, the 9-, 10-, and 11-$M_\odot$ progenitors do not explode in 1D and are not even necessarily easier to explode than higher-mass progenitor stars in 2D. We study the effect of perturbations and of changes to the microphysics and we find that relatively small changes can result in qualitatively different outcomes, even in 1D, for models sufficiently close to the explosion threshold. Finally, we revisit the impact of convection below the protoneutron star (PNS) surface. We analyze, 1D and 2D evolutions of PNSs subject to the same boundary conditions. We find that the impact of PNS convection has been underestimated in previous studies and could result in an increase of the neutrino luminosity by up to factors of two.

The U.S. astronomy decadal surveys have been models for advice to government on how to apportion resources to optimise the scientific return on national investments in facilities and manpower. The U.S. is now gearing up to conduct its 2020 survey and the results are likely to guide international astronomy far into the future. Here, I summarize the current strains in an otherwise world-leading program of ground- and space-based astronomical discovery and some of the issues that will be faced by the participants in this upcoming collective exercise.

We present spectro-photometry spanning 1-5 $\mu$m of 51 Eridani b, a 2-10 M$_\text{Jup}$ planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new $K1$ (1.90-2.19 $\mu$m) and $K2$ (2.10-2.40 $\mu$m) spectra taken with the Gemini Planet Imager as well as an updated $L_P$ (3.76 $\mu$m) and new $M_S$ (4.67 $\mu$m) photometry from the NIRC2 Narrow camera. The new data were combined with $J$ (1.13-1.35 $\mu$m) and $H$ (1.50-1.80 $\mu$m) spectra from the discovery epoch with the goal of better characterizing the planet properties. 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4-T8) and we propose that 51 Eri b might be in the process of undergoing the transition from L-type to T-type. We used two complementary atmosphere model grids including either deep iron/silicate clouds or sulfide/salt clouds in the photosphere, spanning a range of cloud properties, including fully cloudy, cloud free and patchy/intermediate opacity clouds. Model fits suggest that 51 Eri b has an effective temperature ranging between 605-737 K, a solar metallicity, a surface gravity of $\log$(g) = 3.5-4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the SED. From the model atmospheres, we infer a luminosity for the planet of -5.83 to -5.93 ($\log L/L_{\odot}$), leaving 51 Eri b in the unique position as being one of the only directly imaged planet consistent with having formed via cold-start scenario. Comparisons of the planet SED against warm-start models indicates that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 M$_{\oplus}$.

We present new Keck/NIRC2 3$-$5 $\mu$m infrared photometry of the planetary-mass companion to ROXS 42B in $L^\prime$, and for the first time in Brackett-$\alpha$ (Br$\alpha$) and in $M_\mathrm{s}$-band. We combine our data with existing near-infrared photometry and $K$-band (2$-$2.4 $\mu$m) spectroscopy and compare these with models and other directly imaged planetary-mass objects using forward modeling and retrieval methods in order to characterize the atmosphere of ROXS 42B b. ROXS 42B b's 1.25$-$5 $\mu$m spectral energy distribution most closely resembles that of GSC 06214 B and $\kappa$ And b, although it has a slightly bluer $K_{\rm s}$$-$$M_{\rm s}$ color than GSC 06214 B and thus so far lacks evidence for a circumplanetary disk. We cannot formally exclude the possibility that any of the tested dust-free/dusty/cloudy forward models describe atmosphere of ROXS 42B b well. However, models with substantial atmospheric dust/clouds yield temperatures and gravities that are consistent when fit to photometry and spectra separately, whereas dust-free model fits to photometry predict temperatures/gravities inconsistent with ROXS 42B b's $K$-band spectrum and vice-versa. Atmospheric retrieval on the 1$-$5 $\mu$m photometry places a limit on the fractional number density of CO$_2$ of $\log(n_{\rm CO_2})<-2.7$ but provides no other constraints so far. We conclude that ROXS 42B b has mid-IR photometric features that are systematically different from other previously observed planetary-mass and field objects of similar temperature. It remains unclear whether this is in the range of the natural diversity of targets at the very young ($\sim$2 Myr) age of ROXS 42B b, or unique to its early evolution and environment.

Using the Gemini Planet Imager (GPI) located at Gemini South, we measured the near-infrared (1.0-2.4 micron) spectrum of the planetary companion to the nearby, young star $\beta$ Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of $12.9\pm0.2$ $\mathcal{M}_\mathrm{Jup}$, an effective temperature of $1724\pm15$ K, a radius of $1.46\pm0.01$ $\mathcal{R}_\mathrm{Jup}$, and a surface gravity of $\log g = 4.18\pm0.01$ [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of $1700-1800$ K and a surface gravity of $\log g = 3.5$-$4.0$ [dex] depending upon model. These values agree well with other publications and with "hot-start" predictions from planetary evolution models. Further, we find that the spectrum of $\beta$ Pic b best matches a low-surface gravity L2$\pm$1 brown dwarf. Finally comparing the spectrum to field brown dwarfs we find the the spectrum best matches 2MASS J04062677-381210 and 2MASS J03552337+1133437.

Extrasolar planets on eccentric short-period orbits provide a laboratory in which to study radiative and tidal interactions between a planet and its host star under extreme forcing conditions. Studying such systems probes how the planet's atmosphere redistributes the time-varying heat flux from its host and how the host star responds to transient tidal distortion. Here, we report the insights into the planet-star interactions in HAT-P-2's eccentric planetary system gained from the analysis of 350 hr of 4.5 micron observations with the Spitzer Space Telescope. The observations show no sign of orbit-to-orbit variability nor of orbital evolution of the eccentric planetary companion, HAT-P-2 b. The extensive coverage allows us to better differentiate instrumental systematics from the transient heating of HAT-P-2 b's 4.5 micron photosphere and yields the detection of stellar pulsations with an amplitude of approximately 40 ppm. These pulsation modes correspond to exact harmonics of the planet's orbital frequency, indicative of a tidal origin. Transient tidal effects can excite pulsation modes in the envelope of a star, but, to date, such pulsations had only been detected in highly eccentric stellar binaries. Current stellar models are unable to reproduce HAT-P-2's pulsations, suggesting that our understanding of the interactions at play in this system is incomplete.

We provide an example of an analysis to explore the optimization of observations of transiting hot jupiters with JWST to characterize their atmospheres, based on a simple three-parameter forward model. We construct expansive forward model sets for eleven hot jupiters, ten of which are relatively well-characterized, exploring a range of parameters such as equilibrium temperature and metallicity, as well as considering host stars over a wide range in brightness. We compute posterior distributions of our model parameters for each planet with all of the available JWST spectroscopic modes and several programs of combined observations and compute their effectiveness using the metric of estimated mutual information per degree of freedom. From these simulations, clear trends emerge that provide guidelines for designing a JWST observing program. We demonstrate that these guidelines apply over a wide range of planet parameters and target brightnesses for our simple forward model.

We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with non-isotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-$k$ approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci & Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; $k$-distribution tables constructed with $\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$ points per wavenumber bin) may result in errors $\gtrsim 1$-10 % in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct "null-hypothesis" models (chemical equilibrium, radiative equilibrium and solar element abundances) for 6 hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while it consistently under-predicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (ESP; exoclime.net).

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.

We perform the first self-consistent, time-dependent, multi-group calculations in two dimensions (2D) to address the consequences of using the ray-by-ray+ transport simplification in core-collapse supernova simulations. Such a dimensional reduction is employed by many researchers to facilitate their resource-intensive calculations. Our new code (F{\sc{ornax}}) implements multi-D transport, and can, by zeroing out transverse flux terms, emulate the ray-by-ray+ scheme. Using the same microphysics, initial models, resolution, and code, we compare the results of simulating 12-, 15-, 20-, and 25-M$_{\odot}$ progenitor models using these two transport methods. Our findings call into question the wisdom of the pervasive use of the ray-by-ray+ approach. Employing it leads to maximum post-bounce/pre-explosion shock radii that are almost universally larger by tens of kilometers than those derived using the more accurate scheme, typically leaving the post-bounce matter less bound and artificially more "explodable." In fact, for our 25-M$_{\odot}$ progenitor, the ray-by-ray+ model explodes, while the corresponding multi-D transport model does not. Therefore, in two dimensions the combination of ray-by-ray+ with the axial sloshing hydrodynamics that is a feature of 2D supernova dynamics can result in quantitatively, and perhaps qualitatively, incorrect results.

A principal scientific goal of the Gemini Planet Imager (GPI) is obtaining milliarcsecond astrometry to constrain exoplanet orbits. However, astrometry of directly imaged exoplanets is subject to biases, systematic errors, and speckle noise. Here we describe an analytical procedure to forward model the signal of an exoplanet that accounts for both the observing strategy (angular and spectral differential imaging) and the data reduction method (Karhunen-Lo\`eve Image Projection algorithm). We use this forward model to measure the position of an exoplanet in a Bayesian framework employing Gaussian processes and Markov chain Monte Carlo (MCMC) to account for correlated noise. In the case of GPI data on $\beta$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$\sigma$-errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to twelve GPI observations along with previous astrometry from other instruments. The statistical properties of the residuals confirm that BKA is accurate and correctly estimates astrometric errors. Our constraints on the orbit of $\beta$ Pic b firmly rule out the possibility of a transit of the planet at 10-$\sigma$ significance. However, we confirm that the Hill sphere of $\beta$ Pic b will transit, giving us a rare chance to probe the circumplanetary environment of a young, evolving exoplanet. We provide an ephemeris for photometric monitoring of the Hill sphere transit event, which will begin at the start of April in 2017 and finish at the end of January in 2018.

We present HST STIS optical transmission spectroscopy of the cool Saturn-mass exoplanet WASP-39b from 0.29-1.025 micron, along with complementary transit observations from Spitzer IRAC at 3.6 and 4.5 micron. The low density and large atmospheric pressure scale height ofWASP-39b make it particularly amenable to atmospheric characterization using this technique. We detect a Rayleigh scattering slope as well as sodium and potassium absorption features; this is the first exoplanet in which both alkali features are clearly detected with the extended wings predicted by cloud-free atmosphere models. The full transmission spectrum is well matched by a clear, H2-dominated atmosphere or one containing a weak contribution from haze, in good agreement with the preliminary reduction of these data presented in Sing et al. (2016). WASP-39b is predicted to have a pressure-temperature profile comparable to that of HD 189733b and WASP-6b, making it one of the coolest transiting gas giants observed in our HST STIS survey. Despite this similarity, WASP-39b appears to be largely cloud-free while the transmission spectra of HD 189733b and WASP-6b both indicate the presence of high altitude clouds or hazes. These observations further emphasize the surprising diversity of cloudy and cloud-free gas giant planets in short-period orbits and the corresponding challenges associated with developing predictive cloud models for these atmospheres.

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 $\mu$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 $\mu$m, which are consistent with a single blackbody with effective temperature $2372 \pm 60$ K. The measured 3.6 and 4.5 $\mu$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 $\mu$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\sigma$) 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.

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)

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)

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 $\mu$m and 4.5 $\mu$m bands, best match atmospheric model predictions with a dayside temperature inversion and relatively efficient day-night circulation.

We calculate the distance-dependent performance of a few representative terrestrial neutrino detectors in detecting and measuring the properties of the $\nu_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 $\nu_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 $\nu_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 $\nu_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.

Thousands of transiting exoplanets have been discovered, but spectral analysis of their atmospheres has so far been dominated by a small number of exoplanets and data spanning relatively narrow wavelength ranges (such as 1.1 to 1.7 {\mu}m). Recent studies show that some hot-Jupiter exoplanets have much weaker water absorption features in their near-infrared spectra than predicted. The low amplitude of water signatures could be explained by very low water abundances, which may be a sign that water was depleted in the protoplanetary disk at the planet's formation location, but it is unclear whether this level of depletion can actually occur. Alternatively, these weak signals could be the result of obscuration by clouds or hazes, as found in some optical spectra. Here we report results from a comparative study of ten hot Jupiters covering the wavelength range 0.3-5 micrometres, which allows us to resolve both the optical scattering and infrared molecular absorption spectroscopically. Our results reveal a diverse group of hot Jupiters that exhibit a continuum from clear to cloudy atmospheres. We find that the difference between the planetary radius measured at optical and infrared wavelengths is an effective metric for distinguishing different atmosphere types. The difference correlates with the spectral strength of water, so that strong water absorption lines are seen in clear-atmosphere planets and the weakest features are associated with clouds and hazes. This result strongly suggests that primordial water depletion during formation is unlikely and that clouds and hazes are the cause of weaker spectral signatures.

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.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

IW ../submit_V2/abstract.txt ( Row 1 Col 1 6:48 Ctrl-K H for help We report the discovery by the HATSouth network of HATS-7b, a transiting Super-Neptune with a mass of 0.120+/-0.012MJ, a radius of 0.563+/-(0.046,0.034)RJ, and an orbital period of 3.1853days. The host star is a moderately bright (V=13.340+/-0.010mag, K_S=10.976+/-0.026mag) K dwarf star with a mass of 0.849+/-0.027Msun , a radius of 0.815+/-(0.049,-0.035)Rsun, and a metallicity of [Fe/H]=+0.250+/-0.080. The star is photometrically quiet to within the precision of the HATSouth measurements and has low RV jitter. HATS-7b is the second smallest radius planet discovered by a wide-field ground-based transit survey, and one of only a handful of Neptune-size planets with mass and radius determined to 10% precision. Theoretical modeling of HATS-7b yields a hydrogen-helium fraction of 18+/-4% (rock-iron core and H2-He envelope), or 9+/-4% (ice core and H2-He envelope), i.e.it has a composition broadly similar to that of Uranus and Neptune, and very different from that of Saturn, which has 75% of its mass in H2-He. Based on a sample of transiting exoplanets with accurately (<20%) determined parameters, we establish approximate power-law relations for the envelopes of the mass-density distribution of exoplanets. HATS-7b, which, together with the recently discovered HATS-8b, is one of the first two transiting super-Neptunes discovered in the Southern sky, is a prime target for additional follow-up observations with Southern hemisphere facilities to characterize the atmospheres of Super-Neptunes (which we define as objects with mass greater than that of Neptune, and smaller than halfway between that of Neptune and Saturn, i.e. 0.054 MJ<Mp<0.18 MJ).

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 $\mu$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 $\mu$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2\sigma$ upper limit on the nightside flux at 3.6 $\mu$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $\mu$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 $\mu$m, while under-predicting the nightside emission at 4.5 $\mu$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.

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.

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.

We construct models of the structural evolution of super-Earth- and mini-Neptune-type exoplanets with hydrogen-helium envelopes, incorporating radiative cooling and XUV-driven mass loss. We conduct a parameter study of these models, focusing on initial mass, radius, and envelope mass fractions, as well as orbital distance, metallicity, and the specific prescription for mass loss. From these calculations, we investigate how the observed masses and radii of exoplanets today relate to the distribution of their initial conditions. Orbital distance and initial envelope mass fraction are the most important factors determining planetary evolution, particular radius evolution. Initial mass also becomes important below a "turnoff mass," which varies with orbital distance, with mass-radius curves being approximately flat for higher masses. Initial radius is the least important parameter we study, with very little difference between the hot start and cold start limits after an age of 100 Myr. Model sets with no mass loss fail to produce results consistent with observations, but a plausible range of mass loss scenarios is allowed. In addition, we present scenarios for the formation of the Kepler-11 planets. Our best fit to observations Kepler-11b and Kepler-11c involves formation beyond the snow line, after which they moved inward, circularized, and underwent a reduced degree mass loss.

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.

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.

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.

We report Hubble Space Telescope (HST) optical to near-infrared transmission spectroscopy of the hot Jupiter WASP-6b, measured with the Space Telescope Imaging Spectrograph (STIS) and Spitzer's InfraRed Array Camera (IRAC). The resulting spectrum covers the range $0.29-4.5\,\mu$m. We find evidence for modest stellar activity of WASP-6b and take it into account in the transmission spectrum. The overall main characteristic of the spectrum is an increasing radius as a function of decreasing wavelength corresponding to a change of $\Delta (R_p/R_{\ast})=0.0071$ from 0.33 to $4.5\,\mu$m. The spectrum suggests an effective extinction cross-section with a power law of index consistent with Rayleigh scattering, with temperatures of $973\pm144$ K at the planetary terminator. We compare the transmission spectrum with hot-Jupiter atmospheric models including condensate-free and aerosol-dominated models incorporating Mie theory. While none of the clear-atmosphere models is found to be in good agreement with the data, we find that the complete spectrum can be described by models that include significant opacity from aerosols including Fe-poor Mg$_2$SiO$_4$, MgSiO$_3$, KCl and Na$_2$S dust condensates. WASP-6b is the second planet after HD189733b which has equilibrium temperatures near $\sim1200$ K and shows prominent atmospheric scattering in the optical.

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.

We present Hubble Space Telescope optical and near-IR transmission spectra of the transiting hot-Jupiter WASP-31b. The spectrum covers 0.3-1.7 $\mu$m at a resolution $R\sim$70, which we combine with Spitzer photometry to cover the full-optical to IR. The spectrum is dominated by a cloud-deck with a flat transmission spectrum which is apparent at wavelengths $>0.52\mu$m. The cloud deck is present at high altitudes and low pressures, as it covers the majority of the expected optical Na line and near-IR H$_2$O features. While Na I absorption is not clearly identified, the resulting spectrum does show a very strong potassium feature detected at the 4.2-$\sigma$ confidence level. Broadened alkali wings are not detected, indicating pressures below $\sim$10 mbar. The lack of Na and strong K is the first indication of a sub-solar Na/K abundance ratio in a planetary atmosphere (ln[Na/K]$=-3.3\pm2.8$), which could potentially be explained by Na condensation on the planet's night side, or primordial abundance variations. A strong Rayleigh scattering signature is detected at short wavelengths, with a 4-$\sigma$ significant slope. Two distinct aerosol size populations can explain the spectra, with a smaller sub-micron size grain population reaching high altitudes producing a blue Rayleigh scattering signature on top of a larger, lower-lying population responsible for the flat cloud deck at longer wavelengths. We estimate that the atmospheric circulation is sufficiently strong to mix micron size particles upward to the required 1-10 mbar pressures, necessary to explain the cloud deck. These results further confirm the importance of clouds in hot-Jupiters, which can potentially dominate the overall spectra and may alter the abundances of key gaseous species.

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.

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 $\sigma$ 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.

We analyze all existing secondary eclipse time series spectroscopy of hot Jupiter HD189733b acquired with the now defunct Spitzer/IRS instrument. We describe the novel approaches we develop to remove the systematic effects and extract accurate secondary eclipse depths as a function of wavelength in order to construct the emission spectrum of the exoplanet. We compare our results to a previous study by Grillmair et al. that did not examine all data sets available to us. We are able to confirm the detection of a water feature near 6{\mu}m claimed by Grillmair et al. We compare the planetary emission spectrum to three model families -- based on isothermal atmosphere, gray atmosphere, and two realizations of the complex radiative transfer model by Burrows et al., adopted in Grillmair et al.'s study. While we are able to reject the simple isothermal and gray models based on the data at the 97% level just from the IRS data, these rejections hinge on eclipses measured within relatively narrow wavelength range, between 5.5 and 7{\mu}m. This underscores the need for observational studies with broad wavelength coverage and high spectral resolution, in order to obtain robust information on exoplanet atmospheres.

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.

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5~$\mu$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.

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 $\mu$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $\mu$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.

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.

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 $\lambda/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 $\kappa$ And b.

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.

We present new $L^\prime$ (3.8 $\mu m$) and Br-$\alpha$ (4.05 $\mu 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-$\sigma$ 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.

Using the recently installed Gemini Planet Imager (GPI), we have taken the first H-band spectrum of the planetary companion to the nearby young star beta Pictoris. GPI is designed to image and provide low-resolution spectra of Jupiter sized, self-luminous planetary companions around young nearby stars. These observations were taken covering the H-band (1.65 microns). The spectrum has a resolving power of $\sim$ 45 and demonstrates the distinctive triangular shape of a cool substellar object with low surface gravity. Using atmospheric models, we find an effective temperature of $1650 \pm 50$ K and a surface gravity of $\log(g) = 4.0 \pm 0.25$ (cgs units). These values agree well with predictions from planetary evolution models for a gas giant with mass between 10 and 12 $M_{\rm Jup}$ and age between 10 and 20 Myrs.

Many exoplanets have been discovered with radii of 1-4 Earth radii, between that of Earth and Neptune. A number of these are known to have densities consistent with solid compositions, while others are "sub-Neptunes" likely to have significant hydrogen-helium envelopes. Future surveys will no doubt significantly expand these populations. In order to understand how the measured masses and radii of such planets can inform their structures and compositions, we construct models both for solid layered planets and for planets with solid cores and gaseous envelopes, exploring a range of core masses, hydrogen-helium envelope masses, and associated envelope entropies. For planets in the super-Earth/sub-Neptune regime for which both radius and mass are measured, we estimate how each is partitioned into a solid core and gaseous envelope, associating a specific core mass and envelope mass with a given exoplanet. We perform this decomposition for both "Earth-like" rock-iron cores and pure ice cores, and find that the necessary gaseous envelope masses for this important sub-class of exoplanets must range very widely from zero to many Earth masses, even for a given core mass. This result bears importantly on exoplanet formation and envelope evaporation processes.

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.

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.

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.

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.

Super-Earths transiting nearby bright stars are key objects that simultaneously allow for accurate measurements of both their mass and radius, providing essential constraints on their internal composition. We present here the confirmation, based on Spitzer transit observations, that the super-Earth HD 97658 b transits its host star. HD 97658 is a low-mass ($M_*=0.77\pm0.05\,M_{\odot}$) K1 dwarf, as determined from the Hipparcos parallax and stellar evolution modeling. To constrain the planet parameters, we carry out Bayesian global analyses of Keck-HIRES radial velocities, and MOST and Spitzer photometry. HD 97658 b is a massive ($M_P=7.55^{+0.83}_{-0.79} M_{\oplus}$) and large ($R_{P} = 2.247^{+0.098}_{-0.095} R_{\oplus}$ at 4.5 $\mu$m) super-Earth. We investigate the possible internal compositions for HD 97658 b. Our results indicate a large rocky component, by at least 60% by mass, and very little H-He components, at most 2% by mass. We also discuss how future asteroseismic observations can improve the knowledge of the HD 97658 system, in particular by constraining its age. Orbiting a bright host star, HD 97658 b will be a key target for coming space missions TESS, CHEOPS, PLATO, and also JWST, to characterize thoroughly its structure and atmosphere.

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$\mu$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\sigma$ 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.

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.

Unlike hot Jupiters or other gas giants, super-Earths are expected to have a wide variety of compositions, ranging from terrestrial bodies like our own to more gaseous planets like Neptune. Observations of transiting systems, which allow us to directly measure planet masses and radii and constrain atmospheric properties, are key to understanding the compositional diversity of the planets in this mass range. Although Kepler has discovered hundreds of transiting super-Earth candidates over the past four years, the majority of these planets orbit stars that are too far away and too faint to allow for detailed atmospheric characterization and reliable mass estimates. Ground-based transit surveys focus on much brighter stars, but most lack the sensitivity to detect planets in this size range. One way to get around the difficulty of finding these smaller planets in transit is to start by choosing targets that are already known to contain super-Earth sized bodies detected using the radial velocity technique. Here we present results from a Spitzer program to observe six of the most favorable RV detected super-Earth systems, including HD 1461, HD 7924, HD 156668, HIP 57274, and GJ 876. We find no evidence for transits in any of their 4.5 micron flux light curves, and place limits on the allowed transit depths and corresponding planet radii that rule out even the most dense and iron-rich compositions for these objects. We also observed HD 97658, but the observation window was based on a possible ground-based transit detection (Henry et al. 2011) that was later ruled out; thus the window did not include the predicted time for the transit detection recently made by MOST (Dragomir et al. 2013).

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.

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.

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.

We present an optical to near-infrared transmission spectrum of the hot Jupiter HAT-P-1b, based on HST observations, covering the spectral regime from 0.29 to 1.027{\mu}m with STIS, which is coupled with a recent WFC3 transit (1.087 to 1.687{\mu}m). We derive refined physical parameters of the HAT-P-1 system, including an improved orbital ephemeris. The transmission spectrum shows a strong absorption signature shortward of 0.55{\mu}m, with a strong blueward slope into the near-ultraviolet. We detect atmospheric sodium absorption at a 3.3{\sigma} significance level, but find no evidence for the potassium feature. The red data implies a marginally flat spectrum with a tentative absorption enhancement at wavelength longer than ~0.85{\mu}m. The STIS and WFC3 spectra differ significantly in absolute radius level (4.3 +/- 1.6 pressure scale heights), implying strong optical absorption in the atmosphere of HAT-P-1b. The optical to near-infrared difference cannot be explained by stellar activity, as simulta- neous stellar activity monitoring of the G0V HAT-P-1b host star and its identical companion show no significant activity that could explain the result. We compare the complete STIS and WFC3 transmission spectrum with theoretical atmospheric mod- els which include haze, sodium and an extra optical absorber. We find that both an optical absorber and a super-solar sodium to water abundance ratio might be a sce- nario explaining the HAT-P-1b observations. Our results suggest that strong optical absorbers may be a dominant atmospheric feature in some hot Jupiter exoplanets.

We present HST optical transmission spectra of the transiting hot Jupiter WASP-12b, taken with the STIS instrument. From the transmission spectra, we are able to decisively rule out prominent absorption by TiO in the exoplanet's atmosphere. Strong pressure-broadened Na and K absorption signatures are also excluded, as are significant metal-hydride features. We compare our combined broadband spectrum to a wide variety of existing aerosol-free atmospheric models, though none are satisfactory fits. However, we do find that the full transmission spectrum can be described by models which include significant opacity from aerosols: including Rayleigh scattering, Mie scattering, tholin haze, and settling dust profiles. The transmission spectrum follows an effective extinction cross section with a power-law of index alpha, with the slope of the transmission spectrum constraining the quantity alphaT = -3528+/-660 K, where T is the atmospheric temperature. Rayleigh scattering (alpha=-4) is among the best fitting models, though requires low terminator temperatures near 900 K. Sub-micron size aerosol particles can provide equally good fits to the entire transmission spectrum for a wide range of temperatures, and we explore corundum as a plausible dust aerosol. The presence of atmospheric aerosols also helps to explain the modestly bright albedo implied by Spitzer observations, as well as the near black body nature of the emission spectrum. Ti-bearing condensates on the cooler night-side is the most natural explanation for the overall lack of TiO signatures in WASP-12b, indicating the day/night cold-trap is an important effect for very hot Jupiters. These finding indicate that aerosols can play a significant atmospheric role for the entire wide range of hot-Jupiter atmospheres, potentially affecting their overall spectrum and energy balance.(abridged)

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $\beta$ 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 $\beta$ 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. $\kappa$ And b and 1RXJ 1609B -- match $\beta$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $\mu m$ and 5 $\mu m$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $\mu m$) dust grains fail to reproduce the $\beta$ 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 $\mu 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 $\beta$ 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.

Several exoplanets have recently been imaged at wide separations of >10 AU from their parent stars. These span a limited range of ages (<50 Myr) and atmospheric properties, with temperatures of 800--1800 K and very red colors (J - H > 0.5 mag), implying thick cloud covers. Furthermore, substantial model uncertainties exist at these young ages due to the unknown initial conditions at formation, which can lead to an order of magnitude of uncertainty in the modeled planet mass. Here, we report the direct imaging discovery of a Jovian exoplanet around the Sun-like star GJ 504, detected as part of the SEEDS survey. The system is older than all other known directly-imaged planets; as a result, its estimated mass remains in the planetary regime independent of uncertainties related to choices of initial conditions in the exoplanet modeling. Using the most common exoplanet cooling model, and given the system age of 160 [+350, -60] Myr, GJ 504 b has an estimated mass of 4 [+4.5, -1.0] Jupiter masses, among the lowest of directly imaged planets. Its projected separation of 43.5 AU exceeds the typical outer boundary of ~30 AU predicted for the core accretion mechanism. GJ 504 b is also significantly cooler (510 [+30, -20] K) and has a bluer color (J-H = -0.23 mag) than previously imaged exoplanets, suggesting a largely cloud-free atmosphere accessible to spectroscopic characterization. Thus, it has the potential of providing novel insights into the origins of giant planets, as well as their atmospheric properties.

We present Hubble Space Telescope near-infrared transmission spectroscopy of the transiting hot-Jupiter HAT-P-1b. We observed one transit with Wide Field Camera 3 using the G141 low-resolution grism to cover the wavelength range 1.087- 1.678 {\mu}m. These time series observations were taken with the newly available spatial scan mode that increases the duty cycle by nearly a factor of two, thus improving the resulting photometric precision of the data. We measure a planet-to-star radius ratio of Rp/R*=0.11709+/-0.00038 in the white light curve with the centre of transit occurring at 2456114.345+/-0.000133 (JD). We achieve S/N levels per exposure of 1840 (0.061%) at a resolution of {\Delta\lambda}=19.2nm (R~70) in the 1.1173 - 1.6549{\mu}m spectral region, providing the precision necessary to probe the transmission spectrum of the planet at close to the resolution limit of the instrument. We compute the transmission spectrum using both single target and differential photometry with similar results. The resultant transmission spectrum shows a significant absorption above the 5-{\sigma} level matching the 1.4{\mu}m water absorption band. In solar composition models, the water absorption is sensitive to the ~1 mbar pressure levels at the terminator. The detected absorption agrees with that predicted by an 1000 K isothermal model, as well as with that predicted by a planetary-averaged temperature model.

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.

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.

We measure the transmission spectrum of WASP-19b from 3 transits using low-resolution optical spectroscopy from the HST Space Telescope Imaging Spectrograph (STIS). The STIS spectra cover a wavelength range of 0.29-1.03 microns with resolving power R=500. The optical data are combined with archival near-IR data from the HST Wide Field Camera 3 (WFC3) G141 grism, covering the wavelength range 1.087-1.687 micron, with resolving power R=130. We obtain the transmission spectrum from 0.53-1.687 microns with S/N levels between 3000 and 11,000 in 0.1 micron bins. WASP-19 is a very active star, with optical stellar flux varying by a few per cent over time. We correct the transit light curves for the effects of stellar activity using ground-based activity monitoring with the Cerro Tololo Inter-American Observatory (CTIO). While we were not able to construct a transmission spectrum using the blue optical data due to the presence of large occulted star spots, we were able to use the spot crossings to help constrain the mean stellar spot temperature. To search for predicted features in the hot-Jupiter atmosphere, we also define spectral indices for differential radius measurements to specifically search for the presence of TiO and alkali line features. Our measurements rule out TiO features predicted for a planet of WASP-19b's equilibrium temperature (2050 K) in the transmission spectrum at the 2.7-2.9 sigma confidence level, depending on atmospheric model formalism. The WFC3 transmission spectrum shows strong absorption features due to the presence of H2O, which is detected at the 4 sigma confidence level between 1.1 and 1.4 microns. The results indicate that WASP-19b is a planet with no or low levels of TiO and without a high C/O ratio. The lack of observable TiO features are possibly due to rainout, breakdown from stellar activity or the presence of other absorbers in the optical.

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.

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.

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_{\theta} + K_{\phi}). 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.

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$\sigma$ 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.

We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse phase, and a multi-species neutrino leakage scheme after bounce. The latter captures the important effects of deleptonization, neutrino cooling and heating and enables approximate predictions for the neutrino luminosities in the early evolution after core bounce. We consider 12-solar-mass and 40-solar-mass presupernova models and systematically study the effects of (i) rotation, (ii) progenitor structure, and (iii) postbounce neutrino leakage on dynamics, GW, and, neutrino signals. We demonstrate, that the GW signal of rapidly rotating core collapse is practically independent of progenitor mass and precollapse structure. Moreover, we show that the effects of neutrino leakage on the GW signal are strong only in nonrotating or slowly rotating models in which GW emission is not dominated by inner core dynamics. In rapidly rotating cores, core bounce of the centrifugally-deformed inner core excites the fundamental quadrupole pulsation mode of the nascent protoneutron star. The ensuing global oscillations (f~700-800 Hz) lead to pronounced oscillations in the GW signal and correlated strong variations in the rising luminosities of antineutrino and heavy-lepton neutrinos. We find these features in cores that collapse to protoneutron stars with spin periods <~ 2.5 ms and rotational energies sufficient to drive hyper-energetic core-collapse supernova explosions. Hence, joint GW + neutrino observations of a core collapse event could deliver strong evidence for or against rapid core rotation. [abridged]

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.

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 \mu m bands of the Spitzer Space Telescope. The 3.6 and 4.5 \mu 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 \mu 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 \mu 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.

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.

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.

We present a formulation for multigroup radiation hydrodynamics that is correct to order $O(v/c)$ using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically-rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts, one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.

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.

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.

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.

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.

We present new calculations of transit spectra of super-Earths that allow for atmospheres with arbitrary proportions of common molecular species and haze. We test this method with generic spectra, reproducing the expected systematics and absorption features, then apply it to the nearby super-Earth GJ 1214b, which has produced conflicting observational data, leaving the questions of a hydrogen-rich versus hydrogen-poor atmosphere and the water content of the atmosphere ambiguous. We present representative transit spectra for a range of classes of atmosphere models for GJ 1214b. Our analysis supports a hydrogen-rich atmosphere with a cloud or haze layer, although a hydrogen-poor model with less than 10% water is not ruled out. Several classes of models are ruled out, however, including hydrogen-rich atmospheres with no haze, hydrogen-rich atmospheres with a haze of about 0.01-micron tholin particles, and hydrogen-poor atmospheres with major sources of absorption other than water. We propose an observational test to distinguish hydrogen-rich from hydrogen-poor atmospheres. Finally, we provide a library of theoretical transit spectra for super-Earths with a broad range of parameters to facilitate future comparison with anticipated data.

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.

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.

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.

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]

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.

We present results from a suite of axisymmetric, core-collapse supernova simulations in which hydrodynamic recoil from an asymmetric explosion produces large proto-neutron star (PNS) velocities. We use the adaptive-mesh refinement code CASTRO to self-consistently follow core-collapse, the formation of the PNS and its subsequent acceleration. We obtain recoil velocities of up to 620 km/s at ~1 s after bounce. These velocities are consistent with the observed distribution of pulsar kicks and with PNS velocities obtained in other theoretical calculations. Our PNSs are still accelerating at several hundred km/s at the end of our calculations, suggesting that even the highest velocity pulsars may be explained by hydrodynamic recoil in generic, core-collapse supernovae.

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.

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

[Abridged] We report Warm Spitzer full-orbit phase observations of WASP-12b at 3.6 and 4.5 micron. We are able to measure the transit depths, eclipse depths, thermal and ellipsoidal phase variations at both wavelengths. The large amplitude phase variations, combined with the planet's previously-measured day-side spectral energy distribution, is indicative of non-zero Bond albedo and very poor day-night heat redistribution. The transit depths in the mid-infrared indicate that the atmospheric opacity is greater at 3.6 than at 4.5 micron, in disagreement with model predictions, irrespective of C/O ratio. The secondary eclipse depths are consistent with previous studies. We do not detect ellipsoidal variations at 3.6 micron, but our parameter uncertainties -estimated via prayer-bead Monte Carlo- keep this non-detection consistent with model predictions. At 4.5 micron, on the other hand, we detect ellipsoidal variations that are much stronger than predicted. If interpreted as a geometric effect due to the planet's elongated shape, these variations imply a 3:2 ratio for the planet's longest:shortest axes and a relatively bright day-night terminator. If we instead presume that the 4.5 micron ellipsoidal variations are due to uncorrected systematic noise and we fix the amplitude of the variations to zero, the best fit 4.5 micron transit depth becomes commensurate with the 3.6 micron depth, within the uncertainties. The relative transit depths are then consistent with a Solar composition and short scale height at the terminator. Assuming zero ellipsoidal variations also yields a much deeper 4.5 micron eclipse depth, consistent with a Solar composition and modest temperature inversion. We suggest future observations that could distinguish between these two scenarios.

The internal heat loss, or cooling, of a planet determines its structure and evolution. We study the effects of irradiation, metallicity of the atmosphere, heat redistribution, stratospheres, and the depth where the heat redistribution takes place on the atmospheric structure, the core entropy, and subsequently on the cooling of the interior of the planet. We address in a consistent fashion the coupling between the day and the night sides of a planet by means of model atmosphere calculations with heat redistribution. We assume that strong convection leads to the same entropy on the day and night sides and that gravity is the same on both hemispheres. We argue that the core cooling rates from the two hemispheres of a strongly irradiated planet may not be the same and that the difference depends on several important parameters. If the day-night heat redistribution is very efficient or if it takes place at the large optical depth, then the day-side and the night-side cooling may be comparable. However, if the day-night heat transport is not efficient or if it takes place at a shallow optical depth then there can be a large difference between the day- and the night-side cooling and the night side will cool more efficiently. If stellar irradiation becomes stronger, e.g. owing to stellar evolution or migration, cooling from both the day and the night sides is reduced. Enhanced metallicity of the atmosphere would act as an added "blanket" and reduces both the day- and the night-side cooling. However, a stratosphere on the planetary day side can enhance day-side cooling since its opacity acts as a "sunshade" that screens the stellar irradiation. These effects may also influence the well-known gravity darkening and bolometric albedo effects in interacting binaries, especially for strongly irradiated cold components.

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.

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{\mu}m and 4.5{\mu}m for each target. XO-4b exhibits a stronger eclipse depth at 4.5{\mu}m than at 3.6{\mu}m, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6{\mu}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(\omega) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

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.

We describe the development of a flux-limited gray radiation solver for the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically-rectangular variable-sized grids with simultaneous refinement in both space and time. The gray radiation solver is based on a mixed-frame formulation of radiation hydrodynamics. In our approach, the system is split into two parts, one part that couples the radiation and fluid in a hyperbolic subsystem, and another parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem is solved explicitly with a high-order Godunov scheme, whereas the parabolic part is solved implicitly with a first-order backward Euler method.

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 $\beta$ Pic b sculpts debris disk structures and whether a second planet is required to explain the warp/inclined disk.

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.

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.

GJ 758 B is a cold (~600K) companion to a Sun-like star at 29 AU projected separation, which was recently detected with high-contrast imaging. Here we present photometry of the companion in seven photometric bands from Subaru/HiCIAO, Gemini/NIRI and Keck/NIRC2, providing a rich sampling of the spectral energy distribution in the 1-5 micron wavelength range. A clear detection at 1.58 micron combined with an upper limit at 1.69 micron shows methane absorption in the atmosphere of the companion. The mass of the companion remains uncertain, but an updated age estimate indicates that the most likely mass range is ~30-40 Mjup. In addition, we present an updated astrometric analysis that imposes tighter constraints on GJ 758 B's orbit and identifies the proposed second candidate companion, "GJ 758 C", as a background star.

Secondary eclipses are a powerful tool to measure directly the thermal emission from extrasolar planets, and to constrain their type and physical parameters. We started a project to obtain reliable broad-band measurements of the thermal emission of transiting exoplanets. Ground-based high-cadence near-infrared relative photometry was used to obtain sub-millimagnitude precision light curve of a secondary eclipse of WASP-4b -- a 1.12 M_J hot Jupiter on a 1.34 day orbit around G7V star. The data show a clear ~10-\sigma detection of the planet's thermal emission at 2.2 \mu m. The calculated thermal emission corresponds to a fractional eclipse depth of 0.185^{+0.014}_{-0.013}%, with a related brightness temperature in Ks of T_B = 1995 \pm 40 K, centered at T_C = 2455102.61162^{+0.00071}_{-0.00077} HJD. We could set a limit on the eccentricity of e cos \omega=0.0027 \pm 0.0018, compatible with a near-circular orbit. The calculated brightness temperature, as well as the specific models suggest a highly inefficient redistribution of heat from the day-side to the night-side of the planet, and a consequent emission mainly from the day-side. The high-cadence ground-based technique is capable of detecting the faint signal of the secondary eclipse of extrasolar planets, making it a valuable complement to space-based mid-IR observations.

It has been suggested that the observed rotation periods of radio pulsars might be induced by a non-axisymmetric spiral-mode instability in the turbulent region behind the stalled supernova bounce shock, even if the progenitor core was not initially rotating. In this paper, using the three-dimensional AMR code CASTRO with a realistic progenitor and equation of state and a simple neutrino heating and cooling scheme, we present a numerical study of the evolution in 3D of the rotational profile of a supernova core from collapse, through bounce and shock stagnation, to delayed explosion. By the end of our simulation ($\sim$420 ms after core bounce), we do not witness significant spin up of the proto-neutron star core left behind. However, we do see the development before explosion of strong differential rotation in the turbulent gain region between the core and stalled shock. Shells in this region acquire high spin rates that reach $\sim$$150\,$ Hz, but this region contains too little mass and angular momentum to translate, even if left behind, into rapid rotation for the full neutron star. We find also that much of the induced angular momentum is likely to be ejected in the explosion, and moreover that even if the optimal amount of induced angular momentum is retained in the core, the resulting spin period is likely to be quite modest. Nevertheless, induced periods of seconds are possible.

We present new 1--1.25 micron (z and J band) Subaru/IRCS and 2 micron (K band) VLT/NaCo data for HR 8799 and a rereduction of the 3--5 micron MMT/Clio data first presented by Hinz et al. (2010). Our VLT/NaCo data yields a detection of a fourth planet at a projected separation of ~ 15 AU -- "HR 8799e". We also report new, albeit weak detections of HR 8799b at 1.03 microns and 3.3 microns. Empirical comparisons to field brown dwarfs show that at least HR 8799b and HR8799c, and possibly HR 8799d, have near-to-mid IR colors/magnitudes significantly discrepant from the L/T dwarf sequence. Standard cloud deck atmosphere models appropriate for brown dwarfs provide only (marginally) statistically meaningful fits to HR 8799b and c for unphysically small radii. Models with thicker cloud layers not present in brown dwarfs reproduce the planets' SEDs far more accurately and without the need for rescaling the planets' radii. Our preliminary modeling suggests that HR 8799b has log(g) = 4--4.5, Teff = 900K, while HR 8799c, d, and (by inference) e have log(g) = 4--4.5, Teff = 1000--1200K. Combining results from planet evolution models and new dynamical stability limits implies that the masses of HR 8799b, c, d, and e are 6--7 Mj, 7--10 Mj, 7--10 Mj and 7--10 Mj. 'Patchy" cloud prescriptions may provide even better fits to the data and may lower the estimated surface gravities and masses. Finally, contrary to some recent claims, forming the HR 8799 planets by core accretion is still plausible, although such systems are likely rare.

We perform 3+1 general relativistic simulations of rotating core collapse in the context of the collapsar model for long gamma-ray bursts. We employ a realistic progenitor, rotation based on results of stellar evolution calculations, and a simplified equation of state. Our simulations track self-consistently collapse, bounce, the postbounce phase, black hole formation, and the subsequent early hyperaccretion phase. We extract gravitational waves from the spacetime curvature and identify a unique gravitational wave signature associated with the early phase of collapsar formation.

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.

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).

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.

The collapse of a massive star's core, followed by a neutrino-driven, asymmetric supernova explosion, can naturally lead to pulsar recoils and neutron star kicks. Here, we present a two-dimensional, radiation-hydrodynamic simulation in which core collapse leads to significant acceleration of a fully-formed, nascent neutron star (NS) via an induced, neutrino-driven explosion. During the explosion, a ~10% anisotropy in the low-mass, high-velocity ejecta lead to recoil of the high-mass neutron star. At the end of our simulation, the NS has achieved a velocity of ~150 km s$^{-1}$ and is accelerating at ~350 km s$^{-2}$, but has yet to reach the ballistic regime. The recoil is due almost entirely to hydrodynamical processes, with anisotropic neutrino emission contributing less than 2% to the overall kick magnitude. Since the observed distribution of neutron star kick velocities peaks at ~300-400 km s$^{-1}$, recoil due to anisotropic core-collapse supernovae provides a natural, non-exotic mechanism with which to obtain neutron star kicks.

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.

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.

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.

The presence of a close, low-mass companion is thought to play a substantial and perhaps necessary role in shaping post-Asymptotic Giant Branch and Planetary Nebula outflows. During post-main-sequence evolution, radial expansion of the primary star, accompanied by intense winds, can significantly alter the binary orbit via tidal dissipation and mass loss. To investigate this, we couple stellar evolution models (from the zero-age main-sequence through the end of the post-main sequence) to a tidal evolution code. The binary's fate is determined by the initial masses of the primary and the companion, the initial orbit (taken to be circular), and the Reimers mass-loss parameter. For a range of these parameters, we determine whether the orbit expands due to mass loss or decays due to tidal torques. Where a common envelope (CE) phase ensues, we estimate the final orbital separation based on the energy required to unbind the envelope. These calculations predict period gaps for planetary and brown dwarf companions to white dwarfs. The upper end of the gap is the shortest period at which a CE phase is avoided. The lower end is the longest period at which companions survive their CE phase. For binary systems with 1 $M_\odot$ progenitors, we predict no Jupiter-mass companions with periods $\lesssim$270 days. Once engulfed, Jupiter-mass companions do not survive a CE phase. For binary systems consisting of a 1 $M_\odot$ progenitor with a companion 10 times the mass of Jupiter, we predict a period gap between $\sim$0.1 and $\sim$380 days. These results are consistent with both the detection of a $\sim$50 $M_{\rm J}$ brown dwarf in a $\sim$0.003 AU ($\sim$0.08 day) orbit around the white dwarf WD 0137-349 and the tentative detection of a $\sim$2 $M_{\rm J}$ planet in a $\gtrsim$2.7 AU ($\gtrsim$4 year) orbit around the white dwarf GD66.

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.

We explore the dependence on spatial dimension of the viability of the neutrino heating mechanism of core-collapse supernova explosions. We find that the tendency to explode is a monotonically increasing function of dimension, with 3D requiring $\sim$40$-$50\% lower driving neutrino luminosity than 1D and $\sim$15$-$25\% lower driving neutrino luminosity than 2D. Moreover, we find that the delay to explosion for a given neutrino luminosity is always shorter in 3D than 2D, sometimes by many hundreds of milliseconds. The magnitude of this dimensional effect is much larger than the purported magnitude of a variety of other effects, such as nuclear burning, inelastic scattering, or general relativity, which are sometimes invoked to bridge the gap between the current ambiguous and uncertain theoretical situation and the fact of robust supernova explosions. Since real supernovae occur in three dimensions, our finding may be an important step towards unraveling one of the most problematic puzzles in stellar astrophysics. In addition, even though in 3D we do see pre-explosion instabilities and blast asymmetries, unlike the situation in 2D, we do not see an obvious axially-symmetric dipolar shock oscillation. Rather, the free energy available to power instabilites seems to be shared by more and more degrees of freedom as the dimension increases. Hence, the strong dipolar axisymmetry seen in 2D and previously identified as a fundamental characteristic of the shock hydrodynamics may not survive in 3D as a prominent feature.

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.

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.

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.

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.

We present our analysis of the red-channel CoRoT data of CoRoT-2b. A deep secondary eclipse is detected at a level of 1.02+-0.20 x10-4, which suggests that all of the planet-signal detected previously in white light by Alonso et al. (2009) originates from the red channel. CoRoT-2b is the coolest exoplanet that has been detected in the optical so far. In contrast to the other planets, its measured brightness temperature of 2170+-55 K is significantly higher than its maximum hemisphere-averaged effective day-side temperature. However, it is not expected that a hot Jupiter radiates as a black body, and its thermal spectrum can deviate significantly from a Planck curve. We present models of the planet/star flux ratio as function of wavelength, which are calculated for a T/P profile in radiative and hydrostatic equilibrium, using a self-consistent atmosphere code. These are compared with the CoRoT detection, and with measurements at 4.5 and 8 um from the Spitzer Space Telescope from Gillon et al. (2009). We estimate that reflected light contributes only at a 10-20% level to the total optical eclipse depth. The models allow for an 'extra absorber' to be inserted at high altitude in the planet's atmosphere. This produces a thermal inversion layer, recently invoked to explain the photometric reversals and flux enhancements seen in some planets in the infrared. In the 0.5-1.5 um wavelength range, the model-spectra of planets with an extra absorber at high altitude, are relatively suppressed in flux compared to those without such absorber. We therefore argue that, in contrast to the other exoplanets detected in the optical so far, CoRoT-2b may not exhibit a significant thermal inversion in its atmosphere, causing its optical brightness temperature to be boosted above its maximum effective day-side temperature.

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 .

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.

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.

The determination of Saturn's atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful constraint on how they form. The helium abundance will show the extent to which helium has phase separated from hydrogen in the planet's deep interior. Jupiter's striking neon depletion may also be tied to its helium depletion, and must be confirmed or refuted in Saturn. Together with Jupiter's measured atmospheric helium abundance, a consistent evolutionary theory for both planets, including "helium rain" will be possible. We will then be able to calibrate the theory of the evolution of all giant planets, including exoplanets. In addition, high pressure H/He mixtures under giant planet conditions are an important area of condensed matter physics that are beyond the realm of experiment.

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.

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^{\gamma}$, $-1 \leqslant \gamma \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}$.

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{\omega} = 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.

Direct measurements of the spectra of extrasolar giant planets are the keys to determining their physical and chemical nature. The goal of theory is to provide the tools and context with which such data are understood. It is only by putting spectral observations through the sieve of theory that the promise of exoplanet research can be realized. With the new {\em Spitzer} and HST data of transiting "hot Jupiters," we have now dramatically entered the era of remote sensing. We are probing their atmospheric compositions and temperature profiles, are constraining their atmospheric dynamics, and are investigating their phase light curves. Soon, many non-transiting exoplanets with wide separations (analogs of Jupiter) will be imaged and their light curves and spectra measured. In this paper, we present the basic physics, chemistry, and spectroscopy necessary to model the current direct detections and to develop the more sophisticated theories for both close-in and wide-separation extrasolar giant planets that will be needed in the years to come as exoplanet research accelerates into its future.

Using the NICMOS instrument on the Hubble Space Telescope, we have measured the dayside spectrum of HD 209458b between 1.5--2.5 microns. The emergent spectrum is dominated by features due to the presence of methane (CH4) and water vapor (H2O), with smaller contributions from carbon dioxide (CO2). Combining this near-infrared spectrum with existing mid-infrared measurements shows the existence of a temperature inversion and confirms the interpretation of previous photometry measurements. We find a family of plausible solutions for the molecular abundance and detailed temperature profile. Observationally resolving the ambiguity between abundance and temperature requires either (1) improved wavelength coverage or spectral resolution of the dayside emission spectrum, or (2) a transmission spectrum where abundance determinations are less sensitive to the temperature structure.

Core-collapse supernovae are among Nature's most energetic events. They mark the end of massive star evolution and pollute the interstellar medium with the life-enabling ashes of thermonuclear burning. Despite their importance for the evolution of galaxies and life in the universe, the details of the core-collapse supernova explosion mechanism remain in the dark and pose a daunting computational challenge. We outline the multi-dimensional, multi-scale, and multi-physics nature of the core-collapse supernova problem and discuss computational strategies and requirements for its solution. Specifically, we highlight the axisymmetric (2D) radiation-MHD code VULCAN/2D and present results obtained from the first full-2D angle-dependent neutrino radiation-hydrodynamics simulations of the post-core-bounce supernova evolution. We then go on to discuss the new code Zelmani which is based on the open-source HPC Cactus framework and provides a scalable AMR approach for 3D fully general-relativistic modeling of stellar collapse, core-collapse supernovae and black hole formation on current and future massively-parallel HPC systems. We show Zelmani's scaling properties to more than 16,000 compute cores and discuss first 3D general-relativistic core-collapse results.

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.

Some transiting extrasolar giant planets have measured radii larger than predicted by the standard theory. In this paper, we explore the possibility that an earlier episode of tidal heating can explain such radius anomalies and apply the formalism we develop to HD 209458b as an example. We find that for strong enough tides the planet's radius can undergo a transient phase of inflation that temporarily interrupts canonical, monotonic shrinking due to radiative losses. Importantly, an earlier episode of tidal heating can result in a planet with an inflated radius, even though its orbit has nearly circularized. Moreover, we confirm that at late times, and under some circumstances, by raising tides on the star itself a planet can spiral into its host. We note that a 3$\times$ to 10$\times$solar planet atmospheric opacity with no tidal heating is sufficient to explain the observed radius of HD 209458b. However, our model demonstrates that with an earlier phase of episodic tidal heating we can fit the observed radius of HD 209458b even with lower (solar) atmospheric opacities. This work demonstrates that, if a planet is left with an appreciable eccentricity after early inward migration and/or dynamical interaction, coupling radius and orbit evolution in a consistent fashion that includes tidal heating, stellar irradiation, and detailed model atmospheres might offer a generic solution to the inflated radius puzzle for transiting extrasolar giant planets such as WASP-12b, TrES-4, and WASP-6b.

Spitzer Space Telescope infrared observations indicate that several transiting extrasolar giant planets have thermal inversions in their upper atmospheres. Above a relative minimum, the temperature appears to increase with altitude. Such an inversion probably requires a species at high altitude that absorbs a significant amount of incident optical/UV radiation. Some authors have suggested that the strong optical absorbers titanium oxide (TiO) and vanadium oxide (VO) could provide the needed additional opacity. However, if regions of the atmosphere are cold enough for Ti and V to be sequestered into solids they might rain out and be severely depleted. With a model of the vertical distribution of a refractory species in gaseous and condensed form, we address the question of whether enough TiO (or VO) could survive aloft in an irradiated planet's atmosphere to produce a thermal inversion. We find that, without significant macroscopic mixing, a heavy species such as TiO -- especially one that can condense in a cold-trap region -- cannot persist in a planet's upper atmosphere. In our model, the amount of macroscopic mixing that would be required to loft TiO to millibar pressures ranges from ~10^7 to ~10^{11} cm^2/s on HD 209458b, HD 149026b, TrES-4, OGLE-TR-56b, and WASP-12b, depending on the planet and on assumed condensate sizes. Such large values may be problematic for the TiO hypothesis. [Abridged]

The new spectroscopic classes, L and T, are defined by the role of dust clouds in their atmospheres, the former by their presence and the latter by their removal and near absence. Moreover, the M to L and L to T transitions are intimately tied to the condensation and character of silicate and iron grains, and the associated clouds play pivotal roles in the colors and spectra of such brown dwarfs. Spanning the effective temperature range from $\sim$2200 K to $\sim$600 K, these objects are being found in abundance and are a new arena in which condensation chemistry and the optical properties of grains is assuming astronomical importance. In this short paper, I summarize the role played by such refractories in determining the properties of these "stars" and the complexities of their theoretical treatment.

Recent observations of the extrasolar planet HD 189733b did not reveal the presence of water in the emission spectrum of the planet. Yet models of such 'Hot Jupiter' planets predict an abundance of atmospheric water vapour. Validating and constraining these models is crucial for understanding the physics and chemistry of planetary atmospheres in extreme environments. Indications of the presence of water in the atmosphere of HD 189733b have recently been found in transmission spectra, where the planet's atmosphere selectively absorbs the light of the parent star, and in broadband photometry. Here we report on the detection of strong water absorption in a high signal-to-noise, mid-infrared emission spectrum of the planet itself. We find both a strong downturn in the flux ratio below 10 microns and discrete spectral features that are characteristic of strong absorption by water vapour. The differences between these and previous observations are significant and admit the possibility that predicted planetary-scale dynamical weather structures might alter the emission spectrum over time. Models that match the observed spectrum and the broadband photometry suggest that heat distribution from the dayside to the night side is weak. Reconciling this with the high night side temperature will require a better understanding of atmospheric circulation or possible additional energy sources.

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.

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.

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.

We investigate the criteria for successful core-collapse supernova explosions by the neutrino mechanism. We find that a critical-luminosity/mass-accretion-rate condition distinguishes non-exploding from exploding models in hydrodynamic one-dimensional (1D) and two-dimensional (2D) simulations. We present 95 such simulations that parametrically explore the dependence on neutrino luminosity, mass accretion rate, resolution, and dimensionality. While radial oscillations mediate the transition between 1D accretion (non-exploding) and exploding simulations, the non-radial standing accretion shock instability characterizes 2D simulations. We find that it is useful to compare the average dwell time of matter in the gain region with the corresponding heating timescale, but that tracking the residence time distribution function of tracer particles better describes the complex flows in multi-dimensional simulations. Integral quantities such as the net heating rate, heating efficiency, and mass in the gain region decrease with time in non-exploding models, but for 2D exploding models, increase before, during, and after explosion. At the onset of explosion in 2D, the heating efficiency is $\sim$2% to $\sim$5% and the mass in the gain region is $\sim$0.005 M$_{\sun}$ to $\sim$0.01 M$_{\sun}$. Importantly, we find that the critical luminosity for explosions in 2D is $\sim$70% of the critical luminosity required in 1D. This result is not sensitive to resolution or whether the 2D computational domain is a quadrant or the full 180$^{\circ}$. We suggest that the relaxation of the explosion condition in going from 1D to 2D (and to, perhaps, 3D) is of a general character and is not limited by the parametric nature of this study.

To explain their observed radii, we present theoretical radius-age trajectories for the extrasolar giant planets (EGPs) TrES-4, XO-3b, and HAT-P-1b. We factor in variations in atmospheric opacity, the presence of an inner heavy-element core, and possible heating due to orbital tidal dissipation. A small, yet non-zero, degree of core heating is needed to explain the observed radius of TrES-4, unless its atmospheric opacity is significantly larger than a value equivalent to that at 10$\times$solar metallicity with equilibrium molecular abundances. This heating rate is reasonable, and corresponds for an energy dissipation parameter ($Q_p$) of $\sim10^{3.8}$ to an eccentricity of $\sim$0.01, assuming 3$\times$solar atmospheric opacity and a heavy-element core of $M_c = 30$ $M_{\oplus}$. For XO-3b, which has an observed orbital eccentricity of 0.26, we show that tidal heating needs to be taken into account to explain its observed radius. Furthermore, we reexamine the core mass needed for HAT-P-1b in light of new measurements and find that it now generally follows the correlation between stellar metallicity and core mass suggested recently. Given various core heating rates, theoretical grids and fitting formulae for a giant planet's equilibrium radius and equilibration timescale are provided for planet masses $M_p=$ 0.5, 1.0, and 1.5 $M_J$ with $a =$ 0.02-0.06 AU, orbiting a G2V star. When the equilibration timescale is much shorter than that of tidal heating variation, the ``effective age'' of the planet is shortened, resulting in evolutionary trajectories more like those of younger EGPs. Motivated by the work of \citet{jackson08a,jackson08b}, we suggest that this effect could indeed be important in better explaining some observed transit radii.

We show that a consistent fit to observed secondary eclipse data for several strongly irradiated transiting planets demands a temperature inversion (stratosphere) at altitude. Such a thermal inversion significantly influences the planet/star contrast ratios at the secondary eclipse,their wavelength dependences, and, importantly, the day-night flux contrast during a planetary orbit. The presence of the thermal inversion/stratosphere seems to roughly correlate with the stellar flux at the planet. Such temperature inversions might caused by an upper-atmosphere absorber whose exact nature is still uncertain

In this paper, we describe a new hydrodynamics code for 1D and 2D astrophysical simulations, BETHE-hydro, that uses time-dependent, arbitrary, unstructured grids. The core of the hydrodynamics algorithm is an arbitrary Lagrangian-Eulerian (ALE) approach, in which the gradient and divergence operators are made compatible using the support-operator method. We present 1D and 2D gravity solvers that are finite differenced using the support-operator technique, and the resulting system of linear equations are solved using the tridiagonal method for 1D simulations and an iterative multigrid-preconditioned conjugate-gradient method for 2D simulations. Rotational terms are included for 2D calculations using cylindrical coordinates. We document an incompatibility between a subcell pressure algorithm to suppress hourglass motions and the subcell remapping algorithm and present a modified subcell pressure scheme that avoids this problem. Strengths of this code include a straightforward structure, enabling simple inclusion of additional physics packages, the ability to use a general equation of state, and most importantly, the ability to solve self-gravitating hydrodynamic flows on time-dependent, arbitrary grids. In what follows, we describe in detail the numerical techniques employed and, with a large suite of tests, demonstrate that BETHE-hydro finds accurate solutions with 2$^{nd}$-order convergence.

We present VULCAN/2D multi-group flux-limited-diffusion radiation hydrodynamics simulations of binary neutron star (BNS) mergers, using the Shen equation of state, covering ~100 ms, and starting from azimuthal-averaged 2D slices obtained from 3D SPH simulations of Rosswog & Price for 1.4 Msun (baryonic) neutron stars with no initial spins, co-rotating spins, and counter-rotating spins. Snapshots are post-processed at 10 ms intervals with a multi-angle neutrino-transport solver. We find polar-enhanced neutrino luminosities, dominated by $\bar{\nu}_e$ and ``$\nu_\mu$'' neutrinos at peak, although $\nu_e$ emission may be stronger at late times. We obtain typical peak neutrino energies for $\nu_e$, $\bar{\nu}_e$, and ``$\nu_\mu$'' of ~12, ~16, and ~22 MeV. The super-massive neutron star (SMNS) formed from the merger has a cooling timescale of ~1 s. Charge-current neutrino reactions lead to the formation of a thermally-driven bipolar wind with <$\dot{M}$> ~10$^{-3}$ Msun/s, baryon-loading the polar regions, and preventing any production of a GRB prior to black-hole formation. The large budget of rotational free energy suggests magneto-rotational effects could produce a much greater polar mass loss. We estimate that ~10$^{-4}$ Msun of material with electron fraction in the range 0.1-0.2 become unbound during this SMNS phase as a result of neutrino heating. We present a new formalism to compute the $\nu_i\bar{\nu}_i$ annihilation rate based on moments of the neutrino specific intensity computed with our multi-angle solver. Cumulative annihilation rates, which decay as $t^{-1.8}$, decrease over our 100 ms window from a few 10$^{50}$ to ~10$^{49}$ erg/s, equivalent to a few 10$^{54}$ to ~10$^{53}$ $e^-e^+$ pairs per second.

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.

We estimate flux ratios of the extrasolar planet XO-1b to its host star XO-1 at 3.6, 4.5, 5.8 and 8.0 microns with the IRAC on the Spitzer Space Telescope to be 0.00086 +/- 0.00007, 0.00122 +/- 0.00009, 0.00261 +/- 0.00031 and 0.00210 +/- 0.00029, respectively. The fluxes are inconsistent with a canonical cloudless model for the thermal emission from a planet and suggest an atmosphere with a thermal inversion layer and a possible stratospheric absorber. A newly emerging correlation between the presence of a thermal inversion layer in the planetary atmosphere and stellar insolation of the planet (Burrows et al. 2007b) is refined. The sub-stellar point flux from the parent star at XO-1b of ~ 0.49 x 10^9 erg cm^-2 s^-1 sets a new lower limit for the occurrence of a thermal inversion in a planetary atmosphere.

Massive stars end their short lives in spectacular explosions, supernovae, that synthesize new elements and drive galaxy evolution. Throughout history supernovae were discovered chiefly through their delayed optical light, preventing observations in the first moments (hours to days) following the explosion. As a result, the progenitors of some supernovae and the events leading up to their violent demise remain intensely debated. Here we report the serendipitous discovery of a supernova at the time of explosion, marked by an extremely luminous X-ray outburst. We attribute the outburst to the break-out of the supernova shock-wave from the progenitor, and show that the inferred rate of such events agrees with that of all core-collapse supernovae. We forecast that future wide-field X-ray surveys will catch hundreds of supernovae each year in the act of explosion, and thereby enable crucial neutrino and gravitational wave detections that may ultimately unravel the explosion mechanism.

We calculate a new suite of albedo models for close-in extrasolar giant planets and compare with the recent stringent upper limit for HD 209458b of Rowe et al. using MOST. We find that all models without scattering clouds are consistent with this optical limit. We explore the dependence on wavelength and waveband, metallicity, the degree of heat redistribution, and the possible presence of thermal inversions and find a rich diversity of behaviors. Measurements of transiting extrasolar giant planets (EGPs) at short wavelengths by MOST, Kepler, and CoRoT, as well as by proposed dedicated multi-band missions, can complement measurements in the near- and mid-IR using {\it Spitzer} and JWST. Collectively, such measurements can help determine metallicity, compositions, atmospheric temperatures, and the cause of thermal inversions (when they arise) for EGPs with a broad range of radii, masses, degrees of stellar insolation, and ages. With this paper, we reappraise and highlight the diagnostic potential of albedo measurements of hot EGPs shortward of $\sim$1.3 $\mu$m.

We present theoretical atmosphere, spectral, and light-curve models for extrasolar giant planets (EGPs) undergoing strong irradiation for which {\it Spitzer} planet/star contrast ratios or light curves have been published (circa June 2007). These include HD 209458b, HD 189733b, TrES-1, HD 149026b, HD 179949b, and $\upsilon$ And b. By comparing models with data, we find that a number of EGP atmospheres experience thermal inversions and have stratospheres. This is particularly true for HD 209458b, HD 149026b, and $\upsilon$ And b. This finding translates into qualitative changes in the planet/star contrast ratios at secondary eclipse and in close-in EGP orbital light curves. Moreover, the presence of atmospheric water in abundance is fully consistent with all the {\it Spitzer} data for the measured planets. For planets with stratospheres, water absorption features invert into emission features and mid-infrared fluxes can be enhanced by a factor of two. In addition, the character of near-infrared planetary spectra can be radically altered. We derive a correlation between the importance of such stratospheres and the stellar flux on the planet, suggesting that close-in EGPs bifurcate into two groups: those with and without stratospheres. From the finding that TrES-1 shows no signs of a stratosphere, while HD 209458b does, we estimate the magnitude of this stellar flux breakpoint. We find that the heat redistribution parameter, P$_n$, for the family of close-in EGPs assumes values from $\sim$0.1 to $\sim$0.4. This paper provides a broad theoretical context for the future direct characterization of EGPs in tight orbits around their illuminating stars.

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.

Recent stellar evolutionary calculations of low-metallicity massive fast-rotating main-sequence stars yield iron cores at collapse endowed with high angular momentum. It is thought that high angular momentum and black hole formation are critical ingredients of the collapsar model of long-soft gamma-ray bursts (GRBs). Here, we present 2D multi-group, flux-limited-diffusion MHD simulations of the collapse, bounce, and immediate post-bounce phases of a 35-Msun collapsar-candidate model of Woosley & Heger. We find that, provided the magneto-rotational instability (MRI) operates in the differentially-rotating surface layers of the millisecond-period neutron star, a magnetically-driven explosion ensues during the proto-neutron star phase, in the form of a baryon-loaded non-relativistic jet, and that a black hole, central to the collapsar model, does not form. Paradoxically, and although much uncertainty surrounds stellar mass loss, angular momentum transport, magnetic fields, and the MRI, current models of chemically homogeneous evolution at low metallicity yield massive stars with iron cores that may have too much angular momentum to avoid a magnetically-driven, hypernova-like, explosion in the immediate post-bounce phase. We surmise that fast rotation in the iron core may inhibit, rather than enable, collapsar formation, which requires a large angular momentum not in the core but above it. Variations in the angular momentum distribution of massive stars at core collapse might explain both the diversity of Type Ic supernovae/hypernovae and their possible association with a GRB. A corollary might be that, rather than the progenitor mass, the angular momentum distribution, through its effect on magnetic field amplification, distinguishes these outcomes.

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.

We find that a theoretical fit to all the HD 209458b data at secondary eclipse requires that the dayside atmosphere of HD 209458b have a thermal inversion and a stratosphere. This inversion is caused by the capture of optical stellar flux by an absorber of uncertain origin that resides at altitude. One consequence of stratospheric heating and temperature inversion is the flipping of water absorption features into emission features from the near- to the mid-infrared and we see evidence of such a water emission feature in the recent HD 209458b IRAC data of Knutson et al. In addition, an upper-atmosphere optical absorber may help explain both the weaker-than-expected Na D feature seen in transit and the fact that the transit radius at 24 $\mu$m is smaller than the corresponding radius in the optical. Moreover, it may be a factor in why HD 209458b's optical transit radius is as large as it is. We speculate on the nature of this absorber and the planets whose atmospheres may, or may not, be affected by its presence.

We present a systematic study of the spectral consequences of departures from chemical equilibrium in the atmospheres of L and T dwarfs, and for even cooler dwarfs. The temperature/pressure profiles of the non-equilibrium models are fully consistent with the non-equilibrium chemistry. Our grid of non-equilibrium models includes spectra for effective temperatures from 200 K to 1800 K, three surface gravities, four possible values of the coefficient of eddy diffusion in the radiative zone, and three different CO/CH$_4$ chemical reaction prescriptions. We find that the non-equilibrium overabundance of CO translates into flux suppressions in the M ($\sim4-$5 $\mu$m) band of at most $\sim$40% between effective temperatures of 600 and 1800 K. The effect is largest around $T_{\rm eff} \approx 1100$ K. The underabundance of ammonia due to non-equilibrium chemistry translates into flux enhancements of no more than $\sim$20% for the $T_{\rm eff}$ range from 300 to 1800 K, with the largest effects at the lowest values of $T_{\rm eff}$. The magnitude of the departure from chemical equilibrium increases with decreasing gravity, with increasing eddy diffusion coefficient, and with decreasing speed of the CO/CH$_4$ reaction. Though these effects are modest, they lead to better fits with the measured T dwarf spectra. Furthermore, the suppression in the M band due to non-equilibrium enhancements in the CO abundance disappears below $\sim$500 K, and is only partial above $\sim$500 K, preserving the M band flux as a useful diagnostic of cool atmospheres and maintaining its importance for searches for brown dwarfs cooler than T dwarfs.

We present 2D multi-group flux-limited diffusion magnetohydrodynamics (MHD) simulations of the Accretion-Induced Collapse (AIC) of a rapidly-rotating white dwarf. We focus on the dynamical role of MHD processes after the formation of a millisecond-period protoneutron star. We find that including magnetic fields and stresses can lead to a powerful explosion with an energy of a few Bethe, rather than a weak one of at most 0.1 Bethe, with an associated ejecta mass of ~0.1Msun, instead of a few 0.001Msun. The core is spun down by ~30% within 500ms after bounce, and the rotational energy extracted from the core is channeled into magnetic energy that generates a strong magnetically-driven wind, rather than a weak neutrino-driven wind. Baryon loading of the ejecta, while this wind prevails, precludes it from becoming relativistic. This suggests that a GRB is not expected to emerge from such AICs during the early protoneutron star phase, except in the unlikely event that the massive white dwarf has sufficient mass to lead to black hole formation. In addition, we predict both negligible 56Ni-production (that should result in an optically-dark, adiabatically-cooled explosion) and the ejection of 0.1Msun of material with an electron fraction of 0.1-0.2. Such pollution by neutron-rich nuclei puts strong constraints on the possible rate of such AICs. Moreover, being free from ``fallback,'' such highly-magnetized millisecond-period protoneutron stars may later become magnetars, and the magnetically-driven winds may later transition to Poynting-flux-dominated, relativistic winds, eventually detectable as GRBs at cosmological distances. However, the low expected event rate of AICs will constrain them to be, at best, a small subset of GRB and/or magnetar progenitors.

The recently discovered star SCR 1845-6357 is the first late M/T dwarf binary discovered. SCR 1845 is a particular object due to its tight orbit (currently around 4 AU) and its proximity to the Sun (3.85 pc). We present spatially resolved VLT/NACO images and low resolution spectra of SCR 1845 in the J, H and K near-infrared bands. Since the T dwarf companion, SCR 1845B, is so close to the primary SCR 1845A, orbital motion is evident even within a year. Following the orbital motion, the binary's mass can be measured accurately within a decade, making SCR 1845B a key T-dwarf mass-luminosity calibrator. The NIR spectra allow for accurate determination of spectral type and also for rough estimates of the object's physical parameters. The spectral type of SCR 1845B is determined by direct comparison of the flux calibrated JHK spectra with T dwarf standard template spectra and also by NIR spectral indices obtained from synthetic photometry. Constrained values for surface gravity, effective temperature and metallicity are derived by comparison with model spectra. Our data prove that SCR 1845B is a brown dwarf of spectral type T6 that is co-moving with and therefore gravitationally bound to the M8.5 primary. Fitting the NIR spectrum of SCR 1845B to model spectra yields an effective temperature of about 950K and a surface gravity log(g)=5.1 (cgs) assuming solar metallicity. Mass and age of SCR 1845B are in the range 40 to 50 Jupiter masses and 1.8 to 3.1 Gyr.

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.

We calculate the theoretical evolution of the radii of all fourteen of the known transiting extrasolar giant planets (EGPs) for a variety of assumptions concerning atmospheric opacity, dense inner core masses, and possible internal power sources. We incorporate the effects of stellar irradiation and customize such effects for each EGP and star. Looking collectively at the family as a whole, we find that there are in fact two radius anomalies to be explained. Not only are the radii of a subset of the known transiting EGPs larger than expected from previous theory, but many of the other objects are smaller than the default theory would allow. We suggest that the larger EGPs can be explained by invoking enhanced atmospheric opacities that naturally retain internal heat. This explanation might obviate the necessity for an extra internal power source. We explain the smaller radii by the presence in perhaps all the known transiting EGPs of dense cores, such as have been inferred for Saturn and Jupiter. Importantly, we derive a rough correlation between the masses of our "best-fit" cores and the stellar metallicity that seems to buttress the core-accretion model of their formation. Though many caveats and uncertainties remain, the resulting comprehensive theory that incorporates enhanced-opacity atmospheres and dense cores is in reasonable accord with all the current structural data for the known transiting giant planets.

We have created a general methodology for calculating the wavelength-dependent light curves of close-in extrasolar giant planets (EGPs) as they traverse their orbits. Focussing on the transiting EGPs HD189733b, TrES-1, and HD209458b, we calculate planet/star flux ratios during secondary eclipse and compare them with the Spitzer data points obtained so far in the mid-infrared. We introduce a simple parametrization for the redistribution of heat to the planet's nightside, derive constraints on this parameter (P_n), and provide a general set of predictions for planet/star contrast ratios as a function of wavelength, model, and phase. Moreover, we calculate average dayside and nightside atmospheric temperature/pressure profiles for each transiting planet/P_n pair with which existing and anticipated Spitzer data can be used to probe the atmospheric thermal structure of severely irradiated EGPs. We find that the baseline models do a good job of fitting the current secondary eclipse dataset, but that the Spitzer error bars are not yet small enough to discriminate cleanly between all the various possibilities.

We report on the measurement of the 7.5-14.7 micron spectrum for the transiting extrasolar giant planet HD 189733b using the Infrared Spectrograph on the Spitzer Space Telescope. Though the observations comprise only 12 hours of telescope time, the continuum is well measured and has a flux ranging from 0.6 mJy to 1.8 mJy over the wavelength range, or 0.49 +/- 0.02% of the flux of the parent star. The variation in the measured fractional flux is very nearly flat over the entire wavelength range and shows no indication of significant absorption by water or methane, in contrast with the predictions of most atmospheric models. Models with strong day/night differences appear to be disfavored by the data, suggesting that heat redistribution to the night side of the planet is highly efficient.

We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which already includes radiation-hydrodynamics in various approximations and can be used with arbitrarily moving meshes (ALE). The MHD scheme, which consists of cell-centered magnetic field variables, preserves the nodal finite difference representation of $div(\bB)$ by construction, and therefore any initially divergence-free field remains divergence-free through the simulation. In this paper, we describe the new scheme in detail and present comparisons of VULCAN/2D results with those of the code ZEUS/2D for several one-dimensional and two-dimensional test problems. The code now enables two-dimensional simulations of the collapse and explosion of the rotating, magnetic cores of massive stars. Moreover, it can be used to simulate the very wide variety of astrophysical problems for which multi-D radiation-magnetohydrodynamics (RMHD) is relevant.

The observed mass function for all known extrasolar giant planets (EGPs) varies approximately as M^{-1} for mass M between 0.2 Jupiter masses (M_J) and 5 M_J. In order to study evaporation effects for highly-irradiated EGPs in this mass range, we have constructed an observational mass function for a subset of EGPs in the same mass range but with orbital radii <0.07 AU. Surprisingly, the mass function for such highly-irradiated EGPs agrees quantitatively with the M^{-1} law, implying that the mass function for EGPs is preserved despite migration to small orbital radii. Unless there is a remarkable compensation of mass-dependent orbital migration for mass-dependent evaporation, this result places a constraint on orbital migration models and rules out the most extreme mass loss rates in the literature. A theory that predicts more moderate mass loss gives a mass function that is closer to observed statistics but still disagrees for M < 1 M_J.

In this paper, we bring together various of our published and unpublished findings from our recent 2D multi-group, flux-limited radiation hydrodynamic simulations of the collapse and explosion of the cores of massive stars. Aided by 2D and 3D graphical renditions, we motivate the acoustic mechanism of core-collapse supernova explosions and explain, as best we currently can, the phases and phenomena that attend this mechanism. Two major foci of our presentation are the outer shock instability and the inner core g-mode oscillations. The former sets the stage for the latter, which damp by the generation of sound. This sound propagates outward to energize the explosion and is relevant only if the core has not exploded earlier by some other means. Hence, it is a more delayed mechanism than the traditional neutrino mechanism that has been studied for the last twenty years since it was championed by Bethe and Wilson. We discuss protoneutron star convection, accretion-induced-collapse, gravitational wave emissions, pulsar kicks, the angular anisotropy of the neutrino emissions, a subset of numerical issues, and a new code we are designing that should supercede our current supernova code VULCAN/2D. Whatever ideas last from this current generation of numerical results, and whatever the eventual mechanism(s), we conclude that the breaking of spherical symmetry will survive as one of the crucial keys to the supernova puzzle.

We present mid-infrared photometry of 124 white dwarf stars with Spitzer Space Telescope. Objects were observed simultaneously at 4.5 and 8.0um with sensitivities better than 1 mJy. This data can be used to test models of white dwarf atmospheres in a new wavelength regime, as well as to search for planetary companions and debris disks.

We present calculations for the evolution and surviving mass of highly-irradiated extrasolar giant planets (EGPs) at orbital semimajor axes ranging from 0.023 to 0.057 AU using a generalized scaled theory for mass loss, together with new surface-condition grids for hot EGPs and a consistent treatment of tidal truncation. Theoretical estimates for the rate of energy-limited hydrogen escape from giant-planet atmospheres differ by two orders of magnitude, when one holds planetary mass, composition, and irradiation constant. Baraffe et al. (2004, A&A 419, L13-L16) predict the highest rate, based on the theory of Lammer et al. (2003, Astrophys. J. 598, L121-L124). Scaling the theory of Watson et al. (1981, Icarus 48, 150-166) to parameters for a highly-irradiated exoplanet, we find an escape rate ~102 lower than Baraffe's. With the scaled Watson theory we find modest mass loss, occurring early in the history of a hot EGP. In this theory, mass loss including the effect of Roche-lobe overflow becomes significant primarily for masses below a Saturn mass, for semimajor axes = 0.023 AU. This contrasts with the Baraffe model, where hot EGPs are claimed to be remnants of much more massive bodies, originally several times Jupiter and still losing substantial mass fractions at present.

In the context of 2D, axisymmetric, multi-group, radiation/hydrodynamic simulations of core-collapse supernovae over the full 180$^{\circ}$ domain, we present an exploration of the progenitor dependence of the acoustic mechanism of explosion. All progenitor models we have tested with our Newtonian code explode. We investigate the roles of the Standing-Accretion-Shock-Instability (SASI), the excitation of core g-modes, the generation of core acoustic power, the ejection of matter with r-process potential, the wind-like character of the explosion, and the fundamental anisotropy of the blasts. We find that the breaking of spherical symmetry is central to the supernova phenomenon and the blasts, when top-bottom asymmetric, are self-collimating. We see indications that the initial explosion energies are larger for the more massive progenitors, and smaller for the less massive progenitors, and that the neutrino contribution to the explosion energy may be an increasing function of progenitor mass. The degree of explosion asymmetry we obtain is completely consistent with that inferred from the polarization measurements of Type Ic supernovae. Furthermore, we calculate for the first time the magnitude and sign of the net impulse on the core due to anisotropic neutrino emission and suggest that hydrodynamic and neutrino recoils in the context of our asymmetric explosions afford a natural mechanism for observed pulsar proper motions. [abridged]

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.

We present a comprehensive description of the theory and practice of opacity calculations from the infrared to the ultraviolet needed to generate models of the atmospheres of brown dwarfs and extrasolar giant planets. Methods for using existing line lists and spectroscopic databases in disparate formats are presented and plots of the resulting absorptive opacities versus wavelength for the most important molecules and atoms at representative temperature/pressure points are provided. Electronic, ro-vibrational, bound-free, bound-bound, free-free, and collision-induced transitions and monochromatic opacities are derived, discussed, and analyzed. The species addressed include the alkali metals, iron, heavy metal oxides, metal hydrides, $H_2$, $H_2O$, $CH_4$, $CO$, $NH_3$, $H_2S$, $PH_3$, and representative grains. [Abridged]

We present the results of a program to acquire photometry for eighty-six late-M, L, and T dwarfs using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We examine the behavior of these cool dwarfs in various color-color and color-magnitude diagrams composed of near-IR and IRAC data. The T dwarfs exhibit the most distinctive positions in these diagrams. In M_5.8 versus [5.8]-[8.0], the IRAC data for T dwarfs are not monotonic in either magnitude or color, giving the clearest indication yet that the T dwarfs are not a one parameter family in Teff. Because metallicity does not vary enough in the solar neighborhood to act as the second parameter, the most likely candidate then is gravity, which in turn translates to mass. Among objects with similar spectral type, the range of mass suggested by our sample is about a factor of five (~70 M_Jup to ~15 M_Jup), with the less massive objects making up the younger members of the sample. We also find the IRAC 4.5 micron fluxes to be lower than expected, from which we infer a stronger CO fundamental band at ~4.67 microns. This suggests that equilibrium CH_4/CO chemistry underestimates the abundance of CO in T dwarf atmospheres, confirming earlier results based on M-band observations from the ground. In combining IRAC photometry with near-IR JHK photometry and parallax data, we find the combination of Ks, IRAC 3.6 micron, and 4.5 micron bands to provide the best color-color discrimination for a wide range of M, L, and T dwarfs. Also noteworthy is the M_Ks versus Ks-[4.5] relation, which shows a smooth progression over spectral type and splits the M, L, and T types cleanly.

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.

We present the results of a Hubble Space Telescope NICMOS imaging survey of 22 T-type field brown dwarfs. Five are resolved as binary systems with angular separations of 0"05-0"35, and companionship is established on the basis of component F110W-F170M colors (indicative of CH4 absorption) and low probabilities of background contamination. Prior ground-based observations show 2MASS 1553+1532AB to be a common proper motion binary. The properties of these systems - low multiplicity fraction (11[+7][-3]% resolved, as corrected for sample selection baises), close projected separations (a = 1.8-5.0 AU) and near-unity mass ratios - are consistent with previous results for field brown dwarf binaries. Three of the binaries have components that span the poorly-understood transition between L dwarfs and T dwarfs. Spectral decomposition analysis of one of these, SDSS 1021-0304AB, reveals a peculiar flux reversal between its components, as its T5 secondary is ~30% brighter at 1.05 and 1.27 micron than its T1 primary. This system, 2MASS 0518-2828AB and SDSS 1534+1615AB all demonstrate that the J-band brightening observed between late-type L to mid-type T dwarfs is an intrinsic feature of this spectral transition, albeit less pronounced than previously surmised. We also find that the resolved binary fraction of L7 to T3.5 dwarfs is twice that of other L and T dwarfs, an anomaly that can be explained by a relatively rapid evolution of brown dwarfs through the L/T transition, perhaps driven by dynamic (nonequilibrium) depletion of photospheric condensates.

We present a new theory for the gravitational wave signatures of core-collapse supernovae. Previous studies identified axisymmetric rotating core collapse, core bounce, postbounce convection, and anisotropic neutrino emission as the primary processes and phases for the radiation of gravitational waves. Our results, which are based on axisymmetric, Newtonian radiation-hydrodynamics supernova simulations (Burrows et al. 2006), indicate that the dominant emission process of gravitational waves in core-collapse supernovae may be the oscillations of the protoneutron star core. The oscillations are predominantly of g-mode character, are excited hundreds of milliseconds after bounce, and typically last for several hundred milliseconds. Our results suggest that even nonrotating core-collapse supernovae should be visible to current LIGO-class detectors throughout the Galaxy, and depending on progenitor structure, possibly out to Megaparsec distances.

We present 2.5D radiation-hydrodynamics simulations of the accretion-induced collapse (AIC) of white dwarfs, starting from 2D rotational equilibrium configurations of a 1.46-Msun and a 1.92-Msun model. Electron capture leads to the collapse to nuclear densities of these cores within a few tens of milliseconds. The shock generated at bounce moves slowly, but steadily, outwards. Within 50-100ms, the stalled shock breaks out of the white dwarf along the poles. The blast is followed by a neutrino-driven wind that develops within the white dwarf, in a cone of ~40deg opening angle about the poles, with a mass loss rate of 5-8 x 10^{-3} Msun/yr. The ejecta have an entropy on the order of 20-50 k_B/baryon, and an electron fraction distribution that is bimodal. By the end of the simulations, at >600ms after bounce, the explosion energy has reached 3-4 x 10^{49}erg and the total ejecta mass has reached a few times 0.001Msun. We estimate the asymptotic explosion energies to be lower than 10^{50}erg, significantly lower than those inferred for standard core collapse. The AIC of white dwarfs thus represents one instance where a neutrino mechanism leads undoubtedly to a successful, albeit weak, explosion. We document in detail the numerous effects of the fast rotation of the progenitors: The neutron stars are aspherical; the ``nu_mu'' and anti-nu_e neutrino luminosities are reduced compared to the nu_e neutrino luminosity; the deleptonized region has a butterfly shape; the neutrino flux and electron fraction depend strongly upon latitude (a la von Zeipel); and a quasi-Keplerian 0.1-0.5-Msun accretion disk is formed.

Based on multi-dimensional, multi-group, flux-limited-diffusion hydrodynamic simulations of core-collapse supernovae with the VULCAN/2D code, we study the physical conditions within and in the vicinity of the nascent protoneutron star (PNS). Our numerical study follows the evolution of the collapsing envelope of the 11-Msun model of Woosley & Weaver, from ~200 ms before bounce to ~300ms after bounce, on a spatial grid that switches from Cartesian at the PNS center to spherical above a 10km radius. As has been shown previously, we do not see any large-scale overturn of the inner PNS material. Convection, directly connected to the PNS, is found to occur in two distinct regions: between 10 and 20km, coincident with the region of negative lepton gradient, and exterior to the PNS above 50km. Separating these two regions, an interface, with no sizable outward or inward motion, is the site of gravity waves, emerging at ~200-300ms after core bounce, excited by the convection in the outer convective zone. In the PNS, convection is always confined within the neutrinospheric radii for all neutrino energies above just a few MeV. We find that such convective motions do not appreciably enhance the nu_e neutrino luminosity, and that they enhance the anti-nu_e and "nu_mu" luminosities modestly, by ~15% and ~30%, respectively, during the first post-bounce 100--200ms. Moreover, we see no evidence of doubly-diffusive instabilities in the PNS, expected to operate on diffusion timescales of at least a second, much longer than the millisecond timescale associated with PNS convection. PNS convection is thus found to be a secondary feature of the core-collapse phenomenon, rather than a decisive ingredient for a successful explosion.

When extrasolar planets are observed to transit their parent stars, we are granted unprecedented access to their physical properties. It is only for transiting planets that we are permitted direct estimates of the planetary masses and radii, which provide the fundamental constraints on models of their physical structure. In particular, precise determination of the radius may indicate the presence (or absence) of a core of solid material, which in turn would speak to the canonical formation model of gas accretion onto a core of ice and rock embedded in a protoplanetary disk. Furthermore, the radii of planets in close proximity to their stars are affected by tidal effects and the intense stellar radiation. As a result, some of these "hot Jupiters" are significantly larger than Jupiter in radius. Precision follow-up studies of such objects (notably with the space-based platforms of the Hubble and Spitzer Space Telescopes) have enabled direct observation of their transmission spectra and emitted radiation. These data provide the first observational constraints on atmospheric models of these extrasolar gas giants, and permit a direct comparison with the gas giants of the Solar system. Despite significant observational challenges, numerous transit surveys and quick-look radial velocity surveys are active, and promise to deliver an ever-increasing number of these precious objects. The detection of transits of short-period Neptune-sized objects, whose existence was recently uncovered by the radial-velocity surveys, is eagerly anticipated. Ultra-precise photometry enabled by upcoming space missions offers the prospect of the first detection of an extrasolar Earth-like planet in the habitable zone of its parent star, just in time for Protostars and Planets VI.

We present new observations of the white dwarf G 29-38 with the camera (4.5 and 8 microns), photometer (24 microns), and spectrograph (5.5-14 microns) of the Spitzer Space Telescope. This star has an exceptionally large infrared excess amounting to 3% of the bolometric luminosity. The spectral energy distribution has a continuum peak around 4.5 micros and a 9-11 micron emission feature 1.25 times brighter than the continuum. A mixture of amorphous olivine and a small amount of forsterite in an emitting region 1-5 Rsun from the star can reproduce the shape of the 9-11 micron feature. The spectral energy distribution also appears to require amorphous carbon to explain the hot continuum. Our new measurements support the idea that a relatively recent disruption of a comet or asteroid created the cloud.

We present 4.5$\mu$m and 8$\mu$m photometric observations of 18 cool white dwarfs obtained with the Spitzer Space Telescope. Our observations demonstrate that four white dwarfs with T_eff< 6000 K show slightly depressed mid-infrared fluxes relative to white dwarf models. In addition, another white dwarf with a peculiar optical and near-infrared spectral energy distribution (LHS 1126) is found to display significant flux deficits in Spitzer observations. These mid-infrared flux deficits are not predicted by the current white dwarf models including collision induced absorption due to molecular hydrogen. We postulate that either the collision induced absorption calculations are incomplete or there are other unrecognized physical processes occuring in cool white dwarf atmospheres. The spectral energy distribution of LHS 1126 surprisingly fits a Rayleigh-Jeans spectrum in the infrared, mimicking a hot white dwarf with effective temperature well in excess of 10$^5$ K. This implies that the source of this flux deficit is probably not molecular absorption but some other process.

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.

Using a model for refractory clouds, a novel algorithm for handling them, and the latest gas-phase molecular opacities, we have produced a new series of L and T dwarf spectral and atmosphere models as a function of gravity and metallicity, spanning the \teff range from 2200 K to 700 K. The correspondence with observed spectra and infrared colors for early- and mid-L dwarfs and for mid- to late-T dwarfs is good. We find that the width in infrared color-magnitude diagrams of both the T and L dwarf branches is naturally explained by reasonable variations in gravity and, therefore, that gravity is the "second parameter" of the L/T dwarf sequence. We investigate the dependence of theoretical dwarf spectra and color-magnitude diagrams upon various cloud properties, such as particle size and cloud spatial distribution. In the region of the L$\to$T transition, we find that no one cloud-particle-size and gravity combination can be made to fit all the observed data. Furthermore, we note that the new, lower solar oxygen abundances of Allende-Prieto, Lambert, & Asplund (2002) produce better fits to brown dwarf data than do the older values. Finally, we discuss various issues in cloud physics and modeling and speculate on how a better correspondence between theory and observation in the problematic L$\to$T transition region might be achieved.

We present a new mechanism for core-collapse supernova explosions that relies upon acoustic power generated in the inner core as the driver. In our simulation using an 11-solar-mass progenitor, a strong advective-acoustic oscillation a la Foglizzo with a period of ~25-30 milliseconds (ms) arises ~200 ms after bounce. Its growth saturates due to the generation of secondary shocks, and kinks in the resulting shock structure funnel and regulate subsequent accretion onto the inner core. However, this instability is not the primary agent of explosion. Rather, it is the acoustic power generated in the inner turbulent region and most importantly by the excitation and sonic damping of core g-mode oscillations. An l=1 mode with a period of ~3 ms grows to be prominent around ~500 ms after bounce. The accreting protoneutron star is a self-excited oscillator. The associated acoustic power seen in our 11-solar-mass simulation is sufficient to drive the explosion. The angular distribution of the emitted sound is fundamentally aspherical. The sound pulses radiated from the core steepen into shock waves that merge as they propagate into the outer mantle and deposit their energy and momentum with high efficiency. The core oscillation acts like a transducer to convert accretion energy into sound. An advantage of the acoustic mechanism is that acoustic power does not abate until accretion subsides, so that it is available as long as it may be needed to explode the star. [abridged]

We present a method for measuring the physical parameters of the coldest T-type brown dwarfs using low resolution near infrared spectra. By comparing H_2O- and H_2-sensitive spectral ratios between empirical data and theoretical atmosphere models, and calibrating these ratios to measurements for the well-characterized 2-5 Gyr companion brown dwarf Gliese 570D, we derive estimates of the effective temperatures and surface gravities for 13 mid- and late-type field T dwarfs. We also deduce the first quantitative estimate of subsolar metallicity for the peculiar T dwarf 2MASS 0937+2931. Derived temperatures are consistent with prior estimates based on parallax and bolometric luminosity measurements, and examination of possible systematic effects indicate that the results are robust. Two recently discovered late-type T dwarfs, 2MASS 0939-2448 and 2MASS 1114-2618, both appear to be >50 K cooler than the latest-type T dwarf, 2MASS 0415-0935, and are potentially the coldest and least luminous brown dwarfs currently known. We find that, in general, higher surface gravity T dwarfs have lower effective temperatures and luminosities for a given spectral type, explaining previously observed scatter in the Teff/spectral type relation for these objects. Masses, radii and ages are estimated for the T dwarfs in our sample using the evolutionary models of Burrows et al.; we also determine masses and radii independently for eight T dwarfs with measured luminosities. These two determinations are largely consistent, lending support to the validity of evolutionary models at late ages. Our method is well suited to large samples of faint brown dwarfs, and can ultimately be used to directly measure the substellar mass function and formation history in the Galaxy.

We present the discovery of SDSS J042348.57-041403.5 as a closely-separated (0"16) brown dwarf binary, resolved by the Hubble Space Telescope Near Infrared Camera and Multi-Object Spectrometer. Physical association is deduced from the angular proximity of the components and constraints on their common proper motion. SDSS0423-0414AB appears to be composed of two brown dwarfs with spectral types L6+/-1 and T2+/-1. Hence, this system straddles the transition between L dwarfs and T dwarfs, a unique evolutionary phase of brown dwarfs characterized by substantial shifts in spectral morphology over an apparently narrow effective temperature range. Binarity explains a number of unusual properties of SDSS 0423-0414, including its overluminosity and high effective temperature compared to other early-type T dwarfs, and possibly its conflicting spectral classifications (L7.5 in the optical, T0 in the near infrared). The relatively short estimated orbital period of this system (~15-20 yr) and the presence of Li I absorption in its combined light spectrum make it an ideal target for both resolved spectroscopy and dynamical mass measurements. SDSS 0423-0414AB joins a growing list of late-L/early-T dwarf binaries, the high percentage of which (~50%) may provide a natural explanation for observed peculiarities across the L/T transition.

A three-dimensional, Newtonian hydrodynamic technique is used to follow the postbounce phase of a stellar core collapse event. For realistic initial data we have employed post core-bounce snapshots of the iron core of a 20 solar mass star. The models exhibit strong differential rotation but have centrally condensed density stratifications. We demonstrate for the first time that such postbounce cores are subject to a so-called low-T/|W| nonaxisymmetric instability and, in particular, can become dynamically unstable to an m=1 - dominated spiral mode at T/|W| ~ 0.08. We calculate the gravitational wave emission by the instability and find that the emitted waves may be detectable by current and future GW observatories from anywhere in the Milky Way.

In the standard approaches to neutrino transport in the simulation of core-collapse supernovae, one will often start from the classical Boltzmann equation for the neutrino's spatial, temporal, and spectral evolution. For each neutrino species, and its anti-particle, the classical density in phase space, or the associated specific intensity, will be calculated as a function of time. The neutrino radiation is coupled to matter by source and sink terms on the "right-hand-side" of the transport equation and together with the equations of hydrodynamics this set of coupled partial differential equations for classical densities describes, in principle, the evolution of core collapse and explosion. However, with the possibility of neutrino oscillations between species, a purely quantum-physical effect, how to generalize this set of Boltzmann equations for classical quantities to reflect oscillation physics has not been clear. To date, the formalisms developed have retained the character of quantum operator physics involving complex quantities and have not been suitable for easy incorporation into standard supernova codes. In this paper, we derive generalized Boltzmann equations for quasi-classical, real-valued phase-space densities that retain all the standard oscillation phenomenology, including resonant flavor conversion (the MSW effect), neutrino self-interactions, and the interplay between decohering matter coupling and flavor oscillations. With this formalism, any code(s) that can now handle the solution of the classical Boltzmann or transport equation can easily be generalized to include neutrino oscillations in a quantum-physically consistent fashion.

We calculate the planet-star flux-density ratios as a function of wavelength from 0.5 microns to 25 microns for the transiting extrasolar giant planets TrES-1 and HD209458b and compare them with the recent Spitzer/IRAC-MIPS secondary eclipse data in the 4.5, 8.0, and 24 micron bands. With only three data points and generic calibration issues, detailed conclusions are difficult, but inferences regarding atmospheric composition, temperature, and global circulation can be made. Our models reproduce the observations reasonably well, but not perfectly, and we speculate on the theoretical consequences of variations around our baseline models. One preliminary conclusion is that we may be seeing in the data indications that the day side of a close-in extrasolar giant planet is brighter in the mid-infrared than its night side, unlike Jupiter and Saturn. This correspondence will be further tested when the data anticipated in other Spitzer bands are acquired, and we make predictions for what those data may show.

Using the 2D multi-group, flux-limited diffusion version of the code VULCAN/2D, that also incorporates rotation, we have calculated the collapse, bounce, shock formation, and early post-bounce evolutionary phases of a core-collapse supernova for a variety of initial rotation rates. This is the first series of such multi-group calculations undertaken in supernova theory with fully multi-D tools. We find that though rotation generates pole-to-equator angular anisotropies in the neutrino radiation fields, the magnitude of the asymmetries is not as large as previously estimated. Moreover, we find that the radiation field is always more spherically symmetric than the matter distribution, with its plumes and convective eddies. We present the dependence of the angular anisotropy of the neutrino fields on neutrino species, neutrino energy, and initial rotation rate. Only for our most rapidly rotating model do we start to see qualitatively different hydrodynamics, but for the lower rates consistent with the pre-collapse rotational profiles derived in the literature the anisotropies, though interesting, are modest. This does not mean that rotation does not play a key role in supernova dynamics. The decrease in the effective gravity due to the centripetal effect can be quite important. Rather, it means that when a realistic mapping between initial and final rotational profiles and 2D multi-group radiation-hydrodynamics are incorporated into collapse simulations the anisotropy of the radiation fields may be only a secondary, not a pivotal factor, in the supernova mechanism.

Using previous measurements and quantum chemical calculations to derive the molecular properties of the TiH molecule, we obtain new values for its ro-vibrational constants, thermochemical data, spectral line lists, line strengths, and absorption opacities. Furthermore, we calculate the abundance of TiH in M and L dwarf atmospheres and conclude that it is much higher than previously thought. We find that the TiH/TiO ratio increases strongly with decreasing metallicity, and at high temperatures can exceed unity. We suggest that, particularly for subdwarf L and M dwarfs, spectral features of TiH near $\sim$0.52 \mic, 0.94 \mic, and in the $H$ band may be more easily measureable than heretofore thought. The recent possible identification in the L subdwarf 2MASS J0532 of the 0.94 \mic feature of TiH is in keeping with this expectation. We speculate that looking for TiH in other dwarfs and subdwarfs will shed light on the distinctive titanium chemistry of the atmospheres of substellar-mass objects and the dimmest stars.

Astronomy is at times a science of unexpected discovery. When it is, and if we are lucky, new intellectual territories emerge to challenge our views of the cosmos. The recent indirect detections using high-precision Doppler spectroscopy of now more than one hundred giant planets orbiting more than one hundred nearby stars is an example of such rare serendipity. What has been learned has shaken our preconceptions, for none of the planetary systems discovered to date is like our own. However, the key to unlocking a planet's chemical, structural, and evolutionary secrets is the direct detection of the planet's light. I review the embryonic theory of the spectra, atmospheres, and light curves of irradiated giant planets and put this theory into the context of the many proposed astronomical campaigns to image them.

We calculate self-consistent extrasolar giant planet (EGP) phase functions and light curves for orbital distances ranging from 0.2 AU to 15 AU. We explore the dependence on wavelength, cloud condensation, and Keplerian orbital elements. We find that the light curves of EGPs depend strongly on wavelength, the presence of clouds, and cloud particle sizes. Furthermore, the optical and infrared colors of most EGPs are phase-dependent, tending to be reddest at crescent phases in $V-R$ and $R-I$. Assuming circular orbits, we find that at optical wavelengths most EGPs are 3 to 4 times brighter near full phase than near greatest elongation for highly-inclined (i.e., close to edge-on) orbits. Furthermore, we show that the planet/star flux ratios depend strongly on the Keplerian elements of the orbit, particularly inclination and eccentricity. Given a sufficiently eccentric orbit, an EGP's atmosphere may make periodic transitions from cloudy to cloud-free, an effect that may be reflected in the shape and magnitude of the planet's light curve. Such elliptical orbits also introduce an offset between the time of the planet's light curve maximum and the time of full planetary phase, and for some sets of orbital parameters, this light curve maximum can be a steeply increasing function of eccentricity. We investigate the detectability of EGPs by proposed space-based direct-imaging instruments.

We construct models of core-collapse supernovae in one spatial dimension, including rotation, angular momentum transport, and viscous dissipation employing an alpha-prescription. We compare the evolution of a fiducial 11 M_sun non-rotating progenitor with its evolution including a wide range of imposed initial rotation profiles (1.25<P_0<8 s, where P_0 is the initial, approximately solid-body, rotation period of the iron core). This range of P_0 covers the region of parameter space from where rotation begins to modify the dynamics (P_0~8 s) to where angular velocities at collapse approach Keplerian (P_0~1 s). Assuming strict angular momentum conservation, all models in this range leave behind neutron stars with spin periods <10 ms, shorter than those of most radio pulsars, but similar to those expected theoretically for magnetars at birth. A fraction of the gravitational binding energy of collapse is stored in the free energy of differential rotation. This energy source may be tapped by viscous processes, providing a mechanism for energy deposition that is not strongly coupled to the mass accretion rate through the stalled supernova shock. This effect yields qualitatively new dynamics in models of supernovae. We explore several potential mechanisms for viscosity in the core-collapse environment: neutrino viscosity, turbulent viscosity caused by the magnetorotational instability (MRI), and turbulent viscosity by entropy- and composition-gradient-driven convection. We argue that the MRI is the most effective. We find that for rotation periods in the range P_0<~5 s, and a range of viscous stresses, that the post-bounce dynamics is significantly effected by the inclusion of this extra energy deposition mechanism; in several cases we obtain strong supernova explosions.

We review some of the reasons for believing that the generic core-collapse supernova is neutrino-driven, not MHD-jet driven. We include a discussion of the possible role of rotation in supernova blast energetics and morphology, and speculate on the origin of Cas A's and SN1987A's ejecta fields. Two "explosive" phenomena may be associated with most core collapses, the neutrino-driven supernova itself and an underenergetic jet-like ejection that follows. The latter may be a magnetic wind that easily penetrates the debris created by the much more energetic supernova. We predict that many core-collapse supernova remnants should have sub-dominant jet-like features. In Cas A, we associate this sub-dominant collimated wind with the "jet/counter-jet" structure observed. We suggest that the actual Cas A explosion itself is at nearly right angles to this jet, along a rotation axis that coincides with the bulk of the ejecta, the iron lobes, and the putative direction of motion of the point source. It may be that when rotation becomes sufficiently rapid that the strong-neutrino-driven-supernova/weak-jet duality switches to a strong-MHD-jet scenario that might be associated with hypernovae, and in some cases GRBs. Finally, we present a calculation using a new 2D multi-group, flux-limited radiation/hydrodynamics code we have recently developed for the simulation of core-collapse supernovae. We discuss the rotation-induced anisotropy in the neutrino radiation field, neutrino heating, and the neutrino flux vectors and speculate on rotation's possible role in the supernova mechanism and in the overall supernova phenomenon.

The mechanism responsible for the natal kicks of neutron stars continues to be a challenging problem. Indeed, many mechanisms have been suggested, and one hydrodynamic mechanism may require large initial asymmetries in the cores of supernova progenitor stars. Goldreich, Lai, & Sahrling (1997) suggested that unstable g-modes trapped in the iron (Fe) core by the convective burning layers and excited by the $\epsilon$-mechanism may provide the requisite asymmetries. We perform a modal analysis of the last minutes before collapse of published core structures and derive eigenfrequencies and eigenfunctions, including the nonadiabatic effects of growth by nuclear burning and decay by both neutrino and acoustic losses. In general, we find two types of g-modes: inner-core g-modes, which are stabilized by neutrino losses and outer-core g-modes which are trapped near the burning shells and can be unstable. Without exception, we find at least one unstable g-mode for each progenitor in the entire mass range we consider, 11 M$_{\sun}$ to 40 M$_{\sun}$. More importantly, we find that the timescales for growth and decay are an order of magnitude or more longer than the time until the commencement of core collapse. We conclude that the $\epsilon$-mechanism may not have enough time to significantly amplify core g-modes prior to collapse.

I summarize, in the form of an extended abstract, the ongoing efforts at the University of Arizona (and in collaboration) to understand core-collapse supernovae theoretically. Included are short discussions of 1D (SESAME) and 2D (VULCAN/2D) codes and results, as well as discussions of the possible role of rotation. Highlighted are recent developments in multi-dimensional radiation hydrodynamics and the essential physics of the neutrino-driven mechanism.

We calculate radius versus age trajectories for the photometrically-selected transiting extrasolar giant planet, OGLE-TR-56b, and find agreement between theory and observation, without introducing an ad hoc extra source of heat in its core. The fact that the radius of HD209458b seems larger than the radii of the recently discovered OGLE family of extremely close-in transiting planets suggests that HD209458b is anomalous. Nevertheless, our good fit to OGLE-TR-56b bolsters the notion that the generic dependence of transit radii on stellar irradiation, mass, and age is, to within error bars, now quantitatively understood.

Neutrino-matter cross sections and interaction rates are central to the core-collapse supernova phenomenon and, very likely, to the viability of the explosion mechanism itself. In this paper, we describe the major neutrino scattering, absorption, and production processes that together influence the outcome of core collapse and the cooling of protoneutron stars. One focus is on energy redistribution and many-body physics, but our major goal is to provide a useful resource for those interested in supernova neutrino microphysics.

We calculate as a function of orbital distance, mass, and age the theoretical spectra and orbit-averaged planet/star flux ratios for representative wide-separation extrasolar giant planets (EGPs) in the optical, near-infrared, and mid-infrared. Stellar irradiation of the planet's atmosphere and the effects of water and ammonia clouds are incorporated and handled in a consistent fashion. We include predictions for 12 specific known EGPs. In the process, we derive physical diagnostics that can inform the direct EGP detection and remote sensing programs now being planned or proposed. Furthermore, we calculate the effects of irradiation on the spectra of a representative companion brown dwarf as a function of orbital distance.

We have developed a time-dependent, multi-energy-group, and multi-angle (S$_n$) Boltzmann transport scheme for radiation hydrodynamics simulations, in one and two spatial dimensions. The implicit transport is coupled to both 1D (spherically-symmetric) and 2D (axially-symmetric) versions of the explicit Newtonian hydrodynamics code VULCAN. The 2D variant, VULCAN/2D, can be operated in general structured or unstructured grids and though the code can address many problems in astrophysics it was constructed specifically to study the core-collapse supernova problem. Furthermore, VULCAN/2D can simulate the radiation/hydrodynamic evolution of differentially rotating bodies. We summarize the equations solved and methods incorporated into the algorithm and present results of a time-dependent 2D test calculation. A more complete description of the algorithm is postponed to another paper. We highlight a 2D test run that follows for 22 milliseconds the immediate post-bounce evolution of a collapsed core. We present the relationship between the anisotropies of the overturning matter field and the distribution of the corresponding flux vectors, as a function of energy group. This is the first 2D multi-group, multi-angle, time-dependent radiation/hydro calculation ever performed in core collapse studies. Though the transport module of the code is not gray and does not use flux limiters (however, there is a flux-limited variant of VULCAN/2D), it still does not include energy redistribution and most velocity-dependent terms.

There are many interesting topics at the intersection of physics and astrophysics we call Supernova Theory. A small subset of them include the origin of pulsar kicks, gravitational radiation signatures of core bounce, and the possible roles of neutrinos and rotation in the mechanism of explosion. In this brief communication we summarize various recent ideas and calculations that bear on these themes.

We have carried out an extensive set of two-dimensional, axisymmetric, purely-hydrodynamic calculations of rotational stellar core collapse with a realistic, finite-temperature nuclear equation of state and realistic massive star progenitor models. For each of the total number of 72 different simulations we performed, the gravitational wave signature was extracted via the quadrupole formula in the slow-motion, weak-field approximation. We investigate the consequences of variation in the initial ratio of rotational kinetic energy to gravitational potential energy and in the initial degree of differential rotation. Furthermore, we include in our model suite progenitors from recent evolutionary calculations that take into account the effects of rotation and magnetic torques. For each model, we calculate gravitational radiation wave forms, characteristic wave strain spectra, energy spectra, final rotational profiles, and total radiated energy. In addition, we compare our model signals with the anticipated sensitivities of the 1st- and 2nd-generation LIGO detectors coming on line. We find that most of our models are detectable by LIGO from anywhere in the Milky Way.

We explore the spectral and atmospheric properties of brown dwarfs cooler than the latest known T dwarfs. Our focus is on the yet-to-be-discovered free-floating brown dwarfs in the \teff range from $\sim$800 K to $\sim$130 K and with masses from 25 to 1 \mj. This study is in anticipation of the new characterization capabilities enabled by the launch of SIRTF and the eventual launch of JWST. We provide spectra from $\sim$0.4 \mic to 30 \mic, highlight the evolution and mass dependence of the dominant H$_2$O, CH$_4$, and NH$_3$ molecular bands, consider the formation and effects of water-ice clouds, and compare our theoretical flux densities with the sensitivities of the instruments on board SIRTF and JWST. The latter can be used to determine the detection ranges from space of cool brown dwarfs. In the process, we determine the reversal point of the blueward trend in the near-infrared colors with decreasing \teff, the \teffs at which water and ammonia clouds appear, the strengths of gas-phase ammonia and methane bands, the masses and ages of the objects for which the neutral alkali metal lines are muted, and the increasing role as \teff decreases of the mid-infrared fluxes longward of 4 \mic. These changes suggest physical reasons to expect the emergence of at least one new stellar class beyond the T dwarfs. Our spectral models populate, with cooler brown dwarfs having progressively more planet-like features, the theoretical gap between the known T dwarfs and the known giant planets. Such objects likely inhabit the galaxy, but their numbers are as yet unknown.

We show that under certain circumstances the differences between the absorption mean and Planck mean opacities can lead to multiple solutions for an LTE atmospheric structure. Since the absorption and Planck mean opacities are not expected to differ significantly in the usual case of radiative equilibrium, non-irradiated atmospheres, the most interesting situations where the effect may play a role are strongly irradiated stars and planets, and also possibly structures where there is a significant deposition of mechanical energy, such as stellar chromospheres and accretion disks. We have presented an illustrative example of a strongly irradiated giant planet where the bifurcation effect is predicted to occur for a certain range of distances from the star.

Using a full frequency-dependent atmosphere code that can incorporate irradiation by a central primary star, we calculate self-consistent boundary conditions for the evolution of the radius of the transiting planet HD 209458b. Using a well-tested extrasolar giant planet evolutionary code, we then calculate the behavior of this planet's radius with age. The measured radius is in fact a transit radius that resides high in HD 209458b's inflated atmosphere. Using our derived atmospheric and interior structures, we find that irradiation plus the proper interpretation of the transit radius can yield a theoretical radius that is within the measured error bars. We conclude that if HD 209458b's true transit radius is at the lower end of the measured range, an extra source of core heating power is not necessary to explain the transit observations.

We calculate new line lists and opacities for the $F^4\Delta_i-X^4\Delta_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\"onl-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.

We present 6300--10100 {\AA} 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$\alpha$ 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 {\AA} 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.

We present the discovery of the first L-type subdwarf, 2MASS J05325346+8246465. This object exhibits enhanced collision-induced H$_2$ absorption, resulting in blue NIR colors ($J-K_s = 0.26{\pm}0.16$). In addition, strong hydride bands in the red optical and NIR, weak TiO absorption, and an optical/J-band spectral morphology similar to the L7 DENIS 0205$-$1159AB imply a cool, metal-deficient atmosphere. We find that 2MASS 0532+8246 has both a high proper motion, $\mu$ = 2$\farcs60\pm0\farcs$15 yr$^{-1}$, and a substantial radial velocity, $v_{rad} = -195{\pm}11$ km s$^{-1}$, and its probable proximity to the Sun (d = 10--30 pc) is consistent with halo membership. Comparison to subsolar-metallicity evolutionary models strongly suggests that 2MASS 0532+8246 is substellar, with a mass of 0.077 $\lesssim$ M $\lesssim$ 0.085 M$_{\sun}$ for ages 10--15 Gyr and metallicities $Z = 0.1-0.01$ $Z_{\sun}$. The discovery of this object clearly indicates that star formation occurred below the Hydrogen burning mass limit at early times, consistent with prior results indicating a flat or slightly rising mass function for the lowest-mass stellar subdwarfs. Furthermore, 2MASS 0532+8246 serves as a prototype for a new spectral class of subdwarfs, additional examples of which could be found in NIR proper motion surveys.

(Abridged) We present results from dynamical models of core-collapse supernovae in one spatial dimension, employing a newly-developed Boltzmann neutrino radiation transport algorithm, coupled to Lagrangean hydrodynamics and a consistent high-density nuclear equation of state. We focus on shock breakout and its neutrino signature and follow the dynamical evolution of the cores of 11 M_sun, 15 M_sun, and 20 M_sun progenitors through collapse and the first 250 milliseconds after bounce. We examine the effects on the emergent neutrino spectra, light curves, and mix of species of artificial opacity changes, the number of energy groups, the weak magnetism/recoil corrections, nucleon-nucleon bremsstrahlung, neutrino-electron scattering, and the compressibility of nuclear matter. Furthermore, we present the first high-resolution look at the angular distribution of the neutrino radiation field both in the semi-transparent regime and at large radii and explore the accuracy with which our tangent-ray method tracks the free propagation of a pulse of radiation in a near vacuum. Finally, we fold the emergent neutrino spectra with the efficiencies and detection processes for a selection of modern underground neutrino observatories and argue that the prompt electron-neutrino breakout burst from the next galactic supernova is in principle observable and usefully diagnostic of fundamental collapse/supernova behavior. Though we are not in this study focusing on the supernova mechanism per se, our simulations support the theoretical conclusion (already reached by others) that spherical (1D) supernovae do not explode when good physics and transport methods are employed.

We present a comprehensive theory of the spectra and atmospheres of irradiated extrasolar giant planets. We explore the dependences on stellar type, orbital distance, cloud characteristics, planet mass, and surface gravity. Phase-averaged spectra for specific known extrasolar giant planets that span a wide range of the relevant parameters are calculated, plotted, and discussed. The connection between atmospheric composition and emergent spectrum is explored in detail. Furthermore, we calculate the effect of stellar insolation on brown dwarfs. We review a variety of representative observational techniques and programs for their potential for direct detection, in light of our theoretical expectations, and we calculate planet-to-star flux ratios as a function of wavelength. Our results suggest which spectral features are most diagnostic of giant planet atmospheres and reveal the best bands in which to image planets of whatever physical or orbital characteristics.

Using a self-consistent atmosphere code, we construct a new model of the atmosphere of the transiting extrasolar giant planet HD 209458b to investigate the disparity between the observed strength of the sodium absorption feature at 589 nm and the predictions of previous models. For the atmospheric temperature-pressure profile we derive, silicate and iron clouds reside at a pressure of several mbar in the planet's atmosphere. These clouds have significant vertical extent and optical depth due to our slant viewing geometry and lead to increased absorption in bands directly adjacent to the sodium line core. Using a non-LTE sodium ionization model that includes photoionization by stellar UV flux, collisional processes with H_2, and radiative recombination, we show that the ionization depth in the planet's atmosphere reaches ~1/2 mbar at the day/night terminator. Ionization leads to a slight weakening of the sodium feature. We present our baseline model, including ionization and clouds, which falls near the observational error bars. The sensitivity of our conclusions to the derived atmospheric temperature-pressure profile is discussed.

Because the opacity of clouds in substellar mass object (SMO) atmospheres depends on the composition and distribution of particle sizes within the cloud, a credible cloud model is essential for accurately modeling SMO spectra and colors. We present a one--dimensional model of cloud particle formation and subsequent growth based on a consideration of basic cloud microphysics. We apply this microphysical cloud model to a set of synthetic brown dwarf atmospheres spanning a broad range of surface gravities and effective temperatures (g_surf = 1.78 * 10^3 -- 3 * 10^5 cm/s^2 and T_eff = 600 -- 1600 K) to obtain plausible particle sizes for several abundant species (Fe, Mg2SiO4, and Ca2Al2SiO7). At the base of the clouds, where the particles are largest, the particle sizes thus computed range from ~5 microns to over 300 microns in radius over the full range of atmospheric conditions considered. We show that average particle sizes decrease significantly with increasing brown dwarf surface gravity. We also find that brown dwarfs with higher effective temperatures have characteristically larger cloud particles than those with lower effective temperatures. We therefore conclude that it is unrealistic when modeling SMO spectra to apply a single particle size distribution to the entire class of objects.

Neutrino-matter interaction rates are central to the core collapse phenomenon and, perhaps, to the viability of the mechanism of core-collapse supernova explosions. In this paper we catalog and discuss the major neutrino scattering, absorption, and production processes that together influence the outcome of core collapse and the cooling of protoneutron stars. These are the essential inputs into the codes used to simulate the supernova phenomenon and an understanding of these processes is a prerequisite to continuing progress in supernova theory.

We review the status of the current quest to understand the mechanism of core-collapse supernovae, if neutrino-driven. In the process, we discuss the spherical explosion paradigm and its problems, some results from our new suite of collapse calculations performed using a recently-developed 1D implicit, multi-group, Feautrier/tangent-ray, Boltzmann solver coupled to explicit predictor/corrector hydrodynamics, the basic energetics of supernova explosions, and the promise of multi-D radiation/hydro simulations to explain why the cores of massive stars explode.

At the low temperatures achieved in cool brown dwarf and hot giant planet atmospheres, the less refractory neutral alkali metals assume an uncharacteristically prominent role in spectrum formation. In particular, the wings of the Na-D (5890 \AA) and K I (7700 \AA) resonance lines come to define the continuum and dominate the spectrum of T dwarfs from 0.4 to 1.0 \mic. Whereas in standard stellar atmospheres the strengths and shapes of the wings of atomic spectral lines are rarely needed beyond 25 \AA of a line center, in brown dwarfs the far wings of the Na and K resonance lines out to 1000's of \AA detunings are important. Using standard quantum chemical codes and the Unified Franck-Condon model for line profiles in the quasi-static limit, we calculate the interaction potentials and the wing line shapes for the dominant Na and K resonance lines in H$_2$- and helium-rich atmospheres. Our theory has natural absorption profile cutoffs, has no free parameters, and is readily adapted to spectral synthesis calculations for stars, brown dwarfs, and planets with effective temperatures below 2000 Kelvin.

We review some of the characteristics of irradiated extrasolar giant planets (EGPs), in anticipation of their direct detection from the ground and from space. Spectral measurements are the key to unlocking their structural and atmospheric characteristics and to determining the true differences between giant planets and brown dwarfs. In this spirit, the theoretical spectral and atmospheric calculations we summarize here are in support of the many searches for EGPs to be conducted in the coming decade by astronomers from around the world.

In this brief contribution, I touch on a subset of recent efforts in spectral and opacity modeling aimed at improving our understanding of brown dwarfs, L dwarfs, and T dwarfs. I discuss theoretical calculations of the alkali line profiles, newly generated CrH opacities, new evidence for refractory rainout in T dwarfs from optical spectral measurements, and the distinction between brown dwarfs and planets.

In this paper, we calculate new line lists and opacities for the 12 bands of the A$^6\Sigma^{+}$ -- X$^6\Sigma^{+}$ transitions of the CrH molecule. Identified in objects of the new L dwarf spectroscopic class (many of which are brown dwarfs), as well as in sunspots, the CrH molecule plays an important role in the diagnosis of low-temperature atmospheres. As a tentative first application of these opacities, we employ our new theoretical CrH data in an atmospheres code to obtain a CrH/H$_2$ number ratio for the skin of the L5 dwarf 2MASSI J1507038-151648 of $\sim 2-4\times 10^{-9}$, in rough agreement with chemical equilibrium expectations. Since in previous compilations the oscillator strength was off by more than an order of magnitude, this agreement represents a modest advance. However, in order to determine the CrH abundance in an L dwarf atmosphere, silicate clouds need to be incorporated into the model, and cloud modeling is still in a primitive stage of development. Nevertheless, one important step in L dwarf modeling is a reliable CrH opacity and this is what we have here attempted to provide.

We have generated new, self-consistent spectral and atmosphere models for the effective temperature range 600 K to 1300 K thought to encompass the known T dwarfs. For the first time, theoretical models are compared with a {\it family} of measured T dwarf spectra at wavelengths shortward of $\sim$1.0 micron. By defining spectral indices and standard colors in the optical and very near-infrared, we explore the theoretical systematics with \teff, gravity, and metallicity. We conclude that the short- wavelength range is rich in diagnostics that complement those in the near-infrared now used for spectral subtyping. We also conclude that the wings of the Na D and K I (7700\AA) resonance lines and aggressive rainout of heavy metals (with the resulting enhancement of the sodium and potassium abundances at altitude) are required to fit the new data shortward of 1.0 \mic. Furthermore, we find that the water bands weaken with increasing gravity, that modest decreases in metallicity enhance the effect in the optical of the sodium and potassium lines, and that at low \teffs, in a reversal of the normal pattern, optical spectra become bluer with further decreases in \teff. Moreover, we conclude that T dwarf subtype is not a function of \teff alone, but that it is a non-trivial function of gravity and metallicity as well. As do Marley et al. (2001), we see evidence in early T dwarf atmospheres of a residual effect of clouds. With cloudless models, we obtain spectral fits to the two late T dwarfs with known parallaxes, but a residual effect of clouds on the emergent spectra of even late T dwarfs can not yet be discounted.

Use of the highly sensitive Hokupa'a curvature wavefront sensor has allowed for the first time direct adaptive optics (AO) guiding on brown dwarfs and VLM stars (SpT=M7-L2). An initial survey of 9 such objects discovered one 0.15" binary (2MASSJ 1426316+155701). The companion is about half as bright as the primary (Delta K = 0.61+/-0.05$, Delta H = 0.70+/-0.05) and has even redder colors H-K=0.59+/-0.14 than the primary. The blended spectrum of the binary has been previously determined to be M9.0. We modeled a blend of an M8.5 template and a L1-L3 template reproducing a M9.0 spectrum in the case of Delta K = 0.61+/-0.05,Delta H = 0.70\pm0.05$. These spectral types also match the observed H-K colors of each star. Based the previously observed low space motion and $H_{\alpha}$ activity we assign an age of $0.8^{+6.7}_{-0.3} Gyr$. Utilizing this age range and the latest DUSTY models of the Lyon group we assign a photometric distance of $18.8^{+1.44}_{-1.02} pc$ and masses of $M_{A}=0.074^{+0.005}_{-0.011} M_\odot$ and $M_{B}=0.066^{+0.006}_{-0.015} M_\odot$. We therefore estimate a system separation of $2.92_{+0.22}^{-0.16}AU$ and a period of $13.3{+3.18}^{-1.51} yr$ respectively. Hence, 2M1426 is among the smallest separation brown dwarf binaries resolved to date.

We present near-infrared spectra for a sample of T dwarfs, including eleven new discoveries made using the Two Micron All Sky Survey. These objects are distinguished from warmer (L-type) brown dwarfs by the presence of methane absorption bands in the 1--2.5 $\micron$ spectral region. A first attempt at a near-infrared classification scheme for T dwarfs is made, based on the strengths of CH$_4$ and H$_2$O bands and the shapes of the 1.25, 1.6, and 2.1 $\micron$ flux peaks. Subtypes T1 V through T8 V are defined, and spectral indices useful for classification are presented. The subclasses appear to follow a decreasing T$_{eff}$ scale, based on the evolution of CH$_4$ and H$_2$O bands and the properties of L and T dwarfs with known distances. However, we speculate that this scale is not linear with spectral type for cool dwarfs, due to the settling of dust layers below the photosphere and subsequent rapid evolution of spectral morphology around T$_{eff}$ $\sim$ 1300--1500 K. Similarities in near-infrared colors and continuity of spectral features suggest that the gap between the latest L dwarfs and earliest T dwarfs has been nearly bridged. This argument is strengthened by the possible role of CH$_4$ as a minor absorber shaping the K-band spectra of the latest L dwarfs. Finally, we discuss one peculiar T dwarf, 2MASS 0937+2931, which has very blue near-infrared colors (J-K$_s$ = $-0.89\pm$0.24) due to suppression of the 2.1 $\micron$ peak. The feature is likely caused by enhanced collision-induced H$_2$ absorption in a high pressure or low metallicity photosphere.

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 ($\tau_\rho\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.

We present a synthesis of physical effects influencing the observed lightcurve of an extrasolar giant planet (EGP) transiting its host star. The synthesis includes a treatment of Rayleigh scattering, cloud scattering, refraction, and molecular absorption of starlight in the EGP atmosphere. Of these effects, molecular absorption dominates in determining the transit-derived radius $R$. Using a generic model for the atmosphere of EGP HD209458b, we perform a fit to the best available transit lightcurve data, and infer that this planet has a radius at a pressure of 1 bar, $R_1$, equal to 94430 km. We predict that $R$ will be a function of wavelength of observation, with a robust prediction of variations of $\pm 1$% at infrared wavelengths where H$_2$O opacity in the high EGP atmosphere dominates.

Straddling the traditional realms of the planets and the stars, objects below the edge of the main sequence have such unique properties, and are being discovered in such quantities, that one can rightly claim that a new field at the interface of planetary science and and astronomy is being born. In this review, we explore the essential elements of the theory of brown dwarfs and giant planets, as well as of the new spectroscopic classes L and T. To this end, we describe their evolution, spectra, atmospheric compositions, chemistry, physics, and nuclear phases and explain the basic systematics of substellar-mass objects across three orders of magnitude in both mass and age and a factor of 30 in effective temperature. Moreover, we discuss the distinctive features of those extrasolar giant planets that are irradiated by a central primary, in particular their reflection spectra, albedos, and transits. Aspects of the latest theory of Jupiter and Saturn are also presented. Throughout, we highlight the effects of condensates, clouds, molecular abundances, and molecular/atomic opacities in brown dwarf and giant planet atmospheres and summarize the resulting spectral diagnostics. Where possible, the theory is put in its current observational context.

We have observed 4 transits of the planet of HD 209458 using the STIS spectrograph on HST. Summing the recorded counts over wavelength between 582 nm and 638 nm yields a photometric time series with 80 s time sampling and relative precision of about 1.1E-4 per sample. The folded light curve can be fit within observational errors using a model consisting of an opaque circular planet transiting a limb-darkened stellar disk. In this way we estimate the planetary radius R_p = 1.347 +/- 0.060 R_Jup, the orbital inclination i = 86.68 +/- 0.14 degrees, the stellar radius R_* = 1.146 +/- 0.050 R_solar, and one parameter describing the stellar limb darkening. Our estimated radius is smaller than those from earlier studies, but is consistent within measurement errors, and is also consistent with theoretical estimates of the radii of irradiated Jupiter-like planets. Satellites or rings orbiting the planet would, if large enough, be apparent from distortions of the light curve or from irregularities in the transit timings. We find no evidence for either satellites or rings, with upper limits on satellite radius and mass of 1.2 Earth radii and 3 Earth masses, respectively. Opaque rings, if present, must be smaller than 1.8 planetary radii in radial extent. The high level of photometric precision attained in this experiment confirms the feasibility of photometric detection of Earth-sized planets circling Sun-like stars.

I summarize some of the consequences for the optical and very-near-infrared spectra of T dwarfs (in particular) and brown dwarfs (in general) of their possible dominance by the neutral alkali metal lines. As a byproduct of this study, I estimate the true optical color of ``brown'' dwarfs.

We report results from a large Hubble Space Telescope project to observe a significant (~34,000) ensemble of main sequence stars in the globular cluster 47 Tucanae with a goal of defining the frequency of inner-orbit, gas-giant planets. Simulations based on the characteristics of the 8.3 days of time-series data in the F555W and F814W WFPC2 filters show that ~17 planets should be detected by photometric transit signals if the frequency of hot Jupiters found in the solar neighborhood is assumed to hold for 47 Tuc. The experiment provided high-quality data sufficient to detect planets. A full analysis of these WFPC2 data reveals ~75 variables, but NO light curves resulted for which a convincing interpretation as a planet could be made. The planet frequency in 47 Tuc is at least an order of magnitude below that for the solar neighborhood. The cause of the absence of close-in planets in 47 Tuc is not yet known; presumably the low metallicity and/or crowding of 47 Tuc interfered with planet formation, with orbital evolution to close-in positions, or with planet survival.

In preparation for a set of hydrodynamical simulations of core-collapse supernovae and protoneutron star winds, we investigate the rates of production and thermalization of $\nu_\mu$ and $\nu_\tau$ neutrinos in dense nuclear matter. Included are contributions from electron scattering, electron-positron annihilation, nucleon-nucleon bremsstrahlung, and nucleon scattering. We find that nucleon scattering dominates electron scattering as a thermalization process at neutrino energies greater than $\sim$15 MeV. In addition, nucleon-nucleon bremsstrahlung dominates electron-positron annihilation as a production mechanism at low neutrino energies, near and below the $\nu_\mu$ and $\nu_\tau$ neutrinospheres. Furthermore, we have begun a study of steady-state general relativistic protoneutron star winds employing simple neutrino heating and cooling terms. From this analysis we obtain acceleration profiles as well as asymptotic lepton fractions and baryon entropies essential in assessing the wind as a potential site for $r$-process nucleosynthesis.

A Keck~II low resolution spectrum shortward of ome-micron is presented for SDSS 1624+0029, the first field methane or T dwarf discovered in the Sloan Digital Sky Survey. Significant flux is detected down to the spectrum's short wavelength limit of 6200\AA. The spectrum exhibits a broad absorption feature centered at 7700\AA, which we interpret as the K~I 7665/7699 resonance doublet. The observed flux declines shortward of 7000\AA, due most likely to the red wing of the Na~I doublet. Both Cs~I doublet lines are detected more strongly than in an earlier red spectrum. Neither Li~I absorption nor H$\alpha$ emission are detected. An exploratory model fit to the spectrum suggests that the shape of the red spectrum can be primarily accounted for by the broad wings of the K~I and Na~I doublets. This behavior is consistent with the argument proffered by Burrows, Marley and Sharp that strong alkali absorption is principally responsible for depressing T dwarf spectra shortward of 1$\mu$m. In particular, there seems no compelling reason at this time to introduce dust or an additional opacity source in the atmosphere of the SDSS object. The width of the K~I and strengths of the Cs~I lines also indicate that the Sloan object is warmer than Gl~229B.

We have developed an implicit, multi-group, time-dependent, spherical neutrino transport code based on the Feautrier variables, the tangent-ray method, and accelerated ${\bf \Lambda}$ iteration. The code achieves high angular resolution, is good to O($v/c$), is equivalent to a Boltzmann solver (without gravitational redshifts), and solves the transport equation at all optical depths with precision. In this paper, we present our formulation of the relevant numerics and microphysics and explore protoneutron star atmospheres for snapshot post-bounce models. Our major focus is on spectra, neutrino-matter heating rates, Eddington factors, angular distributions, and phase-space occupancies. In addition, we investigate the influence on neutrino spectra and heating of final-state electron blocking, stimulated absorption, velocity terms in the transport equation, neutrino-nucleon scattering asymmetry, and weak magnetism and recoil effects. Furthermore, we compare the emergent spectra and heating rates obtained using full transport with those obtained using representative flux-limited transport formulations to gauge their accuracy and viability. Finally, we derive useful formulae for the neutrino source strength due to nucleon-nucleon bremsstrahlung and determine bremsstrahlung's influence on the emergent $\nu_{\mu}$ and $\nu_{\tau}$ neutrino spectra.

We generate theoretical albedo and reflection spectra for a full range of extrasolar giant planet (EGP) models, from Jovian to 51-Pegasi class objects. Our albedo modeling utilizes the latest atomic and molecular cross sections, a Mie theory treatment of extinction by condensates, a variety of particle size distributions, and an extension of the Feautrier radiative transfer method which allows for a general treatment of the scattering phase function. We find that due to qualitative similarities in the compositions and spectra of objects within each of five broad effective temperature ranges, it is natural to establish five representative EGP albedo classes: a ``Jovian'' class (T$_{\rm eff} \lesssim 150$ K; Class I) with tropospheric ammonia clouds, a ``water cloud'' class (T$_{\rm eff} \sim 250$ K; Class II) primarily affected by condensed H$_2$O, a ``clear'' class (T$_{\rm eff} \gtrsim 350$ K; Class III) which lacks clouds, and two high-temperature classes: Class IV (900 K $\lesssim$ T$_{\rm{eff}}$ $\lesssim$ 1500 K) for which alkali metal absorption predominates, and Class V (T$_{\rm{eff}}$ $\gtrsim$ 1500 K and/or low surface gravity ($\lesssim$ 10$^3$ cm s$^{-2}$)) for which a high silicate layer shields a significant fraction of the incident radiation from alkali metal and molecular absorption. The resonance lines of sodium and potassium are expected to be salient features in the reflection spectra of Class III, IV, and V objects. We derive Bond albedos and effective temperatures for the full set of known EGPs and explore the possible effects of non-equilibrium condensed products of photolysis above or within principal cloud decks. As in Jupiter, such species can lower the UV/blue albedo substantially, even if present in relatively small mixing ratios.

The recent discovery that the close-in extrasolar giant planet, HD209458b, transits its star has provided a first-of-its-kind measurement of the planet's radius and mass. In addition, there is a provocative detection of the light reflected off of the giant planet, $\tau$ Boo b. Including the effects of stellar irradiation, we estimate the general behavior of radius/age trajectories for such planets and interpret the large measured radii of HD209458b and $\tau$ Boo b in that context. We find that HD209458b must be a hydrogen-rich gas giant. Furthermore, the large radius of close-in gas giant is not due to the thermal expansion of its atmosphere, but to the high residual entropy that remains throughout its bulk by dint of its early proximity to a luminous primary. The large stellar flux does not inflate the planet, but retards its otherwise inexorable contraction from a more extended configuration at birth. This implies either that such a planet was formed near its current orbital distance or that it migrated in from larger distances ($\geq$0.5 A.U.), no later than a few times $10^7$ years of birth.

We identify the pressure--broadened red wings of the saturated potassium resonance lines at 7700 \AA as the source of anomalous absorption seen in the near-infrared spectra of Gliese 229B and, by extension, of methane dwarfs in general. This conclusion is supported by the recent work of Tsuji {\it et al.} 1999, though unlike them we find that dust need not be invoked to explain the spectra of methane dwarfs shortward of 1 micron. We find that a combination of enhanced alkali abundances due to rainout and a more realistic non-Lorentzian theory of resonant line shapes may be all that is needed to properly account for these spectra from 0.5 to 1.0 microns. The WFPC2 $I$ measurement of Gliese 229B is also consistent with this theory. Furthermore, a combination of the blue wings of this K I resonance doublet, the red wings of the Na D lines at 5890 \AA, and, perhaps, the Li I line at 6708 \AA can explain in a natural way the observed WFPC2 $R$ band flux of Gliese 229B. Hence, we conclude that the neutral alkali metals play a central role in the near-infrared and optical spectra of methane dwarfs and that their lines have the potential to provide crucial diagnostics of brown dwarfs. We speculate on the systematics of the near-infrared and optical spectra of methane dwarfs, for a given mass and composition, that stems from the progressive burial with decreasing \teff of the alkali metal atoms to larger pressures and depths.

We investigate the rates of production and thermalization of $\nu_\mu$ and $\nu_\tau$ neutrinos at temperatures and densities relevant to core-collapse supernovae and protoneutron stars. Included are contributions from electron scattering, electron-positron annihilation, nucleon-nucleon bremsstrahlung, and nucleon scattering. For the scattering processes, in order to incorporate the full scattering kinematics at arbitrary degeneracy, the structure function formalism developed by Reddy et al. (1998) and Burrows and Sawyer (1998) is employed. Furthermore, we derive formulae for the total and differential rates of nucleon-nucleon bremsstrahlung for arbitrary nucleon degeneracy in asymmetric matter. We find that electron scattering dominates nucleon scattering as a thermalization process at low neutrino energies ($\epsilon_\nu\lesssim 10$ MeV), but that nucleon scattering is always faster than or comparable to electron scattering above $\epsilon_\nu\simeq10$ MeV. In addition, for $\rho\gtrsim 10^{13}$ g cm$^{-3}$, $T\lesssim14$ MeV, and neutrino energies $\lesssim60$ MeV, nucleon-nucleon bremsstrahlung always dominates electron-positron annihilation as a production mechanism for $\nu_\mu$ and $\nu_\tau$ neutrinos.

One method of discriminating between the many Type Ia progenitor scenarios is by searching for contaminating hydrogen and helium stripped from the companion star. We present several high-resolution 2-D numerical simulations of the impact of a Type Ia supernova explosion with hydrogen-rich main sequence, subgiant, and red giant companions. The binary parameters were chosen to represent several classes of single-degenerate Type Ia progenitor models that have been suggested in the literature. We use realistic stellar models and supernova debris profiles to represent each binary system. For each scenario, we explore the hydrodynamics of the supernova-secondary interaction, calculate the amount of stellar material stripped from the secondary and the kick delivered by the impact, and construct the velocity and solid angle distributions of the stripped material. After the impact, the main sequence star is puffed up, much like a pre-main sequence star. The subgiant companion behaves similarly. The secondary will not be contaminated by much supernova debris from the initial impact, but it may accrete low-velocity iron-group elements (or oxygen and silicon if the ejecta is radially mixed) at late times. A He pre-white dwarf will be left behind after a supernova explosion with a red giant companion.

Analysis of initial observations from near-infrared sky surveys has shown that the resulting photometric catalogues, combined with far-red optical data, provide an extremely effective method of finding isolated, very low-temperature objects in the general field. Follow-up observations have already identified more than 25 sources with temperatures cooler than the latest M dwarfs. A comparison with detailed model predictions (Burrows & Sharp) indicates that these L dwarfs have effective temperatures between ~2000\pm100 K and 1500\pm100 K, while the available trigonometric parallax data place their luminosities at between 10^{-3.5} and 10^{-4.3} L_solar. Those properties, together with the detection of lithium in one-third of the objects, are consistent with the majority having substellar masses. The mass function cannot be derived directly, since only near-infrared photometry and spectral types are available for most sources, but we can incorporate VLM/brown dwarf models in simulations of the Solar Neighbourhood population and constrain Psi(M) by comparing the predicted L-dwarf surface densities and temperature distributions against observations from the DENIS and 2MASS surveys. The data, although sparse, can be represented by a power-law mass function, Psi(M) ~ M^{-alpha}, with 1 < alpha < 2. Current results favour a value nearer the lower limit. If alpha = 1.3, then the local space density of 0.075 > M/M_solar > 0.01 brown dwarfs is 0.10 systems pc^{-3}. In that case brown dwarfs are twice as common as main-sequence stars, but contribute no more than ~15% of the total mass of the disk.

We explore the effects on nucleosynthesis in Type II supernovae of various parameters (mass cut, neutron excess, explosion energy, progenitor mass) in order to explain the observed trends of the iron-peak element abundance ratios ([Cr/Fe], [Mn/Fe], [Co/Fe] and [Ni/Fe]) in halo stars as a function of metallicity for the range $ -4 \le$ [Fe/H] $\le -2.5$. [Cr/Fe] and [Mn/Fe] decrease with decreasing [Fe/H], while [Co/Fe] behaves the opposite way and increases. We show that such a behavior can be explained by a variation of mass cuts in Type II supernovae as a function of progenitor mass, which provides a changing mix of nucleosynthesis from an alpha-rich freeze-out of Si-burning and incomplete Si-burning. This explanation is consistent with the amount of ejected $^{56}$Ni determined from modeling the early light curves of individual supernovae. We also suggest that the ratio [H/Fe] of halo stars is mainly determined by the mass of interstellar hydrogen mixed with the ejecta of a single supernova which is larger for larger explosion energy and the larger Str\"omgren radius of the progenitor.

We present evolutionary sequences for Jupiter and Saturn, based on new nongray model atmospheres, which take into account the evolution of the solar luminosity and partitioning of dense components to deeper layers. The results are used to set limits on the extent to which possible interior phase separation of hydrogen and helium may have progressed in the two planets. When combined with static models constrained by the gravity field, our evolutionary calculations constrain the helium mass fraction in Jupiter to be between 0.20 and 0.27, relative to total hydrogen and helium. This is in agreement with the Galileo determination. The helium mass fraction in Saturn's atmosphere lies between 0.11 and 0.25, higher than the Voyager determination. Based on the discrepancy between the Galileo and Voyager results for Jupiter, and our models, we predict that Cassini measurements will yield a higher atmospheric helium mass fraction for Saturn relative to the Voyager value.

We summarize and extend recent work on the theory of extrasolar giant planets (EGPs) and brown dwarfs, paying particular attention to Gliese 229 B, the albedos of EGPs, the compositions of substellar atmospheres, the connections with the giant planets in the solar system, cloud physics, and non-gray spectral synthesis. The role of condensates in altering the optical spectrum of Gliese 229 B is explored, as are the systematics of the reflection spectra from extrasolar giant planets near their primaries. In addition, we discuss the role of convection and disequilibrium chemistry in explaining the anomalous detection of CO in Gliese 229 B. Throughout, we highlight the distinctive chemistry that defines this new class of objects and set goals for future study.

We calculate detailed chemical abundance profiles for a variety of brown dwarf and extrasolar giant planet atmosphere models, focusing in particular on Gliese 229B, and derive the systematics of the changes in the dominant reservoirs of the major elements with altitude and temperature. We assume an Anders and Grevesse (1989) solar composition of 27 chemical elements and track 330 gas--phase species, including the monatomic forms of the elements, as well as about 120 condensates. We address the issue of the formation and composition of clouds in the cool atmospheres of substellar objects and explore the rain out and depletion of refractories. We conclude that the opacity of clouds of low--temperature ($\le$900 K), small--radius condensibles (specific chlorides and sulfides), may be responsible for the steep spectrum of Gliese 229B observed in the near infrared below 1 \mic. Furthermore, we assemble a temperature sequence of chemical transitions in substellar atmospheres that may be used to anchor and define a sequence of spectral types for substellar objects with T$_{eff}$s from $\sim$2200 K to $\sim$100 K.

Recent observations of supernovae, supernova remnants, and radio pulsars suggest that there are correlations between pulsar kicks and spins, infrared and gamma-ray line profiles, supernova polarizations, and ejecta debris fields. A framework is emerging in which explosion asymmetries play a central role. The new perspective meshes recent multi-dimensional theoretical investigations of the explosion mechanism with trends in $^{56}$Ni yields and explosion kinetic energies. These trends imply that the mass above which black holes form after collapse is $\sim$30 M$_{\odot}$ and that supernova explosion energies may vary by as much as a factor of four. In addition, new neutrino-matter opacity calculations reveal that the inner cores of protoneutron stars are more transparent than hitherto suspected. This may have consequences for the delayed neutrino-driven mechanism of explosion itself.

Including nucleon--nucleon correlations due to both Fermi statistics and nuclear forces, we have developed a general formalism for calculating the charged--current neutrino--nucleon absorption rates in nuclear matter. We find that at one half nuclear density many--body effects alone suppress the rates by a factor of two and that the suppression factors increase to $\sim$5 at $4\times10^{14}$ g cm$^{-3}$. The associated increase in the neutrino--matter mean--free--paths parallels that found for neutral--current interactions and opens up interesting possibilities in the context of the delayed supernova mechanism and protoneutron star cooling.

Giant planets in circumstellar disks can migrate inward from their initial (formation) positions. Radial migration is caused by inward torques between the planet and the disk; by outward torques between the planet and the spinning star; and by outward torques due to Roche lobe overflow and consequent mass loss from the planet. We present self-consistent numerical considerations of the problem of migrating giant planets. Summing torques on planets for various physical parameters, we find that Jupiter-mass planets can stably arrive and survive at small heliocentric distances, thus reproducing observed properties of some of the recently discovered extra-solar planets. Inward migration timescales can be approximately equal to or less than disk lifetimes and star spindown timescales. Therefore, the range of fates of massive planets is broad, and generally comprises three classes: (I) planets which migrate inward too rapidly and lose all their mass; (II) planets which migrate inward, lose some but not all of their mass, and survive in very small orbits; and (III) planets which do not lose any mass. Some planets in Class III do not migrate very far from their formation locations. Our results show that there is a wide range of possible fates for Jupiter-mass planets for both final heliocentric distance and final mass.

Including nucleon-nucleon correlations due to both Fermi statistics and nuclear forces, we have developed a general formalism for calculating the neutral-current neutrino-nucleon opacities in nuclear matter. We derive corrections to the dynamic structure factors due to both density and spin correlations and find that neutrino-nucleon cross sections are suppressed by large factors around and above nuclear density. In addition, we find that the spectrum of energy transfers in neutrino scattering is considerably broadened by the interactions in the medium. An identifiable component of this broadening comes from the absorption and emission of quanta of collective modes akin to the Gamow-Teller and Giant Dipole resonances in nuclei (zero-sound; spin waves), with \v{C}erenkov kinematics. Under the assumption that both the charged-current and the neutral-current cross sections are decreased by many-body effects, we calculate a set of ad hoc protoneutron star cooling models to gauge the potential importance of the new opacities to the supernova itself. We find that after many hundreds of milliseconds to seconds the driving neutrino luminosities might be increased by from 10% to 100%. However, the actual consequences, if any, of these new neutrino opacities remain to be determined.

We use a Monte Carlo binary synthesis code to model the formation and evolution of neutron star systems including high-mass X-ray binaries, low-mass X-ray binaries, double neutron star systems and radio pulsars. Our focus is on the signature imprinted on such systems due to natal kicks to neutron stars over and above that imparted by orbital motions. The code incorporates the effect of the galactic potential (including rotation) on the velocities of these systems. A comparison between our models and the observations leads us to infer mean natal kicks between 400-500 km/s. Moreover, to be consistent with all the data, we require a bimodal kick distribution with one peak in the distribution near 0 km/s and the other above 600 km/s.

Since October of 1995, employing Doppler spectroscopy astronomers have discovered as many as 20 giant planets and brown dwarfs, including companions to $\tau$ Boo, 51 Peg, $\upsilon$ And, 55 Cnc, $\rho$ CrB, 70 Vir, 16 Cyg B, and 47 UMa. These discoveries have excited the planetary science community, astronomers, and the public at large. Within hours of the announcement of the planet 51 Peg b, the first direct detection of an unimpeachable brown dwarf, Gl229 B, was also announced. Gl229 B is a watershed since it has methane spectral features and a surface temperature below 1000 Kelvin, characteristics unique to objects with substellar masses. During the last two years, building upon our previous experience in the modeling of brown dwarfs and M stars, we published theoretical studies of the evolution and spectra of extrasolar giant planets. We have recently upgraded our capabilities and now generate non-gray spectral and color models of both giant planets and brown dwarfs. This theory will soon encompass objects whose effective temperatures range from 100 K to 4000 K and whose masses span three orders of magnitude. The evolutionary, spectral, and color calculations upon which we have embarked are in direct support of the searches now being planned in earnest with the HST (WFPC2, NICMOS), the IRTF, the MMT 6.5-meter upgrade, the LBT, Keck's I and II, ISO, UKIRT, NGST, the VLT, COROT, DENIS, 2MASS, and SIRTF.

Though it is agreed that the post-bounce, pre-explosion cores of massive stars are unstable to transonic, Rayleigh-Taylor-like instabilities, the role convective motions may play in igniting the supernova explosion itself is not yet clear. Whatever that role, the morphology and nucleosynthetic yields of the explosion, the recoil of the protoneutron star, the spatial distribution of fresh $^{56}$Ni, and the gravitational wave and neutrino signatures of the event will be affected in interesting ways. I review a subset of the issues raised by the new paradigm and some of the technical obstacles that remain. I focus on neutrino transfer, the role of progenitor structure, kick mechanisms, and the location and potential consequence of various hydrodynamic instabilities.

We have developed a new non-gray theory of the evolution, spectra, and colors of extrasolar giant planets (EGPs) and brown dwarfs that reveals their exotic nature and uniqueness. We have discovered that the fluxes of such objects for T$_{eff}$s from 1300 K to 100 K can be spectacularly higher in the near infrared bands than black body values and that their infrared colors are anomalously blue. As a consequence, EGPs and brown dwarfs reside in hitherto unoccupied realms of the H--R diagram and may be more easily found with current and planned telescopes than previously imagined.

We present the results of a new series of non-gray calculations of the atmospheres, spectra, colors, and evolution of extrasolar giant planets (EGPs) and brown dwarfs for effective temperatures below 1300 K. This theory encompasses most of the mass/age parameter space occupied by substellar objects and is the first spectral study down to 100 K. These calculations are in aid of the multitude of searches being conducted or planned around the world for giant planets and brown dwarfs and reveal the exotic nature of the class. Generically, absorption by H_2 at longer wavelengths and H_2O opacity windows at shorter wavelengths conspire to redistribute flux blueward. Below 1200 K, methane is the dominant carbon bearing molecule and is a universal diagnostic feature of EGP and brown dwarf spectra. We find that the primary bands in which to search are $Z$ (\sim 1.05 \mic), $J$ (\sim 1.2 \mic), $H$ (\sim 1.6 \mic), $K$ (\sim 2.2 \mic), $M$ (\sim 5 \mic), and $N$ (\sim 10 \mic), that enhancements of the emergent flux over blackbody values, in particular in the near infrared, can be by many orders of magnitude, and that the infrared colors of EGPs and brown dwarfs are much bluer than previously believed. In particular, relative to $J$ and $H$, the $K$ band flux is reduced by CH_4 and H_2 absorption. Furthermore, we derive that for T_{eff}s below 1200 K most or all true metals are sequestered below the photosphere, that an interior radiative zone is a generic feature of substellar objects, and that clouds of H_2O and NH_3 are formed for T_{eff}s below \sim 400 K and \sim 200 K, respectively. This study is done for solar-metallicity objects in isolation and does not include the effects of stellar insolation. Nevertheless, it is an attempt to bridge the gap between the planetary and stellar realms and to develop a non-gray theory of objects from 0.3 \mj (``saturn'') to 70 \mj ($\sim$0.07 \mo). We find that the detection ranges for brown dwarf/EGP discovery of both ground-- and space-based telescopes are larger than previously estimated.

With the discovery of the companions of 51 Peg, 55 Cnc, $\tau$ Boo, $\upsilon$ And, 70 Vir, 47 UMa, and Gl229, evolutionary and spectral models of gas giants and/or brown dwarfs with masses from 0.3 through 60 times that of Jupiter assume a new and central role in the emerging field of extrasolar planetary studies. In this contribution, we describe the structural, spectral, and evolutionary characteristics of such exotic objects, as determined by our recent theoretical calculations. These calculations can be used to establish direct search strategies via SIRTF, ISO, and HST (NICMOS), and via various ground-based adaptive optics and interferometric platforms planned for the near future.

Electron-degenerate, pressure-ionized hydrogen (usually referred to as metallic hydrogen) is the principal constituent of brown dwarfs, the long-sought objects which lie in the mass range between the lowest-mass stars (about eighty times the mass of Jupiter) and the giant planets. The thermodynamics and transport properties of metallic hydrogen are important for understanding the properties of these objects, which, unlike stars, continually and slowly cool from initial nondegenerate (gaseous) states. Within the last year, a brown dwarf (Gliese 229 B) has been detected and its spectrum observed and analyzed, and several examples of extrasolar giant planets have been discovered. The brown dwarf appears to have a mass of about forty to fifty Jupiter masses and is now too cool to be fusing hydrogen or deuterium, although we predict that it will have consumed all of its primordial deuterium. This paper reviews the current understanding of the interrelationship between its interior properties and its observed spectrum, and also discusses the current status of research on the structure of giant planets, both in our solar system and elsewhere.

There has been a new infusion of ideas in the study of the mechanism and early character of core--collapse supernovae. However, despite recent conceptual and computational progress, fundamental questions remain. Some are summarize herein.

The recent discoveries of giant planets around nearby stars have galvanized the planetary science community, astronomers, and the public at large. Since {\it direct} detection is now feasible, and is suggested by the recent acquisition of Gl229 B, it is crucial for the future of extrasolar planet searches that the fluxes, evolution, and physical structure of objects from Saturn's mass to 15 Juipter masses be theoretically investigated. We discuss our first attempts to explore the characteristics of extrasolar giant planets (EGPs), in aid of both NASA's and ESA's recent plans to search for such planets around nearby stars.

New insights into the mechanism and character of core--collapse supernova explosions are transforming the approach of theorists to their subject. The universal realization that the direct hydrodynamic mechanism does not work and that a variety of hydrodynamic instabilities can influence the viability of theoretical explosions has ushered in a new era in supernova modeling. In this paper, I discuss the important physical and technical issues that remain. I review the neutrino--driven mechanism, the possible roles of Rayleigh--Taylor instabilities, questions in neutrino transport, and the various observational constraints within which theorists must operate. However, a consensus has yet to be achieved among active workers concerning many important details and some essential phenomenology. This synopsis is meant to accomplish two things: 1) to focus attention on the interesting problems whose resolution will bring needed progress, and 2) to assess the current status of the theoretical art.

Theoretical spectra and evolutionary models that span the giant planet--brown dwarf continuum have been computed based on the recent discovery of the brown dwarf, Gliese 229 B. A flux enhancement in the 4--5 micron window is a universal feature from Jovian planets to brown dwarfs. We confirm the existence of methane and water in Gl 229 B's spectrum and find its mass to be 30 to 55 Jovian masses. Although these calculations focus on Gliese 229 B, they are also meant to guide future searches for extra-solar giant planets and brown dwarfs.

New data imply that the average velocity of radio pulsars is large \cite{hla93}. Under the assumption that these data imply that a pulsar is born with an ``intrinsic'' kick, we investigate whether such kicks can be a consequence of asymmetrical stellar collapse and explosion. We calculate the gravitational wave (GW) signature of such asymmetries due to anisotropic neutrino radiation and mass motions. We predict that any recoils imparted to the neutron star at birth will result in a gravitational wave strain, h$^{TT}_{zz}$, that does not go to zero with time. Hence, there may be ``memory'' \cite{bt87} in the gravitational waveform from a protoneutron star that is correlated with its recoil and neutrino emissions.

Using Doppler spectroscopy to detect the reflex motion of the nearby star, 51 Pegasi, Mayor \& Queloz (1995) claim to have discovered a giant planet in a 0.05-A.U., 4.23-day orbit. They estimate its mass to be in the range 0.5 \mj to 2 \mj, but are not able to determine its nature or origin. Including the effects of the severe stellar insolation implied, we extend the theory of giant planets we have recently developed to encompass those at very small orbital distances. Our calculations can be used to help formulate search strategies for luminous planets in tight orbits around other nearby stars. We calculate the radii and luminosities of such giants for a variety of compositions (H/He, He, H$_2$O, and olivine), evolutionary tracks for solar- composition gas giants, and the geometry of the Hayashi forbidden zone in the gas-giant mass regime. We show that such planets are stable and estimate the magnitude of classical Jeans evaporation and of photodissociation and loss due to EUV radiation. Even over the lifetime of the primary, the present companion would not have lost a large fraction of its mass. In addition, we demonstrate that for the mass range quoted, such planets are well within their Roche lobes. We show that the strong composition-dependence of the model radii and distinctive spectral signatures provide clear diagnostics that might reveal 51 Peg B's nature, should interferometric or adaptive-optics techniques ever succeed in photometrically separating planet from star.

We present a broad suite of models of extrasolar giant planets (EGP's), ranging in mass from 0.3 to 15 Jupiter masses. The models predict luminosity (both reflected and emitted) as a function of age, mass, deuterium abundance and distance from parent stars of various spectral type. We also explore the effects of helium mass fraction, rotation rate and the presence of a rock-ice core. The models incorporate the most accurate available equation of state for the interior, including a new theory for the enhancement of deuterium fusion by electron screening which is potentially important in these low mass objects. The results of our calculations reveal the enormous sensitivity of EGP's to the presence of the parent star, particularly for G and earlier spectral types. They also show a strong sensitivity of the flux contrast in the mid-infrared between parent star and EGP to the mass and age of the EGP's. We interpret our results in terms of search strategies for ground- and space-based observatories in place or anticipated in the near future.

We investigate in this paper the core-collapse supernova explosion mechanism in both one and two dimensions. We verify the usefulness of neutrino-driven overturn (``convection'') between the shock and the neutrinosphere in igniting the supernova explosion. The 2-D simulation of the core of a 15$\mdot$ star that we present here indicates that the breaking of spherical symmetry may be central to the explosion itself and that a multitude of bent and broken fingers is a common feature of the ejecta. In 2-D, the steady-state shock radius is larger than its value in 1-D by 30\%--100\%. The critical condition for explosion is thereby relaxed. Since supernovae are driven by neutrino heating, they are coronal phenomena, akin to winds, though initially bounded by an accretion tamp. Neutrino energy is pumped into the supernova during the shock's propagation through the inner many thousands of kilometers and not instantaneously. Curiously, just after the explosion is triggered, the matter that will eventually be ejected is still bound. The high-speed fingers that emerge from the core seem a natural explanation for the nickel bullets seen in 87A and the shrapnel inferred in some supernova remnants. Within 100 ms of the explosion, a strong, neutrino-driven wind is blowing outward from the protoneutron star.

We have performed Monte Carlo calculations of gamma-ray transport in models of the clumpy debris cloud of the LMC supernova, SN1987A, to study the influence of composition mixing and heterogeneity on its emergent gamma-ray and X-ray fluxes. In particular, we have focused on the problematic Ginga band (16 -- 28 keV) flux at day 600, whose measured value was an order of magnitude higher than predicted by previous theory. We find that the hydrogen of the envelope could not have been intimately mixed with the heavy elements of the core and that the hydrogen/helium volume filling factor interior to 4000 km/s must have been large ($\ge 40\%$). A physical separation of the scattering region and the regions occupied by the high-Z elements is required. The 600-day models that best fit both the line data at 847 keV and 1238 keV and the measured Ginga band fluxes suggest that as much as 50\% of the explosively produced $^{56}$Ni stayed interior to 1000 km s$^{-1}$ and 2 M$_\odot$. The $^{56}$Ni may have been more centrally-concentrated than in the benchmark models. $^{56}$Ni filling factors greater than 60\% are not preferred, since such models are too good at absorbing photons below 100 keV. Furthermore, a total envelope mass between 10 $\mdot$ and 15 $\mdot$ is favored.

It is now believed that pulsars comprise the fastest population of stars in the galaxy. With inferred mean, root-mean-square, and maximum 3-D pulsar speeds of $\sim$300-500 km/s, $\sim$500 km/s, and $\sim$2000 km/s, respectively, the question of the origin of such singular proper motions becomes acute. What mechanism can account for speeds that range from zero to twice the galactic escape velocity? We speculate that a major vector component of a neutron star's proper motion comes from the hydrodynamic recoil of the nascent neutron star during the supernova explosion in which it is born. Recently, theorists have shown that asymmetries and instabilities are a natural aspect of supernova dynamics. In this paper, we highlight two phenomena: 1) the ``Brownian-like'' stochastic motion of the core in response to the convective ``boiling'' of the mantle of the protoneutron star during the post-bounce, pre-explosion accretion phase, and 2) the asymmetrical bounce and explosion of an aspherically collapsing Chandrasekhar core.In principle,either phenomenon can leave the young neutron star with a speed of hundreds of kilometers per second.

Recent novel observations have probed the baryonic fraction of the galactic dark matter that has eluded astronomers for decades. Late in 1993, the MACHO and EROS collaborations announced in this journal the detection of transient and achromatic brightenings of a handful of stars in the Large Magellanic Cloud that are best interpreted as gravitational microlensing by low-mass foreground objects (MACHOS). This tantalized astronomers, for it implied that the population of cool, compact objects these lenses represent could be the elusive dark matter of our galactic halo. A year later in 1994, Sackett et al. reported the discovery of a red halo in the galaxy NGC 5907 that seems to follow the inferred radial distribution of its dark matter. This suggested that dwarf stars could constitute its missing component. Since NGC 5907 is similar to the Milky Way in type and radius, some surmised that the solution of the galactic dark matter problem was an abundance of ordinary low-mass stars. Now Bahcall et al., using the Wide-Field Camera of the recently repaired Hubble Space Telescope, have dashed this hope.

Interest among astronomers in the detection of extra-solar planets is accelerating with the growing realization that it may soon be technically feasible. The ongoing renaissance in telescope construction and the anticipated launches of new space platforms are encouraging many scientists to review and improve the means by which planets can be discovered. The direct detection of the light from a distant planet would be the most compelling means of discovery and to gauge the feasibility of various search strategies, astronomers have traditionally used the current Jupiter as a benchmark planet. However, in principle, extra-solar giant planets (EGPs) can have a wide range of masses and, hence, can be significantly brighter than Jupiter. Furthermore, the maximum mass a planet can have is not known a priori, and observations will be needed to determine it. We predict the optical and infrared fluxes of EGPs with masses from 0.3 through 15 Jupiter masses and ages from 10$^7$ through $5 \times 10^{9}$ years that searches in the next few years may reveal.