We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.

Crater ellipticity determination is a complex and time consuming task that so far has evaded successful automation. We train a state of the art computer vision algorithm to identify craters in Lunar digital elevation maps and retrieve their sizes and 2D shapes. The computational backbone of the model is MaskRCNN, an "instance segmentation" general framework that detects craters in an image while simultaneously producing a mask for each crater that traces its outer rim. Our post-processing pipeline then finds the closest fitting ellipse to these masks, allowing us to retrieve the crater ellipticities. Our model is able to correctly identify 87\% of known craters in the longitude range we hid from the network during training and validation (test set), while predicting thousands of additional craters not present in our training data. Manual validation of a subset of these "new" craters indicates that a majority of them are real, which we take as an indicator of the strength of our model in learning to identify craters, despite incomplete training data. The crater size, ellipticity, and depth distributions predicted by our model are consistent with human-generated results. The model allows us to perform a large scale search for differences in crater diameter and shape distributions between the lunar highlands and maria, and we exclude any such differences with a high statistical significance. The predicted test set catalogue and trained model are available here: https://github.com/malidib/Craters_MaskRCNN/.

We present calculations of the occurrence rate of small close-in planets around low mass dwarf stars using the known planet populations from the $Kepler$ and $K2$ missions. Applying completeness corrections clearly reveals the radius valley in the maximum a-posteriori occurrence rates as a function of orbital separation and planet radius. We measure the slope of the valley to be $r_{p,\text{valley}} \propto F^{-0.060\pm 0.025}$ which bears the opposite sign from that measured around Sun-like stars thus suggesting that thermally driven atmospheric mass loss may not dominate the evolution of planets in the low stellar mass regime or that we are witnessing the emergence of a separate channel of planet formation. The latter notion is supported by the relative occurrence of rocky to non-rocky planets increasing from $0.5\pm 0.1$ around mid-K dwarfs to $8.5\pm 4.6$ around mid-M dwarfs. Furthermore, the center of the radius valley at $1.54\pm 0.16$ R$_{\oplus}$ is shown to shift to smaller sizes with decreasing stellar mass in agreement with physical models of photoevaporation, core-powered mass loss, and gas-poor formation. Although current measurements are insufficient to robustly identify the dominant formation pathway of the radius valley, such inferences may be obtained by $TESS$ with $\mathcal{O}(85,000)$ mid-to-late M dwarfs observed with 2-minute cadence. The measurements presented herein also precisely designate the subset of planetary orbital periods and radii that should be targeted in radial velocity surveys to resolve the rocky to non-rocky transition around low mass stars.

Global Circulation Models (GCMs) of atmospheric flows are now routinely used to interpret observational data on Hot Jupiters. Localized "equatorial $\beta$-plane" simulations by Fromang et al. (2016) have revealed that a barotropic (horizontal shear) instability of the equatorial jet appears at horizontal resolutions beyond those typically achieved in global models; this instability could limit wind speeds and lead to increased atmospheric variability. To address this possibility, we adapt the computationally efficient, pseudo-spectral PlaSim GCM, originally designed for Earth studies, to model Hot Jupiter atmospheric flows and validate it on the Heng et al. (2011) reference benchmark. We then present high resolution global models of HD209458b, with horizontal resolutions of T85 (128x256) and T127 (192x384). The barotropic instability phenomenology found in $\beta$-plane simulations is not reproduced in these global models, despite comparably high resolutions. Nevertheless, high resolution models do exhibit additional flow variability on long timescales (of order 100 planet days or more), which is absent from the lower resolution models. It manifests as a breakdown of north-south symmetry of the equatorial wind. From post-processing the atmospheric flows at various resolutions (assuming a cloud-free situation), we show that the stronger flow variability achieved at high resolution does not translate into noticeably stronger dayside infrared flux variability. More generally, our results suggest that high horizontal resolutions are not required to capture the key features of hot Jupiter atmospheric flows.

A large number of studies have responded to the growing body of confirmed terrestrial habitable zone exoplanets by presenting models of various possible climates. However, the impact of the partial pressure of background gases such as N$_2$ has been poorly-explored, despite the abundance of N$_2$ in Earth's atmosphere and the lack of constraints on its typical abundance in terrestrial planet atmospheres. We use PlaSim, a fast 3D climate model, to simulate many hundreds of climates with varying N$_2$ partial pressures, insolations, and surface characteristics to identify the impact of the background gas partial pressure on the climate. We find that the climate's response is nonlinear and highly sensitive to the background gas partial pressure. We identify pressure broadening of CO$_2$ and H$_2$O absorption lines, amplification of warming or cooling by the water vapor greenhouse positive feedback, heat transport efficiency, and cooling through Rayleigh scattering as the dominant competing mechanisms that determine the equilibrium climate for a given N$_2$ partial pressure. Finally, we show that different amounts of N$_2$ should have a significant effect on broadband reflected light observations of terrestrial exoplanets.

Machine learning (ML) is a standard approach for estimating the redshifts of galaxies when only photometric information is available. ML photo-z solutions have traditionally ignored the morphological information available in galaxy images or partly included it in the form of hand-crafted features, with mixed results. We train a morphology-aware photometric redshift machine using modern deep learning tools. It uses a custom architecture that jointly trains on galaxy fluxes, colors and images. Galaxy-integrated quantities are fed to a Multi-Layer Perceptron (MLP) branch while images are fed to a convolutional (convnet) branch that can learn relevant morphological features. This split MLP-convnet architecture, which aims to disentangle strong photometric features from comparatively weak morphological ones, proves important for strong performance: a regular convnet-only architecture, while exposed to all available photometric information in images, delivers comparatively poor performance. We present a cross-validated MLP-convnet model trained on 130,000 SDSS-DR12 galaxies that outperforms a hyperoptimized Gradient Boosting solution (hyperopt+XGBoost), as well as the equivalent MLP-only architecture, on the redshift bias metric. The 4-fold cross-validated MLP-convnet model achieves a bias $\delta z / (1+z) =-0.70 \pm 1 \times 10^{-3} $, approaching the performance of a reference ANNZ2 ensemble of 100 distinct models trained on a comparable dataset. The relative performance of the morphology-aware and morphology-blind models indicates that galaxy morphology does improve ML-based photometric redshift estimation.

Habitable planets are commonly imagined to be temperate planets like Earth, with areas of open ocean and warm land. In contrast, planets in snowball states, where oceans are entirely ice-covered, are believed to be inhospitable. However, we show using a general circulation model that terrestrial planets in the inner habitable zone are able to support large unfrozen areas of land while in a snowball state. Due to their lower albedo, these unfrozen regions reach summer temperatures in excess of 10 $^\circ$Celsius. Such conditions permit CO$_2$ weathering, suggesting that continental weathering can provide a mechanism for trapping planets in stable snowball states. The presence of land areas with warm temperatures and liquid surface water motivates a more-nuanced understanding of habitability during these snowball events.

L 98-59 (TIC 307210830, TOI-175) is a nearby M3 dwarf around which TESS revealed three terrestrial-sized transiting planets (0.80, 1.35, 1.57 Earth radii) in a compact configuration with orbital periods shorter than 7.5 days. Here we aim to measure the masses of the known transiting planets in this system using precise radial velocity (RV) measurements taken with the HARPS spectrograph. We consider both trained and untrained Gaussian process regression models of stellar activity to simultaneously model the RV data with the planetary signals. Our RV analysis is then supplemented with dynamical simulations to provide strong constraints on the planets' orbital eccentricities by requiring long-term stability. We measure the planet masses of the two outermost planets to be $2.46\pm 0.31$ and $2.26\pm 0.50$ Earth masses which confirms their bulk terrestrial compositions. We are able to place an upper limit on the mass of the smallest, innermost planet of $<0.98$ Earth masses with 95% confidence. Our RV plus dynamical stability analysis places strong constraints on the orbital eccentricities and reveals that each planet's orbit likely has $e<0.1$ to ensure a dynamically stable system. The L 98-59 compact system of three likely rocky planets offers a unique laboratory for studies of planet formation, dynamical stability, and comparative atmospheric planetology. Continued RV monitoring will help refine the characterization of the innermost planet and potentially reveal additional planets in the system at wider separations.

Turbulent vertical transport driven by double-diffusive shear instabilities is identified as likely important in hot exoplanet atmospheres. In hot Jupiter atmospheres, the resulting vertical mixing appears sufficient to alleviate the nightside cold trap, thus facilitating the maintenance of nocturnal clouds on these planets. The strong level of vertical mixing expected near hot Jupiter thermal photospheres will impact their atmospheric chemistry and even their vertical structures where cloud radiative feedback proves important.

Crater counting on the Moon and other bodies is crucial to constrain the dynamical history of the Solar System. This has traditionally been done by visual inspection of images, thus limiting the scope, efficiency, and/or accuracy of retrieval. In this paper we demonstrate the viability of using convolutional neural networks (CNNs) to determine the positions and sizes of craters from Lunar digital elevation maps (DEMs). We recover 92% of craters from the human-generated test set and almost double the total number of crater detections. Of these new craters, 15% are smaller in diameter than the minimum crater size in the ground-truth dataset. Our median fractional longitude, latitude and radius errors are 11% or less, representing good agreement with the human-generated datasets. From a manual inspection of 361 new craters we estimate the false positive rate of new craters to be 11%. Moreover, our Moon-trained CNN performs well when tested on DEM images of Mercury, detecting a large fraction of craters in each map. Our results suggest that deep learning will be a useful tool for rapidly and automatically extracting craters on various Solar System bodies. We make our code and data publicly available at https://github.com/silburt/DeepMoon.git and https://doi.org/10.5281/zenodo.1133969 .

In an earlier campaign to characterize the mass of the transiting temperate super-Earth K2-18b with HARPS, a second, non-transiting planet was posited to exist in the system at $\sim 9$ days. Further radial velocity follow-up with the CARMENES spectrograph visible channel revealed a much weaker signal at 9 days which also appeared to vary chromatically and temporally leading to the conclusion that the origin of the 9 day signal was more likely to be related to stellar activity than to being planetary. Here we conduct a detailed re-analysis of all available RV time-series, including a set of 31 previously unpublished HARPS measurements, to investigate the effects of time-sampling and of simultaneous modelling of planetary + activity signals on the existence and origin of the curious 9 day signal. We conclude that the 9 day signal is real and was initially seen to be suppressed in the CARMENES data due to a small number of anomalous measurements, although the exact cause of these anomalies remains unknown. Investigation of the signal's evolution in time, with wavelength, and detailed model comparison reveals that the 9 day signal is most likely planetary in nature. By this analysis, we reconcile the conflicting HARPS and CARMENES results and measure precise and self-consistent planet masses of $m_{p,b} = 8.63 \pm 1.35$ and $m_{p,c}\sin{i_c}=5.62 \pm 0.84$ M$_{\oplus}$. This work, along with the previously published RV papers on the K2-18 planetary system, highlight the importance of understanding one's time-sampling and of simultaneous planet + stochastic activity modelling, particularly when searching for sub-Neptune-sized planets with radial velocities.

GJ1132 is a nearby red dwarf known to host a transiting Earth-size planet. After its initial detection, we pursued an intense follow-up with the HARPS velocimeter. We now confirm the detection of GJ1132b with radial velocities only. We refined its orbital parameters and, in particular, its mass ($m_b = 1.66\pm0.23 M_\oplus$), density ($\rho_b = 6.3\pm1.3$ g.cm$^{-3}$) and eccentricity ($e_b < 0.22 $; 95\%). We also detect at least one more planet in the system. GJ1132c is a super-Earth with period $P_c = 8.93\pm0.01$ days and minimum mass $m_c \sin i_c = 2.64\pm0.44~M_\oplus$. Receiving about 1.9 times more flux than Earth in our solar system, its equilibrium temperature is that of a temperate planet ($T_{eq}=230-300$ K for albedos $A=0.75-0.00$) and places GJ1132c near the inner edge of the so-called habitable zone. Despite an a priori favourable orientation for the system, $Spitzer$ observations reject most transit configurations, leaving a posterior probability $<1\%$ that GJ1132c transits. GJ1132(d) is a third signal with period $P_d = 177\pm5$ days attributed to either a planet candidate with minimum mass $m_d \sin i_d = 8.4^{+1.7}_{-2.5}~M_\oplus$ or stellar activity. (abridged)

The SPIRou near infrared spectro-polarimeter is destined to begin science operations at the Canada-France-Hawaii Telescope in mid-2018. One of the instrument's primary science goals is to discover the closest exoplanets to the Solar System by conducting a 3-5 year long radial velocity survey of nearby M dwarfs at an expected precision of $\sim 1$ m s$^{-1}$; the SPIRou Legacy Survey-Planet Search (SLS-PS). In this study we conduct a detailed Monte-Carlo simulation of the SLS-PS using our current understanding of the occurrence rate of M dwarf planetary systems and physical models of stellar activity. From simultaneous modelling of planetary signals and activity, we predict the population of planets detected in the SLS-PS. With our fiducial survey strategy and expected instrument performance over a nominal survey length of $\sim 3$ years, we expect SPIRou to detect $85.3^{+29.3}_{-12.4}$ planets including $20.0^{+16.8}_{-7.2}$ habitable zone planets and $8.1^{+7.6}_{-3.2}$ Earth-like planets from a sample of 100 M1-M8.5 dwarfs out to 11 pc. By studying mid-to-late M dwarfs previously inaccessible to existing optical velocimeters, SPIRou will put meaningful constraints on the occurrence rate of planets around those stars including the value of $\eta_{\oplus}$ at an expected level of precision of $\lesssim 45$%. We also predict a subset of $46.7^{+16.0}_{-6.0}$ planets may be accessible with dedicated high-contrast imagers on the next generation of ELTs including $4.9^{+4.7}_{-2.0}$ potentially imagable Earth-like planets. Lastly, we compare the results of our fiducial survey strategy to other foreseeable survey versions to quantify which strategy is optimized to reach the SLS-PS science goals. The results of our simulations are made available to the community on github.

It has recently been proposed that Earth-like planets in the outer regions of the habitable zone experience unstable climates, repeatedly cycling between glaciated and deglaciated climatic states (Menou 2015). While this result has been confirmed and also extended to explain early Mars climate records (Haqq-Misra et al. 2016; Batalha et al. 2016), all existing work relies on highly idealized low-dimensional climate models. Here, we confirm that the phenomenology of climate cycles remains in 3D Earth climate models with considerably more degrees of freedom. To circumvent the computational barrier of integrating climate on Gyr timescales, we design a hybrid 0D-3D integrator which uses a general circulation model (GCM) as a short relaxation step along a long evolutionary climate sequence. We find that GCM climate cycles are qualitatively consistent with reported low-dimensional results. This establishes on a firmer ground the notion that outer habitable zone planets may be preferentially found in transiently glaciated states.

