The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility designed to study a range of astrophysically relevant plasma processes as well as novel geometries that mimic astrophysical systems. A multi-cusp magnetic bucket constructed from strong samarium cobalt permanent magnets now confines a 10 m$^3$, fully ionized, magnetic-field free plasma in a spherical geometry. Plasma parameters of $ T_{e}\approx5$ to $20$ eV and $n_{e}\approx10^{11}$ to $5\times10^{12}$ cm$^{-3}$ provide an ideal testbed for a range of astrophysical experiments including self-exciting dynamos, collisionless magnetic reconnection, jet stability, stellar winds, and more. This article describes the capabilities of WiPAL along with several experiments, in both operating and planning stages, that illustrate the range of possibilities for future users.

The Dark Energy Camera is a new imager with a 2.2-degree diameter field of view mounted at the prime focus of the Victor M. Blanco 4-meter telescope on Cerro Tololo near La Serena, Chile. The camera was designed and constructed by the Dark Energy Survey Collaboration, and meets or exceeds the stringent requirements designed for the wide-field and supernova surveys for which the collaboration uses it. The camera consists of a five element optical corrector, seven filters, a shutter with a 60 cm aperture, and a CCD focal plane of 250 micron thick fully-depleted CCDs cooled inside a vacuum Dewar. The 570 Mpixel focal plane comprises 62 2kx4k CCDs for imaging and 12 2kx2k CCDs for guiding and focus. The CCDs have 15 microns x15 microns pixels with a plate scale of 0.263 arc sec per pixel. A hexapod system provides state-of-the-art focus and alignment capability. The camera is read out in 20 seconds with 6-9 electrons readout noise. This paper provides a technical description of the camera's engineering, construction, installation, and current status.

We perform numerical optimization of the axisymmetric flows in a sphere to minimize the critical magnetic Reynolds number Rm_cr required for dynamo onset. The optimization is done for the class of laminar incompressible flows of von Karman type satisfying the steady-state Navier-Stokes equation. Such flows are determined by equatorially antisymmetric profiles of driving azimuthal (toroidal) velocity specified at the spherical boundary. The model is relevant to the Madison plasma dynamo experiment (MPDX), whose spherical boundary is capable of differential driving of plasma in the azimuthal direction. We show that the dynamo onset in this system depends strongly on details of the driving velocity profile and the fluid Reynolds number Re. It is found that the overall lowest Rm_cr~200 is achieved at Re~240 for the flow, which is hydrodynamically marginally stable. We also show that the optimized flows can sustain dynamos only in the range Rm_cr<Rm<Rm_cr2, where Rm_cr2 is the second critical magnetic Reynolds number, above which the dynamo is quenched. Samples of the optimized flows and the corresponding dynamo fields are presented.

Several processes can cause the shape of an extrasolar giant planet's shadow, as viewed in transit, to depart from circular. In addition to rotational effects, cloud formation, non-homogenous haze production and movement, and dynamical effects (winds) could also be important. When such a planet transits its host star as seen from Earth, the asphericity will introduce a deviation in the transit lightcurve relative to the transit of a perfectly spherical (or perfectly oblate) planet. We develop a theoretical framework to interpret planet shapes. We then generate predictions for transiting planet shapes based on a published theoretical dynamical model of HD189733b. Using these shape models we show that planet shapes are unlikely to introduce detectable lightcurve deviations (those >~1e-5 of the host star), but that the shapes may lead to astrophysical sources of systematic error when measuring planetary oblateness, transit time, and impact parameter.

We present global, three-dimensional numerical simulations of the atmospheric circulation on HD 209458b and HD 189733b and calculate the infrared spectra and light curves predicted by these simulations, which we compare with available observations. Radiative heating/cooling is parameterized with a simplified Newtonian relaxation scheme. Our simulations develop day-night temperature contrasts that vary strongly with pressure. At low pressure (<10 mbar), air flows from the substellar point toward the antistellar point, both along the equator and over the poles. At deeper levels, the flow develops an eastward equatorial jet with speeds of 3-4 km/sec, with weaker westward flows at high latitudes. This basic flow pattern is robust to variations in model resolution, gravity, radiative time constant, and initial temperature structure. Nightside spectra show deep absorption bands of H2O, CO, and/or CH4, whereas on the dayside these absorption bands flatten out or even flip into emission. This results from the strong effect of dynamics on the vertical temperature-pressure structure; the temperature decreases strongly with altitude on the nightside but becomes almost isothermal on the dayside. In Spitzer bandpasses, our predicted planet-to-star flux ratios vary by a factor of 2-10 with orbital phase, depending on the wavelength and chemistry. For HD 189733b, where a detailed 8-micron light curve has been obtained, we correctly produce the observed phase offset of the flux maximum, but we do not explain the flux minimum and we overpredict the total flux variation. This discrepancy likely results from the simplifications inherent in the Newtonian relaxation scheme and provides motivation for incorporating realistic radiative transfer in future studies.

