Context - Exoplanetary upper atmospheres are low density environments where radiative processes can compete with collisional ones and introduce non-local thermodynamic equilibrium (NLTE) effects into transmission spectra. Aims - We develop a NLTE radiative transfer framework capable of modelling exoplanetary transmission spectra over a wide range of planetary properties. Methods - We adapt the NLTE spectral synthesis code Cloudy to produce an atmospheric structure and atomic transmission spectrum in both NLTE and local thermodynamic equilibrium (LTE) for the hot Jupiter HD209458b, given a published T-P profile and assuming solar metallicity. Selected spectral features, including H$\alpha$, Na I D, He I $\lambda$10830, Fe I & II ultra-violet (UV) bands, and C, O and Si UV lines, are compared with literature observations and models where available. The strength of NLTE effects are measured for individual spectral lines to identify which features are most strongly affected. Results - The developed modelling framework computing NLTE synthetic spectra reproduces literature results for the He I $\lambda$10830 triplet, the Na I D lines, and the forest of Fe I lines in the optical. Individual spectral lines in the NLTE spectrum exhibit up to 40 % stronger absorption relative to the LTE spectrum.

We present two prescriptions for broadband (~77 - 252 GHz), millimeter-wave antireflection coatings for cryogenic, sintered polycrystalline aluminum oxide optics: one for large-format (700 mm diameter) planar and plano-convex elements, the other for densely packed arrays of quasi-optical elements, in our case 5 mm diameter half-spheres (called "lenslets"). The coatings comprise three layers of commercially-available, polytetrafluoroethylene-based, dielectric sheet material. The lenslet coating is molded to fit the 150 mm diameter arrays directly while the large-diameter lenses are coated using a tiled approach. We review the fabrication processes for both prescriptions then discuss laboratory measurements of their transmittance and reflectance. In addition, we present the inferred refractive indices and loss tangents for the coating materials and the aluminum oxide substrate. We find that at 150 GHz and 300 K the large-format coating sample achieves (97 +/- 2)% transmittance and the lenslet coating sample achieves (94 +/- 3)% transmittance.

The inflated transiting hot Jupiter HD 209458b is one of the best studied objects since the beginning of exoplanet characterization. Transmission observations of this system between the mid infrared and the far ultraviolet have revealed the signature of atomic, molecular, and possibly aerosol species in the lower atmosphere of the planet, as well as escaping hydrogen and metals in the upper atmosphere. From a re-analysis of near-ultraviolet (NUV) transmission observations of HD 209458b, we detect ionized iron (Fe II) absorption in a 100 A-wide range around 2370 A, lying beyond the planetary Roche lobe. However, we do not detect absorption of equally strong Fe II lines expected to be around 2600 A. Further, we find no evidence for absorption by neutral magnesium (Mg I), ionized magnesium (Mg II), nor neutral iron (Fe I). These results avoid the conflict with theoretical models previously found by Vidal-Madjar et al. (2013), which detected Mg I but did not detect Mg II from this same data set. Our results indicate that hydrodynamic escape is strong enough to carry atoms as heavy as iron beyond the planetary Roche lobe, even for planets less irradiated than the extreme ultra-hot-Jupiters such as WASP-12b and KELT-9b. The detection of iron and non-detection of magnesium in the upper atmosphere of HD 209458b can be explained by a model in which the lower atmosphere forms (hence, sequesters) primarily magnesium-bearing condensates, rather than iron condensates. This is suggested by current microphysical models. The inextricable synergy between upper- and lower-atmosphere properties highlights the value of combining observations that probe both regions.

We conduct an in-depth analysis of statistical flow properties calculated from the reference high-resolution Saturn simulation obtained by global climate modelling in Part II. In the steady state of this reference simulation, strongly energetic, zonally dominated, large-scale structures emerge, which scale with the Rhines scale. Spectral analysis reveals a strong anisotropy in the kinetic energy spectra, consistent with the zonostrophic turbulent flow regime. By computing spectral energy and enstrophy fluxes we confirm the existence of a double cascade scenario related to 2D-turbulent theory. To diagnose the relevant 3D dynamical mechanisms in Saturn's turbulent atmosphere, we run a set of four simulations using an idealized version of our Global Climate Model devoid of radiative transfer, with a well-defined Taylor-Green forcing and over several rotation rates (4, 1, 0.5, and 0.25 times Saturn's rotation rate). This allows us to identify dynamics in three distinctive inertial ranges: (1) a ``residual-dominated'' range, in which non-axisymmetric structures dominate with a -5/3 spectral slope; (2) a ``zonostrophic inertial'' range, dominated by axisymmetric jets and characterized by the pile-up of strong zonal modes with a steeper, nearly -3, spectral slope; and (3) a ``large-scale'' range, beyond Rhines' typical length scale, in which the reference Saturn simulation and our idealized simulations differ. In the latter range, the dynamics is dominated by long-lived zonal modes 2 and 3 when a Saturn-like seasonal forcing is considered (reference simulation), and a steep energetic decrease with the idealized Taylor-Green forcing. Finally, instantaneous spectral fluxes show the coexistence of upscale and downscale enstrophy/energy transfers at large scales, specific to the regime of zonostrophic turbulence in a 3D atmosphere.