The TRAPPIST-1, Proxima Centauri, and LHS 1140 systems are the most exciting prospects for future follow-up observations of potentially inhabited planets. All orbit nearby M-stars and are likely tidally locked in 1:1 spin-orbit states, which motivates the consideration of the effects that tidal locking might have on planetary habitability. On Earth, periods of global glaciation (snowballs) may have been essential for habitability and remote signs of life (biosignatures) because they are correlated with increases in the complexity of life and in the atmospheric oxygen concentration. In this paper we investigate the snowball bifurcation (sudden onset of global glaciation) on tidally locked planets using both an energy balance model and an intermediate-complexity global climate model. We show that tidally locked planets are unlikely to exhibit a snowball bifurcation as a direct result of the spatial pattern of insolation they receive. Instead they will smoothly transition from partial to complete ice coverage and back. A major implication of this work is that tidally locked planets with an active carbon cycle should not be found in a snowball state. Moreover, this work implies that tidally locked planets near the outer edge of the habitable zone with low CO2 outgassing fluxes will equilibrate with a small unglaciated substellar region rather than cycling between warm and snowball states. More work is needed to determine how the lack of a snowball bifurcation might affect the development of life on a tidally locked planet.

The bright M dwarf K2-18 at 34 pc is known to host a transiting super-Earth-sized planet orbiting within the star's habitable zone; K2-18b. Given the superlative nature of this system for studying an exoplanetary atmosphere receiving similar levels of insolation as the Earth, we aim to characterize the planet's mass which is required to interpret atmospheric properties and infer the planet's bulk composition. We obtain precision radial velocity measurements with the HARPS spectrograph and couple those measurements with the K2 photometry to jointly model the observed radial velocity variation with planetary signals and a radial velocity jitter model based on Gaussian process regression. We measure the mass of K2-18b to be $8.0 \pm 1.9$ M$_{\oplus}$ with a bulk density of $3.7 \pm 0.9$ g/cm$^3$ which may correspond to a predominantly rocky planet with a significant gaseous envelope or an ocean planet with a water mass fraction $\gtrsim 50$%. We also find strong evidence for a second, warm super-Earth K2-18c at $\sim 9$ days with a semi-major axis 2.4 times smaller than the transiting K2-18b. After re-analyzing the available light curves of K2-18 we conclude that K2-18c is not detected in transit and therefore likely has an orbit that is non-coplanar with K2-18b. A suite of dynamical integrations with varying simulated orbital eccentricities of the two planets are used to further constrain each planet's eccentricity posterior from which we measure $e_b < 0.43$ and $e_c < 0.47$ at 99% confidence. The discovery of the inner planet K2-18c further emphasizes the prevalence of multi-planet systems around M dwarfs. The characterization of the density of K2-18b reveals that the planet likely has a thick gaseous envelope which along with its proximity to the Solar system makes the K2-18 planetary system an interesting target for the atmospheric study of an exoplanet receiving Earth-like insolation.

The requirement that planetary systems be dynamically stable is often used to vet new discoveries or set limits on unconstrained masses or orbital elements. This is typically carried out via computationally expensive N-body simulations. We show that characterizing the complicated and multi-dimensional stability boundary of tightly packed systems is amenable to machine learning methods. We find that training an XGBoost machine learning algorithm on physically motivated features yields an accurate classifier of stability in packed systems. On the stability timescale investigated ($10^7$ orbits), it is 3 orders of magnitude faster than direct N-body simulations. Optimized machine learning classifiers for dynamical stability may thus prove useful across the discipline, e.g., to characterize the exoplanet sample discovered by the upcoming Transiting Exoplanet Survey Satellite (TESS). This proof of concept motivates investing computational resources to train algorithms capable of predicting stability over longer timescales and over broader regions of phase space.

M-dwarfs are known to commonly host high-multiplicity planetary systems. Therefore M-dwarf planetary systems with a known transiting planet are expected to contain additional small planets ($r_p \le 4$ R$_{\oplus}$, $m_p \lesssim 20$ M$_{\oplus}$) that are not seen in transit. In this study we investigate the effort required to detect such planets using precision velocimetry around the sizable subset of M-dwarfs which are slowly rotating ($P_{\mathrm{rot}} \gtrsim 40$ days) and hence more likely to be inactive. We focus on the test case of GJ 1132. Specifically, we perform a suite of Monte-Carlo simulations of the star's radial velocity signal featuring astrophysical contributions from stellar jitter due to rotationally modulated active regions and keplarian signals from the known transiting planet and hypothetical additional planets not seen in transit. We then compute the detection completeness of non-transiting planets around GJ 1132 and consequently estimate the number of RV measurements required to detect those planets. We show that with 1 m s$^{-1}$ RV precision per measurement, only $\sim 50$ measurements are required to achieve a 50\% detection completeness to all non-transiting planets in the system and to planets which are potentially habitable. Throughout we advocate the use of Gaussian process regression as an effective tool for mitigating the effects of stellar jitter including stars with high activity. Given that GJ 1132 is representative of a large population of slowly rotating M-dwarfs, we conclude with a discussion of how our results may be extended to other systems with known transiting planets such as those which will be discovered with TESS.

Massive giant planets, such as the ones being discovered by direct imaging surveys, likely experience the majority of their growth through a circumplanetary disc. We argue that the entropy of accreted material is determined by boundary layer processes, unlike the "cold-" or "hot-start" hypotheses usually invoked in the core accretion and direct collapse scenarios. A simple planetary evolution model illustrates how a wide range of radius and luminosity tracks become possible, depending on details of the accretion process. Specifically, the proto-planet evolves towards "hot-start" tracks if the scale-height of the boundary layer is $\gtrsim0.24$, a value not much larger than the scale-height of the circumplanetary disc. Understanding the luminosity and radii of young giant planets will thus require detailed models of circumplanetary accretion.

A recent ALMA image revealed several concentric gaps in the protoplanetary disk surrounding the young star HL Tau. We consider the hypothesis that these gaps are carved by planets, and present a general framework for understanding the dynamical stability of such systems over typical disk lifetimes, providing estimates for the maximum planetary masses. We collect these easily evaluated constraints into a workflow that can help guide the design and interpretation of new observational campaigns and numerical simulations of gap opening in such systems. We argue that the locations of resonances should be significantly shifted in massive disks like HL Tau, and that theoretical uncertainties in the exact offset, together with observational errors, imply a large uncertainty in the dynamical state and stability in such disks. This presents an important barrier to using systems like HL Tau as a proxy for the initial conditions following planet formation. An important observational avenue to breaking this degeneracy is to search for eccentric gaps, which could implicate resonantly interacting planets. Unfortunately, massive disks like HL Tau should induce swift pericenter precession that would smear out any such eccentric features of planetary origin. This motivates pushing toward more typical, less massive disks. For a nominal non-resonant model of the HL Tau system with five planets, we find a maximum mass for the outer three bodies of approximately 2 Neptune masses. In a resonant configuration, these planets can reach at least the mass of Saturn. The inner two planets' masses are unconstrained by dynamical stability arguments.

Planets in the habitable zone of lower-mass stars are often assumed to be in a state of tidally synchronized rotation, which would considerably affect their putative habitability. Although thermal tides cause Venus to rotate retrogradely, simple scaling arguments tend to attribute this peculiarity to the massive Venusian atmosphere. Using a global climate model, we show that even a relatively thin atmosphere can drive terrestrial planets' rotation away from synchronicity. We derive a more realistic atmospheric tide model that predicts four asynchronous equilibrium spin states, two being stable, when the amplitude of the thermal tide exceeds a threshold that is met for habitable Earth-like planets with a 1-bar atmosphere around stars more massive than 0.5-0.7Msun. Thus, many recently discovered terrestrial planets could exhibit asynchronous spin-orbit rotation, even with a thin atmosphere.

The carbon-silicate cycle regulates the atmospheric $CO_2$ content of terrestrial planets on geological timescales through a balance between the rates of $CO_2$ volcanic outgassing and planetary intake from rock weathering. It is thought to act as an efficient climatic thermostat on Earth and, by extension, on other habitable planets. If, however, the weathering rate increases with the atmospheric $CO_2$ content, as expected on planets lacking land vascular plants, the carbon-silicate cycle feedback can become severely limited. Here we show that Earth-like planets receiving less sunlight than current Earth may no longer possess a stable warm climate but instead repeatedly cycle between unstable glaciated and deglaciated climatic states. This has implications for the search for life on exoplanets in the habitable zone of nearby stars.

Although tidally-locked habitable planets orbiting nearby M-dwarf stars are among the best astronomical targets to search for extrasolar life, they may also be deficient in volatiles and water. Climate models for this class of planets show atmospheric transport of water from the dayside to the nightside, where it is precipitated as snow and trapped as ice. Since ice only slowly flows back to the dayside upon accumulation, the resulting hydrological cycle can trap a large amount of water in the form of nightside ice. Using ice sheet dynamical and thermodynamical constraints, I illustrate how planets with less than about a quarter the Earth's oceans could trap most of their surface water on the nightside. This would leave their dayside, where habitable conditions are met, potentially dry. The amount and distribution of residual liquid water on the dayside depend on a variety of geophysical factors, including the efficiency of rock weathering at regulating atmospheric CO2 as dayside ocean basins dry-up. Water-trapped worlds with dry daysides may offer similar advantages as land planets for habitability, by contrast with worlds where more abundant water freely flows around the globe.

We present the first three-dimensional circulation models for extrasolar gas giant atmospheres with geometrically and energetically consistent treatments of magnetic drag and ohmic dissipation. Atmospheric resistivities are continuously updated and calculated directly from the flow structure, strongly coupling the magnetic effects with the circulation pattern. We model the hot Jupiters HD 189733b (Teq \approx 1200 K) and HD 209458b (Teq \approx 1500 K) and test planetary magnetic field strengths from 0 to 30 G. We find that even at B = 3 G the atmospheric structure and circulation of HD 209458b are strongly influenced by magnetic effects, while the cooler HD 189733b remains largely unaffected, even in the case of B = 30 G and super-solar metallicities. Our models of HD 209458b indicate that magnetic effects can substantially slow down atmospheric winds, change circulation and temperature patterns, and alter observable properties. These models establish that longitudinal and latitudinal hot spot offsets, day-night flux contrasts, and planetary radius inflation are interrelated diagnostics of the magnetic induction process occurring in the atmospheres of hot Jupiters and other similarly forced exoplanets. Most of the ohmic heating occurs high in the atmosphere and on the day side of the planet, while the heating at depth is strongly dependent on the internal heat flux assumed for the planet, with more heating when the deep atmosphere is hot. We compare the ohmic power at depth in our models, and estimates of the ohmic dissipation in the bulk interior (from general scaling laws), to evolutionary models that constrain the amount of heating necessary to explain the inflated radius of HD 209458b. Our results suggest that deep ohmic heating can successfully inflate the radius of HD 209458b for planetary magnetic field strengths of B \geq 3 - 10 G.

The atmospheres of hot Jupiters and other strongly-forced exoplanets are susceptible to a thermal instability in the presence of ohmic dissipation, weak magnetic drag and strong winds. The instability occurs in radiatively-dominated atmospheric regions when the ohmic dissipation rate increases with temperature faster than the radiative (cooling) rate. The instability domain covers a specific range of atmospheric pressures and temperatures, typically P ~ 3-300 mbar and T ~ 1500-2500K for hot Jupiters, which makes it a candidate mechanism to explain the dayside thermal "inversions" inferred for a number of such exoplanets. The instability is suppressed by high levels of non-thermal photoionization, in possible agreement with a recently established observational trend. We highlight several shortcomings of the instability treatment presented here. Understanding the emergence and outcome of the instability, which should result in locally hotter atmospheres with stronger levels of drag, will require global non-linear atmospheric models with adequate MHD prescriptions.

Due to the gravitational influence of density fluctuations driven by magneto-rotational instability in the gas disk, planetesimals and protoplanets undergo diffusive radial migration as well as changes in other orbital properties. The magnitude of the effect on particle orbits can have important consequences for planet formation scenarios. We use the local-shearing-box approximation to simulate an ideal, isothermal, magnetized gas disk with vertical density stratification and simultaneously evolve numerous massless particles moving under the gravitational field of the gas and the host star. We measure the evolution of the particle orbital properties, including mean radius, eccentricity, inclination, and velocity dispersion, and its dependence on the disk properties and the particle initial conditions. Although the results converge with resolution for fixed box dimensions, we find the response of the particles to the gravity of the turbulent gas correlates with the horizontal box size, up to 16 disk scale heights. This correlation indicates that caution should be exercised when interpreting local-shearing-box models involving gravitational physics of magneto-rotational turbulence. Based on heuristic arguments, nevertheless, the criterion L_h / R ~ O(1), where L_h is the horizontal box size and R is the distance to the host star, is proposed to possibly circumvent this conundrum. If this criterion holds, we can still conclude that magneto-rotational turbulence seems likely to be ineffective at driving either diffusive migration or collisional erosion under most circumstances.

We present a new version of our code for modeling the atmospheric circulation on gaseous exoplanets, now employing a "double-gray" radiative transfer scheme, which self-consistently solves for fluxes and heating throughout the atmosphere, including the emerging (observable) infrared flux. We separate the radiation into infrared and optical components, each with its own absorption coefficient, and solve standard two-stream radiative transfer equations. We use a constant optical absorption coefficient, while the infrared coefficient can scale as a powerlaw with pressure. Here we describe our new code in detail and demonstrate its utility by presenting a generic hot Jupiter model. We discuss issues related to modeling the deepest pressures of the atmosphere and describe our use of the diffusion approximation for radiative fluxes at high optical depths. In addition, we present new models using a simple form for magnetic drag on the atmosphere. We calculate emitted thermal phase curves and find that our drag-free model has the brightest region of the atmosphere offset by ~12 degrees from the substellar point and a minimum flux that is 17% of the maximum, while the model with the strongest magnetic drag has an offset of only ~2 degrees and a ratio of 13%. Finally, we calculate rates of numerical loss of kinetic energy at ~15% for every model except for our strong-drag model, where there is no measurable loss; we speculate that this is due to the much decreased wind speeds in that model.

