Optical Albedo Theory of Strongly-Irradiated Giant Planets: The Case of HD 209458b

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.Comment: 6 pages, 1 table, 1 color figure; accepted to the Astrophysical Journa

Theoretical Spectral Models of the Planet HD 209458b with a Thermal Inversion and Water Emission Bands

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.Comment: Accepted to the Astrophysical Journal Letters on August 28, 2007, six pages in emulateapj forma

Three years of global carbon monoxide from SCIAMACHY: comparison with MOPITT and first results related to the detection of enhanced CO over cities

Carbon monoxide (CO) is an important atmospheric constituent affecting air quality and climate. SCIAMACHY on ENVISAT is currently the only satellite instrument that can measure the vertical column of CO with nearly equal sensitivity at all altitudes down to the Earth's surface because of its near-infrared nadir observations of reflected solar radiation. Here we present three years' (2003–2005) of SCIAMACHY CO columns consistently retrieved with the latest version of our retrieval algorithm (WFMDv0.6). We describe the retrieval method and discuss the multi-year global CO data set focusing on a comparison with the operational CO column data product of MOPITT. We found reasonable to good agreement (~20%) with MOPITT, with the best agreement for 2004. We present detailed results for various regions (Europe, Middle East, India, China) and discuss to what extent enhanced levels of CO can be detected over populated areas including individual cities. The expected CO signal from cities is close to or even below the detection limit of individual measurements. We show that cities can be identified when averaging long time series

Microscopic calculation of neutrino mean free path inside hot neutron matter

We calculate the neutrino mean free path and the Equation of State of pure neutron matter at finite temperature within a selfconsistent scheme based on the Brueckner--Hartree--Fock approximation. We employ the nucleon-nucleon part of the recent realistic baryon-baryon interaction (model NSC97e) constructed by the Nijmegen group. The temperatures considered range from 10 to 80 MeV. We report on the calculation of the mean field, the residual interaction and the neutrino mean free path including short and long range correlations given by the Brueckner--Hartree--Fock plus Random Phase Approximation (BHF+RPA) framework. This is the first fully consistent calculation in hot neutron matter dedicated to neutrino mean free path. We compare systematically our results to those obtain with the D1P Gogny effective interaction, which is independent of the temperature. The main differences between the present calculation and those with nuclear effective interactions come from the RPA corrections to BHF (a factor of about 8) while the temperature lack of consistency accounts for a factor of about 2

Conservative formulations of general relativistic kinetic theory

Experience with core-collapse supernova simulations shows that accurate accounting of total particle number and 4-momentum can be a challenge for computational radiative transfer. This accurate accounting would be facilitated by the use of particle number and 4-momentum transport equations that allow transparent conversion between volume and surface integrals in both configuration and momentum space. Such conservative formulations of general relativistic kinetic theory in multiple spatial dimensions are presented in this paper, and their relevance to core-collapse supernova simulations is described.Comment: 48 page

Crucial Physical Dependencies of the Core-Collapse Supernova Mechanism

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.Comment: 25 pages; 10 figure

Evolutionary models for low-mass stars and brown dwarfs: uncertainties and limits at very young ages

We analyse pre-Main Sequence evolutionary tracks for low mass stars with masses m \le 1.4 \msol based on the Baraffe et al. (1998) input physics. We also extend the recent Chabrier et al. (2000) evolutionary models based on dusty atmosphere to young brown dwarfs down to one mass of Jupiter. We analyse current theoretical uncertainties due to molecular line lists, convection and initial conditions. Simple tests on initial conditions show the high uncertainties of models at ages \simle 1 Myr. We find a significant sensitivity of atmosphere profiles to the treatment of convection at low gravity and \te < 4000 K, whereas it vanishes as gravity increases. This effect adds another source of uncertainty on evolutionary tracks at very early phases. We show that at low surface gravity (\log g \simle 3.5,) the common picture of vertical Hayashi lines with constant \te is oversimplified. The effect of a variation of initial deuterium abundance is studied. We compare our models with evolutionary tracks available in the literature and discuss the main differences. We finally analyse to which extent current observations of young systems provide a good test for pre-Main Sequence tracks.Comment: 12 pages, Latex file, uses aa.cls, accepted for publication in A&

A New Algorithm for Supernova Neutrino Transport and Some Applications

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.Comment: 58 pages, single-spaced LaTeX, 23 figures, revised title, also available at http://jupiter.as.arizona.edu/~burrows/papers, accepted for publication in the Ap.

Anisotropies in the Neutrino Fluxes and Heating Profiles in Two-dimensional, Time-dependent, Multi-group Radiation Hydrodynamics Simulations of Rotating Core-Collapse Supernovae

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.Comment: Includes 11 low-resolution color figures, accepted to the Astrophysical Journal (June 10, 2005; V. 626); high-resolution figures and movies available from the authors upon reques

Operator method in solving non-linear equations of the Hartree-Fock type

The operator method is used to construct the solutions of the problem of the polaron in the strong coupling limit and of the helium atom on the basis of the Hartree-Fock equation. $E_0=-0.1085128052\alpha^2$ is obtained for the polaron ground-state energy. Energies for 2s- and 3s-states are also calculated. The other excited states are briefly discussed.Comment: 7 page