We present three-dimensional numerical simulations of the atmospheric dynamics of close-orbiting planets such as HD209458b. Our simulations show that winds of several km/sec and day-night temperature differences reaching 500-1000 K are possible at and above the photosphere. The circulation takes the form of a broad superrotating (eastward) equatorial jet. At 0.1-1 bar, the superrotation blows the hottest regions of the atmosphere downwind by 60 degrees of longitude, but at lower pressures the temperature pattern tracks the stellar illumination. We predict factors of several variation in the infrared flux received at Earth throughout an orbital cycle; if the photosphere is deep enough (>50-100 mbar pressure), the peak infrared emission should lead the time of secondary eclipse by 10 hours or more. Dynamics plays a key role in shaping the spectrum, clouds, chemistry, and long-term planetary evolution.

"Hot Jupiter" extrasolar planets are expected to be tidally locked because they are close (<0.05 astronomical units, where 1 AU is the average Sun-Earth distance) to their parent stars, resulting in permanent daysides and nightsides. By observing systems where the planet and star periodically eclipse each other, several groups have been able to estimate the temperatures of the daysides of these planets. A key question is whether the atmosphere is able to transport the energy incident upon the dayside to the nightside, which will determine the temperature at different points on the planet's surface. Here we report observations of HD 189733, the closest of these eclipsing planetary systems, over half an orbital period, from which we can construct a 'map' of the distribution of temperatures. We detected the increase in brightness as the dayside of the planet rotated into view. We estimate a minimum brightness temperature of 973 +/- 33 K and a maximum brightness temperature of 1212 +/- 11 K at a wavelength of 8 microns, indicating that energy from the irradiated dayside is efficiently redistributed throughout the atmosphere, in contrast to a recent claim for another hot Jupiter. Our data indicate that the peak hemisphere-integrated brightness occurs 16$\pm$6 degrees before opposition, corresponding to a hot spot shifted east of the substellar point. The secondary eclipse (when the planet moves behind the star) occurs 120 +/- 24 s later than predicted, which may indicate a slightly eccentric orbit.

We explore the infrared spectrum of a three-dimensional dynamical model of planet HD209458b as a function of orbital phase. The dynamical model predicts day-side atmospheric pressure-temperature profiles that are much more isothermal at pressures less than 1 bar than one-dimensional radiative-convective models have found. The resulting day-side thermal spectra are very similar to a blackbody, and only weak water absorption features are seen at short wavelengths. The dayside emission is consequently in significantly better agreement with ground-based and space-based secondary eclipse data than any previous models, which predict strong flux peaks and deep absorption features. At other orbital phases, absorption due to carbon monoxide and methane is also predicted. We compute the spectra under two treatments of atmospheric chemistry: one uses the predictions of equilibrium chemistry, and the other uses non-equilibrium chemistry, which ties the timescales of methane and carbon monoxide chemistry to dynamical timescales. As a function of orbital phase, we predict planet-to-star flux ratios for standard infrared bands and all Spitzer Space Telescope bands. In Spitzer bands, we predict 2-fold to 15-fold variation in planetary flux as a function of orbital phase with equilibrium chemistry, and 2-fold to 4-fold variation with non-equilibrium chemistry. Variation is generally more pronounced in bands from 3-10 $\mu$m than at longer wavelengths. The orbital phase of maximum thermal emission in infrared bands is 15--45 orbital degrees before the time of secondary eclipse. Changes in flux as a function of orbital phase for HD209458b should be observable with Spitzer, given the previously acheived observational error bars.