The exoplanet GJ1214b presents an interesting example of compositional degeneracy for low-mass planets. Its atmosphere may be composed of water, super-solar or solar metallicity material. We present atmospheric circulation models of GJ1214b for these three compositions, with explicit grey radiative transfer and an optional treatment of MHD bottom drag. All models develop strong, superrotating zonal winds (~ 1-2 km/s). The degree of eastward heat advection, which can be inferred from secondary eclipse and thermal phase curve measurements, varies greatly between the models. These differences are understood as resulting from variations in the radiative times at the thermal photosphere, caused by separate molecular weight and opacity effects. Our GJ1214b models illustrate how atmospheric circulation can be used as a probe of composition for similar tidally-locked exoplanets in the mini-Neptune/waterworld class.

Pulsar timing arrays (PTAs) are expected to detect gravitational waves (GWs) from individual low-redshift (z<1.5) compact supermassive (M>10^9 Msun) black hole (SMBH) binaries with orbital periods of approx. 0.1 - 10 yrs. Identifying the electromagnetic (EM) counterparts of these sources would provide confirmation of putative direct detections of GWs, present a rare opportunity to study the environments of compact SMBH binaries, and could enable the use of these sources as standard sirens for cosmology. Here we consider the feasibility of such an EM identification. We show that because the host galaxies of resolved PTA sources are expected to be exceptionally massive and rare, it should be possible to find unique hosts of resolved sources out to redshift z=0.2. At higher redshifts, the PTA error boxes are larger, and may contain as many as 100 massive-galaxy interlopers. The number of candidates, however, remains tractable for follow-up searches in upcoming wide-field EM surveys. We develop a toy model to characterize the dynamics and the thermal emission from a geometrically thin, gaseous disc accreting onto a PTA-source SMBH binary. Our model predicts that at optical and infrared frequencies, the source should appear similar to a typical luminous active galactic nucleus (AGN). However, owing to the evacuation of the accretion flow by the binary's tidal torques, the source might have an unusually low soft X-ray luminosity and weak UV and broad optical emission lines, as compared to an AGN powered by a single SMBH with the same total mass. For sources near z=1, the decrement in the rest-frame UV should be observable as an extremely red optical color. These properties would make the PTA sources stand out among optically luminous AGN, and could allow their unique identification.

We present scaling laws for advection, radiation, magnetic drag and ohmic dissipation in the atmospheres of hot giant exoplanets. In the limit of weak thermal ionization, ohmic dissipation increases with the planetary equilibrium temperature (T_eq >~ 1000 K) faster than the insolation power does, eventually reaching values >~ 1% of the insolation power, which may be sufficient to inflate the radii of hot Jupiters. At higher T_eq values still, magnetic drag rapidly brakes the atmospheric winds, which reduces the associated ohmic dissipation power. For example, for a planetary field strength B=10G, the fiducial scaling laws indicate that ohmic dissipation exceeds 1% of the insolation power over the equilibrium temperature range T_eq ~ 1300-2000 K, with a peak contribution at T_eq ~ 1600 K. Evidence for magnetically dragged winds at the planetary thermal photosphere could emerge in the form of reduced longitudinal offsets for the dayside infrared hotspot. This suggests the possibility of an anticorrelation between the amount of hotspot offset and the degree of radius inflation, linking the atmospheric and interior properties of hot giant exoplanets in an observationally testable way. While providing a useful framework to explore the magnetic scenario, the scaling laws also reveal strong parameter dependencies, in particular with respect to the unknown planetary magnetic field strength.

Hot super-Earths likely possess minimal atmospheres established through vapor saturation equilibrium with the ground. We solve the hydrodynamics of these tenuous atmospheres at the surface of Corot-7b, Kepler 10b and 55 Cnc-e, including idealized treatments of magnetic drag and ohmic dissipation. We find that atmospheric pressures remain close to their local saturation values in all cases. Despite the emergence of strongly supersonic winds which carry sublimating mass away from the substellar point, the atmospheres do not extend much beyond the day-night terminators. Ground temperatures, which determine the planetary thermal (infrared) signature, are largely unaffected by exchanges with the atmosphere and thus follow the effective irradiation pattern. Atmospheric temperatures, however, which control cloud condensation and thus albedo properties, can deviate substantially from the irradiation pattern. Magnetic drag and ohmic dissipation can also strongly impact the atmospheric behavior, depending on atmospheric composition and the planetary magnetic field strength. We conclude that hot super-Earths could exhibit interesting signatures in reflection (and possibly in emission) which would trace a combination of their ground, atmospheric and magnetic properties.

In contrast to the Earth, where frictional heating is typically negligible, we show that drag mechanisms could act as an important heat source in the strongly-forced atmospheres of some exoplanets, with the potential to alter the circulation. We modify the standard formalism of the atmospheric energy cycle to explicitly track the loss of kinetic energy and the associated frictional (re)heating, for application to exoplanets such as the asymmetrically heated "hot Jupiters" and gas giants on highly eccentric orbits. We establish that an understanding of the dominant drag mechanisms and their dependence on local atmospheric conditions is critical for accurate modeling, not just in their ability to limit wind speeds, but also because they could possibly change the energetics of the circulation enough to alter the nature of the flow. We discuss possible sources of drag and estimate the strength necessary to significantly influence the atmospheric energetics. As we show, the frictional heating depends on the magnitude of kinetic energy dissipation as well as its spatial variation, so that the more localized a drag mechanism is, the weaker it can be and still affect the circulation. We also use the derived formalism to estimate the rate of numerical loss of kinetic energy in a few previously published hot Jupiter models with and without magnetic drag and find it to be surprisingly large, at 5-10% of the incident stellar irradiation.

The complexity of atmospheric modelling and its inherent non-linearity, together with the limited amount of data of exoplanets available, motivate model intercomparisons and benchmark tests. In the geophysical community, the Held-Suarez test is a standard benchmark for comparing dynamical core simulations of the Earth's atmosphere with different solvers, based on statistically-averaged flow quantities. In the present study, we perform analogues of the Held-Suarez test for tidally-locked exoplanets with the GFDL-Princeton Flexible Modeling System (FMS) by subjecting both the spectral and finite difference dynamical cores to a suite of tests, including the standard benchmark for Earth, a hypothetical tidally-locked Earth, a "shallow" hot Jupiter model and a "deep" model of HD 209458b. We find qualitative and quantitative agreement between the solvers for the Earth, tidally-locked Earth and shallow hot Jupiter benchmarks, but the agreement is less than satisfactory for the deep model of HD 209458b. Further investigation reveals that closer agreement may be attained by arbitrarily adjusting the values of the horizontal dissipation parameters in the two solvers, but it remains the case that the magnitude of the horizontal dissipation is not easily specified from first principles. Irrespective of radiative transfer or chemical composition considerations, our study points to limitations in our ability to accurately model hot Jupiter atmospheres with meteorological solvers at the level of ten percent for the temperature field and several tens of percent for the velocity field. Direct wind measurements should thus be particularly constraining for the models. Our suite of benchmark tests also provides a reference point for researchers wishing to adapt their codes to study the atmospheric circulation regimes of tidally-locked Earths/Neptunes/Jupiters.

Hot Jupiter atmospheres exhibit fast, weakly-ionized winds. The interaction of these winds with the planetary magnetic field generates drag on the winds and leads to ohmic dissipation of the induced electric currents. We study the magnitude of ohmic dissipation in representative, three-dimensional atmospheric circulation models of the hot Jupiter HD 209458b. We find that ohmic dissipation can reach or exceed 1% of the stellar insolation power in the deepest atmospheric layers, in models with and without dragged winds. Such power, dissipated in the deep atmosphere, appears sufficient to slow down planetary contraction and explain the typically inflated radii of hot Jupiters. This atmospheric scenario does not require a top insulating layer or radial currents that penetrate deep in the planetary interior. Circulation in the deepest atmospheric layers may actually be driven by spatially non-uniform ohmic dissipation. A consistent treatment of magnetic drag and ohmic dissipation is required to further elucidate the consequences of magnetic effects for the atmospheres and the contracting interiors of hot Jupiters.

In the outer regions of the habitable zone, the risk of transitioning into a globally frozen "snowball" state poses a threat to the habitability of planets with the capacity to host water-based life. We use a one-dimensional energy balance climate model (EBM) to examine how obliquity, spin rate, orbital eccentricity, and ocean coverage might influence the onset of such a snowball state. For an exoplanet, these parameters may be strikingly different from the values observed for Earth. Since, for constant semimajor axis, the annual mean stellar irradiation scales with (1-e^2)^(-1/2), one might expect the greatest habitable semimajor axis (for fixed atmospheric composition) to scale as (1-e^2)^(-1/4). We find that this standard ansatz provides a reasonable lower bound on the outer boundary of the habitable zone, but the influence of obliquity and ocean fraction can be profound in the context of planets on eccentric orbits. For planets with eccentricity 0.5, our EBM suggests that the greatest habitable semimajor axis can vary by more than 0.8 AU (78%!) depending on obliquity, with higher obliquity worlds generally more stable against snowball transitions. One might also expect that the long winter at an eccentric planet's apoastron would render it more susceptible to global freezing. Our models suggest that this is not a significant risk for Earth-like planets around Sun-like stars since such planets are buffered by the thermal inertia provided by oceans covering at least 10% of their surface. Since planets on eccentric orbits spend much of their year particularly far from the star, such worlds might turn out to be especially good targets for direct observations with missions such as TPF-Darwin. Nevertheless, the extreme temperature variations achieved on highly eccentric exo-Earths raise questions about the adaptability of life to marginally or transiently habitable conditions.

White dwarfs inspiraling into black holes of mass $\MBH\simgt 10^5M_\odot$ are detectable sources of gravitational waves in the LISA band. In many of these events, the white dwarf begins to lose mass during the main observational phase of the inspiral. The mass loss starts gently and can last for thousands of orbits. The white dwarf matter overflows the Roche lobe through the $L_1$ point at each pericenter passage and the mass loss repeats periodically. The process occurs very close to the black hole and the released gas can accrete, creating a bright source of radiation with luminosity close to the Eddington limit, $L\sim 10^{43}$~erg~s$^{-1}$. This class of inspirals offers a promising scenario for dual detections of gravitational waves and electromagnetic radiation.

We report the detection of 24 micron variations from the planet-hosting upsilon Andromedae system consistent with the orbital periodicity of the system's innermost planet, upsilon And b. We find a peak-to-valley phase curve amplitude of 0.00130 times the mean system flux. Using a simple model with two hemispheres of constant surface brightness and assuming a planetary radius of 1.3 Jupiter radii gives a planetary temperature contrast of >900 K and an orbital inclination of >28 degrees. We further report the largest phase offset yet observed for an extrasolar planet: the flux maximum occurs ~80 degrees before phase 0.5. Such a large phase offset is difficult to reconcile with most current atmospheric circulation models. We improve on earlier observations of this system in several important ways: (1) observations of a flux calibrator star demonstrate the MIPS detector is stable to 10^-4 on long timescales, (2) we note that the background light varies systematically due to spacecraft operations, precluding use of this background as a flux calibrator (stellar flux measured above the background is not similarly affected), and (3) we calibrate for flux variability correlated with motion of the star on the MIPS detector. A reanalysis of our earlier observations of this system is consistent with our new result.

The coalescence of a supermassive black hole binary (SMBHB) is thought to be accompanied by an electromagnetic (EM) afterglow, produced by the viscous infall of the surrounding circumbinary gas disk after the merger. It has been proposed that once the merger has been detected in gravitational waves (GWs) by LISA, follow-up EM searches for this afterglow can help identify the EM counterpart of the LISA source. Here we study whether the afterglows may be sufficiently bright and numerous to be detectable in EM surveys alone. The viscous afterglow, which lasts for years to decades for SMBHBs in LISA's sensitivity window, is characterized by rapid increases in both the bolometric luminosity and in the spectral hardness of the source. If quasar activity is triggered by the same major galaxy mergers that produce SMBHBs, then the afterglow could be interpreted as a signature of the birth of a quasar. Using an idealized model for the post-merger viscous spreading of the circumbinary disk and the resulting light curve, and using the observed luminosity function of quasars as a proxy for the SMBHB merger rate, we delineate the survey requirements for identifying such birthing quasars. If circumbinary disks have a high disk surface density and viscosity, an all-sky soft X-ray survey with a sensitivity of ~<3x10^-14 erg s^-1 cm^-2 and a time resolution of ~months could identify dozens of birthing quasars with sustained brightening rates of >10%/yr. If >1% of the X-ray emission is reprocessed into optical frequencies, birthing quasars could also be identified in optical transient surveys such as the LSST. Distinguishing a birthing quasar from other variable sources may be facilitated by the monotonic hardening of its spectrum, but will likely remain challenging. This reinforces the notion that joint EM-plus-GW observations offer the best prospects for identifying the EM signatures of SMBHB mergers.

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 present a three dimensional hot Jupiter model, extending from 200 bar to 1 mbar, using the Intermediate General Circulation Model from the University of Reading. Our horizontal spectral resolution is T31 (equivalent to a grid of 48x96), with 33 logarithmically spaced vertical levels. A simplified (Newtonian) scheme is employed for the radiative forcing. We adopt a physical set up nearly identical to the model of HD 209458b by Cooper & Showman (2005,2006) to facilitate a direct model inter-comparison. Our results are broadly consistent with theirs but significant differences also emerge. The atmospheric flow is characterized by a super-rotating equatorial jet, transonic wind speeds, and eastward advection of heat away from the dayside. We identify a dynamically-induced temperature inversion ("stratosphere") on the planetary dayside and find that temperatures at the planetary limb differ systematically from local radiative equilibrium values, a potential source of bias for transit spectroscopic interpretations. While our model atmosphere is quasi-identical to that of Cooper & Showman (2005,2006) and we solve the same meteorological equations, we use different algorithmic methods, spectral-implicit vs. grid-explicit, which are known to yield fully consistent results in the Earth modeling context. The model discrepancies identified here indicate that one or both numerical methods do not faithfully capture all of the atmospheric dynamics at work in the hot Jupiter context. We highlight the emergence of a shock-like feature in our model, much like that reported recently by Showman et al. (2009), and suggest that improved representations of energy conservation may be needed in hot Jupiter atmospheric models, as emphasized by Goodman (2009).

Hot Jupiters, with atmospheric temperatures T ~ 1000 K, have residual thermal ionization levels sufficient for the interaction of the ions with the planetary magnetic field to result in a sizable magnetic drag on the (neutral) atmospheric winds. We evaluate the magnitude of magnetic drag in a representative three-dimensional atmospheric model of the hot Jupiter HD 209458b and find that it is a plausible mechanism to limit wind speeds in this class of atmospheres. Magnetic drag has a strong geometrical dependence, both meridionally and from the day to the night side (in the upper atmosphere), which could have interesting consequences for the atmospheric flow pattern. By extension, close-in eccentric planets with transiently heated atmospheres will experience time-variable levels of magnetic drag. A robust treatment of magnetic drag in circulation models for hot atmospheres may require iterated solutions to the magnetic induction and Saha equations as the hydrodynamical flow is evolved.