With a day-side temperature in excess of 4500K, comparable to a mid-K-type star, KELT-9b is the hottest planet known. Its extreme temperature makes KELT-9b a particularly interesting test bed for investigating the nature and diversity of gas giant planets. We observed the transit of KELT-9b at high spectral resolution (R$\sim$94,600) with the CARMENES instrument on the Calar Alto 3.5-m telescope. Using these data, we detect for the first time ionized calcium (CaII triplet) absorption in the atmosphere of KELT-9b; this is the second time that CaII has been observed in a hot Jupiter. Our observations also reveal prominent H$\alpha$ absorption, confirming the presence of an extended hydrogen envelope around KELT-9b. We compare our detections with an atmospheric model and find that all four lines form between atmospheric temperatures of 6100 K and 8000 K and that the CaII lines form at pressures between 10 and 50 nbar while the H$\alpha$ line forms at a lower pressure ($\sim$6 nbar), higher up in the atmosphere. The altitude that the core of H$\alpha$ line forms is found to be $\sim$1.4 R$_{p}$, well within the planetary Roche lobe ($\sim$1.9 R$_{p}$). Therefore, rather than probing the escaping upper atmosphere directly, the H$\alpha$ line and the other observed Balmer and metal lines serve as atmospheric thermometers enabling us to probe the planet's temperature profile, thus energy budget.

The cosmic microwave background (CMB) encodes information about the content and evolution of the universe. The presence of light, weakly interacting particles impacts the expansion history of the early universe, which alters the temperature and polarization anisotropies of the CMB. In this way, current measurements of the CMB place interesting constraints on the neutrino energy density and mass, as well as on the abundance of other possible light relativistic particle species. We present the status of an on-going 1500 sq. deg. survey with the SPT-3G receiver, a new mm-wavelength camera on the 10-m diameter South Pole Telescope (SPT). The SPT-3G camera consists of 16,000 superconducting transition edge sensors, a 10x increase over the previous generation camera, which allows it to map the CMB with an unprecedented combination of sensitivity and angular resolution. We highlight projected constraints on the abundance of sterile neutrinos and the sum of the neutrino masses for the SPT-3G survey, which could help determine the neutrino mass hierarchy.

SPT-3G is a polarization-sensitive receiver, installed on the South Pole Telescope, that measures the anisotropy of the cosmic microwave background (CMB) from degree to arcminute scales. The receiver consists of ten 150~mm-diameter detector wafers, containing a total of 16,000 transition-edge sensor (TES) bolometers observing at 95, 150, and 220 GHz. During the 2018-2019 austral summer, one of these detector wafers was replaced by a new wafer fabricated with Al-Mn TESs instead of the Ti/Au design originally deployed for SPT-3G. We present the results of in-lab characterization and on-sky performance of this Al-Mn wafer, including electrical and thermal properties, optical efficiency measurements, and noise-equivalent temperature. In addition, we discuss and account for several calibration-related systematic errors that affect measurements made using frequency-domain multiplexing readout electronics.

Frequency-domain multiplexing (fMux) is an established technique for the readout of large arrays of transition edge sensor (TES) bolometers. Each TES in a multiplexing module has a unique AC voltage bias that is selected by a resonant filter. This scheme enables the operation and readout of multiple bolometers on a single pair of wires, reducing thermal loading onto sub-Kelvin stages. The current receiver on the South Pole Telescope, SPT-3G, uses a 68x fMux system to operate its large-format camera of $\sim$16,000 TES bolometers. We present here the successful implementation and performance of the SPT-3G readout as measured on-sky. Characterization of the noise reveals a median pair-differenced 1/f knee frequency of 33 mHz, indicating that low-frequency noise in the readout will not limit SPT-3G's measurements of sky power on large angular scales. Measurements also show that the median readout white noise level in each of the SPT-3G observing bands is below the expectation for photon noise, demonstrating that SPT-3G is operating in the photon-noise-dominated regime.