We present a method that employs the secondary eclipse light curves of transiting extrasolar planets to probe the spatial variation of their thermal emission. This technique permits an observer to resolve the surface of the planet without the need to spatially resolve its central star. We evaluate the feasibility of this technique for the HD 209458 system [..]. We consider two representations of the planetary thermal emission; a simple model parameterized by a sinusoidal dependence on longitude and latitude, as well as the results of a three-dimensional dynamical simulation of the planetary atmosphere previously published by Cooper & Showman. We find that observations of the secondary eclipse light curve are most sensitive to a longitudinal offset in the geometric and photometric centroids of the hemisphere of the planet visible near opposition. To quantify this signal, we define a new parameter, the ``uniform time offset,'' which measures the time lag between the observed secondary eclipse and that predicted by a planet with a uniform surface flux distribution. We compare the predicted amplitude of this parameter for HD 209458 with the precision with which it could be measured with IRAC. We find that IRAC observations at 3.6um a single secondary eclipse should permit sufficient precision to confirm or reject the Cooper & Showman model of the surface flux distribution for this planet. We quantify the signal-to-noise ratio for this offset in the remaining IRAC bands (4.5um, 5.8um, and 8.0um), and find that a modest improvement in photometric precision (as might be realized through observations of several eclipse events) should permit a similarly robust detection.

Chemical equilibrium considerations suggest that, assuming solar elemental abundances, carbon on HD 209458b is sequestered primarily as carbon monoxide (CO) and methane (CH4). The relative mole fractions of CO(g) and CH4(g) in chemical equilibrium are expected to vary greatly according to variations in local temperature and pressure. We show, however, that in the p = 1--1000 mbar range, chemical equilibrium does not hold. To explore disequilibrium effects, we couple the chemical kinetics of CO and CH4 to a three-dimensional numerical model of HD 209458b's atmospheric circulation. These simulations show that vigorous dynamics caused by uneven heating of this tidally locked planet homogenize the CO and CH4 concentrations at p < 1 bar, even in the presence of lateral temperature variations of ~500--1000 K. In the 1--1000 mbar pressure range, we find that over 98% of the carbon is in CO. This is true even in cool regions where CH4 is much more stable thermodynamically. Our work shows furthermore that planets 300--500 K cooler than HD 209458b can also have abundant CO in their upper layers due to disequilibrium effects. We demonstrate several interesting observational consequences of these results.

We calculate the meteorology of the close-in transiting extrasolar planet HD 209458b using a global, three-dimensional atmospheric circulation model. Dynamics are driven by perpetual irradiation of one hemisphere of this tidally locked planet. The simulation predicts global temperature contrasts of ~500 K at the photosphere and the development of a steady superrotating jet. The jet extends from the equator to mid-latitudes and from the top model layer at 1 mbar down to 10 bars at the base of the heated region. Wind velocities near the equator exceed 4 km/s at 300 mbar. The hottest regions of the atmosphere are blown downstream from the substellar point by 60 degrees of longitude. We predict from these results a factor of ~2 ratio between the maximum and minimum observed radiation from the planet over a full orbital period, with peak infrared emission preceding the time of the secondary eclipse by ~14 hours.

We have posted the source code for our cloud model for public use as a tool for the intercomparison of planetary radiation transport models attempting to incorporate the physics of cloud condensation.

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

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

The variability observed in many complex gamma-ray bursts (GRBs) is inconsistent with causally connected variations in a single, symmetric, relativistic shell interacting with the ambient material ("external shocks"). Rather, the symmetry of the shell must be broken on an angular scale much smaller than Gamma^{-1} where Gamma is the bulk Lorentz factor for the shell. The observed variability in the external shock models arises from the number of causally connected regions that (randomly) become active. We define the "surface filling factor" to be the ratio of the area of causally connected regions that become active to the observable area of the shell. From the observed variability in 52 BATSE bursts, we estimate the surface filling factor to be typically 0.005 although some values are near unity. We find that the surface filling factor is about 0.1 Delta T/T in both the constant Gamma phase (which probably produces the GRB) and the deaccelerating phase (which probably produces the x-ray afterglows). Here, \Delta T is a typical time scale of variability and T is the time since the initial signal. We analyze the 2 hr flare seen by ASCA 36 hr after the GRB and conclude that the surface filling factor must be small (0.001) in the x-ray afterglow phase as well. Explanations for low surface filling factor can either require more or less energy (by a factor of about 1000) compared to that expected for a symmetric shell.