The Laser Interferometer Space Antenna (LISA) will detect gravitational wave signals from coalescing pairs of massive black holes in the total mass range (10^5 - 10^7)/Msol out to cosmological distances. Identifying and monitoring the electromagnetic counterparts of these events would enable cosmological studies and offer new probes of gas physics around well-characterized massive black holes. Milosavljevic & Phinney (2005) proposed that a circumbinary disk around a binary of mass ~10^6 Msol will emit an accretion-powered X-ray afterglow approximately one decade after the gravitational wave event. We revisit this scenario by using Green's function solutions to calculate the temporal viscous evolution and the corresponding electromagnetic signature of the circumbinary disk. Our calculations suggest that an electromagnetic counterpart may become observable as a rapidly brightening source soon after the merger, i.e. several years earlier than previously thought. The afterglow can reach super-Eddington luminosities without violating the local Eddington flux limit. It is emitted in the soft X-ray by the innermost circumbinary disk, but it may be partially reprocessed at optical and infrared frequencies. We also find that the spreading disk becomes increasingly geometrically thick close to the central object as it evolves, indicating that the innermost flow could become advective and radiatively inefficient, and generate a powerful outflow. We conclude that the mergers of massive black holes detected by LISA offer unique opportunities for monitoring on humanly tractable timescales the viscous evolution of accretion flows and the emergence of outflows around massive black holes with precisely known masses, spins and orientations.

Radiative transfer in planetary atmospheres is usually treated in the static limit, i.e., neglecting atmospheric motions. We argue that hot Jupiter atmospheres, with possibly fast (sonic) wind speeds, may require a more strongly coupled treatment, formally in the regime of radiation-hydrodynamics. To lowest order in v/c, relativistic Doppler shifts distort line profiles along optical paths with finite wind velocity gradients. This leads to flow-dependent deviations in the effective emission and absorption properties of the atmospheric medium. Evaluating the overall impact of these distortions on the radiative structure of a dynamic atmosphere is non-trivial. We present transmissivity and systematic equivalent width excess calculations which suggest possibly important consequences for radiation transport in hot Jupiter atmospheres. If winds are fast and bulk Doppler shifts are indeed important for the global radiative balance, accurate modeling and reliable data interpretation for hot Jupiter atmospheres may prove challenging: it would involve anisotropic and dynamic radiative transfer in a coupled radiation-hydrodynamical flow. On the bright side, it would also imply that the emergent properties of hot Jupiter atmospheres are more direct tracers of their atmospheric flows than is the case for Solar System planets. Radiation-hydrodynamics may also influence radiative transfer in other classes of hot exoplanetary atmospheres with fast winds.

We survey the basic principles of atmospheric dynamics relevant to explaining existing and future observations of exoplanets, both gas giant and terrestrial. Given the paucity of data on exoplanet atmospheres, our approach is to emphasize fundamental principles and insights gained from Solar-System studies that are likely to be generalizable to exoplanets. We begin by presenting the hierarchy of basic equations used in atmospheric dynamics, including the Navier-Stokes, primitive, shallow-water, and two-dimensional nondivergent models. We then survey key concepts in atmospheric dynamics, including the importance of planetary rotation, the concept of balance, and scaling arguments to show how turbulent interactions generally produce large-scale east-west banding on rotating planets. We next turn to issues specific to giant planets, including their expected interior and atmospheric thermal structures, the implications for their wind patterns, and mechanisms to pump their east-west jets. Hot Jupiter atmospheric dynamics are given particular attention, as these close-in planets have been the subject of most of the concrete developments in the study of exoplanetary atmospheres. We then turn to the basic elements of circulation on terrestrial planets as inferred from Solar-System studies, including Hadley cells, jet streams, processes that govern the large-scale horizontal temperature contrasts, and climate, and we discuss how these insights may apply to terrestrial exoplanets. Although exoplanets surely possess a greater diversity of circulation regimes than seen on the planets in our Solar System, our guiding philosophy is that the multi-decade study of Solar-System planets reviewed here provides a foundation upon which our understanding of more exotic exoplanetary meteorology must build.

Sharp density features in protoplanetary discs, for instance at the edge of a magnetically dead zone, have recently been proposed as effective barriers to slow down or even stop the problematically fast migration of planetary cores into their central star. Density features on a radial scale approaching the disc vertical scale height might not exist, however, since they could be Rayleigh (or more generally Solberg-Hoiland) unstable. Stability must be checked explicitly in one-dimensional viscous accretion disc models because these instabilities are artificially eliminated in the process of reducing the full set of axisymmetric equations. The disc thermodynamics, via the entropy stratification, and its vertical structure also influence stability when sharp density features are present. We propose the concept of Rayleigh adjustment for viscous disc models: any density feature that violates Rayleigh stability (or its generalization) should be diffused radially by hydrodynamical turbulence on a dynamical time-scale, approaching marginal stability in a quasi-continuous manner.

The dynamics of planetesimals and planetary cores may be strongly influenced by density perturbations driven by magneto-rotational turbulence in their natal protoplanetary gas disks. Using the local shearing box approximation, we perform numerical simulations of planetesimals moving as massless particles in a turbulent, magnetized, unstratified gas disk. Our fiducial disk model shows turbulent accretion characterized by a Shakura-Sunyaev viscosity parameter of $\alpha \sim 10^{-2}$, with root-mean-square density perturbations of $\sim$10%. We measure the statistical evolution of particle orbital properties in our simulations including mean radius, eccentricity, and velocity dispersion. We confirm random walk growth in time of all three properties, the first time that this has been done with direct orbital integration in a local model. We find that the growth rate increases with the box size used at least up to boxes of eight scale heights in horizontal size. However, even our largest boxes show velocity dispersions sufficiently low that collisional destruction of planetesimals should be unimportant in the inner disk throughout its lifetime. Our direct integrations agree with earlier torque measurements showing that type I migration dominates over diffusive migration by stochastic torques for most objects in the planetary core and terrestrial planet mass range. Diffusive migration remains important for objects in the mass range of kilometer-sized planetesimals. Discrepancies in the derived magnitude of turbulence between local and global simulations of magneto-rotationally unstable disks remains an open issue, with important consequences for planet formation scenarios.

Using a one-dimensional height integrated model, we calculate the evolution of an unequal mass binary black hole with a coplanar gas disk that contains a gap due to the presence of the secondary black hole. Viscous evolution of the outer circumbinary disk initially hardens the binary, while the inner disk drains onto the primary (central) black hole. As long as the inner disk remains cool and thin at low $\dot{M}_{\rm ext}$ (rather than becoming hot and geometrically thick), the mass of the inner disk reaches an asymptotic mass typically $\sim 10^{-3}-10^{-4}\Msun$. Once the semimajor axis shrinks below a critical value, angular momentum losses from gravitational waves dominate over viscous transport in hardening the binary. The inner disk then no longer responds viscously to the inspiraling black holes. Instead, tidal interactions with the secondary rapidly drive the inner disk into the primary. Tidal and viscous dissipation in the inner disk lead to a late time brightening in luminosity $L\propto t_{\rm minus}^{5/4}$, where $t_{\rm minus}$ is the time prior to the final merger. This late time brightening peaks $\sim 1$ day prior to the final merger at $\sim 0.1 L_{\rm Edd}$. This behavior is relatively robust because of self regulation in the coupled viscous-gravitational evolution of such binary systems. It constitutes a unique electromagnetic signature of a binary supermassive black hole merger and may allow the host galaxy to be identified if used in conjunction with the Laser Interferometric Space Antenna (LISA) localization.

Remote observing of exoplanetary atmospheres is now possible, offering us access to circulation regimes unlike any of the familiar Solar System cases. Atmospheric circulation models are being developed to study these new regimes but model validations and intercomparisons are needed to establish their consistency and accuracy. To this end, we present a simple Earth-like validation of the pseudo-spectral solver of meteorological equations called IGCM (Intermediate General Circulation Model), based on Newtonian relaxation to a prescribed latitudinal profile of equilibrium temperatures. We then describe a straightforward and idealized model extension to the atmospheric flow on a hot Jupiter with the same IGCM solver. This shallow, three-dimensional hot Jupiter model is based on Newtonian relaxation to a permanent day-night pattern of equilibrium temperatures and the absence of surface drag. The baroclinic regime of the Earth's lower atmosphere is contrasted with the more barotropic regime of the simulated hot Jupiter flow. For plausible conditions at the 0.1-1 bar pressure level on HD 209458b, the simulated flow is characterized by unsteadiness, subsonic wind speeds, a zonally-perturbed superrotating equatorial jet and large scale polar vortices. Violation of the Rayleigh-Kuo inflexion point criterion on the flanks of the accelerating equatorial jet indicates that barotropic (horizontal shear) instabilities may be important dynamical features of the simulated flow. Similarities and differences with previously published simulated hot Jupiter flows are briefly noted.

Supermassive black hole binaries (SMBHBs) produced in galaxy mergers are thought to complete their coalescence, below separations of r_GW=10^{-3} (M_BH/10^8 M_sun)^{3/4} pc, as their orbit decays due to the emission of gravitational waves (GWs). It may be possible to identify such GW-driven inspirals statistically in an electromagnetic (EM) survey for variable sources. A GW-driven binary spends a characteristic time T_GW at each orbital separation r_orb < r_GW that scales with the corresponding orbital time t_orb as T_GW = (const) t_orb^{8/3}. If the coalescing binary produces variations in the EM emission on this timescale, then it could be identified as a variable source with a characteristic period t_var = t_orb. The incidence rate of sources with similar inferred BH masses, showing near-periodic variability on the time-scale t_var, would then be proportional to t_var^{8/3}. Luminosity variations corresponding to a fraction f_Edd<0.01 of the Eddington luminosity would have been missed in current surveys. However, if the binary inspirals are associated with quasars, we show that a dedicated survey could detect the population of SMBHBs with a range of periods around tens of weeks. The discovery of a population of periodic sources whose abundance obeys N_var = (const) t_var^{8/3} would confirm (i) that the orbital decay is indeed driven by GWs, and (ii) that circumbinary gas is present at small orbital radii and is being perturbed by the BHs. Deviations from the t_var^{8/3} power-law could constrain the structure of the circumbinary gas disk and viscosity-driven orbital decay. We discuss constraints from existing data, and quantify the sensitivity and sky coverage that could yield a detection in future surveys.

Supermassive black hole binaries (SMBHBs) in galactic nuclei are thought to be a common by-product of major galaxy mergers. We use simple disk models for the circumbinary gas and for the binary-disk interaction to follow the orbital decay of SMBHBs with a range of total masses (M) and mass ratios (q), through physically distinct regions of the disk, until gravitational waves (GWs) take over their evolution. Prior to the GW-driven phase, the viscous decay is in the stalled "secondary-dominated" regime. SMBHBs spend a non-negligible fraction of 10^7 years at orbital periods t_var between a day and a year. A dedicated optical or X-ray survey could identify coalescing SMBHBs statistically, as a population of periodically variable quasars, whose abundance N_var is proportional to t_var^alpha, in a range of periods t_var around tens of weeks. SMBHBs with M < 10^7 M_sun, with 0.5 < alpha < 1.5, would probe the physics of viscous orbital decay, whereas the detection of a population of higher-mass binaries, with alpha=8/3, would confirm that their decay is driven by GWs. The lowest mass SMBHBs (M < 10^{5-6} M_sun) enter the GW-driven regime at short orbital periods, in the frequency band of the Laser Interferometric Space Antenna (LISA). While viscous processes are strongly sub-dominant in the last few years of coalescence, they could reduce the amplitude of any unresolved background of near-stationary LISA sources. We discuss constraints on the SMBHB population available from existing data, and the sensitivity and sky coverage requirements for a detection in future surveys. SMBHBs may also be identified from velocity shifts in their spectra; we discuss the expected abundance of SMBHBs as a function of their orbital velocity.

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

Recent imaging campaigns indicate the likely existence of massive planets (~ 1-10 MJ) on ~1000 year orbits about a few percent of stars. Such objects are not easily explained in most current planet formation models. In this Letter we use ensembles of 100 N-body simulations to evaluate the potential for planet scattering during relaxation of dynamically active systems to explain the population of giant planets with projected separations up to a few 100 AU. We find that such a mechanism could indeed be at play, and that statistical samples of long period planets could place interesting constraints on early stage planet formation scenarios. Results from direct imaging and microlensing surveys are complementary probes of this dynamical relaxation process.

Supermassive black hole binaries (SMBHBs) with masses in the range 10^4-10^7 M_sun/(1+z), produced in galaxy mergers, are thought to complete their coalescence due to the emission of gravitational waves (GWs). The anticipated detection of the GWs by the LISA will constitute a milestone for fundamental physics and astrophysics. While the GW signatures themselves will provide a treasure trove of information, if the source can be securely identified in electromagnetic (EM) bands, this would open up entirely new scientific opportunities, to probe fundamental physics, astrophysics, and cosmology. We discuss several ideas, involving wide-field telescopes, that may be useful in locating electromagnetic counterparts to SMBHBs detected by LISA. In particular, the binary may produce a variable electromagnetic flux, such as a roughly periodic signal due to the orbital motion prior to coalescence, or a prompt transient signal caused by shocks in the circumbinary disk when the SMBHB recoils and "shakes" the disk. We discuss whether these time-variable EM signatures may be detectable, and how they can help in identifying a unique counterpart within the localization errors provided by LISA. We also discuss a possibility of identifying a population of coalescing SMBHBs statistically, in a deep optical survey for periodically variable sources, before LISA detects the GWs directly. The discovery of such sources would confirm that gas is present in the vicinity and is being perturbed by the SMBHB - serving as a proof of concept for eventually finding actual LISA counterparts.