Many previous authors have attempted to find explanations for Westeros's climate, characterized by a generally moderate, Earth-like climate punctuated by extremely long and cold winters, separated by thousands of years. One explanation that has been proposed is that the planet orbits in a Sitnikov configuration, where two equal-mass stars (or a star and a black hole) orbit each other on slightly eccentric orbits, and the planet moves along a line through the barycenter perpendicular to the primaries' orbital plane (Freistetter & Gr\"utzbauch 2018). We modify an intermediate-complexity GCM to include the effects of such an orbit and integrate it for thousands of years to determine whether such an orbit can a) be habitable and b) explain the climatic variations observed by the inhabitants of Westeros, in both double-star and star-black hole configurations. While configurations with low primary eccentricity and initial conditions that permit only small excursions from the ecliptic plane are habitable, these orbits are too stable to explain Westerosi climate. We find that while orbits with more bounded chaos are able to produce rare anomalously long and cold winters similar to Westeros's Long Night, huge variations in incident stellar flux on normal orbital timescales should render these planets uninhabitable. We note that the presence of an orbital megastructure, either around the planet or the barycenter, could block some of the sunlight during crossings of the primaries' orbital plane and preserve Westeros's habitability. While we find that bounded chaotic Sitnikov orbits are a viable explanation for Westeros's Long Night, we propose that chaotic variations of the planet's axial tilt or semimajor axis, potentially due to torques from nearby planets or stars, may be a more realistic explanation than Sitnikov orbits.

During the austral summer of 2016-17, the third-generation camera, SPT-3G, was installed on the South Pole Telescope, increasing the detector count in the focal plane by an order of magnitude relative to the previous generation. Designed to map the polarization of the cosmic microwave background, SPT-3G contains ten 6-in-hexagonal modules of detectors, each with 269 trichroic and dual-polarization pixels, read out using 68x frequency-domain multiplexing. Here we discuss design, assembly, and layout of the modules, as well as early performance characterization of the first-year array, including yield and detector properties.

Extremely irradiated, close-in planets to early-type stars might be prone to strong atmospheric escape. We review the literature showing that X-ray-to-optical measurements indicate that for intermediate-mass stars (IMS) cooler than $\approx$8250 K, the X-ray and EUV (XUV) fluxes are on average significantly higher than those of solar-like stars, while for hotter IMS, because of the lack of surface convection, it is the opposite. We construct spectral energy distributions for prototypical IMS, comparing them to solar. The XUV fluxes relevant for upper planet atmospheric heating are highest for the cooler IMS and lowest for the hotter IMS, while the UV fluxes increase with increasing stellar temperature. We quantify the influence of this characteristic of the stellar fluxes on the mass loss of close-in planets by simulating the atmospheres of planets orbiting EUV-bright (WASP-33) and EUV-faint (KELT-9) A-type stars. For KELT-9b, we find that atmospheric expansion caused by heating due to absorption of the stellar UV and optical light drives mass-loss rates of $\approx$10$^{11}$ g s$^{-1}$, while heating caused by absorption of the stellar XUV radiation leads to mass-loss rates of $\approx$10$^{10}$ g s$^{-1}$, thus underestimating mass loss. For WASP-33b, the high XUV stellar fluxes lead to mass-loss rates of $\approx$10$^{11}$ g s$^{-1}$. Even higher mass-loss rates are possible for less massive planets orbiting EUV-bright IMS. We argue that it is the weak XUV stellar emission, combined with a relatively high planetary mass, which limit planetary mass-loss rates, to allow the prolonged existence of KELT-9-like systems.

The SPT-3G receiver was commissioned in early 2017 on the 10-meter South Pole Telescope (SPT) to map anisotropies in the cosmic microwave background (CMB). New optics, detector, and readout technologies have yielded a multichroic, high-resolution, low-noise camera with impressive throughput and sensitivity, offering the potential to improve our understanding of inflationary physics, astroparticle physics, and growth of structure. We highlight several key features and design principles of the new receiver, and summarize its performance to date.

The third-generation instrument for the 10-meter South Pole Telescope, SPT-3G, was first installed in January 2017. In addition to completely new cryostats, secondary telescope optics, and readout electronics, the number of detectors in the focal plane has increased by an order of magnitude from previous instruments to ~16,000. The SPT-3G focal plane consists of ten detector modules, each with an array of 269 trichroic, polarization-sensitive pixels on a six-inch silicon wafer. Within each pixel is a broadband, dual-polarization sinuous antenna; the signal from each orthogonal linear polarization is divided into three frequency bands centered at 95, 150, and 220 GHz by in-line lumped element filters and transmitted via superconducting microstrip to Ti/Au transition-edge sensor (TES) bolometers. Properties of the TES film, microstrip filters, and bolometer island must be tightly controlled to achieve optimal performance. For the second year of SPT-3G operation, we have replaced all ten wafers in the focal plane with new detector arrays tuned to increase mapping speed and improve overall performance. Here we discuss the TES superconducting transition temperature and normal resistance, detector saturation power, bandpasses, optical efficiency, and full array yield for the 2018 focal plane.

The South Pole Telescope (SPT) is a millimeter-wavelength telescope designed for high-precision measurements of the cosmic microwave background (CMB). The SPT measures both the temperature and polarization of the CMB with a large aperture, resulting in high resolution maps sensitive to signals across a wide range of angular scales on the sky. With these data, the SPT has the potential to make a broad range of cosmological measurements. These include constraining the effect of massive neutrinos on large-scale structure formation as well as cleaning galactic and cosmological foregrounds from CMB polarization data in future searches for inflationary gravitational waves. The SPT began observing in January 2017 with a new receiver (SPT-3G) containing $\sim$16,000 polarization-sensitive transition-edge sensor bolometers. Several key technology developments have enabled this large-format focal plane, including advances in detectors, readout electronics, and large millimeter-wavelength optics. We discuss the implementation of these technologies in the SPT-3G receiver as well as the challenges they presented. In late 2017 the implementations of all three of these technologies were modified to optimize total performance. Here, we present the current instrument status of the SPT-3G receiver.

The desire for higher sensitivity has driven ground-based cosmic microwave background (CMB) experiments to employ ever larger focal planes, which in turn require larger reimaging optics. Practical limits to the maximum size of these optics motivates the development of quasi-optically-coupled (lenslet-coupled), multi-chroic detectors. These detectors can be sensitive across a broader bandwidth compared to waveguide-coupled detectors. However, the increase in bandwidth comes at a cost: the lenses (up to $\sim$700 mm diameter) and lenslets ($\sim$5 mm diameter, hemispherical lenses on the focal plane) used in these systems are made from high-refractive-index materials (such as silicon or amorphous aluminum oxide) that reflect nearly a third of the incident radiation. In order to maximize the faint CMB signal that reaches the detectors, the lenses and lenslets must be coated with an anti-reflective (AR) material. The AR coating must maximize radiation transmission in scientifically interesting bands and be cryogenically stable. Such a coating was developed for the third generation camera, SPT-3G, of the South Pole Telescope (SPT) experiment, but the materials and techniques used in the development are general to AR coatings for mm-wave optics. The three-layer polytetrafluoroethylene-based AR coating is broadband, inexpensive, and can be manufactured with simple tools. The coating is field tested; AR coated focal plane elements were deployed in the 2016-2017 austral summer and AR coated reimaging optics were deployed in 2017-2018.

The third-generation South Pole Telescope camera is designed to measure the cosmic microwave background across three frequency bands (95, 150 and 220 GHz) with ~16,000 transition-edge sensor (TES) bolometers. Each multichroic pixel on a detector wafer has a broadband sinuous antenna that couples power to six TESs, one for each of the three observing bands and both polarization directions, via lumped element filters. Ten detector wafers populate the focal plane, which is coupled to the sky via a large-aperture optical system. Here we present the frequency band characterization with Fourier transform spectroscopy, measurements of optical time constants, beam properties, and optical and polarization efficiencies of the focal plane. The detectors have frequency bands consistent with our simulations, and have high average optical efficiency which is 86%, 77% and 66% for the 95, 150 and 220 GHz detectors. The time constants of the detectors are mostly between 0.5 ms and 5 ms. The beam is round with the correct size, and the polarization efficiency is more than 90% for most of the bolometers

We present results from WIYN pODI imaging of Lacerta I (And XXXI), a satellite dwarf galaxy discovered in the outskirts of the Andromeda galaxy (M31) in Pan-STARRS1 survey data. Our deep, wide-field $g,i$ photometry reaches $\sim$3 magnitudes fainter than the photometry in the Pan-STARRS1 discovery paper and allows us to trace the stellar population of Lac I beyond two half-light radii from the galaxy center. We measure a Tip of the Red Giant Branch (TRGB) distance for Lac I of $(m-M)_0=24.44\pm0.11$ mag (773$\pm$40 kpc, or 264$\pm$6 kpc from M31), which is consistent with the Pan-STARRS1 distance. We use a maximum-likelihood technique to derive structural properties for the galaxy, and find a half-light radius ($r_h$) of 3.24$\pm$0.21 arcmin (728$\pm$47 pc), ellipticity ($\epsilon$) of 0.44$\pm$0.03, total magnitude $M_V$ $=$ $-$11.4$\pm$0.3, and central surface brightness $\mu_{V,0}$ $=$ 24.8$\pm$0.3 mag arcsec$^{-2}$. We find no HI emission in archival data and set a limit on Lac I's neutral gas mass-to-light ratio of $M_{HI}/L_V$ $<$ 0.06 $M_{sun}$/$L_{sun}$, confirming Lac I as a gas-poor dwarf spheroidal galaxy. Photometric metallicities derived from Red Giant Branch stars within 2 $r_h$ yield a median [Fe/H] of $-$1.68$\pm$0.03, which is more metal-rich than the spectroscopically-derived value from Martin et al. (2014). Combining our measured magnitude with this higher metallicity estimate places Lac I closer to its expected position on the luminosity-metallicity relation for dwarf galaxies.

We present an investigation of the globular cluster population synthesis method of McWilliam & Bernstein (2008), focusing on the impact of NLTE modeling effects and CMD discretization. Johnson-Cousins-Bessel U-B, B-V, V-I, and J-K colors are produced for 96 synthetic integrated light spectra with two different discretization prescriptions and three degrees of NLTE treatment. These color values are used to compare NLTE and LTE derived population ages. Relative contributions of different spectral types to the integrated light spectra for different wavebands are measured. Integrated light NLTE spectra are shown to be more luminous in the UV and optical than LTE spectra, but show stronger absorption features in the IR. The main features showing discrepancies between NLTE and LTE integrated light spectra may be attributed to light metals, primarily Fe I, Ca I, and Ti I, as well as TiO molecular bands. Main Sequence stars are shown to have negligible NLTE effects at IR wavelengths compared to more evolved stars. Photometric color values are shown to vary at the millimagnitude level as a function of CMD discretization. Finer CMD sampling for the upper main sequence and turnoff, base of the red giant branch, and the horizontal branch minimizes this variation. Differences in ages derived from LTE and NLTE IL spectra are found to range from 0.55 to 2.54 Gyr, comparable to the uncertainty in GC ages derived from color indices with observational uncertainties of 0.01 magnitudes, the limiting precision of the Harris catalog (Harris 1996).

We describe the SpArc science gateway for spectral data obtained during the period from 1975 through 1995 at the Kitt Peak National Observatory using the Fourier Transform Spectrometer (FTS) in operation at the Mayall 4-m telescope. SpArc is hosted by Indiana University Bloomington and is available for public access. The archive includes nearly 10,000 individual spectra of more than 800 different astronomical sources including stars, nebulae, galaxies, and Solar System objects. We briefly describe the FTS instrument itself, and summarize the conversion of the original interferograms into spectral data and the process for recovering the data into FITS files. The architecture of the archive is discussed, and the process for retrieving data from the archive is introduced. Sample use cases showing typical FTS spectra are presented.

Using 2D smoothed particle hydrodynamics, we investigate the distribution of wait times between strong shocks in a turbulent, self-gravitating accretion disc. We show the resulting distributions do not depend strongly on the cooling time or resolution of the disc and that they are consistent with the predictions of earlier work (Young & Clarke 2015; Cossins et al. 2009, 2010). We use the distribution of wait times between shocks to estimate the likelihood of stochastic fragmentation by gradual contraction of shear-resistant clumps on the cooling time scale. We conclude that the stochastic fragmentation mechanism (Paardekooper 2012) cannot change the radius at which fragmentation is possible by more than ~20%, restricting direct gravitational collapse as a mechanism for giant planet formation to the outer regions of protoplanetary discs.

We investigate the dependence of the complete system of 22 Lick indices on overall metallicity scaled from solar abundances, [M/H], from the solar value, 0.0, down to the extremely-metal-poor (XMP) value of -6.0, for late-type giant stars (MK luminosity class III, log(g)=2.0) of MK spectral class late-K to late-F (3750 < Teff < 6500 K) of the type that are detected as "fossils" of early galaxy formation in the Galactic halo and in extra-galactic structures. Our investigation is based on synthetic index values, I, derived from atmospheric models and synthetic spectra computed with PHOENIX in LTE and Non-LTE (NLTE), where the synthetic spectra have been convolved to the spectral resolution, R, of both IDS and SDSS (and LAMOST) spectroscopy. We identify nine indices, that we designate "Lick-XMP", that remain both detectable and significantly [M/H]-dependent down to [M/H] values of at least ~-5.0, and down to [M/H] ~ -6.0 in five cases, while also remaining well-behaved . For these nine, we study the dependence of I on NLTE effects, and on spectral resolution. For our LTE I values for spectra of SDSS resolution, we present the fitted polynomial coefficients, C_n, from multi-variate linear regression for I with terms up to third order in the independent variable pairs (Teff, [M/H]), and (V-K, [M/H]), and compare them to the fitted C_n values of Worthey et al. (1994) at IDS spectral resolution.

We perform a series of 2D smoothed particle hydrodynamics (SPH) simulations of gas accretion onto binaries via a circumbinary disc, for a range of gas temperatures and binary mass ratios ($q$). We show that increasing the gas temperature increases the accretion rate onto the primary for all values of the binary mass ratio: for example, for $q=0.1$ and a fixed binary separation, an increase of normalised sound speed by a factor of $5$ (from our "cold" to "hot" simulations) changes the fraction of the accreted gas that flows on to the primary from $ 10\%$ to $\sim40\%$. We present a simple parametrisation for the average accretion rate of each binary component accurate to within a few percent and argue that this parametrisation (rather than those in the literature based on warmer simulations) is relevant to supermassive black hole accretion and all but the widest stellar binaries. We present trajectories for the growth of $q$ during circumbinary disc accretion and argue that the period distribution of stellar "twin" binaries is strong evidence for the importance of circumbinary accretion. We also show that our parametrisation of binary accretion increases the minimum mass ratio needed for spin alignment of supermassive black holes to $q \sim 0.4$, with potentially important implications for the magnitude of velocity kicks acquired during black hole mergers.

We propose a framework for understanding the fragmentation criterion for self-gravitating discs which, in contrast to studies that emphasise the `gravoturbulent' nature of such discs, instead focuses on the properties of their quasi-regular spiral structures. Within this framework there are two evolutionary paths to fragmentation: i) collapse on the free-fall time, which requires that the ratio of cooling time to dynamical time ($\beta$) $< 3$ and ii) quasistatic collapse on the cooling time at a rate that is sufficiently fast that fragments are compact enough to withstand disruption when they encounter spiral features in the disc. We perform 2D grid simulations which demonstrate numerically converged fragmentation at $\beta < 3$ (in good agreement with Paardekooper et al. (2011) and others) and argue that this is a consequence of the fact that such simulations smooth the gravitational force on the scale $H$, the scale height of the disc. Such simulations thus only allow fragmentation via route i) above since they suppress the quasistatic contraction of fragments on scales $< H$; the inability of fragments to contract to significantly smaller scales then renders them susceptible to disruption at the next spiral arm encounter. On the other hand, 3D simulations indeed show fragmentation at higher $\beta$ via route ii). We derive an analytic prediction of fragmentation by route ii) when $\beta \lesssim 12$, based on the requirement that fragments must contract sufficiently to withstand disruption by spiral arms. We also discuss the necessary numerical requirements on both grid based and SPH codes if they are to model fragmentation via route ii).

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

We explore an unresolved controversy in the literature about the accuracy of Smoothed Particle Hydrodynamics (SPH) in modeling the accretion of gas onto a binary system, a problem with important applications to the evolution of proto-binaries as well as accreting binary super massive black holes. It has previously been suggested that SPH fails to model the flow of loosely bound material from the secondary to primary Roche lobe and that its general prediction that accretion drives mass ratios upwards is numerically flawed. Here we show with 2D SPH that this flow from secondary to primary Roche lobe is a sensitive function of gas temperature and that this largely explains the conflicting claims in the literature which have hitherto been based on either 'cold' SPH simulations or 'hot' grid based calculations. We present simulations of a specimen 'cold' and 'hot' accretion scenario which are numerically converged and evolved into a steady state. Our analysis of the conservation of the Jacobi integral of accreting particles indicates that our results are not strongly compromised by numerical dissipation. We also explore the low resolution limit and find that simulations where the ratio of SPH smoothing length to disc scale height at the edge of the circumsecondary is less than 1 accurately capture binary accretion rates.

Spectra for 2D stars in the 1.5D approximation are created from synthetic spectra of 1D non-local thermodynamic equilibrium (NLTE) spherical model atmospheres produced by the PHOENIX code. The 1.5D stars have the spatially averaged Rayleigh-Jeans flux of a K3-4 III star, while varying the temperature difference between the two 1D component models ($\Delta T_{\mathrm{1.5D}}$), and the relative surface area covered. Synthetic observable quantities from the 1.5D stars are fitted with quantities from NLTE and local thermodynamic equilibrium (LTE) 1D models to assess the errors in inferred $T_{\mathrm{eff}}$ values from assuming horizontal homogeneity and LTE. Five different quantities are fit to determine the $T_{\mathrm{eff}}$ of the 1.5D stars: UBVRI photometric colors, absolute surface flux SEDs, relative SEDs, continuum normalized spectra, and TiO band profiles. In all cases except the TiO band profiles, the inferred $T_{\mathrm{eff}}$ value increases with increasing $\Delta T_{\mathrm{1.5D}}$. In all cases, the inferred $T_{\mathrm{eff}}$ value from fitting 1D LTE quantities is higher than from fitting 1D NLTE quantities and is approximately constant as a function of $\Delta T_{\mathrm{1.5D}}$ within each case. The difference between LTE and NLTE for the TiO bands is caused indirectly by the NLTE temperature structure of the upper atmosphere, as the bands are computed in LTE. We conclude that the difference between $T_{\mathrm{eff}}$ values derived from NLTE and LTE modelling is relatively insensitive to the degree of the horizontal inhomogeneity of the star being modeled, and largely depends on the observable quantity being fit.

We present results from ground-based optical imaging of a low-mass dwarf galaxy discovered by the ALFALFA 21-cm HI survey. Broadband (BVR) data obtained with the WIYN 3.5-m telescope at Kitt Peak National Observatory (KPNO) are used to construct color-magnitude diagrams of the galaxy's stellar population down to V_0 ~ 25. We also use narrowband H-alpha imaging from the KPNO 2.1-m telescope to identify an HII region in the galaxy. We use these data to constrain the distance to the galaxy to be between 1.5 and 2.0 Mpc. This places Leo P within the Local Volume but beyond the Local Group. Its properties are extreme: it is the lowest-mass system known that contains significant amounts of gas and is currently forming stars.

The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.

We present 11 candidate late-type companions to nearby stars identified with data from the Wide-field Infrared Survey Explorer (WISE) and the Two Micron All-Sky Survey (2MASS). Eight of the candidates are likely to be companions based on their common proper motions with the primaries. The remaining three objects are rejected as companions, one of which is a free-floating T7 dwarf. Spectral types are available for five of the companions, which consist of M2V, M8.5V, L5, T8, and T8. Based on their photometry, the unclassified companions are probably two mid-M dwarfs and one late-M/early-L dwarf. One of the T8 companions, WISE J142320.84+011638.0, has already been reported by Pinfield and coworkers. The other T8 companion, ULAS J095047.28+011734.3, was discovered by Burningham and coworkers through the United Kingdom Infrared Telescope Infrared Deep Sky Survey, but its companionship has not been previously recognized in the literature. The L5 companion, 2MASS J17430860+8526594, is a new member of a class of L dwarfs that exhibit unusually blue near-IR colors. Among the possible mechanisms that have been previously proposed for the peculiar colors of these L dwarfs, low metallicity does not appear to be a viable explanation for 2MASS J17430860+8526594 since our spectrum of the primary suggests that its metallicity is not significantly subsolar.

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.

We present results from a study of the globular cluster (GC) systems of four spiral and S0 galaxies imaged as part of an ongoing wide-field survey of the GC systems of giant galaxies. The target galaxies -- the SB0 galaxy NGC1023, the SBb galaxy NGC1055, and an isolated pair comprised of the Sbc galaxy NGC7339 and the S0 galaxy NGC7332 -- were observed in BVR filters with the WIYN 3.5-m telescope and Minimosaic camera. For two of the galaxies, we combined the WIYN imaging with previously-published data from the Hubble Space Telescope and the Keck Observatory to help characterize the GC distribution in the central few kiloparsecs. We determine the radial distribution (surface density of GCs versus projected radius) of each galaxy's GC system and use it to calculate the total number of GCs (N_GC). We find N_GC = 490+/-30, 210+/-40, 175+/-15, and 75+/-10 for NGC1023, NGC1055, NGC7332, and NGC7339, respectively. We also calculate the GC specific frequency (N_GC normalized by host galaxy luminosity or mass) and find values typical of those of the other spiral and E/S0 galaxies in the survey. The two lenticular galaxies have sufficient numbers of GC candidates for us to perform statistical tests for bimodality in the GC color distributions. We find evidence at a high confidence level (>95%) for two populations in the B-R distribution of the GC system of NGC1023. We find weaker evidence for bimodality (>81% confidence) in the GC color distribution of NGC7332. Finally, we identify eight GC candidates that may be associated with the Magellanic dwarf galaxy NGC1023A, a satellite of NGC1023.

Using the 4 Ms Chandra Deep Field-South (CDF-S) survey, we have identified a sample of 6845 X-ray undetected galaxies that dominates the unresolved ~ 20-25% of the 6-8 keV cosmic X-ray background (XRB). This sample was constructed by applying mass and color cuts to sources from a parent catalog based on GOODS-South HST z-band imaging of the central 6'-radius area of the 4 Ms CDF-S. The stacked 6-8 keV detection is significant at the 3.9 sigma level, but the stacked emission was not detected in the 4-6 keV band which indicates the existence of an underlying population of highly obscured active galactic nuclei (AGNs). Further examinations of these 6845 galaxies indicate that the galaxies on the top of the blue cloud and with redshifts of 1 < z < 3, magnitudes of 25 < z_850 < 28, and stellar masses of 2E8 < M_star/M_sun < 2E9 make the majority contributions to the unresolved 6-8 keV XRB. Such a population is seemingly surprising given that the majority of the X-ray detected AGNs reside in massive (> ~1E10 M_sun) galaxies. We discuss constraints upon this underlying AGN population, supporting evidence for relatively low-mass galaxies hosting highly obscured AGNs, and prospects for further boosting the stacked signal.

We describe a simple method for estimating the vertical column density in Smoothed Particle Hydrodynamics (SPH) simulations of discs. As in the method of Stamatellos et al. (2007), the column density is estimated using pre-computed local quantities and is then used to estimate the radiative cooling rate. The cooling rate is a quantity of considerable importance, for example, in assessing the probability of disc fragmentation. Our method has three steps: (i) the column density from the particle to the mid plane is estimated using the vertical component of the gravitational acceleration, (ii) the "total surface density" from the mid plane to the surface of the disc is calculated, (iii) the column density from each particle to the surface is calculated from the difference between (i) and (ii). This method is shown to greatly improve the accuracy of column density estimates in disc geometry compared with the method of Stamatellos. On the other hand, although the accuracy of our method is still acceptable in the case of high density fragments formed within discs, we find that the Stamatellos method performs better than our method in this regime. Thus, a hybrid method (where the method is switched in regions of large over-density) may be optimal.

We present near infra-red light curves of supernova (SN) 2011fe in M101, including 34 epochs in H band starting fourteen days before maximum brightness in the B-band. The light curve data were obtained with the WIYN High-Resolution Infrared Camera (WHIRC). When the data are calibrated using templates of other Type Ia SNe, we derive an apparent H-band magnitude at the epoch of B-band maximum of 10.85 \pm 0.04. This implies a distance modulus for M101 that ranges from 28.86 to 29.17 mag, depending on which absolute calibration for Type Ia SNe is used.

The 4 Ms Chandra Deep Field-South (CDF-S) and other deep X-ray surveys have been highly effective at selecting active galactic nuclei (AGN). However, cosmologically distant low-luminosity AGN (LLAGN) have remained a challenge to identify due to significant contribution from the host galaxy. We identify long-term X-ray variability (~month-years, observed frame) in 20 of 92 CDF-S galaxies spanning redshifts z~0.08-1.02 that do not meet other AGN selection criteria. We show that the observed variability cannot be explained by X-ray binary populations or ultraluminous X-ray sources, so the variability is most likely caused by accretion onto a supermassive black hole. The variable galaxies are not heavily obscured in general, with a stacked effective power-law photon index of Gamma_stack~1.93+/-0.13, and are therefore likely LLAGN. The LLAGN tend to lie a factor of ~6-80 below the extrapolated linear variability-luminosity relation measured for luminous AGN. This may be explained by their lower accretion rates. Variability-independent black-hole mass and accretion-rate estimates for variable galaxies show that they sample a significantly different black-hole mass-accretion rate space, with masses a factor of 2.4 lower and accretion rates a factor of 22.5 lower than variable luminous AGN at the same redshift. We find that an empirical model based on a universal broken power-law PSD function, where the break frequency depends on SMBH mass and accretion rate, roughly reproduces the shape, but not the normalization, of the variability-luminosity trends measured for variable galaxies and more luminous AGN.

We use the 4Ms CDF-S survey to place direct X-ray constraints on the ubiquity of z~2 heavily obscured AGNs in K<22 BzK galaxies. Forty seven of the 222 BzK galaxies in the central region of the CDF-S are detected at X-ray energies, 11 of which have hard X-ray spectral slopes (Gamma<1) indicating the presence of heavily obscured AGN activity. The other 36 X-ray detected BzK galaxies appear to be relatively unobscured AGNs and starburst galaxies; we use X-ray variability analyses over a rest-frame baseline of ~3 years to further confirm the presence of AGN activity in many of these systems. The majority (7 out of 11) of the heavily obscured AGNs have excess IR emission over that expected from star formation (termed "IR-excess galaxies"). However, we find that X-ray detected heavily obscured AGNs only comprise ~25% of the IR-excess galaxy population, which is otherwise composed of relatively unobscured AGNs and starburst galaxies. We find that the typical X-ray spectrum of the heavily obscured AGNs is better characterized by a pure reflection model than an absorbed power-law model, suggesting extreme Compton-thick absorption in some systems. We verify this result by producing a composite rest-frame 2-20 keV spectrum, which has a similar shape as a reflection-dominated X-ray spectrum and reveals an emission feature at rest-frame energy ~6.4 keV, likely to be due to Fe K. These heavily obscured AGNs are likely to be the distant analogs of the reflection-dominated AGNs recently identified at z~0 with >10 keV observatories. On the basis of these analyses we estimate the space density for typical (intrinsic X-ray luminosities of L_X>1E43 erg/s) heavily obscured and Compton-thick AGNs at z~2. Our space-density constraints are conservative lower limits but they are already consistent with the range of predictions from X-ray background models.

We present results from a WIYN 3.5m telescope imaging study of the globular cluster (GC) systems of the edge-on spiral galaxies NGC891 and NGC4013. We used the 10' x 10' Minimosaic Imager to observe the galaxies in BVR filters to projected radii of ~20 kpc from the galaxy centers. We combined the WIYN data with archival and published data from WFPC2 and ACS on the Hubble Space Telescope to assess the contamination level of the WIYN GC candidate sample and to follow the GC systems further in toward the galaxies' centers. We constructed radial distributions for the GC systems using both the WIYN and HST data. The GC systems of NGC891 and NGC4013 extend to 9+/-3 kpc and 14+/-5 kpc, respectively, before falling off to undetectable levels in our images. We use the radial distributions to calculate global values for the total number (N_GC) and specific frequencies (S_N and T) of GCs. NGC4013 has N_GC = 140+/-20, S_N = 1.0+/-0.2 and T = 1.9+/-0.5; our N_GC value is ~40% smaller than a previous determination from the literature. The HST data were especially useful for NGC891, because the GC system is concentrated toward the plane of the galaxy and was only weakly detected in our WIYN images. Although NGC891 is thought to resemble the Milky Way in its overall properties, it has only half as many GCs, with N_GC = 70+/-20, S_N = 0.3+/-0.1 and T = 0.6+/-0.3. We also calculate the galaxy-mass-normalized number of blue (metal-poor) GCs in NGC891 and NGC4013 and find that they fall along a general trend of increasing specific frequency of blue GCs with increasing galaxy mass. Given currently available resources, the optimal method for studying the global properties of extragalactic GC systems is to combine HST data with wide-field, ground-based imaging with good resolution. The results here demonstrate the advantage gained by using both methods when possible.

We develop a new approach to the well-studied anti-correlation between the optical-to-X-ray spectral index, alpha_ox, and the monochromatic optical luminosity, l_opt. By cross-correlating the SDSS DR5 quasar catalog with the XMM-Newton archive, we create a sample of 327 quasars with X-ray S/N > 6, where both optical and X-ray spectra are available. This allows alpha_ox to be defined at arbitrary frequencies, rather than the standard 2500 Angstroms and 2 keV. We find that while the choice of optical wavelength does not strongly influence the alpha_ox-l_opt relation, the slope of the relation does depend on the choice of X-ray energy. The slope of the relation becomes steeper when alpha_ox is defined at low (~ 1 keV) X-ray energies. This change is significant when compared to the slope predicted by a decrease in the baseline over which alpha_ox is defined. The slopes are also marginally flatter than predicted at high (~ 10 keV) X-ray energies. Partial correlation tests show that while the primary driver of alpha_ox is l_opt, the Eddington ratio correlates strongly with alpha_ox when l_opt is taken into account, so accretion rate may help explain these results. We combine the alpha_ox-l_opt and Gamma -L_bol/L_Edd relations to naturally explain two results: 1) the existence of the Gamma-l_x relation as reported in Young et al. (2009) and 2) the lack of a Gamma-l_opt relation. The consistency of the optical/X-ray correlations establishes a more complete framework for understanding the relation between quasar emission mechanisms. We also discuss two correlations with the hard X-ray bolometric correction, which we show correlates with both alpha_ox and Eddington ratio. This confirms that an increase in accretion rate correlates with a decrease in the fraction of up-scattered disk photons.

We present a catalog of 792 DR5 SDSS quasars with optical spectra that have been observed serendipitously in the X-rays with the \emph{XMM-Newton}. These quasars cover a redshift range of z = 0.11 - 5.41 and a magnitude range of i = 15.3-20.7. Substantial numbers of radio-loud (70) and broad absorption line (51) quasars exist within this sample. Significant X-ray detections at >2sigma account for 87% of the sample (685 quasars), and 473 quasars are detected at >6sigma, sufficient to allow X-ray spectral fits. For detected sources, ~ 60% have X-ray fluxes between F(2-10 keV) = 1 - 10 x 10^-14 ergs cm^-2 s^-1. We fit a single power-law, a fixed power-law with intrinsic absorption left free to vary, and an absorbed power-law model to all quasars with X-ray S/N > 6, resulting in a weighted mean photon index Gamma = 1.91 +/- 0.08, with an intrinsic dispersion sigma_Gamma = 0.38. For the 55 sources (11.6%) that prefer intrinsic absorption, we find a weighted mean N_H = 1.5 +/- 0.3 x 10^21 cm^-2. We find that Gamma correlates significantly with optical color, Delta(g-i), the optical-to-X-ray spectral index (alpha_ox) and the X-ray luminosity. While the first two correlations can be explained as artefacts of undetected intrinsic absorption, the correlation between Gamma and X-ray luminosity appears to be a real physical correlation, indicating a pivot in the X-ray slope.

Using pulsewidth data for 872 isolated radio pulsars we test the hypothesis that pulsars evolve through a progressive narrowing of the emission cone combined with progressive alignment of the spin and magnetic axes. The new data provide strong evidence for the alignment over a time-scale of about 1 Myr with a log standard deviation of around 0.8 across the observed population. This time-scale is shorter than the time-scale of about 10 Myr found by previous authors, but the log standard deviation is larger. The results are inconsistent with models based on magnetic field decay alone or monotonic counter-alignment to orthogonal rotation. The best fits are obtained for a braking index parameter n_gamma approximately equal to 2.3, consistent the mean of the six measured values, but based on a much larger sample of young pulsars. The least-squares fitted models are used to predict the mean inclination angle between the spin and magnetic axes as a function of log characteristic age. Comparing these predictions to existing estimates it is found that the model in which pulsars are born with a random angle of inclination gives the best fit to the data. Plots of the mean beaming fraction as a function of characteristic age are presented using the best-fitting model parameters.

We present a correlation between the 2-10 keV spectral slope Gamma and the Eddington ratio L/L_{EDD} in a sample of ~400 Sloan Digital Sky Survey quasars with available hard X-ray spectra from XMM-Newton serendipitous observations. We find that the Gamma-L/L_{EDD} correlation is strongest in objects with black hole (BH) masses determined from the Hbeta line, and weaker (but still present) for those based on Mg II. An empirical non-linear correction of the Mg II-based masses, obtained by comparing the mass estimates in SDSS quasars having both Hbeta and Mg II measurements, significantly increases the strength of the correlation. No correlation is found among objects with BH masses derived from C IV, confirming that this line is not a reliable indicator of the BH mass. No significant correlation is found with the bolometric luminosity, while a Gamma-M_BH relation is present, though with a lower statistical significance than between Gamma and L/L_{EDD}. Our results imply a physical link between the accretion efficiency in the (cold) accretion disc of AGNs and the physical status of the (hot) corona responsible for the X-ray emission.

We have cross-correlated the SDSS DR3 Schneider et al. (2005) quasar catalog with the XMM-Newton archive. Color and redshift selections (g - r > 0.5 and 0.9 z < 2.1) result in a sample of 17 red, moderate redshift quasars. The redshift selection minimizes possible contamination due to host galaxy emission and Lyalpha forest absorption. Both optical and X-ray information are required to distinguish between the two likely remaining causes of the red colors: 1) dust-reddening and 2) an intrinsically red continuum. We find that 7 of 17 quasars can be classified as probable `intrinsically red' objects. These 7 quasars have unusually broad MgII emission lines (<FWHM>=10,500 km s^{-1}), moderately flat, but unabsorbed X-ray spectra <Gamma>=1.66+/-0.08), and low accretion rates (mdot/mdot_{Edd}} ~ 0.01). We suggest low accretion rates as a possible physical explanation for quasars with intrinsically red optical continua. We find that 8 of 17 quasars can be classified as dust-reddened. Three of these have upper-limits on the absorption column from X-ray spectral fits of N_H = 3-13 x 10^{22} cm^2, while the other five quasars must be absorbed by at least N_H = 10^{23} cm^2 in order to be consistent with a comparably selected alpha_{ox}-l_{uv} distribution. Two objects in the sample are unclassified.

It is argued that the polar gap and flux tube in the pulsar magnetosphere act as a resonant cavity/waveguide system which is excited by oscillations in the primary beam current and accelerating potential. The modes will be converted, probably scattered, to produce radio beams in the frequency range of those observed.