The obliquity of a terrestrial planet is an important clue about its formation and critical to its climate. Previous studies using simulated photometry of Earth show that continuous observations over most of a planet's orbit can be inverted to infer obliquity. However, few studies of more general planets with arbitrary albedo markings have been made and, in particular, a simple theoretical understanding of why it is possible to extract obliquity from light curves is missing. Reflected light seen by a distant observer is the product of a planet's albedo map, its host star's illumination, and the visibility of different regions. It is useful to treat the product of illumination and visibility as the kernel of a convolution. Time-resolved photometry constrains both the albedo map and the kernel, the latter of which sweeps over the planet due to rotational and orbital motion. The kernel's movement distinguishes prograde from retrograde rotation for planets with non-zero obliquity on inclined orbits. We demonstrate that the kernel's longitudinal width and mean latitude are distinct functions of obliquity and axial orientation. Notably, we find that a planet's spin axis affects the kernel -- and hence time-resolved photometry -- even if this planet is East-West uniform or spinning rapidly, or if it is North-South uniform. We find that perfect knowledge of the kernel at 2-4 orbital phases is usually sufficient to uniquely determine a planet's spin axis. Surprisingly, we predict that East-West albedo contrast is more useful for constraining obliquity than North-South contrast.

Future space telescopes will directly image extrasolar planets at visible wavelengths. Time-resolved reflected light from an exoplanet encodes information about atmospheric and surface inhomogeneities. Previous research has shown that the light curve of an exoplanet can be inverted to obtain a low-resolution map of the planet, as well as constraints on its spin orientation. Estimating the uncertainty on 2D albedo maps has so far remained elusive. Here we present exocartographer, a flexible open-source Bayesian framework for solving the exo-cartography inverse problem. The map is parameterized with equal-area HEALPix pixels. For a fiducial map resolution of 192 pixels, a four-parameter Gaussian process describing the spatial scale of albedo variations, and two unknown planetary spin parameters, exocartographer explores a 198-dimensional parameter space. To test the code, we produce a light curve for a cloudless Earth in a face-on orbit with a 90$^\circ$ obliquity. We produce synthetic white light observations of the planet: 5 epochs of observations throughout the planet's orbit, each consisting of 24 hourly observations with a photometric uncertainty of $1\%$ (120 data). We retrieve an albedo map and$-$for the first time$-$its uncertainties, along with spin constraints. The albedo map is recognizably of Earth, with typical uncertainty of $30\%$. The retrieved characteristic length scale is 88$\pm 7 ^\circ$, or 9800 km. The obliquity is recovered with a $1-\sigma$ uncertainty of $0.8^\circ$. Despite the uncertainty in the retrieved albedo map, we robustly identify a high albedo region (the Sahara desert) and a large low-albedo region (the Pacific Ocean).

The disk-integrated reflected brightness of an exoplanet changes as a function of time due to orbital and rotational motion coupled with an inhomogeneous albedo map. We have previously derived analytic reflected lightcurves for spherical harmonic albedo maps in the special case of a synchronously-rotating planet on an edge-on orbit (Cowan, Fuentes & Haggard 2013). In this paper, we present analytic reflected lightcurves for the general case of a planet on an inclined orbit, with arbitrary spin period and non-zero obliquity. We do so for two different albedo basis maps: bright points ($\delta$-maps), and spherical harmonics ($Y_l^m$-maps). In particular, we use Wigner $D$-matrices to express an harmonic lightcurve for an arbitrary viewing geometry as a non-linear combination of harmonic lightcurves for the simpler edge-on, synchronously rotating geometry. These solutions will enable future exploration of the degeneracies and information content of reflected lightcurves, as well as fast calculation of lightcurves for mapping exoplanets based on time-resolved photometry. To these ends we make available Exoplanet Analytic Reflected Lightcurves (EARL), a simple open-source code that allows rapid computation of reflected lightcurves.

Photometry of short-period planetary systems allows astronomers to monitor exoplanets, their host stars, and their mutual interactions. In addition to the transits of a planet in front of its star and the eclipses of the planet by its star, researchers have reported flux variations at the orbital frequency and its harmonics: planetary reflection and/or emission and Doppler beaming of starlight produce one peak per orbit, while ellipsoidal variations of a tidally distorted star and/or planet produce two maxima per orbit. Researchers have also reported significant photometric variability at three times the orbital frequency, sometimes much greater than the predictions of tidal theory. The reflected phase variations of a homogeneous planet contains power at even orbital harmonics-important for studies of ellipsoidal variations-but cannot contain odd orbital harmonics. We show that odd harmonics can, however, be produced by an edge-on planet with a time-variable map, or an inclined planet with a North-South (N-S) asymmetric map. Either of these scenarios entail weather: short-period planets are expected to have zero obliquity and hence N-S symmetric stellar forcing. North-South asymmetry in a giant planet would therefore suggest stochastic localized features, such as weather. However, we find that previous claims of large-amplitude odd modes in Kepler photometry are artifacts of removing planetary transits rather than modeling them. The only reliable claims of odd harmonics remain HAT-P-7b and Kepler-13Ab, for which the third mode amplitude is 6-8% of the planetary flux. Although time-variable albedo maps could in principle explain these odd harmonics, upper-limits on the infrared variability of hot Jupiters make this hypothesis unlikely. We urge theorists to study the effects of close-in planets on stellar atmospheres, as this remains the only plausible hypothesis.

We show that there is a classical metric satisfying the Einstein equations outside a finite spacetime region where matter collapses into a black hole and then emerges from a white hole. We compute this metric explicitly. We show how quantum theory determines the (long) time for the process to happen. A black hole can thus quantum-tunnel into a white hole. For this to happen, quantum gravity should affect the metric also in a small region outside the horizon: we show that contrary to what is commonly assumed, this is not forbidden by causality or by the semiclassical approximation, because quantum effects can pile up over a long time. This scenario alters radically the discussion on the black hole information puzzle.

[Abridged] Distant stars and planets will remain spatially unresolved for the foreseeable future. It is nonetheless possible to infer aspects of their brightness markings and viewing geometries by analyzing disk-integrated rotational and orbital brightness variations. We compute the harmonic lightcurves, F_l^m(t), resulting from spherical harmonic maps of intensity or albedo, Y_l^m(theta,phi), where l and m are the total and longitudinal order. Notably, odd m>1 are present in an inclined lightcurve, but not seen by an equatorial observer. We therefore suggest that the Fourier spectrum of a thermal lightcurve may be sufficient to determine the orbital inclination of non-transiting short-period planets, the rotational inclination of stars and brown dwarfs, and the obliquity of directly imaged planets. In the best-case scenario of a nearly edge-on geometry, measuring the m=3 mode of a star's rotational lightcurve to within a factor of two provides an inclination estimate good to +/- 6 degrees, assuming stars have randomly distributed spots. Alternatively, if stars have brightness maps perfectly symmetric about the equator, their lightcurves will have no m=3 power, regardless of orientation. In general, inclination estimates will remain qualitative until detailed hydrodynamic simulations and/or occultation maps can be used as a calibrator. We further derive harmonic reflected lightcurves for tidally-locked planets; these are higher-order versions of the well-known Lambert phase curve. We show that a non-uniform planet may have an apparent albedo 25% lower than its intrinsic albedo, even if it exhibits precisely Lambertian phase variations. Lastly, we provide low-order analytic expressions for harmonic lightcurves that can be used for fitting observed photometry; as a general rule, edge-on solutions cannot simply be scaled by sin(i) to mimic inclined lightcurves.