Hyper-accreting discs occur in compact-object mergers and collapsars, and may power gamma-ray bursts (GRBs). We calculate the microscopic viscosity and resistivity of plasma in these discs, and discuss the implications for their global structure and evolution. In the neutrino-cooled innermost regions, the viscosity is provided mainly by mildly degenerate electrons, while the resistivity is modified from the Spitzer value due to the effects of both relativity and degeneracy. The plasma behaves as an almost ideal MHD fluid. Among the non-ideal MHD effects the Hall term is relatively the most important, while the magnetic Prandtl number, Pr (the ratio of viscosity to resistivity), is typically larger than unity: 10 < Pr < 6000. The outer radiatively inefficient regions also display high Pr. Numerical simulations of the magneto-rotational instability indicate that the saturation level and angular momentum transport efficiency may be greatly enhanced at high Pr. If this behaviour persists in the presence of a strong Hall effect we would expect that hyper-accreting discs should be strongly magnetised and highly variable. The expulsion of magnetic field that cannot be dissipated at small scales may also favour a magnetic outflow. We note the similaries between the Prandtl number in hyper-accreting discs and X-ray binary discs, which suggests that a comparison between late-time activity in GRBs and X-ray binary accretion states may be fruitful. Our results imply that the behavior of high Prandtl number MHD flows needs to be considered in studies of hyper-accreting discs.

Extrasolar terrestrial planets with the potential to host life might have large obliquities or be subject to strong obliquity variations. We revisit the habitability of oblique planets with an energy balance climate model (EBM) allowing for dynamical transitions to ice-covered snowball states as a result of ice-albedo feedback. Despite the great simplicity of our EBM, it captures reasonably well the seasonal cycle of global energetic fluxes at Earth's surface. It also performs satisfactorily against a full-physics climate model of a highly oblique Earth-like planet, in an unusual regime of circulation dominated by heat transport from the poles to the equator. Climates on oblique terrestrial planets can violate global radiative balance through much of their seasonal cycle, which limits the usefulness of simple radiative equilibrium arguments. High obliquity planets have severe climates, with large amplitude seasonal variations, but they are not necessarily more prone to global snowball transitions than low obliquity planets. We find that terrestrial planets with massive CO2 atmospheres, typically expected in the outer regions of habitable zones, can also be subject to such dynamical snowball transitions. Some of the snowball climates investigated for CO2-rich atmospheres experience partial atmospheric collapse. Since long-term CO2 atmospheric build-up acts as a climatic thermostat for habitable planets, partial CO2 collapse could limit the habitability of such planets. A terrestrial planet's habitability may thus depend sensitively on its short-term climatic stability.

Turbulence is ubiquitous in Solar System planetary atmospheres. In hot Jupiter atmospheres, the combination of moderately slow rotation and thick pressure scale height may result in dynamical weather structures with unusually large, planetary-size scales. Using equivalent-barotropic, turbulent circulation models, we illustrate how such structures can generate a variety of features in the thermal phase curves of hot Jupiters, including phase shifts and deviations from periodicity. Such features may have been spotted in the recent infrared phase curve of HD 189733b. Despite inherent difficulties with the interpretation of disk-integrated quantities, phase curves promise to offer unique constraints on the nature of the circulation regime present on hot Jupiters.

The notion that microparsec-scale black holes can be used to probe gigaparsec-scale physics may seem counterintuitive, at first. Yet, the gravitational observatory LISA will detect cosmologically-distant coalescing pairs of massive black holes, accurately measure their luminosity distance and help identify an electromagnetic counterpart or a host galaxy. A wide variety of new black hole studies and a gravitational version of Hubble's diagram become possible if host galaxies are successfully identified. Furthermore, if dark energy is a manifestation of large-scale modified gravity, deviations from general relativistic expectations could become apparent in a gravitational signal propagated over cosmological scales, especially when compared to the electromagnetic signal from a same source. Finally, since inspirals of white dwarfs into massive black holes at cosmological distances may permit pre-merger localizations, we suggest that careful monitoring of these events and any associated electromagnetic counterpart could lead to high-precision cosmological measurements with LISA.

According to the standard liquid-water definition, the Earth is only partially habitable. We reconsider planetary habitability in the framework of energy-balance models, the simplest seasonal models in physical climatology, to assess the spatial and temporal habitability of Earth-like planets. We quantify the degree of climatic habitability of our models with several metrics of fractional habitability. Previous evaluations of habitable zones may have omitted important climatic conditions by focusing on close Solar System analogies. For example, we find that model pseudo-Earths with different rotation rates or different land-ocean fractions have fractional habitabilities that differ significantly from that of the Earth itself. Furthermore, the stability of a planet's climate against albedo-feedback snowball events strongly impacts its habitability. Therefore, issues of climate dynamics may be central in assessing the habitability of discovered terrestrial exoplanets, especially if astronomical forcing conditions are different from the moderate Solar System cases.

Electromagnetic (EM) counterparts to supermassive black hole binary mergers observed by LISA can be localized to within the field of view of astronomical instruments ~10 deg^2 hours to weeks prior to coalescence. The temporal coincidence of any prompt EM counterpart with a gravitationally-timed merger may offer the best chance of identifying a unique host galaxy. We discuss the challenges posed by searches for prompt EM counterparts and propose novel observational strategies to address them. In particular, we discuss the size and shape evolution of the LISA localization error ellipses on the sky, and quantify the requirements for dedicated EM surveys of the area prior to coalescence. A triggered EM counterpart search campaign will require monitoring a several-square degree area. It could aim for variability at the 24-27 mag level in optical bands, for example, which corresponds to 1-10% of the Eddington luminosity of the prime LISA sources of 10^6-10^7 Msun BHs at z=1-2, on time-scales of minutes to hours, the orbital time-scale of the binary in the last 2-4 weeks. A cross-correlation of the period of any variable EM signal with the quasi-periodic gravitational waveform over 10-1000 cycles may aid the detection. Alternatively, EM searches can detect a transient signal accompanying the coalescence. We highlight the measurement of differences in the arrival times of photons and gravitons from the same cosmological source as a valuable independent test of the massive character of gravity, and of possible violations of Lorentz invariance in the gravity sector.

We present results from a set of over 300 pseudospectral simulations of atmospheric circulation on extrasolar giant planets with circular orbits. The simulations are of high enough resolution (up to 341 total and sectoral modes) to resolve small-scale eddies and waves, required for reasonable physical accuracy. In this work, we focus on the global circulation pattern that emerges in a shallow, ``equivalent-barotropic'', turbulent atmosphere on both tidally synchronized and unsynchronized planets. A full exploration of the large physical and numerical parameter-space is performed to identify robust features of the circulation. For some validation, the model is first applied to Solar System giant planets. For extrasolar giant planets with physical parameters similar to HD209458b--a presumably synchronized extrasolar giant planet, representative in many dynamical respects--the circulation is characterized by the following features: 1) a coherent polar vortex that revolves around the pole in each hemisphere; 2) a low number--typically two or three--of slowly-varying, broad zonal (east-west) jets that form when the maximum jet speed is comparable to, or somewhat stronger than, those observed on the planets in the Solar System; and, 3) motion-associated temperature field, whose detectability and variability depend on the strength of the net heating rate and the global root mean square wind speed in the atmosphere. In many ways, the global circulation is Earth-like, rather than Jupiter-like. However, if extrasolar giant planets rotate faster and are not close-in (therefore not synchronized), their circulations become more Jupiter-like, for Jupiter-like rotation rates.

Hot Jupiters are new laboratories for the physics of giant planet atmospheres. Subject to unusual forcing conditions, the circulation regime on these planets may be unlike anything known in the Solar System. Characterizing the atmospheric circulation of hot Jupiters is necessary for reliable interpretation of the multifaceted data currently being collected on these planets. We discuss several fundamental concepts of atmospheric dynamics that are likely central to obtaining a solid understanding of these fascinating atmospheres. A particular effort is made to compare the various modeling approaches employed so far to address this challenging problem.

A gravitational observatory such as LISA will detect coalescing pairs of massive black holes, accurately measure their luminosity distance and help identify a host galaxy or an electromagnetic counterpart. If dark energy is a manifestation of modified gravity on large scales, gravitational waves from cosmologically-distant spacetime sirens are direct probes of this new physics. For example, a gravitational Hubble diagram based on black hole pair luminosity distances and host galaxy redshifts could reveal a large distance extra-dimensional leakage of gravity. Various additional signatures may be expected in a gravitational signal propagated over cosmological scales.

Solar active regions, produced by the emergence of tubes of strong magnetic field in the photosphere, are restricted to within 35 degrees of the solar equator. The nature of the dynamo processes that create and renew these fields, and are therefore responsible for solar magnetic phenomena, are not well understood. We analyze the magneto-rotational stability of the solar tachocline for general field geometry. This thin region of strong radial and latitudinal differential rotation, between the radiative and convective zones, is unstable at latitudes above 37 degrees, yet is stable closer to the equator. We propose that small-scale magneto-rotational turbulence prevents coherent magnetic dynamo action in the tachocline except in the vicinity of the equator, thus explaining the latitudinal restriction of active regions. Tying the magnetic dynamo to the tachocline elucidates the physical conditions and processes relevant to solar magnetism.

Migration of protoplanets in their gaseous host disks may be largely responsible for the observed orbital distribution of extrasolar planets. Recent simulations have shown that the magnetorotational turbulence thought to drive accretion in protoplanetary disks can affect migration by turning it into an orbital random walk. However, these simulations neglected the disk's ionization structure. Low ionization fraction near the midplane of the disk can decouple the magnetic field from the gas, forming a dead zone with reduced or no turbulence. Here, to understand the effect of dead zones on protoplanetary migration, we perform numerical simulations of a small region of a stratified disk with magnetorotational turbulence confined to thin active layers above and below the midplane. Turbulence in the active layers exerts decreased, but still measurable, gravitational torques on a protoplanet located at the disk midplane. We find a decrease of two orders of magnitude in the diffusion coefficient for dead zones with dead-to-active surface density ratios approaching realistic values in protoplanetary disks. This torque arises primarily from density fluctuations within a distance of one scale height of the protoplanet. Turbulent torques have correlation times of only $\sim 0.3$ orbital periods and apparently time-stationary distributions. These properties are encouraging signs that stochastic methods can be used to determine the orbital evolution of populations of protoplanets under turbulent migration. Our results indicate that dead zones may be dynamically distinct regions for protoplanetary migration.

We study the saturation near threshold of the axisymmetric magnetorotational instability (MRI) of a viscous, resistive, incompressible fluid in a thin-gap Taylor-Couette configuration. A vertical magnetic field, Keplerian shear and no-slip, conducting radial boundary conditions are adopted. The weakly non-linear theory leads to a real Ginzburg-Landau equation for the disturbance amplitude, like in our previous idealized analysis. For small magnetic Prandtl number (P \ll 1), the saturation amplitude scales as P^{2/3} while the magnitude of angular momentum transport scales as P^{4/3}. The difference with the previous scalings (~P^{1/2} and P respectively) is attributed to the emergence of radial boundary layers. Away from those, steady-state non-linear saturation is achieved through a modest reduction in the destabilizing shear. These results will be useful to understand MRI laboratory experiments and associated numerical simulations.

The continuous improvement in localization errors (sky position and distance) in real time as LISA observes the gradual inspiral of a supermassive black hole (SMBH) binary can be of great help in identifying any prompt electromagnetic counterpart associated with the merger. We develop a new method, based on a Fourier decomposition of the time-dependent, LISA-modulated gravitational-wave signal, to study this intricate problem. The method is faster than standard Monte Carlo simulations by orders of magnitude. By surveying the parameter space of potential LISA sources, we find that counterparts to SMBH binary mergers with total mass M~10^5-10^7 M_Sun and redshifts z<~3 can be localized to within the field of view of astronomical instruments (~deg^2) typically hours to weeks prior to coalescence. This will allow targeted searches for variable electromagnetic counterparts as the merger proceeds, as well as monitoring of the most energetic coalescence phase. A rich set of astrophysical and cosmological applications would emerge from the identification of electromagnetic counterparts to these gravitational-wave standard sirens.

Physical conditions in the atmospheres of tidally-locked, slowly-rotating hot Jupiters correspond to dynamical circulation regimes with Rhines scales and Rossby deformation radii comparable to the planetary radii. Consequently, the large spatial scales of moving atmospheric structures could generate significant photospheric variability. Here, we estimate the level of thermal infrared variability expected in successive secondary eclipse depths, according to hot Jupiter turbulent ``shallow-layer'' models. The variability, at the few percent level or more in models with strong enough winds, is within the reach of Spitzer measurements. Eclipse depth variability is thus a valuable tool to constrain the circulation regime and global wind speeds in hot Jupiter atmospheres.

Recent Spitzer infrared measurements of hot Jupiter eclipses suggest that eclipse mapping techniques could be used to spatially resolve the day-side photospheric emission of these planets using partial occultations. As a first step in this direction, we simulate ingress/egress lightcurves for the three brightest known eclipsing hot Jupiters and evaluate the degree to which parameterized photospheric emission models can be distinguished from each other with repeated, noisy eclipse measurements. We find that the photometric accuracy of Spitzer is insufficient to use this tool effectively. On the other hand, the level of photospheric details that could be probed with a few JWST eclipse measurements could greatly inform hot Jupiter atmospheric modeling efforts. A JWST program focused on non-parametric eclipse map inversions for hot Jupiters should be actively considered.

We show by means of a perturbative weakly nonlinear analysis that the axisymmetric magnetorotational instability (MRI) of a viscous, resistive, incompressible rotating shear flow in a thin channel gives rise to a real Ginzburg-Landau equation for the disturbance amplitude. For small magnetic Prandtl number (${\cal P}_{\rm m}$), the saturation amplitude is $\propto \sqrt{{\cal P}_{\rm m}}$ and the resulting momentum transport scales as ${\cal R}^{-1}$, where $\cal R$ is the {\em hydrodynamic} Reynolds number. Simplifying assumptions, such as linear shear base flow, mathematically expedient boundary conditions and continuous spectrum of the vertical linear modes, are used to facilitate this analysis. The asymptotic results are shown to comply with numerical calculations using a spectral code. They suggest that the transport due to the nonlinearly developed MRI may be very small in experimental setups with ${\cal P}_{\rm m} \ll 1$.

The star upsilon Andromeda is orbited by three known planets, the innermost of which has an orbital period of 4.617 days and a mass at least 0.69 that of Jupiter. This planet is close enough to its host star that the radiation it absorbs overwhelms its internal heat losses. Here we present the 24 micron light curve of this system, obtained with the Spitzer Space Telescope. It shows a clear variation in phase with the orbital motion of the innermost planet. This is the first demonstration that such planets possess distinct hot substellar (day) and cold antistellar (night) faces.

Angular momentum transport must have occurred in the Sun's radiative zone to explain its current solid body rotation. We survey the stability of the early Sun's radiative zone with respect to diffusive rotational instabilities, for a variety of plausible past configurations. We find that the (faster rotating) early Sun was prone to rotational instabilities even if only weak levels of radial differential rotation were present, while the current Sun is not. Stability domains are determined by approximate balance between dynamical and diffusive timescales, allowing generalizations to other stellar contexts. Depending on the strength and geometry of the weak magnetic field present, the fastest growing unstable mode can be hydrodynamic or magneto-hydrodynamic (MHD) in nature. Our results suggest that diffusive MHD modes may be more efficient at transporting angular momentum than their hydrodynamic (``Goldreich-Schubert-Fricke'') counterparts because the minimum spatial scale required for magnetic tension to be destabilizing limits the otherwise very small scales favored by double-diffusive instabilities. Diffusive magneto-rotational instabilities are thus attractive candidates for angular momentum transport in the early Sun's radiative zone.

Motivated by the low-collisionality of gas accreted onto black holes in Sgr A* and other nearby galactic nuclei, we study a family of 2D advective accretion solutions with thermal conduction. While we only impose global inflow, the accretion flow spontaneously develops bipolar outflows. The role of conduction is key in providing the extra degree of freedom (latitudinal energy transport) necessary to launch these rotating thermal outflows. The sign of the Bernoulli constant does not discriminate between inflowing and outflowing regions. Our parameter survey covers mass outflow rates from ~ 0 to 13% of the net inflow rate, outflow velocities from ~0 to 11% of the local Keplerian velocity and outflow opening angles from ~ 0 to 60 degs. As the magnitude of conduction is increased, outflows can adopt a conical geometry, pure inflow solutions emerge, and the limit of 2D non-rotating Bondi-like solutions is eventually reached. These results confirm that radiatively-inefficient, hot accretion flows have a hydrodynamical propensity to generate bipolar thermal outflows.

Torque fluctuations due to magnetorotational turbulence in proto-planetary disks may greatly influence the migration patterns and survival probabilities of nascent planets. Provided that the turbulence is a stationary stochastic process with finite amplitude and correlation time, the resulting diffusive migration can be described with a Fokker-Planck equation, which we reduce to an advection-diffusion equation. We calibrate the coefficients with existing turbulent-disk simulations and mean-migration estimates, and solve the equation both analytically and numerically. Diffusion tends to dominate over advection for planets of low-mass and those in the outer regions of proto-planetary disks, whether they are described by the Minimum Mass Solar Nebula (MMSN) or by T-Tauri alpha disks. Diffusion systematically reduces the lifetime of most planets, yet it allows a declining fraction of them to survive for extended periods of time at large radii. Mean planet lifetimes can even be formally infinite (e.g. in an infinite steady MMSN), though median lifetimes are always finite. Surviving planets may linger near specific radii where the combined effects of advection and diffusion are minimized, or at large radii, depending on model specifics. The stochastic nature of migration in turbulent disks challenges deterministic planet formation scenarios and suggests instead that a wide variety of planetary outcomes are possible from similar initial conditions. This would contribute to the diversity of (extrasolar) planetary systems.

Cataclysmic Variables (CVs) should be present in large numbers in Globular Clusters (GCs). Numerous low-luminosity X-ray sources identified over the past few years as candidate CVs in GCs support this notion. Yet, very few "cataclysms," the characteristic feature of this class of objects in the field, have been observed in GCs. We address this discrepancy here, within the framework of the standard Disk Instability Model for CV outbursts. We argue that the paucity of outbursts in GCs is probably not a direct consequence of the donors' low metallicities. We present diagnostics based on outburst properties allowing tests of the hypothesis that rare cataclysms are entirely due to lower mass transfer rates in GCs relative to the field, and we argue against this explanation. Instead, we propose that a combination of low mass transfer rates (>~ 10^14-15 g/s) and moderately strong white dwarf magnetic moments (>~ 10^30 G cm^3) stabilize CV disks in GCs and thus prevent most of them from experiencing frequent outbursts. If it is so, rare cataclysms in GCs would signal important evolutionary differences between field and cluster CVs.

The gravitational waves (GW) emitted during the coalescence of supermassive black holes (SMBHs) in the mass range 10^4-10^7 M_sun will be detectable out to high redshifts with LISA. We calculate the size and orientation of the three-dimensional error ellipse in solid angle and redshift within which the LISA event could be localized using the GW signatures alone. We take into account uncertainties in LISA's measurements of the luminosity distance and direction to the source, in the background cosmology, in weak gravitational lensing magnification due to inhomogeneities along the line of sight, and potential source peculiar velocities. We find that weak lensing errors exceed other sources of uncertainties by nearly an order of magnitude. Under the plausible assumption that BH mergers are accompanied by gas accretion leading to Eddington-limited quasar activity, we then compute the number of quasars that would be found in a typical LISA error volume, as a function of BH mass and redshift. We find that low redshifts offer the best opportunities to identify quasar counterparts to cosmological standard sirens, and that the LISA error volume will typically contain a single near-Eddington quasar at z=1. This will allow a straightforward test of the hypothesis that BH mergers are accompanied by bright quasar activity and, if the hypothesis proves correct, will guarantee the identification of a unique quasar counterpart. This would yield unprecedented tests of the physics of SMBH accretion, and offer an alternative method to precisely constrain cosmological parameters [abridged].

Observations of the hot gas surrounding Sgr A* and a few other nearby galactic nuclei imply electron and proton mean free paths comparable to the gas capture radius: hot accretion likely proceeds under weakly-collisional conditions in these systems. As a result, thermal conduction, rather than convection, may be important on all scales and affect the global flow properties. The self-similar ADAF solution of Narayan & Yi (1994) is generalized to include a saturated form of thermal conduction, as is appropriate for the weakly-collisional regime of interest. Conduction provides extra heating and yet it reduces the free-free radiative efficiency of the accretion flow (by potentially large factors). These idealized solutions suggest that thermal conduction may be an important physical ingredient to understand hot accretion onto dim accreting black holes. Conduction could also play a role in reducing the rate at which black holes capture ambient gas and in providing an evaporation mechanism for an underlying cold thin disk.

We discuss atmosphere models of HD209458b in light of the recent day-side flux measurement of HD209458b's secondary eclipse by Spitzer-MIPS at 24 microns. In addition, we present a revised secondary eclipse IRTF upper limit at 2.2 microns which places a stringent constraint on the adjacent H2O absorption band depths. These two measurements are complementary because they are both shaped by H2O absorption and because the former is on the Wien tail of the planet's thermal emission spectrum and the latter is near the thermal emission peak. A wide range of models fit the observational data, confirming our basic understanding of hot Jupiter atmospheric physics. Although a range of models are viable, some models at the hot and cold end of the plausible temperature range can be ruled out. One class of previously unconsidered hot Jupiter atmospheric models that fit the data are those with C/O >~ 1 (as Jupiter may have), which have a significant paucity of H2O compared to solar abundance models with C/O = 0.5. The models indicate that HD209458b is in a situation intermediate between pure in situ reradiation and very efficient redistribution of heat; one which will require a careful treatment of atmospheric circulation. We discuss how future wavelength-dependent and phase-dependent observations will further constrain the atmospheric circulation regime. In the shorter term, additional planned measurements for HD209458b, especially Spitzer IRAC photometry, should lift many of the model degeneracies. Multiwavelength IR observations constrain the atmospheric structure and circulation properties of hot Jupiters and thus open a new chapter in quantitative extrasolar planetology.

Space-born gravitational-wave interferometers such as {\it LISA} will detect the gravitational wave (GW) signal from the inspiral, plunge and ringdown phases of massive black hole binary mergers at cosmological distances. From the inspiral waves, we will be able to measure the masses of the binaries' members; from the ringdown waves, we will be able to measure the mass of the final merged remnant. A subset of detected events allow the identification of both the inspiral and the ringdown waveforms in a given source, and thus allow a measurement of the total mass-energy lost to GWs over the coalescence, $M_{\rm GW}$. We define ``golden'' binary mergers to be those with measurement errors likely to be small enough for a physically useful determination of $M_{\rm GW}$. A detailed sensitivity study, combined with simple black hole population models, suggests that a few golden binary mergers may be detected during a 3-year {\it LISA} mission lifetime. Any such mass deficit measurement would constitute a robust and valuable observational test of strong-field relativistic gravity. An extension of this concept to include spin measurements may allow a direct empirical test of the black hole area theorem.

We present spectroscopic observations of the short-period cataclysmic variable SW Ursa Majoris, obtained by the Far Ultraviolet Spectroscopic Explorer (FUSE) satellite while the system was in quiescence. The data include the resonance lines of O VI at 1031.91 and 1037.61 A. These lines are present in emission, and they exhibit both narrow (~ 150 km/s) and broad (~ 2000 km/s) components. The narrow O VI emission lines exhibit unusual double-peaked and redshifted profiles. We attribute the source of this emission to a cooling flow onto the surface of the white dwarf primary. The broad O VI emission most likely originates in a thin, photoionized surface layer on the accretion disk. We searched for emission from H_2 at 1050 and 1100 A, motivated by the expectation that the bulk of the quiescent accretion disk is in the form of cool, molecular gas. If H_2 is present, then our limits on the fluxes of the H_2 lines are consistent with the presence of a surface layer of atomic H that shields the interior of the disk. These results may indicate that accretion operates primarily in the surface layers of the disk in SW UMa. We also investigate the far-UV continuum of SW UMa and place an upper limit of 15,000 K on the effective temperature of the white dwarf.

Recent improved determinations of the mass density rho_BH of supermassive black holes (SMBHs) in the local universe have allowed accurate comparisons of rho_BH with the amount of light received from past quasar activity. These comparisons support the notion that local SMBHs are ``dead quasars'' and yield a value epsilon >~ 0.1 for the average radiative efficiency of cosmic SMBH accretion. BH coalescences may represent an important component of the quasar mass assembly and yet not produce any observable electromagnetic signature. Therefore, ignoring gravitational wave (GW) emission during such coalescences, which reduces the amount of mass locked into remnant BHs, results in an overestimate of epsilon. Here, we put constraints on the magnitude of this bias. We calculate the cumulative mass loss to GWs experienced by a representative population of BHs during repeated cosmological mergers, using loss prescriptions based on detailed general relativistic calculations. Despite the possibly large number of mergers in the assembly history of each individual SMBH, we find that near--equal mass mergers are rare, and therefore the cumulative loss is likely to be modest, amounting at most to an increase by 20 percent of the inferred epsilon value. Thus, recent estimates of epsilon >~ 0.1 appear robust. The space interferometer LISA should provide empirical constraints on the dark side of quasar evolution, by measuring the masses and rates of coalescence of massive BHs to cosmological distances.

Two spectacular results emerging from recent studies of nearby dead quasars and distant active quasars are (i) the existence of tight relations between the masses of black holes (BHs) and the properties of their host galaxies (spheroid luminosity or velocity dispersion, galaxy mass), and (ii) a consistency between the local mass density in BHs and that expected by summing up the light received from distant active quasars. These results are partly shaped by successive galactic mergers and BH coalescences, since mergers redistribute the population of BHs in galaxies and BH binary coalescences reduce the mass density in BHs through losses to gravitational waves. Here, we isolate and quantify these effects by following the cosmological merger history of a population of massive BHs representing the quasar population between z = 3 and z =0. Our results suggest that the relation between BH mass and host galaxy properties, and inferences on the global efficiency of BH accretion during active quasar phases, could be influenced by the cumulative effect of repeated mergers.

The atmospheres of weakly-magnetized neutron stars expand hydrostatically and rotate differentially during thermonuclear X-ray bursts. Differential rotation is probably related to the frequency drifts of millisecond burst oscillations exhibited by about a dozen nuclear-powered X-ray pulsars. Here, we analyze the linear stability of this differential rotation with respect to local, axisymmetric, multi-diffusive MHD perturbations, at various heights in the neutron star atmosphere. Unstable magneto-rotational modes are identified from within to well above the burning layers. Properties of the fastest growing modes depend sensitively on the local magnetic field geometry. Linear estimates suggest that momentum transport due to magneto-rotational instabilities can affect atmospheric rotation profiles on time-scales relevant to burst oscillations. This transport would likely strengthen the coherence of burst oscillations and contribute to their drifts.

We present detailed estimates of ''type-I'' migration rates for low-mass proto-planets embedded in steady-state T-Tauri alpha-disks, based on Lindblad torque calculations ignoring feedback on the disk. Differences in migration rates for several plausible background disk models are explored and we contrast results obtained using the standard two dimensional formalism of spiral density wave theory with those obtained from a simple treatment of three-dimensional effects. Opacity transitions in the disk result in sudden radial variations of the migration rates. Regions with minimal migration rates may be preferred sites of gravitational interactions between proto-planets. Three-dimensional torques are significantly weaker than two-dimensional ones and they are sensitive to the surface density profile of the background disk. We find that migration times in excess of runaway envelope accretion times or T-Tauri disk lifetimes are possible for Earth-mass proto-planets in some background disk models, even at sub-AU distances. We conclude that an understanding of the background disk structure and ''viscosity'', as well as a proper treatment of three-dimensional effects in torque calculations, are necessary to obtain reliable estimates of ``type-I'' migration rates.

We study the local stability of stratified, differentially-rotating fluids to axisymmetric perturbations in the presence of a weak magnetic field and of finite resistivity, viscosity and heat conductivity. This is a generalization of the Goldreich-Schubert-Fricke (GSF) double-diffusive analysis to the magnetized and resistive, triple-diffusive case. Our fifth-order dispersion relation admits a novel branch which describes a magnetized version of multi-diffusive modes. We derive necessary conditions for axisymmetric stability in the inviscid and perfect-conductor (double-diffusive) limits. In each case, rotation must be constant on cylinders and angular velocity must not decrease with distance from the rotation axis for stability, irrespective of the relative strength of viscous, resistive and heat diffusion. Therefore, in both double-diffusive limits, solid body rotation marginally satisfies our stability criteria. The role of weak magnetic fields is essential to reach these conclusions. The triple-diffusive situation is more complex, and its stability criteria are not easily stated. Numerical analysis of our general dispersion relation confirms our analytic double-diffusive criteria, but also shows that an unstable double-diffusive situation can be significantly stabilized by the addition of a third, ostensibly weaker, diffusion process. We describe a numerical application to the Sun's upper radiative zone and establish that it would be subject to unstable multi-diffusive modes if moderate or strong radial gradients of angular velocity were present.

We report Chandra observations of the dwarf nova WX Hyi in quiescence. The X-ray spectrum displays strong and narrow emission lines of N, O, Mg, Ne, Si, S and Fe. The various ionization states implied by the lines suggest that the emission is produced within a flow spanning a wide temperature range, from T ~ 10^6 K to T >~ 10^8 K. Line diagnostics indicate that most of the radiation originates from a very dense region, with n ~ 10^{13}-10^{14} cm^{-3}. The Chandra data allow the first tests of specific models proposed in the literature for the X-ray emission in quiescent dwarf novae. We have computed the spectra for a set of models ranging from hot boundary layers, to hot settling flows solutions, to X-ray emitting coronae. WX Hyi differs from other dwarf novae observed at minimum in having much stronger low temperature lines, which prove difficult to fit with existing models, and possibly a very strong, broad O VII line, perhaps produced in a wind moving at a few x 10^3 km/s. The accretion rate inferred from the X-rays is lower than the value inferred from the UV. The presence of high-velocity mass ejection could account for this discrepancy while at the same time explaining the presence of the broad O VII line. If this interpretation is correct, it would provide the first detection of a wind from a dwarf nova in quiescence.

The local axisymmetric stability of hydrodynamical and magnetized, nearly-Keplerian gaseous accretion disks around non-rotating black holes is examined in the vicinity of the classical marginally-stable orbit (at radii ~ R_ms). An approximate Paczynski-Wiita pseudo-Newtonian potential is used. Hydrodynamical disks are linearly unstable inside a radius which differs slightly from the classical R_ms value because of finite pressure and radial stratification effects. Linear stresses associated with unstable hydrodynamical modes vanish exactly at the radius of marginal stability and are generally positive inside of that radius. When a magnetic field is introduced, however, the concept of radius of marginal stability becomes largely irrelevant because there are linearly unstable magneto-rotational modes everywhere. Associated linear stresses are positive and continuous across the region of hydrodynamical marginal stability, even for large-scale "hydro-like" modes subject only to weak magnetic tension. This conclusion is valid for arbitrarily thin disks (in ideal MHD) and it does not require a large-scale "radially-connecting" magnetic field. Results on hydrodynamical diskoseismic modes trapped in deep relativistic potential wells should be revised to account for the short-wavelength Alfven-like behavior of inertio-gravity waves in magnetized disks.

Shallow-water numerical simulations show that the atmospheric circulation of the close-in extrasolar giant planet (EGP) HD 209458b is characterized by moving circumpolar vortices and few bands/jets (in contrast with ~10 bands/jets and absence of polar vortices on cloud-top Jupiter and Saturn). The large spatial scales of moving circulation structures on HD 209458b may generate detectable variability of the planet's atmospheric signatures. In this Letter, we generalize these results to other close-in EGPs, by noting that shallow-water dynamics is essentially specified by the values of the Rossby (Ro) and Burger (Bu) dimensionless numbers. The range of likely values of Ro (~ 0.01 - 10) and Bu (~ 1 - 200) for the atmospheric flow of known close-in EGPs indicates that their circulation should be qualitatively similar to that of HD 209458b. This results mostly from the slow rotation of these tidally-synchronized planets.

High-resolution atmospheric flow simulations of the tidally-locked extrasolar giant planet, HD 209458b, show large-scale spatio-temporal variability. This is in contrast to the simple, permanent day/night (i.e., hot/cold) picture. The planet's global circulation is characterized by a polar vortex in motion around each pole and a banded structure corresponding to ~3 broad zonal (east-west) jets. For very strong jets, the circulation-induced temperature difference between moving hot and cold regions can reach up to ~1000 K, suggesting that atmospheric variability could be observed in the planet's spectral and photometric signatures.

Idealized models are used to illustrate the potential of the Laser Interferometer Space Antenna (LISA) as a probe of the largely unknown population of cosmologically-distant Massive Black Holes (MBHs) and as a tool to measure their masses with unprecedented accuracy. The models suggest that LISA will most efficiently probe a MBH population of lower mass than the one found in bright quasars and nearby galactic nuclei. The mass spectrum of these MBHs could constrain formation scenarios for high-redshift, low-mass galaxies.

Habitability is usually defined as the requirement for a terrestrial planet's atmosphere to sustain liquid water. This definition can be complemented by the dynamical requirement that other planets in the system do not gravitationally perturb terrestrial planets outside of their habitable zone, the orbital region allowing the existence of liquid water. We quantify the dynamical habitability of 85 known extrasolar planetary systems via simulations of their orbital dynamics in the presence of potentially habitable terrestrial planets. When requiring that habitable planets remain strictly within their habitable zone at all time, the perturbing influence of giant planets extends beyond the traditional Hill sphere for close encounters: terrestrial planet excursions outside of the habitable zone are also caused by secular eccentricity variations and, in some cases, strong mean-motion resonances. Our results indicate that more than half the known extrasolar planetary systems (mostly those with distant, eccentric giant planets) are unlikely to harbor habitable terrestrial planets. About 1/4 of the systems (mostly those with close-in giant planets), including 1/3 of the potential targets for the Terrestrial Planet Finder, appear as dynamically habitable as our own Solar System. The influence of yet undetected giant planets in these systems could compromise their dynamical habitability. Some habitable terrestrial planets in our simulations have substantial eccentricities (e > 0.1) which may lead to large seasonal climate variations and thus affect their habitability.

We report on XMM-Newton observations of GRO J1655-40 and GRS 1009-45, which are two black hole X-ray transients currently in their quiescent phase. GRO J1655-40 was detected with a 0.5 - 10 keV luminosity of 5.9 10^{31} erg/s. This luminosity is comparable to a previous Chandra measurement, but ten times lower than the 1996 ASCA value, most likely obtained when the source was not yet in a true quiescent state. Unfortunately, XMM-Newton failed to detect GRS 1009-45. A stringent upper limit of 8.9 10^{30} erg/s was derived by combining data from the EPIC-MOS and PN cameras. The X-ray spectrum of GRO J1655-40 is very hard as it can be fitted with a power law model of photon index ~ 1.3 +/- 0.4. Similarly hard spectra have been observed from other systems; these rule out coronal emission from the secondary or disk flares as the origin of the observed X-rays. On the other hand, our observations are consistent with the predictions of the disc instability model in the case that the accretion flow forms an advection dominated accretion flow (ADAF) at distances less than a fraction ~ 0.1 - 0.3) of the circularization radius. This distance corresponds to the greatest extent of the ADAF that is thought to be possible.

We model the 1990 giant X-ray flare of the quiescent galaxy NGC 5905 as the tidal disruption of a star by a supermassive black hole. From the observed rapid decline of the luminosity, over a timescale of a few years, we argue that the flare was powered by the fallback of debris rather than subsequent accretion via a thin disk. The fallback model allows constraints to be set on the black hole mass and the mass of debris. The latter must be very much less than a solar mass to explain the very low luminosity of the flare. The observations can be explained either as the partial stripping of the outer layers of a low-mass main sequence star or as the disruption of a brown dwarf or a giant planet. We find that the X-ray emission in the flare must have originated within a small patch rather than over the entire torus of circularized material surrounding the black hole. We suggest that the patch corresponds to the ``bright spot'' where the stream of returning debris impacts the torus. Interestingly, although the peak luminosity of the flare was highly sub-Eddington, the peak flux from the bright spot was close to the Eddington limit. We speculate on the implications of this result for observations of other flare events.

We discuss the optical and radio properties of 30,000 FIRST sources positionally associated with an SDSS source in 1230 deg$^2$ of sky. The majority (83%) of the FIRST sources identified with an SDSS source brighter than r=21 are optically resolved. We estimate an upper limit of 5% for the fraction of quasars with broad-band optical colors indistinguishable from those of stars. The distribution of quasars in the radio flux -- optical flux plane supports the existence of the "quasar radio-dichotomy"; 8% of all quasars with i<18.5 are radio-loud and this fraction seems independent of redshift and optical luminosity. The radio-loud quasars have a redder median color by 0.08 mag, and a 3 times larger fraction of objects with red colors. FIRST galaxies represent 5% of all SDSS galaxies with r<17.5, and 1% for r<20, and are dominated by red galaxies. Magnitude and redshift limited samples show that radio galaxies have a different optical luminosity distribution than non-radio galaxies selected by the same criteria; when galaxies are further separated by their colors, this result remains valid for both blue and red galaxies. The distributions of radio-to-optical flux ratio are similar for blue and red galaxies in redshift-limited samples; this similarity implies that the difference in their luminosity functions, and resulting selection effects, are the dominant cause for the preponderance of red radio galaxies in flux-limited samples. We confirm that the AGN-to-starburst galaxy number ratio increases with radio flux, and find that radio emission from AGNs is more concentrated than radio emission from starburst galaxies (abridged).

We present dynamically consistent solutions for hot accretion onto unmagnetized, rotating white dwarfs (WDs) in five quiescent dwarf novae. The measured WD rotation rates (and other system parameters) in RX And, SS Cyg, U Gem, VW Hyi and WZ Sge imply spindown of the WD by an extended hot flow emitting most of its X-rays in the vicinity of the stellar surface. In general, energy advection is absent and the flow is stable to convection and hydrodynamical outflows. In rapidly rotating systems, the X-ray luminosity provides only an upper limit on the quiescent accretion rate because of substantial stellar spindown luminosity. We suggest that the presence of hot flows in quiescent dwarf novae may limit the long-term WD rotation rates to significantly sub-Keplerian values.

We show that accretion disks in several compact X-ray binaries with hydrogen-depleted donors are likely subject to a thermal ionization instability, unless they are strongly irradiated. These disks are particularly interesting in that their MHD-turbulent properties in the neutral phase may be quite different from those of standard, hydrogen-rich disks.

Recent work on the modes of accretion onto black holes (BHs) in Soft X-Ray Transients (SXTs) is reviewed, with an emphasis on uncertainties affecting models of accretion during quiescence (inner hot flow, outer thin disk). Various interpretations of the quiescent X-ray luminosity difference between systems containing neutron stars (NSs) and systems containing BH candidates are also summarized. A new scenario, which does not require BH candidates to possess an event horizon, is presented here. This scenario may be ruled out in the future, from detailed X-ray spectroscopic diagnostics or from the absence of type I X-ray bursts in systems containing BH candidates.

The tidal disruption of a star by a supermassive black hole is expected to lead to a short bright flare followed by an extended period of low-level emission. Existing models of the late-time accretion of the stellar debris via a thin disk imply that the Galactic Center source Sgr A* should be brighter than currently observed by several orders of magnitude. A similar problem exists for M31 and M32. This problem could be avoided if thin disk accretion transitions to low radiative efficiency accretion (e.g., ADAF, CDAF) at low accretion rates (via ``evaporation'' of the thin disk). Alternatively, we show that the outer portions of a thin disk may become neutral and unable to sustain MHD turbulence on a timescale less than the time between tidal disruption events; this may dramatically shorten the duration of the late accretion phase.

We positionally match sources observed by the Sloan Digital Sky Survey (SDSS), the Two Micron All Sky Survey (2MASS), and the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. Practically all 2MASS sources are matched to an SDSS source within 2 arcsec; ~11% of them are optically resolved galaxies and the rest are dominated by stars. About 1/3 of FIRST sources are matched to an SDSS source within 2 arcsec; ~80% of these are galaxies and the rest are dominated by quasars. Based on these results, we project that by the completion of these surveys the matched samples will include about 10^7 stars and 10^6 galaxies observed by both SDSS and 2MASS, and about 250,000 galaxies and 50,000 quasars observed by both SDSS and FIRST. Here we present a preliminary analysis of the optical, infrared and radio properties for the extragalactic sources from the matched samples. In particular, we find that the fraction of quasars with stellar colors missed by the SDSS spectroscopic survey is probably not larger than ~10%, and that the optical colors of radio-loud quasars are ~0.05 mag. redder (with 4-sigma significance) than the colors of radio-quiet quasars.

Our understanding of the structure and properties of accretion disks in quiescent Dwarf Novae and X-ray Transients is very limited. Observations during quiescence challenge some of the standard Disk Instability Model (DIM) predictions. Our ignorance of the nature of viscosity may be the main source of uncertainties in quiescent disk models. A ``layered accretion'' alternative to the DIM, in which accretion proceeds via X-ray ionized surface layers, illustrates the magnitude of these uncertainties. A detection of molecular hydrogen in quiescent disks may provide direct constraints on the nature of viscosity.

We present 13 Broad Absorption Line (BAL) quasars, including 12 new objects, which were identified in the Sloan Digital Sky Survey (SDSS) and matched within 2'' to sources in the FIRST radio survey catalog. The surface density of this sample of radio-detected BAL quasars is 4.5 +- 1.2 per 100 deg^2, i.e. approximately 4 times larger than previously found by the shallower FIRST Bright Quasar Survey (FBQS). A majority of these radio-detected BAL quasars are moderately radio-loud objects. The fraction of BAL quasars in the entire radio quasar sample, 4.8 +- 1.3 %, is comparable to the fraction of BAL quasars among the SDSS optical quasar sample (ignoring selection effects). We estimate that the true fraction of BAL quasars (mostly HiBALs) in the radio sample is 9.2 +- 2.6 % once selection effects are accounted for. We caution that the absorption troughs of 4 of the 13 radio-detected quasars considered do not strictly satisfy the standard BALnicity criterion. One or possibly two of the new radio-detected BAL quasars are of the rare ``iron LoBAL'' type. BAL quasars are generally redder than the median SDSS quasar at the same redshift.

We show that thin accretion disks made of Carbon or Oxygen are subject to the same thermal ionization instability as Hydrogen and Helium disks. We argue that the instability applies to disks of any metal content. The relevance of the instability to supernova fallback disks probably means that their power-law evolution breaks down when they first become neutral. We construct simple analytical models for the viscous evolution of fallback disks to show that it is possible for these disks to become neutral when they are still young (ages of a few 10^3 to 10^4 years), compact in size (a few 10^9 cm to 10^11 cm) and generally accreting at sub-Eddington rates (Mdot ~ a few 10^14 - 10^18 g/s). Based on recent results on the nature of viscosity in the disks of close binaries, we argue that this time may also correspond to the end of the disk activity period. Indeed, in the absence of a significant source of viscosity in the neutral phase, the entire disk will likely turn to dust and become passive. We discuss various applications of the evolutionary model, including anomalous X-ray pulsars and young radio pulsars. Our analysis indicates that metal-rich fallback disks around newly-born neutron stars and black holes become neutral generally inside the tidal truncation radius (Roche limit) for planets, at \~10^11 cm. Consequently, the efficiency of the planetary formation process in this context will mostly depend on the ability of the resulting disk of rocks to spread via collisions beyond the Roche limit. It appears easier for the merger product of a doubly degenerate binary, whether it is a massive white dwarf or a neutron star, to harbor planets because it can spread beyond the Roche limit before becoming neutral.[Abridged]

We present a model for the spindown of young radio pulsars in which the neutron star loses rotational energy not only by emitting magnetic dipole radiation but also by torquing a surrounding accretion disk produced by supernova fallback. The braking index predicted in our model is in general less than n=3 (the value for pure dipole magnetic radiation), in agreement with the reported values of n<3 for five young radio pulsars. With an accuracy of 30% or better, our model reproduces the age, braking index and third frequency derivative of the Crab pulsar for a disk accretion rate in the range 3 x 10^16 - 10^17 g/s.

The ubiquity of supermassive black holes (SMBHs) at the centers of nearby luminous galaxies can arise from the multiple mergers experienced by dark matter halos in hierarchical structure formation models, even if only a small fraction of these galaxies harbor SMBHs at high redshifts. We illustrate this possibility using cosmological Monte Carlo simulations of the merger history of dark matter halos and their associated SMBHs. In our most extreme models, in order to populate nearly every bright galaxy with a SMBH at z=0, only a few percent of the halos with virial temperatures above 10^4 K are required to harbor a SMBH at high redshift. This possibility must be included in studies of the luminosity function and the clustering properties of quasars. We predict the number of SMBH merger events that are detectable by the gravitational wave experiment LISA, as a function of redshift, out to z=5. Although the event rates can be significantly reduced in scenarios with rare SMBHs, a minimum of \~10 detectable merger events per year is predicted if SMBH binaries coalesce efficiently. The observed distribution of events with redshift could yield valuable information on the SMBH formation process. If SMBH binaries do not coalesce, we find that at least several SMBH slingshot ejections probably occurred from z=5 to the present in each galaxy more massive than ~10^11 Msun at z=0. Although our results are sensitive to the minimum cooling mass assumed for the formation of SMBHs, we expect the qualitative predictions of our models to be robust.

We use the observed optical-UV and X-ray emission spectrum of Cen X-4 during quiescence to constrain models for the accretion flow in this system. We argue that the optical-UV emission is not due to an optically-thick quiescent accretion disk, nor due to synchrotron emission from an Advection-Dominated Accretion Flow (ADAF). Emission from the bright spot could account for the observed optical-UV component if the mass transfer rate in Cen X-4 is >~ 2.10^16 g/s. Although the presence of an ADAF around the neutron star leads to Compton upscattering of the soft X-ray photons radiated from the stellar surface, we find that this process alone cannot account for the power law component seen in the quiescent X-ray spectrum of Cen X-4 and other X-ray transients containing neutron stars; this result is independent of whether the source of soft photons is incandescent thermal emission or accretion-powered emission. We conclude that, in models which invoke the presence of an ADAF and a propeller effect for the quiescence of X-ray transients containing neutron stars, the intrinsic emission from the ADAF must contribute very little to the optical-UV and X-ray emission observed. If these ADAF+propeller models are correct, the X-ray power law component observed must arise from regions where the gas impacts the neutron star surface. Variability studies could greatly help clarify the role of the various emission mechanisms involved.

We report on optical spectroscopic observations of the X-ray Transient XTE J1118+480 covering the period from April 7, 2000 to July 4, 2000. The spectrum is characterized by weak, broad, double-peaked Balmer and He lines on top of a blue continuum of slope p ~ 1/3, as expected for an optically-thick accretion disk. The weak Bowen blend seen in our spectra may indicate a low intrinsic metallicity for the source. The presence of a partial S-wave pattern in the HeII lambda-4686 line appears consistent with the reported photometric orbital period Porb = 4.1 hr for XTE J1118+480. By using a combination of Doppler mapping and various theoretical arguments, we constrain plausible orbital parameters for the system: a mass ratio 0.02 <~ q <~ 0.1, an inclination i >~ 70 deg for a neutron star primary, or 30 deg <~ i <~ 50 deg for a black hole primary with a mass between 4 and 10 Msun. CaII absorption features observed at very high resolution constrain the interstellar hydrogen absorption column log[N_HI (cm^{-2})] ~ 20.45 +/- 0.2 and the identification of three absorbing clouds indicate a distance to the source <~ 1 kpc, assuming the line-of-sight to XTE J1118+480 has average high-latitude properties. These results are discussed in the context of previous multiwavelength observations of this unusual system.

We consider the combined role of the thermal ionization, magneto-rotational and gravitational instabilities in thin accretion disks around supermassive black holes. We find that in the portions of the disk unstable to the ionization instability, the gas remains well coupled to the magnetic field even on the cold, neutral branch of the thermal limit cycle. This suggests that the ionization instability is not a significant source of large amplitude time-dependent accretion in AGN. We also argue that, for accretion rates greater than 10^{-2} solar masses per year, the gravitationally unstable and magneto-rotationally unstable regions of the accretion disk overlap; for lower accretion rates they may not. Some low-luminosity AGN, e.g. NGC 4258, may thus be in a transient phase in which mass is building up in a non-accreting gravitationally and magneto-rotationally stable ``dead zone.'' We comment on possible implications of these findings.

Rather than accreting via a disk, some White Dwarfs (WDs) in quiescent Dwarf Novae (DN) could accrete via an Advection-Dominated Accretion Flow (ADAF) possibly responsible for the X-ray Bremsstrahlung emission observed. Such a hot accretion flow is also expected to produce characteristic thermal line emission. Using SS Cyg as a specific example, we show that knowing a priori the inclination and the WD mass in quiescent DN makes X-ray line diagnostics powerful probes of the flow structure in these systems. Current X-ray instruments can discriminate, from their width, between lines emitted from a flow with a Keplerian rotation rate and those emitted at a substantially sub-Keplerian rate. This could be used to observationally test the property of energy advection, which is at the origin of the partial radial pressure support by the hot gas and the sub-Keplerian rotation rate in an ADAF.

The self-sustained turbulence which develops in magnetized accretion disks is suppressed in the weakly-ionized, quiescent disks of close binary stars. Because accretion still proceeds during quiescence, another viscosity mechanism operates in these systems. An anticorrelation of the recurrence times of SU UMa dwarf novae with their mass ratio supports spiral waves or shock-waves tidally induced by the companion star as the main process responsible for accretion in the quiescent disks. Other weakly-ionized gaseous disks in systems lacking a massive companion have to rely on yet another transport mechanism or they could be essentially passive.

The boundary layers of weakly-magnetized white dwarfs (WDs) accreting at rates <=10^16 g/s are radially extended, hot, optically-thin, and they advect some of their internally-dissipated energy (Narayan & Popham 1993). Motivated by this, I construct here idealized spectral models of an Advection-Dominated Accretion Flow (ADAF) around a WD, for application to quiescent Dwarf Novae (DN). The Bremsstrahlung cooling of the gas in the ADAF, with temperatures ranging from a few keV to a few tens of keV, can account for the X-ray emission properties of quiescent DN. If the energy advected by the flow is thermalized in the WD atmosphere, the resulting emission from the entire stellar surface (blackbody of temperature T_eff ~ 5 eV) outshines the X-ray luminosity substantially. This extreme-UV component provides a flux in the 0.055-0.28 keV band which is sufficient to power the strong HeII lambda4686 emission lines of quiescent DN by photoionization of the disk material. Reprocessing of the ADAF X-ray emission by a cold outer thin disk could also lead to an observable iron Kalpha fluorescence emission line, which can be used to probe the geometry of the accretion flow. Existing observational data indicate that the presence of ADAFs in quiescent DN is not ubiquitous, while future observations, in particular with the X-ray satellites Chandra and XMM-Newton, have the potential to detect signatures of the hot flow in promising candidates.

We construct time-dependent models of accretion discs around black holes and neutron stars. We investigate the effect of evaporating the disc inner regions during quiescence on the predictions of the Disc Instability Model (DIM) for these systems. We do not include irradiation of the disc in the models. Removing the inner, most unstable parts of the accretion disc increases the predicted recurrence times. However, DIMs with values of the viscosity parameter alpha_hot ~ 0.1 and alpha_cold ~ 0.02 (typical of applications of the DIM to standard dwarf nova outbursts) fail to reproduce the long recurrence times of SXTs (unless we resort to fine-tuning of the parameters) independent of the evaporation strength. We show that models with evaporation and a smaller value of alpha_cold (~ 0.005) do reproduce the long recurrence times and the accretion rates at the level of the Eddington rate observed in outburst. The large difference between the values of alpha_hot and alpha_cold, if confirmed when disc irradiation is included, suggests that several viscosity mechanisms operate in these accretion discs. For some choices of parameters our models predict reflares during the decline from outburst. They are a physical property of the model and result from a heating front forming in the wake of an initial cooling front and subsequent, multiple front reflections. The reflares disappear in low-alpha models where front reflection can not occur.

We consider a class of models for the redshift evolution (between $0\lsim z \lsim 4$) of the observed optical and X-ray quasar luminosity functions (LFs), with the following assumptions: (i) the mass-function of dark matter halos follows the Press-Schechter theory, (ii) the black hole (BH) mass scales linearly with the halo mass, (iii) quasars have a constant universal lifetime, and (iv) a thin accretion disk provides the optical luminosity of quasars, while the X-ray/optical flux ratio is calibrated from a sample of observed quasars. The mass accretion rate $\dot{M}$ onto quasar BHs is a free parameter of the models, that we constrain using the observed LFs. The accretion rate $\dot M$ inferred from either the optical or X-ray data under these assumptions generally decreases as a function of cosmic time from $z \simeq 4$ to $z \simeq 0$. We find that a comparable accretion rate is inferred from the X-ray and optical LF only if the X-ray/optical flux ratio decreases with BH mass. Near $z\simeq 0$, $\dot M$ drops to substantially sub-Eddington values at which advection-dominated accretion flows (ADAFs) exist. Such a decline of $\dot M$, possibly followed by a transition to radiatively inefficient ADAFs, could explain both the absence of bright quasars in the local universe and the faintness of accreting BHs at the centers of nearby galaxies. We argue that a decline of the accretion rate of the quasar population is indeed expected in cosmological structure formation models.

We consider the X-ray luminosity difference between neutron star and black hole soft X-ray transients (NS and BH SXTs) in quiescence. The current observational data suggest that BH SXTs are significantly fainter than NS SXTs. The luminosities of quiescent BH SXTs are consistent with the predictions of binary evolution models for the mass transfer rate if (1) accretion occurs via an ADAF in these systems and (2) the accreting compact objects have event horizons. The luminosities of quiescent NS SXTs are not consistent with the predictions of ADAF models when combined with binary evolution models, unless most of the mass accreted in the ADAF is prevented from reaching the neutron star surface. We consider the possibility that mass accretion is reduced in quiescent NS SXTs because of an efficient propeller and develop a model of the propeller effect that accounts for the observed luminosities. We argue that modest winds from ADAFs are consistent with the observations while strong winds are probably not.

We show that the outcome of disc instability models is strongly influenced by boundary conditions such as the position of the inner and outer disc edges. We discuss other sources of uncertainties, such as the tidal torque, and we conclude that disc illumination, disk size variations and a proper prescription for the tidal torque must be included in models if one wishes to extract meaningful physical information on e.g. viscosity from the comparison of predicted and observed lightcurves.

We use high-resolution time-dependent numerical simulations of accretion discs around white dwarfs to study the structure and properties of transition fronts in the context of the thermal-viscous disc instability model. The thermal structure of cooling and heating fronts is dominated by radiative cooling and viscous heating, respectively, except in a very narrow precursor region in heating fronts where advection and radial transport of energy dominate. Cooling fronts are much broader than heating fronts, but the widths of both types of fronts scale with the local vertical scale height of the disc. We confirm that during a fair fraction of the propagation time of a cooling front, the structure of the inner disc is close to self-similar. The speed of heating fronts is ~ a few km/s, while the speed of cooling fronts is ~ a fraction of a km/s. We show that direct measurements of the speed of transition fronts probably cannot discriminate between various prescriptions proposed for the viscosity parameter alpha. A natural prediction of the disc instability model is that fronts decelerate as they propagate in the disc, independent of the prescription for alpha. Observation of this effect would confirm that dwarf nova outbursts are driven by the thermal-viscous instability. Most of our results also apply to low mass X-ray binaries in which the accreting object is a neutron star or a black hole.

We study the thermal and viscous stability of accretion flows in Low Mass Black Hole Binaries (LMBHBs). We consider a model in which an inner advection-dominated accretion flow (ADAF) is surrounded by a geometrically thin accretion disk, the transition between the two zones occurring at a radius R_tr. In all the known LMBHBs, R_tr appears to be such that the outer disks could suffer from a global thermal-viscous instability. This instability is likely to cause the transient behavior of these systems. However, in most cases, if R_tr were slightly larger than the estimated values, the systems would be globally stable. This suggests that a population of faint persistent LMBHBs with globally stable outer disks could be present in the Galaxy. Such LMBHBs would be hard to detect because they would lack large amplitude outbursts, and because their ADAF zones would have very low radiative efficiencies, making the systems very dim. We present model spectra of such systems covering the optical and X-ray bands.

We have developed 1D time-dependent numerical models of accretion discs, using an adaptive grid technique and an implicit numerical scheme, in which the disc size is allowed to vary with time. The code fully resolves the cooling and heating fronts propagating in the disc. We show that models in which the radius of the outer edge of the disc is fixed produce incorrect results, from which probably incorrect conclusions about the viscosity law have been inferred. In particular we show that outside-in outbursts are possible when a standard bimodal behaviour of the Shakura-Sunyaev viscosity parameter alpha is used. We also discuss to what extent insufficient grid resolutions have limited the predictive power of previous models. We find that the global properties (magnitudes, etc. ...) of transient discs can be addressed by codes using a high, but reasonable, number of fixed grid points. However, the study of the detailed physical properties of the transition fronts generally requires resolutions which are out of reach of fixed grid codes. It appears that most time-dependent models of accretion discs published in the literature have been limited by resolution effects, improper outer boundary conditions, or both.

Astronomers have discovered many potential black holes in X-ray binaries and galactic nuclei. These black holes are usually identified by the fact that they are too massive to be neutron stars. Until recently, however, there was no convincing evidence that the objects identified as black hole candidates actually have event horizons. This has changed with extensive applications of a class of accretion models for describing the flow of gas onto compact objects; for these solutions, called advection-dominated accretion flows (ADAFs), the black hole nature of the accreting star, specifically its event horizon, plays an important role. We review the evidence that, at low luminosities, accreting black holes in both X-ray binaries and galactic nuclei contain ADAFs rather than the standard thin accretion disk.

We show that dwarf nova disks in quiescence have rather low magnetic Reynolds number, of order 10^3. Numerical simulations of magnetized accretion disks suggest that under these conditions magnetohydrodynamic turbulence and the associated angular momentum transport is sharply reduced. This could be the physical origin of episodic accretion in dwarf nova disks. If so, the standard disk instability model needs to be revised.