We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

We report mass measurements of neutron-rich Ga isotopes $^{80-85}$Ga with TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The measurements determine the masses of $^{80-83}$Ga in good agreement with previous measurements. The masses of $^{84}$Ga and $^{85}$Ga were measured for the first time. Uncertainties between $25-48$ keV were reached. The new mass values reduce the nuclear uncertainties associated with the production of A $\approx$ 84 isotopes by the \emph{r}-process for astrophysical conditions that might be consistent with a binary neutron star (BNS) merger producing a blue kilonova. Our nucleosynthesis simulations confirm that BNS merger may contribute to the first abundance peak under moderate neutron-rich conditions with electron fractions $Y_e=0.35-0.38$.

We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$\alpha$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

A new bound dark energy, BDE, cosmology has been proposed where the dark energy is the binding energy between light meson fields that condense a few tens of years after the big bang. It is reported that the correct dark energy density emerges using particle physics without fine tuning. This alone makes the BDE cosmology worthy of further investigation. This work looks at the late time BDE predictions of the evolution of cosmological parameters and the values of fundamental constants to determine whether the cosmology's predictions are consistent with observation. The work considers the time period between a scale factor of 0.1 and 1.0. A model BDE cosmology is considered with current day values of the cosmological parameters well within the observational limits. The calculations use three different values of the current day dark energy equation of state close to minus one. All three cases produce evolutions of the cosmological parameters and fundamental constants consistent with the observational constraints. Analytic relations between the BDE and cosmological parameters are developed to insure a consistent set of parameters.

This paper uses the beta function formalism to extend the analysis of quintessence cosmological parameters to the logarithmic and exponential dark energy potentials. The previous paper (Thompson 2018) demonstrated the formalism using power and inverse power potentials. The essentially identical evolution of the Hubble parameter for all of the quintessence cases and LambdaCDM is attributed to the flatness of the quintessence dark energy potentials in the dark energy dominated era. The Hubble parameter is therefore incapable of discriminating between static and dynamic dark energy. Unlike the other three potentials considered in the two papers the logarithmic dark energy potential requires a numerical integration in the formula for the superpotential rather than being an analytic function. The dark energy equation of state and the fundamental constants continue to be good discriminators between static and dynamical dark energy. A new analysis of quintessence with all four of the potentials relative the swampland conjectures indicates that the conjecture on the change in the scalar field is satisfied but that the conjecture on the change of the potential is not.

This investigation explores using the beta function formalism to calculate analytic solutions for the observable parameters in rolling scalar field cosmologies. The beta function in this case is the derivative of the scalar $\phi$ with respect to the natural log of the scale factor $a$, $\beta(\phi)=\frac{d \phi}{d \ln(a)}$. Once the beta function is specified, modulo a boundary condition, the evolution of the scalar $\phi$ as a function of the scale factor is completely determined. A rolling scalar field cosmology is defined by its action which can contain a range of physically motivated dark energy potentials. The beta function is chosen so that the associated "beta potential" is an accurate, but not exact, representation of the appropriate dark energy model potential. The basic concept is that the action with the beta potential is so similar to the action with the model potential that solutions using the beta action are accurate representations of solutions using the model action. The beta function provides an extra equation to calculate analytic functions of the cosmologies parameters as a function of the scale factor that are that are not calculable using only the model action. As an example this investigation uses a quintessence cosmology to demonstrate the method for power and inverse power law dark energy potentials. An interesting result of the investigation is that the Hubble parameter H is almost completely insensitive to the power of the potentials and that $\Lambda$CDM is part of the family of quintessence cosmology power law potentials with a power of zero.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

Non-parametric morphology measures are a powerful tool for identifying galaxy mergers at low redshifts. We employ cosmological zoom simulations using Gizmo with the Mufasa feedback scheme, post-processed using 3D dust radiative transfer into mock observations, to study whether common morphological measures Gini G, M20, concentration C, and asymmetry A are effective at identifying major galaxy mergers at z ~ 2 - 4, i.e. "Cosmic Noon". Our zoom suite covers galaxies with 10^8.6 < M_* < 10^11 M_sun at z ~ 2, and broadly reproduces key global galaxy observations. Our primary result is that these morphological measures are unable to robustly pick out galaxies currently undergoing mergers during Cosmic Noon, typically performing no better than a random guess. This improves only marginally if we consider whether galaxies have undergone a merger within the last Gyr. When also considering minor mergers, galaxies display no trend of moving towards the merger regime with increasing merger ratio. From z = 4 -> 2, galaxies move from the non-merger towards the merger regime in all statistics, but this is primarily an effect of mass: Above a given noise level, higher mass galaxies display a more complex outer morphology induced by their clustered environment. We conclude that during Cosmic Noon, these morphological statistics are of limited value in identifying galaxy mergers.

We utilize the hydrodynamic and N-body code {\small GIZMO} coupled with our newly developed sub-grid Population~III (Pop~III) Legacy model, designed specifically for cosmological volume simulations, to study the baseline metal enrichment from Pop~III star formation at $z>7$. In this idealized numerical experiment, we only consider Pop~III star formation. We find that our model Pop~III star formation rate density (SFRD), which peaks at $\sim 10^{-3}\ {\rm M_\odot yr^{-1} Mpc^{-1}}$ near $z\sim10$, agrees well with previous numerical studies and is consistent with the observed estimates for Pop~II SFRDs. The mean Pop~III metallicity rises smoothly from $z=25-7$, but does not reach the critical metallicity value, $Z_{\rm crit}=10^{-4}\ Z_\odot$, required for the Pop~III to Pop~II transition in star formation mode until $z\simeq7$. This suggests that, while individual halos can suppress in-situ Pop~III star formation, the external enrichment is insufficient to globally terminate Pop~III star formation. The maximum enrichment from Pop~III star formation in star forming dark matter halos is $Z\sim10^{-2}\ Z_\odot$, whereas the minimum found in externally enriched haloes is $Z\gtrsim10^{-7}\ Z_\odot$. Finally, mock observations of our simulated IGM enriched with Pop~III metals produce equivalent widths similar to observations of an extremely metal poor damped Lyman alpha (DLA) system at $z=7.04$, which is thought to be enriched by Pop~III star formation only.

We present new photometric data of 11 hot Jupiter transiting exoplanets (CoRoT-12b, HAT-P-5b, HAT-P-12b, HAT-P-33b, HAT-P-37b, WASP-2b, WASP-24b, WASP-60b, WASP-80b, WASP-103b, XO-3b) in order to update their planetary parameters and to constrain information about their atmospheres. These observations of CoRoT-12b, HAT-P-37b and WASP-60b are the first follow-up data since their discovery. Additionally, the first near-UV transits of WASP-80b and WASP-103b are presented. We compare the results of our analysis with previous work to search for transit timing variations (TTVs) and a wavelength dependence in the transit depth. TTVs may be evidence of a third body in the system and variations in planetary radius with wavelength can help constrain the properties of the exoplanet's atmosphere. For WASP-103b and XO-3b, we find a possible variation in the transit depths that may be evidence of scattering in their atmospheres. The B-band transit depth of HAT-P-37b is found to be smaller than its near-IR transit depth and such a variation may indicate TiO/VO absorption. These variations are detected from 2-4.6$\sigma$, so follow-up observations are needed to confirm these results. Additionally, a flat spectrum across optical wavelengths is found for 5 of the planets (HAT-P-5b, HAT-P-12b, WASP-2b, WASP-24b, WASP-80b), suggestive that clouds may be present in their atmospheres. We calculate a refined orbital period and ephemeris for all the targets, which will help with future observations. No TTVs are seen in our analysis with the exception of WASP-80b and follow-up observations are needed to confirm this possible detection.

Of the almost 40 star forming galaxies at z>~5 (not counting QSOs) observed in [CII] to date, nearly half are either very faint in [CII], or not detected at all, and fall well below expectations based on locally derived relations between star formation rate (SFR) and [CII] luminosity. Combining cosmological zoom simulations of galaxies with SIGAME (SImulator of GAlaxy Millimeter/submillimeter Emission) we have modeled the multi-phased interstellar medium (ISM) and its emission in [CII], [OI] and [OIII], from 30 main sequence galaxies at z~6 with star formation rates ~3-23Msun/yr, stellar masses ~(0.7-8)x10^9Msun, and metallicities ~(0.1-0.4)xZsun. The simulations are able to reproduce the aforementioned [CII]-faintness at z>5, match two of the three existing z>~5 detections of [OIII], and are furthermore roughly consistent with the [OI] and [OIII] luminosity relations with SFR observed for local starburst galaxies. We find that the [CII] emission is dominated by the diffuse ionized gas phase and molecular clouds, which on average contribute ~66% and ~27%, respectively. The molecular gas, which constitutes only ~10% of the total gas mass is thus a more efficient emitter of [CII] than the ionized gas making up ~85% of the total gas mass. A principal component analysis shows that the [CII] luminosity correlates with the star formation activity as well as average metallicity. The low metallicities of our simulations together with their low molecular gas mass fractions can account for their [CII]-faintness, and we suggest these factors may also be responsible for the [CII]-faint normal galaxies observed at these early epochs.

We examine the growth and evolution of quenched galaxies in the Mufasa cosmological hydrodynamic simulations that include an evolving halo mass-based quenching prescription, with galaxy colours computed accounting for line-of-sight extinction to individual star particles. Mufasa reproduces the observed present-day red sequence reasonably well, including its slope, amplitude, and scatter. In Mufasa, the red sequence slope is driven entirely by the steep stellar mass-stellar metallicity relation, which independently agrees with observations. High-mass star-forming galaxies blend smoothly onto the red sequence, indicating the lack of a well-defined green valley at M*>10^10.5 Mo. The most massive galaxies quench the earliest and then grow very little in mass via dry merging; they attain their high masses at earlier epochs when cold inflows more effectively penetrate hot halos. To higher redshifts, the red sequence becomes increasingly contaminated with massive dusty star-forming galaxies; UVJ selection subtly but effectively separates these populations. We then examine the evolution of the mass functions of central and satellite galaxies split into passive and star-forming via UVJ. Massive quenched systems show good agreement with observations out to z~2, despite not including a rapid early quenching mode associated with mergers. However, low-mass quenched galaxies are far too numerous at z<1 in Mufasa, indicating that Mufasa strongly over-quenches satellites. A challenge for hydrodynamic simulations is to devise a quenching model that produces enough early massive quenched galaxies and keeps them quenched to z=0, while not being so strong as to over-quench satellites; Mufasa's current scheme fails at the latter.

We utilise a series of high-resolution cosmological zoom simulations of galaxy formation to investigate the relationship between the ultraviolet (UV) slope, beta, and the ratio of the infrared luminosity to UV luminosity (IRX) in the spectral energy distributions (SEDs) of galaxies. We employ dust radiative transfer calculations in which the SEDs of the stars in galaxies propagate through the dusty interstellar medium. Our main goals are to understand the origin of, and scatter in the IRX-beta relation; to assess the efficacy of simplified stellar population synthesis screen models in capturing the essential physics in the IRX-beta relation; and to understand systematic deviations from the canonical local IRX-beta relations in particular populations of high-redshift galaxies. Our main results follow. Galaxies that have young stellar populations with relatively cospatial UV and IR emitting regions and a Milky Way-like extinction curve fall on or near the standard Meurer relation. This behaviour is well captured by simplified screen models. Scatter in the IRX-beta relation is dominated by three major effects: (i) older stellar populations drive galaxies below the relations defined for local starbursts due to a reddening of their intrinsic UV SEDs; (ii) complex geometries in high-z heavily star forming galaxies drive galaxies toward blue UV slopes owing to optically thin UV sightlines; (iii) shallow extinction curves drive galaxies downward in the IRX-beta plane due to lowered NUV/FUV extinction ratios. We use these features of the UV slopes of galaxies to derive a fitting relation that reasonably collapses the scatter back toward the canonical local relation. Finally, we use these results to develop an understanding for the location of two particularly enigmatic populations of galaxies in the IRX-beta plane: z~2-4 dusty star forming galaxies, and z>5 star forming galaxies.

The observed constraints on the variability of the proton to electron mass ratio $\mu$ and the fine structure constant $\alpha$ are used to establish constraints on the variability of the Quantum Chromodynamic Scale and a combination of the Higgs Vacuum Expectation Value and the Yukawa couplings. Further model dependent assumptions provide constraints on the Higgs VEV and the Yukawa couplings separately. A primary conclusion is that limits on the variability of dimensionless fundamental constants such as $\mu$ and $\alpha$ provide important constraints on the parameter space of new physics and cosmologies.

Many physical theories beyond the Standard Model predict time variations of basic physics parameters. Direct measurement of the time variations of these parameters is very difficult or impossible to achieve. By contrast, measurements of fundamental constants are relatively easy to achieve, both in the laboratory and by astronomical spectra of atoms and molecules in the early universe. In this work measurements of the proton to electron mass ratio $\mu$ and the fine structure constant $\alpha$ are combined to place mildly model dependent limits on the fractional variation of the Quantum Chromodynamic Scale and the sum of the fractional variations of the Higgs Vacuum Expectation Value and the Yukawa couplings on time scales of more than half the age of the universe. The addition of another model parameter allows the fractional variation of the Higgs VEV and the Yukawa couplings to be computed separately. Limits on their variation are found at the level of less than $5 \times 10^{-5}$ over the past seven gigayears. A model dependent relation between the expected fractional variation of $\alpha$ relative to $\mu$ tightens the limits to $10^{-7}$ over the same time span. Limits on the present day rate of change of the constants and parameters are then calculated using slow roll quintessence. A primary result of this work is that studies of the dimensionless fundamental constants such as $\alpha$ and $\mu$, whose values depend on the values of the physics parameters, are excellent monitors of the limits on the time variation of these parameters.

We examine galaxy star formation rates (SFRs), metallicities, and gas contents predicted by the MUFASA cosmological hydrodynamic simulations, which employ meshless hydrodynamics and novel feedback prescriptions that yield a good match to observed galaxy stellar mass assembly. We combine 50, 25, and 12.5 Mpc/h boxes with a quarter billion particles each to show that MUFASA broadly reproduces a wide range of relevant observations, including SFR and specific SFR functions, the mass-metallicity relation, HI and H2 fractions, HI (21 cm) and CO luminosity functions, and cosmic gas density evolution. There are mild but significant discrepancies, such as too many high-SFR galaxies, overly metal-rich and HI-poor galaxies at M*>10^{10} Mo, and sSFRs that are too low at z~1-2. The HI mass function increases by x2 out to z~1 then steepens to higher redshifts, while the CO luminosity function computed using the Narayanan et al. conversion factor shows a rapid increase of CO-bright galaxies out to z~2 in accord with data. Omega_HI and Omega_H2 both scale roughly as (1+z)^0.7 out to z~3, comparable to the rise in HI and H2 fractions. MUFASA galaxies with high SFR at a given M* have lower metallicities and higher HI and H2 fractions, following observed trends; we make quantitative predictions for how fluctuations in the baryon cycle drive correlated scatter around galaxy scaling relations. Most of these trends are well converged with numerical resolution. These successes highlight MUFASA as a viable platform to study many facets of cosmological galaxy evolution.

gPhoton is a new database product and software package that enables analysis of GALEX ultraviolet data at the photon level. The project's stand-alone, pure-Python calibration pipeline reproduces the functionality of the original mission pipeline to reduce raw spacecraft data to lists of time-tagged, sky-projected photons, which are then hosted in a publicly available database by the Mikulski Archive at Space Telescope (MAST). This database contains approximately 130 terabytes of data describing approximately 1.1 trillion sky-projected events with a timestamp resolution of five milliseconds. A handful of Python and command line modules serve as a front-end to interact with the database and to generate calibrated light curves and images from the photon-level data at user-defined temporal and spatial scales. The gPhoton software and source code are in active development and publicly available under a permissive license. We describe the motivation, design, and implementation of the calibration pipeline, database, and tools, with emphasis on divergence from prior work, as well as challenges created by the large data volume. We summarize the astrometric and photometric performance of gPhoton relative to the original mission pipeline. For a brief example of short time domain science capabilities enabled by gPhoton, we show new flares from the known M dwarf flare star CR Draconis. The gPhoton software has permanent object identifiers with the ASCL (ascl:1603.004) and DOI (doi:10.17909/T9CC7G). This paper describes the software as of version v1.27.2.

Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection.

We present the MUFASA suite of cosmological hydrodynamic simulations, which employs the GIZMO meshless finite mass (MFM) code including H2-based star formation, nine-element chemical evolution, two-phase kinetic outflows following scalings from the Feedback in Realistic Environments zoom simulations, and evolving halo mass-based quenching. Our fiducial (50 Mpc/h)^3 volume is evolved to z=0 with a quarter billion particles, The predicted galaxy stellar mass functions (GSMF) reproduce observations from z=4-0 to <1.2sigma in cosmic variance, providing an unprecedented match to this key diagnostic. The cosmic star formation history and stellar mass growth show general agreement with data, with a strong archaeological downsizing trend such that dwarf galaxies form the majority of their stars after z~1. We run 25 Mpc/h and 12.5 Mpc/h volumes to z=2 with identical feedback prescriptions, the latter resolving all hydrogen-cooling halos, and the three runs display fair resolution convergence. The specific star formation rates broadly agree with data at z=0, but are underpredicted at z~2 by a factor of three, re-emphasizing a longstanding puzzle in galaxy evolution models. We compare runs using MFM and two flavours of Smoothed Particle Hydrodynamics, and show that the GSMF is sensitive to hydrodynamics methodology at the ~x2 level, which is sub-dominant to choices for parameterising feedback.

The sources that drove cosmological reionization left clues regarding their identity in the slope and inhomogeneity of the ultraviolet ionizing background (UVB): Bright quasars (QSOs) generate a hard UVB with predominantly large-scale fluctuations while Population II stars generate a softer one with smaller-scale fluctuations. Metal absorbers probe the UVB's slope because different ions are sensitive to different energies. Likewise, they probe spatial fluctuations because they originate in regions where a galaxy-driven UVB is harder and more intense. We take a first step towards studying the reionization-epoch UVB's slope and inhomogeneity by comparing observations of 12 metal absorbers at $z\sim6$ versus predictions from a cosmological hydrodynamic simulation using three different UVBs: a soft, spatially-inhomogeneous "galaxies+QSOs" UVB; a homogeneous "galaxies+QSOs" UVB (Haardt & Madau 2012); and a QSOs-only model. All UVBs reproduce the observed column density distributions of CII, SiIV, and CIV reasonably well although high-column, high-ionization absorbers are underproduced, reflecting numerical limitations. With upper limits treated as detections, only a soft, fluctuating UVB reproduces both the observed SiIV/CIV and CII/CIV distributions. The QSOs-only UVB overpredicts both CIV/CII and CIV/SiIV, indicating that it is too hard. The Haardt & Madau (2012) UVB underpredicts CIV/SiIV, suggesting that it lacks amplifications near galaxies. Hence current observations prefer a soft, fluctuating UVB as expected from a predominantly Population II background although they cannot rule out a harder one. Future observations probing a factor of two deeper in metal column density will distinguish between the soft, fluctuating and QSOs-only UVBs.

Transits of exoplanets observed in the near-UV have been used to study the scattering properties of their atmospheres and possible star-planet interactions. We observed the primary transits of 15 exoplanets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-16b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-12b, WASP-33b, WASP-36b, WASP-44b, WASP-48b, and WASP-77Ab) in the near-UV and several optical photometric bands to update their planetary parameters, ephemerides, search for a wavelength dependence in their transit depths to constrain their atmospheres, and determine if asymmetries are visible in their light curves. Here we present the first ground-based near-UV light curves for 12 of the targets (CoRoT-1b, GJ436b, HAT-P-1b, HAT-P-13b, HAT-P-22b, TrES-2b, TrES-4b, WASP-1b, WASP-33b, WASP-36b, WASP-48b, and WASP-77Ab). We find that none of the near-UV transits exhibit any non-spherical asymmetries, this result is consistent with recent theoretical predictions by Ben-Jaffel et al. and Turner et al. The multi-wavelength photometry indicates a constant transit depth from near-UV to optical wavelengths in 10 targets (suggestive of clouds), and a varying transit depth with wavelength in 5 targets (hinting at Rayleigh or aerosol scattering in their atmospheres). We also present the first detection of a smaller near-UV transit depth than that measured in the optical in WASP-1b and a possible opacity source that can cause such radius variations is currently unknown. WASP-36b also exhibits a smaller near-UV transit depth at 2.6$\sigma$. Further observations are encouraged to confirm the transit depth variations seen in this study.

Astronomical observations have a unique ability to determine the laws of physics at distant times in the universe. They, therefore, have particular relevance in answering the basic question as to whether the laws of physics are invariant with time. The dimesionless fundamental constants, such as the proton to electron mass ratio and the fine structure constant are key elements in the investigation. If they vary with time then the answer is clearly that the laws of physics are not invariant with time and significant new physics must be developed to describe the universe. Limits on their variance, on the other hand, constrains the parameter space available to new physics that requires a variation with time of basic physical law. There are now observational constraints on the time variation of the proton to electron mass ratio mu at the 1.E-7 level. Constraints on the variation of the fine structure constant alpha are less rigorous, 1E-5, but are imposed at higher redshift. The implications of these limits on new cosmologies that require rolling scalar fields has already had its first investigations. Here we address the implications on basic particle physics. The proton to electron mass ratio is obviously dependent on the particle physics parameters that set the mass of the proton and the electron. To first order the ratio is dependent on a combination of the Quantum Chromodynamic scale, the Yukawa couplings, and the Higgs Vacuum Expectation Value. Here that relationship is quantitative defined for the first time. When coupled with previous determinations of the relation of the fine structure constant to the same parameters two constraints exist on the fractional variation of these parameters with time. A third independent constraint involving only the three parameters could set the stage for constraints on their individual fractional variation.

We present S\'IGAME simulations of the [CII]157.7$\mu$ fine structure line emission from cosmological smoothed particle hydrodynamics (SPH) simulations of seven main sequence galaxies at z=2. Using sub-grid physics prescriptions the gas in our simulations is modeled as a multi-phased interstellar medium (ISM) comprised of molecular gas residing in giant molecular clouds, an atomic gas phase associated with photo-dissociation regions (PDRs) at the cloud surfaces, and a diffuse, ionized gas phase. Adopting logotropic cloud density profiles and accounting for heating by the local FUV radiation field and cosmic rays by scaling both with local star formation rate (SFR) volume density, we calculate the [CII] emission using a photon escape probability formalism. The [CII] emission peaks in the central $\lesssim$1 kpc of our galaxies as do the SFR radial profiles, with most [CII] ($\gtrsim$70%) originating in the molecular gas phase, whereas further out ($\gtrsim$2 kpc), the atomic/PDR gas dominates ($\gtrsim$90%) the [CII] emission, no longer tracing on-going star formation. Throughout, the ionized gas contribution is negligible ($\lesssim$3%). The [CII] luminosity vs. SFR ([CII]-SFR) relationship, integrated as well as spatially resolved (on scales of 1 kpc), delineated by our simulated galaxies is in good agreement with the corresponding relations observed locally and at high redshifts. In our simulations, the molecular gas dominates the [CII] budget at SFR$\gtrsim$20 M$_{\odot}$yr$^{-1}$ (SFR density $\gtrsim$0.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$), while atomic/PDR gas takes over at lower SFRs, suggesting a picture in which [CII] predominantly traces the molecular gas in high-density/pressure regions where star formation is on-going, and otherwise reveals the atomic/PDR gas phase.

Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence).

We investigate the time evolution of non-parametric morphological quantities and their relationship to major mergers between $4\geq z \geq 2$ in high-resolution cosmological zoom simulations of disk galaxies that implement kinetic wind feedback, $H_2$-based star formation, and minimal ISM pressurisation. We show that the resulting galaxies broadly match basic observed physical properties of $z\sim 2$ objects. We measure the galaxies' concentrations ($C$), asymmetries ($A$), and $Gini$ ($G$) and $M_{20}$ coefficients, and correlate these with major merger events identified from the mass growth history. We find that high values of asymmetry provide the best indicator for identifying major mergers of $>1:4$ mass ratio within our sample, with $Gini$-$M_{20}\,$ merger classification only as effective for face-on systems and much less effective for edge-on or randomly-oriented galaxies. The canonical asymmetry cut of $A\geq0.35$, however, is only able to correctly identify major mergers $\sim 10\%$ of the time, while a higher cut of $A\geq 0.8$ more efficiently picks out mergers at this epoch. We further examine the temporal correlation between morphological statistics and mergers, and show that for randomly-oriented galaxies, half the galaxies with $A\geq0.8$ undergo a merger within $\pm0.2\,{\rm Gyr}$, whereas $Gini$-$M_{20}\,$ identification only identifies about a third correctly. The fraction improves further using $A\geq 1.5$, but about the half the mergers are missed by this stringent cut.

The Bullet Cluster has provided some of the best evidence for the $\Lambda$ cold dark matter ($\Lambda\rm{CDM}$) model via direct empirical proof of the existence of collisionless dark matter, while posing a serious challenge owing to the unusually high inferred pairwise velocities of its progenitor clusters. Here we investigate the probability of finding such a high-velocity pair in large-volume N-body simulations, particularly focusing on differences between halo finding algorithms. We find that algorithms that do not account for the kinematics of infalling groups yield vastly different statistics and probabilities. When employing the ROCKSTAR halo finder that considers particle velocities, we find numerous Bullet-like pair candidates that closely match not only the high pairwise velocity, but also the mass, mass ratio, separation distance, and collision angle of the initial conditions that have been shown to produce the Bullet Cluster in non-cosmological hydrodynamic simulations. The probability of finding a high pairwise velocity pair among haloes with $M_{\rm halo}\geq10^{14} M_{\odot}$ is $4.6\times 10^{-4}$ using ROCKSTAR, while it is $\approx 34\times$ lower using a friends-of-friends (FOF) based approach as in previous studies. This is because the typical spatial extent of Bullet progenitors is such that FOF tends to group them into a single halo despite clearly distinct kinematics. Further requiring an appropriately high average mass among the two progenitors, we find the comoving number density of potential Bullet-like candidates to be on the order of $\approx10^{-10} \,{\rm Mpc}^{-3}$. Our findings suggest that $\Lambda\rm{CDM}$ straightforwardly produces massive, high relative velocity halo pairs analogous to Bullet Cluster progenitors, and hence the Bullet Cluster does not present a challenge to the $\Lambda\rm{CDM}$ model.

Observations suggest that CII was more abundant than CIV in the intergalactic medium towards the end of the hydrogen reionization epoch. This transition provides a unique opportunity to study the enrichment history of intergalactic gas and the growth of the ionizing background (UVB) at early times. We study how carbon absorption evolves from z=10-5 using a cosmological hydrodynamic simulation that includes a self-consistent multifrequency UVB as well as a well-constrained model for galactic outflows to disperse metals. Our predicted UVB is within 2-4 times that of Haardt & Madau (2012), which is fair agreement given the uncertainties. Nonetheless, we use a calibration in post-processing to account for Lyman-alpha forest measurements while preserving the predicted spectral slope and inhomogeneity. The UVB fluctuates spatially in such a way that it always exceeds the volume average in regions where metals are found. This implies both that a spatially-uniform UVB is a poor approximation and that metal absorption is not sensitive to the epoch when HII regions overlap globally even at column densites of 10^{12} cm^{-2}. We find, consistent with observations, that the CII mass fraction drops to low redshift while CIV rises owing the combined effects of a growing UVB and continued addition of carbon in low-density regions. This is mimicked in absorption statistics, which broadly agree with observations at z=6-3 while predicting that the absorber column density distributions rise steeply to the lowest observable columns. Our model reproduces the large observed scatter in the number of low-ionization absorbers per sightline, implying that the scatter does not indicate a partially-neutral Universe at z=6.

Iris is an extensible application that provides astronomers with a user-friendly interface capable of ingesting broad-band data from many different sources in order to build, explore, and model spectral energy distributions (SEDs). Iris takes advantage of the standards defined by the International Virtual Observatory Alliance, but hides the technicalities of such standards by implementing different layers of abstraction on top of them. Such intermediate layers provide hooks that users and developers can exploit in order to extend the capabilities provided by Iris. For instance, custom Python models can be combined in arbitrary ways with the Iris built-in models or with other custom functions. As such, Iris offers a platform for the development and integration of SED data, services, and applications, either from the user's system or from the web. In this paper we describe the built-in features provided by Iris for building and analyzing SEDs. We also explore in some detail the Iris framework and software development kit, showing how astronomers and software developers can plug their code into an integrated SED analysis environment.

It is important to explore the diversity of characteristics of low-mass, low-density planets to understand the nature and evolution of this class of planets. We present a homogeneous analysis of 12 new and 9 previously published broadband photometric observations of the Uranus-sized extrasolar planet GJ 3470b, which belongs to the growing sample of sub-Jovian bodies orbiting M dwarfs. The consistency of our analysis explains some of the discrepancies between previously published results and provides updated constraints on the planetary parameters. Our data are also consistent with previous transit observations of this system. We also provide new spectroscopic measurements of GJ 3470 from 0.33 to 2.42 $\mu$$m$ to aid our analysis. We find $R_{\star}$ = 0.48$\pm$0.04 $R_{\odot}$, $M_{\star}$ = 0.51$\pm$0.06 $M_{\odot}$, and $T_{\rm eff}$ = 3652$\pm$50 K for GJ 3470, along with a rotation period of $20.70\pm{0.15}$ d and an R-band amplitude of 0.01 mag, which is small enough that current transit measurements should not be strongly affected by stellar variability. We also present the most precise orbital ephemeris for this system: T$_{o}$ = 2455983.70472$\pm$0.00021 BJD$_{TDB}$, P = 3.3366487$^{+0.0000043}_{-0.0000033}$ d, and we see no evidence for transit timing variations greater than 1 minute. Our reported planet to star radius ratio is 0.07642$\pm$0.00037. The physical parameters of this planet are $R_{p}$ = 3.88$\pm$0.32 $R_{\oplus}$, and $M_{p}$ = 13.73$\pm$1.61 $M_{\oplus}$. Because of our revised stellar parameters, the planetary radius we present is smaller than previously reported values. We also perform a second analysis of the transmission spectrum of the entire ensemble of transit observations to date, supporting the existence of a H$_{2}$ dominated atmosphere exhibiting a strong Rayleigh scattering slope.

We present the first large, blind and uniformly selected survey for molecular hydrogen (H_2) in damped Lyman alpha systems (DLAs) with moderate-to-high resolution spectra. 86 DLAs were searched for absorption in the many Lyman and Werner H_2 transitions, with approximately 79% completeness for H_2 column densities above N(H_2)=10^17.5 cm^-2 for an assumed Doppler broadening parameter b=2 km s^-1. Only a single strong H_2 absorber was found -- a system detected previously in VLT/UVES spectra. Given our distribution of N(H_2) upper limits, this ~1% detection rate is smaller than expected from previous surveys at 99.8% confidence. Assuming the N(H_2) distribution shape from previous surveys, our detection rate implies a covering factor of ~1% for N(H_2) >= 10^17.5 cm^-2 gas in DLAs (<6% at 95% confidence). We obtained new Magellan/MagE spectra for 53 DLAs; 8 km s^-1 -resolution spectra were available for 27 DLAs. MagE's moderate resolution (~71 km s^-1) yields weaker N(H_2) upper limits and makes them dependent on the assumed Doppler parameter. For example, half the (relevant) previous H_2 detections have N(H_2) >= 10^18.1 cm^-2, a factor of just 3 higher than our median upper limit. Nevertheless, several tests suggest our upper limits are accurate, and they would need to be increased by 1.8 dex to bring our detection rate within 95% confidence of previous surveys.

Recent observations have indicated the existence of dust in high-redshift galaxies, however, the dust properties in them are still unknown. Here we present theoretical constraints on dust properties in Lyman break galaxies (LBGs) at z=3 by post-processing a cosmological smoothed particle hydrodynamics simulation with radiative transfer calculations. We calculate the dust extinction in 2800 dark matter halos using the metallicity information of individual gas particles in our simulation. We use only bright galaxies with rest-frame UV magnitude M_1700 < -20 mag, and study the dust size, dust-to-metal mass ratio, and dust composition. From the comparison of calculated color excess between B and V-band (i.e., E(B-V)) and the observations, we constrain the typical dust size, and show that the best-fitting dust grain size is ~ 0.05 micron, which is consistent with the results of theoretical dust models for Type-II supernova. Our simulation with the dust extinction effect can naturally reproduce the observed rest-frame UV luminosity function of LBGs at z=3 without assuming an ad hoc constant extinction value. In addition, in order to reproduce the observed mean E(B-V), we find that the dust-to-metal mass ratio needs to be similar to that of the local galaxies, and that the graphite dust is dominant or at least occupy half of dust mass.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

The values of the fundamental constants such as $\mu = m_P/m_e$, the proton to electron mass ratio and $\alpha$, the fine structure constant, are sensitive to the product $\sqrt{\zeta_x^2(w+1)}$ where $\zeta_x$ is a coupling constant between a rolling scalar field responsible for the acceleration of the expansion of the universe and the electromagnetic field with x standing for either $\mu$ or $\alpha$. The dark energy equation of state $w$ can assume values different than $-1$ in cosmologies where the acceleration of the expansion is due to a scalar field. In this case the value of both $\mu$ and $\alpha$ changes with time. The values of the fundamental constants, therefore, monitor the equation of state and are a valuable tool for determining $w$ as a function of redshift. In fact the rolling of the fundamental constants is one of the few definitive discriminators between acceleration due to a cosmological constant and acceleration due to a quintessence rolling scalar field. $w$ is often given in parameterized form for comparison with observations. In this manuscript the predicted evolution of $\mu$, is calculated for a range of parameterized equation of state models and compared to the observational constraints on $\Delta \mu / \mu$. We find that the current limits on $\Delta \mu / \mu$ place significant constraints on linear equation of state models and on thawing models where $w$ deviates from $-1$ at late times. They also constrain non-dynamical models that have a constant $w$ not equal to $-1$. These constraints are an important compliment to geometric tests of $w$ in that geometric tests are sensitive to the evolution of the universe before the epoch of observation while fundamental constants are sensitive to the evolution of the universe after the observational epoch. Abstract truncated.

WFIRST is the highest priority space mission of the Decadal review, however, it is unlikely to begin in this decade primarily due to a anticipated NASA budget that is unlikely to have sufficient resources to fund such a mission. For this reason we present a lower cost mission that accomplishes all of the WFIRST science as described in the Design Reference Mission 1 with a probe class design. This is effort is motivated by a desire to begin WFIRST in a timely manner and within a budget that can fit within the assets available to NASA on a realistic basis. The design utilizes dichroics to form four focal planes all having the same field of view to use the majority of available photons from a 1.2 meter telescope.

It has been shown observationally that star formation (SF) correlates tightly with the presence of molecular hydrogen (H$_2$). Therefore it would be important to investigate its implication on galaxy formation in a cosmological context. In the present work, we track the H$_2$ mass fraction within our cosmological smoothed particle hydrodynamics (SPH) code GADGET-3 using an equilibrium analytic model by Krumholz et al. This model allows us to regulate the star formation in our simulation by the local abundance of H$_2$ rather than the total cold gas density, and naturally introduce the dependence of star formation on metallicity. We investigate implications of the equilibrium H$_2$-based SF model on galaxy population properties, such as the stellar-to-halo mass ratio (SHMR), baryon fraction, cosmic star formation rate density (SFRD), galaxy specific SFR, galaxy stellar mass functions (GSMF), and Kennicutt-Schmidt (KS) relationship. The advantage of our work over the previous ones is having a large sample of simulated galaxies in a cosmological volume from high-redshift to $z=0$. We find that low-mass halos with M$_{DM}<10^{10.5}$ M$_\odot$ are less efficient in producing stars in the H$_2$-based SF model at $z\geq6$, which brings the simulations to a better agreement with observational estimates of SHMR and GSMF at the low-mass end. This is particularly evident by a reduction in the number of low-mass galaxies at M$_\star\leq10^{8}$ M$_\odot$ in the GSMF. The overall SFRD is also reduced at high-$z$ in the H$_2$ run, which results in slightly higher SFRD at low-redshift due to more abundant gas available for star formation at later times. This new H$_2$ model is able to reproduce the empirical KS relationship at $z=0$ naturally without the need for setting its normalization by hand, and overall it seems to have more advantages than the previous pressure-based SF model.

We examine how HI and metal absorption lines within low-redshift galaxy halos trace the dynamical state of circumgalactic gas, using cosmological hydrodynamic simulations that include a well-vetted heuristic model for galactic outflows. We categorize inflowing, outflowing, and ambient gas based on its history and fate as tracked in our simulation. Following our earlier work showing that the ionisation level of absorbers was a primary factor in determining the physical conditions of absorbing gas, we show here that it is also a governing factor for its dynamical state. Low-ionisation metal absorbers (e.g. MgII) tend to arise in gas that will fall onto galaxies within several Gyr, while high-ionisation metal absorbers (e.g. OVI) generally trace material that was deposited by outflows many Gyr ago. Inflowing gas is dominated by enriched material that was previously ejected in an outflow, hence accretion at low redshifts is typically substantially enriched. Recycling wind material is preferentially found closer to galaxies, and is more dominant in lower-mass halos since high-mass halos have more hot gas that is able to support itself against infall. Low-mass halos also tend to re-eject more of their accreted material, owing to our outflow prescription that employs higher mass loading factors for lower-mass galaxies. Typical HI absorbers trace unenriched ambient material that is not participating in the baryon cycle, but stronger HI absorbers arise in cool, enriched inflowing gas. Instantaneous radial velocity measures of absorbers are generally poor at distinguishing between inflowing and outflowing gas, except in the case of very recent outflows. These results suggest that probing halo gas using a range of absorbers can provide detailed information about the amount and physical conditions of material that is participating in the baryon cycle.

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

We present the chemical abundance measurements of the first large, medium-resolution, uniformly selected damped Lyman alpha system (DLA) survey. The sample contains 99 DLAs towards 89 quasars selected from the SDSS DR5 DLA sample in a uniform way. We analyze the metallicities and kinematic diagnostics, including the velocity width of 90% of the optical depth, \Delta v_{90}, and the equivalent widths of the SiII 1526, CIV 1548 and MgII 2796 transitions. To avoid strong line-saturation effects on the metallicities measured in medium-resolution spectra (FWHM~71 km/s), we derived metallicities from metal transitions which absorbed at most 35% of the quasar continuum flux. We find the evolution in cosmic mean metallicity of the sample, <Z> = (-0.04 +- 0.013)z -(1.06 +- 0.36), consistent with no evolution over the redshift range z ~ [2.2, 4.4], but note that the majority of our sample falls at z ~ [2.2, 3.5]. The apparent lack of metallicity evolution with redshift is also seen in a lack of evolution in the median \Delta v_{90} and SiII 1526 equivalent width values. While this result may seem to conflict with other large surveys that have detected significant metallicity evolution, such as Rafelski et al. 2012 who found <Z> = (-0.22 +- 0.03)z - (0.65 +- 0.09) over z ~ [0, 5], several tests show that these surveys are not inconsistent with our new result. However, over the smaller redshift range covered by our uniformly-selected sample, the true evolution of the cosmic mean metallicity in DLAs may be somewhat flatter than the Rafelski et al. 2012 estimate.

We present the results of a numerical study examining the effect of H_2-based star formation (SF) model on the rest-frame UV luminosity function and star formation rate function (SFRF) of z>=6 galaxies, and the implications for reionization. Using cosmological hydrodynamical simulations outfitted with an H_2-SF model, we find good agreement with our previous results (non-H_2 SF model) and observations at Muv<=-18. However at Muv>-18, we find that the LF deviates from both our previous work and current observational extrapolations, producing significantly fewer low-luminosity galaxies and exhibiting additional turnover at the faint end. We constrain the redshift evolution of this turnover point using a modified Schechter function that includes additional terms to quantify the turnover magnitude (Muv^t) and subsequent slope (Beta). We find that Muv^t evolves from Muv^t=-17.33 (at z=8) to -15.38 (z=6), while Beta becomes shallower by \Delta Beta=0.22 during the same epoch. This occurs in an Muv range which will be observable by JWST. By integrating the SFRF, we determine that even though H_2-SF model significantly reduces the number density of low luminosity galaxies at Muv>-18, it does not suppress the total SFR density enough to affect the capability of SF to maintain reionization.

New PKS1830-211 radio frequency observations of methanol at a redshift of 0.88582 have established the most stringent limits on changes in the proton to electron mass ratio mu to date. The observations place the limit of (delta mu)/mu </= (0.0 +/- 1.0) x 10^{-7} which is approximately a factor of four lower than the previous lowest limit at a redshift of 0.6742. This stringent limit at a look back time of roughly half the age of the universe has profound implications for rolling scalar field cosmologies and the new physics that they require. Many of these cosmologies invoke a scalar field phi that is also coupled to the electromagnetic field causing the values of the fundamental constants, mu and the fine structure constant alpha to roll with time. If the lowest expected value of the coupling to mu, zeta_{mu}$ is invoked the new limit requires a limit on the dark energy equation of state parameter w such that w+1 </= 0.001 at a redshift of 0.88582. This eliminates almost all of the expected parameter space for such cosmologies and new physics that have a coupling to the electromagnetic field. In these cases the limit requires that w must be extremely close to -1 for the last half of the age of the universe or that the coupling of the rolling scalar field to mu and the electromagnetic field be significantly below or at the limit of its expected range. The new observations solidify the role of fundamental constants in providing probes of the possible cosmologies and new physics to explain the acceleration of the expansion of the universe.

We discuss scientific, technical and programmatic issues related to the use of an NRO 2.4m telescope for the WFIRST initiative of the 2010 Decadal Survey. We show that this implementation of WFIRST, which we call "NEW WFIRST," would achieve the goals of the NWNH Decadal Survey for the WFIRST core programs of Dark Energy and Microlensing Planet Finding, with the crucial benefit of deeper and/or wider near-IR surveys for GO science and a potentially Hubble-like Guest Observer program. NEW WFIRST could also include a coronagraphic imager for direct detection of dust disks and planets around neighboring stars, a high-priority science and technology precursor for future ambitious programs to image Earth-like planets around neighboring stars.

Many cosmological models invoke rolling scalar fields to account for the observed acceleration of the expansion of the universe. These theories generally include a potential V(phi) which is a function of the scalar field phi. Although V(phi) can be represented by a very diverse set of functions, recent work has shown the under some conditions, such as the slow roll conditions, the equation of state parameter w is either independent of the form of V(phi) or is part of family of solutions with only a few parameters. In realistic models of this type the scalar field couples to other sectors of the model leading to possibly observable changes in the fundamental constants such as the fine structure constant alpha and the proton to electron mass ratio mu. This paper explores the limits this puts on the validity of various cosmologies that invoke rolling scalar fields. We find that the limit on the variation of mu puts significant constraints on the product of a cosmological parameter w+1 times a new physics parameter zeta_mu^2, the coupling constant between mu and the rolling scalar field. Even when the cosmologies are restricted to very slow roll conditions either the value of zeta_mu must be at the lower end of or less than its expected values or the value of w+1 must be restricted to values vanishingly close to 0. This implies that either the rolling scalar field is very weakly coupled with the electromagnetic field, small zeta_mu, very weakly coupled with gravity, w+1 ~ 0 or both. These results stress that adherence to the measured invariance in mu is a very significant test of the validity of any proposed cosmology and any new physics it requires. The limits on the variation of mu also produces a significant tension with the reported changes in the value of alpha.

We present a data model describing the structure of spectrophotometric datasets with spectral and temporal coordinates and associated metadata. This data model may be used to represent spectra, time series data, segments of SED (Spectral Energy Distributions) and other spectral or temporal associations.

The Simple Spectral Access (SSA) Protocol (SSAP) defines a uniform interface to remotely discover and access one dimensional spectra. SSA is a member of an integrated family of data access interfaces altogether comprising the Data Access Layer (DAL) of the IVOA. SSA is based on a more general data model capable of describing most tabular spectrophotometric data, including time series and spectral energy distributions (SEDs) as well as 1-D spectra; however the scope of the SSA interface as specified in this document is limited to simple 1-D spectra, including simple aggregations of 1-D spectra. The form of the SSA interface is simple: clients first query the global resource registry to find services of interest and then issue a data discovery query to selected services to determine what relevant data is available from each service; the candidate datasets available are described uniformly in a VOTable format document which is returned in response to the query. Finally, the client may retrieve selected datasets for analysis. Spectrum datasets returned by an SSA spectrum service may be either precomputed, archival datasets, or they may be virtual data which is computed on the fly to respond to a client request. Spectrum datasets may conform to a standard data model defined by SSA, or may be native spectra with custom project-defined content. Spectra may be returned in any of a number of standard data formats. Spectral data is generally stored externally to the VO in a format specific to each spectral data collection; currently there is no standard way to represent astronomical spectra, and virtually every project does it differently. Hence spectra may be actively mediated to the standard SSA-defined data model at access time by the service, so that client analysis programs do not have to be familiar with the idiosyncratic details of each data collection to be accessed.

Changes in the values of the fundamental constants mu, the proton to electron mass ratio, and alpha, the fine structure constant due to rolling scalar fields have been discussed both in the context of cosmology and in new physics such as Super Symmetry (SUSY) models. This article examines the changes in these fundamental constants in a particular example of such fields, freezing and thawing slow roll quintessence. Constraints are placed on the product of a cosmological quantity, w, the equation of state parameter, and the square of the coupling constants for mu and alpha with the field, zeta_x, x = mu,alpha, using the existing observational limits on the values of Delta x/x. Various examples of slow rolling quintessence models are used to further quantify the constraints. Some of the examples appear to be rejected by the existing data which strongly suggests that conformation to the values of the fundamental constants in the early universe is a standard test that should be applied to any cosmological model or suggested new physics.

The data archives from space and ground-based telescopes present a vast opportunity for the astronomical community. We describe a classification encoding system for stellar spectra designed for archival databases that organizes the spectral data by "spectral classes." These classes are encoded into a digital format of the form TT.tt.LL.PPPP, where TT and tt refer to spectral type and subtype, LL to luminosity class, and PPPP to possible spectral peculiarities. Archive centers may wish to utilize this system to quantify classes of formerly arbitrary spectral classification strings found in classification catalogs corresponding to datasets of pointed spectroscopic observations in their holdings. The encoding system will also allow users to request archived data based on spectral class ranges, thereby streamlining an otherwise tedious data discovery process. Material in Appendix A is "normative" (part of the defined standard). Appendices B and C are "informative," meant to show how one data provider (MAST) has opted to handle some practical details.

We present the results of a numerical study comparing photometric and physical properties of simulated z=6-9 galaxies to the observations taken by the WFC3 instrument aboard the Hubble Space Telescope. Using cosmological hydrodynamical simulations we find good agreement with observations in color-color space at all studied redshifts. We also find good agreement between observations and our Schechter luminosity function fit in the observable range, Muv<= -18, provided that a moderate dust extinction effect exists for massive galaxies. However beyond what currently can be observed, simulations predict a very large number of low-mass galaxies and evolving steep faint-end slopes from alpha_L = -2.15 at z=6 to alpha_L = -2.64 at z=9, with a dependence of |alpha_L| \propto (1+z)^0.59. During the same epoch, the normalization phi* increases and the characteristic magnitude Muv* becomes moderately brighter with decreasing redshift. We find similar trends for galaxy stellar mass function with evolving low-mass end slope from alpha_M = - 2.26 at z=6 to alpha_M = -2.87 at z=9, with a dependence of |alpha_M| \propto (1+z)^0.65. Together with our recent result on the high escape fraction of ionizing photons for low-mass galaxies, our results suggest that the low-mass galaxies are important contributor of ionizing photons for the reionisation of the Universe at z>=6.

The existence of 1E0657-56 poses a challenge to the concordance Lambda cold dark matter model. Here we investigate the velocity distribution of dark matter halo pairs in large N-body simulations with differing box sizes (250Mpc/h-2 Gpc/h) and resolutions. We examine statistics such as the halo masses, pairwise halo velocities (v_{12}), and pair separation distances. We then compare our results to the initial conditions (ICs) required to reproduce the observational properties of 1E0657-56 in non-cosmological hydrodynamical simulations. We find that the high velocity tail of the v_{12} distribution extends to greater velocities as we increase the box size. We also find that the number of high-v_{12} pairs increases as we increase the particle count and resolution with a fixed box size, however, this increase is mostly due to lower mass halos which do not match the observed masses of 1E0657-56. We find that the redshift evolution is not very strong for the v_{12} distribution function between z=0.0 and z~0.5. We identify some pairs whose v_{12} resemble the required ICs, however, even the best candidates have either wrong halo mass ratios, or too large separations. Our simulations suggest that it is very difficult to produce such ICs at z=0.0, 0.296,& 0.489 in comoving volumes as large as (2Gpc/h)^3. Based on the extrapolation of our cumulative v_{12} function, we find that one needs a simulation with a comoving box size of (4.48Gpc/h)^3 and 2240^3 DM particles in order to produce at least one pair of halos that resembles the required v_{12} and observed masses of 1E0657-56. We find that the probability of finding a halo pair with v_{12}>=3000km/s and masses >=10^{14}Msun to be 2.76x10^{-8} at z=0.489. We conclude that either 1E0657-56 is incompatible with the concordance LCDM universe, or the ICs suggested by the non-cosmological simulations must be revised to give a lower value of v_{12}.

We present integral field spectroscopy in the near infrared of the nearby dwarf starburst galaxies NGC 3125-A and of the low metallicity dwarf galaxy SBS 0335-052. The use of adaptive optics in the observations produces sub-arcsecond angular resolution. We pinpoint the star forming cores of both galaxies, identify relevant ISM components such as dust, photo ionized gas, shock excited gas and molecular gas. We relate these components to the large scale star formation process of the galaxies. In particular we find the emission of the near infrared lines of H2 and especially [FeII] does not coincide with the HII region in NGC 3125. We have the first clear detection of [FeII] in SBS 0335-052.

Kepler's primary mission is a search for earth-size exoplanets in the habitable zone of late-type stars using the transit method. To effectively accomplish this mission, Kepler orbits the Sun and stares nearly continuously at one field-of-view which was carefully selected to provide an appropriate density of target stars. The data transmission rates, operational cycles, and target management requirements implied by this mission design have been optimized and integrated into a comprehensive plan for science operations. The commissioning phase completed all critical tasks and accomplished all objectives within a week of the pre-launch plan. Since starting science, the nominal data collection timeline has been interrupted by two safemode events, several losses of fine point, and some small pointing adjustments. The most important anomalies are understood and mitigated, so Kepler's technical performance metrics have improved significantly over this period and the prognosis for mission success is excellent. The Kepler data archive is established and hosting data for the science team, guest observers, and public. The first data sets to become publicly available include the monthly full-frame images, dropped targets, and individual sources as they are published. Data are released through the archive on a quarterly basis; the Kepler Results Catalog will be released annually starting in 2011.

In an effort to resolve the discrepancy between two measurements of the fundamental constant mu, the proton to electron mass ratio, at early times in the universe we reanalyze the same data used in the earlier studies. Our analysis of the molecular hydrogen absorption lines in archival VLT/UVES spectra of the damped Lyman alpha systems in the QSOs Q0347-383 and Q0405-443 yields a combined measurement of a (Delta mu)/mu value of (-7 +/- 8) x 10^{-6}, consistent with no change in the value of mu over a time span of 11.5 gigayears. Here we define (Delta mu) as (mu_z - mu_0) where mu_z is the value of mu at a redshift of z and mu_0 is the present day value. Our null result is consistent with the recent measurements of King et al. 2009, (Delta mu)/u = (2.6 +/- 3.0) x 10^{-6}, and inconsistent with the positive detection of a change in mu by Reinhold et al. 2006. Both of the previous studies and this study are based on the same data but with differing analysis methods. Improvements in the wavelength calibration over the UVES pipeline calibration is a key element in both of the null results. This leads to the conclusion that the fundamental constant mu is unchanged to an accuracy of 10^{-5} over the last 80% of the age of the universe, well into the matter dominated epoch. This limit provides constraints on models of dark energy that invoke rolling scalar fields and also limits the parameter space of Super Symmetric or string theory models of physics. New instruments, both planned and under construction, will provide opportunities to greatly improve the accuracy of these measurements.

This paper addresses alternative procedures to the ESO supplied pipeline procedures for the reduction of UVES spectra of two quasar spectra to determine the value of the fundamental constant mu = Mp/Me at early times in the universe. The procedures utilize intermediate product images and spectra produced by the pipeline with alternative wavelength calibration and spectrum addition methods. Spectroscopic studies that require extreme wavelength precision need customized wavelength calibration procedures beyond that usually supplied by the standard data reduction pipelines. An example of such studies is the measurement of the values of the fundamental constants at early times in the universe. This article describes a wavelength calibration procedure for the UV-Visual Echelle Spectrometer on the Very Large Telescope, however, it can be extended to other spectrometers as well. The procedure described here provides relative wavelength precision of better than 3E-7 for the long-slit Thorium-Argon calibration lamp exposures. The gain in precision over the pipeline wavelength calibration is almost entirely due to a more exclusive selection of Th/Ar calibration lines.

As one of the lowest metallicity star forming galaxies, with a nucleus of several super star clusters, SBS0335-052E is the subject of substantial current study. We present new insights on this galaxy based on new and archival high spatial resolution NICMOS and ACS images. We provide new measurements and limits on the size of several of the SSCs. The images have sufficient resolution to divide the star formation into compact regions and newly discovered extended regions, indicating a bi-modal form of star formation. The star formation regions are dated via the equivalent width of the Pa alpha emission and we find that two of the extended regions of star formation are less than 10 million years old. Our previous finding that stellar winds confine the photo-ionizing flux to small regions around individual stars is consistent with the new observations. This may allow planet formation in what would traditionally be considered a harsh environment and has implications for the number of planets around globular cluster stars. In addition the images pinpoint the regions of H2 emission as located in, but not at the center of the two star forming super star clusters, S1 and S2.

This paper addresses the nature of the near infrared background. We investigate whether there is an excess background at 1.4 microns, what is the source of the near infrared background and whether that background after the subtraction of all known sources contains the signature of high redshift objects (Z > 10). Based on NICMOS observations in the Hubble Ultra Deep Field and the Northern Hubble Deep Field we find that there is no excess in the background at 1.4 microns and that the claimed excess is due to inaccurate models of the zodiacal background. We find that the near infrared background is now spatially resolved and is dominated by galaxies in the redshift range between 0.5 and 1.5. We find no signature than can be attributed to high redshift sources after subtraction of all known sources either in the residual background or in the fluctuations of the residual background. We show that the color of the fluctuations from both NICMOS and Spitzer observations are consistent with low redshift objects and inconsistent with objects at redshifts greater than 10. It is most likely that the residual fluctuation power after source subtraction is due to the outer regions of low redshift galaxies that are below the source detection limit and therefore not removed during the source subtraction.

The Palomar Testbed Interferometer (PTI) archive of observations between 1998 and 2005 is examined for objects appropriate for calibration of optical long-baseline interferometer observations - stars that are predictably point-like and single. Approximately 1,400 nights of data on 1,800 objects were examined for this investigation. We compare those observations to an intensively studied object that is a suitable calibrator, HD217014, and statistically compare each candidate calibrator to that object by computing both a Mahalanobis distance and a Principal Component Analysis. Our hypothesis is that the frequency distribution of visibility data associated with calibrator stars differs from non-calibrator stars such as binary stars. Spectroscopic binaries resolved by PTI, objects known to be unsuitable for calibrator use, are similarly tested to establish detection limits of this approach. From this investigation, we find more than 350 observed stars suitable for use as calibrators (with an additional $\approx 140$ being rejected), corresponding to $\gtrsim 95%$ sky coverage for PTI. This approach is noteworthy in that it rigorously establishes calibration sources through a traceable, empirical methodology, leveraging the predictions of spectral energy distribution modeling but also verifying it with the rich body of PTI's on-sky observations.

We have resolved the classical nova V1663 Aql using long-baseline near-IR interferometry covering the period from 5--18 days after peak brightness. We directly measure the shape and size of the fireball, which we find to be asymmetric. In addition we measure an apparent expansion rate of 0.21 +/- 0.03 mas/day. Assuming a linear expansion model we infer a time of initial outburst approximately 4 days prior to peak brightness. When combined with published spectroscopic expansion velocities our angular expansion rate implies a distance of 8.9 +/- 3.6 kpc. This distance measurement is independent of, but consistent with, determinations made using widely available photometric relations for novae.

This letter extends our previous fluctuation analysis of the near infrared background at 1.6 microns to the 1.1 micron (F110W) image of the Hubble Ultra Deep field. When all detectable sources are removed the ratio of fluctuation power in the two images is consistent with the ratio expected for faint, z<8, sources, and is inconsistent with the expected ratio for galaxies with z>8. We also use numerically redshifted model galaxy spectral energy distributions for 50 and 10 million year old galaxies to predict the expected fluctuation power at 3.6 microns and 4.5 microns to compare with recent Spitzer observations. The predicted fluctuation power for galaxies at z = 0-12 matches the observed Spitzer fluctuation power while the predicted power for z>13 galaxies is much higher than the observed values. As was found in the 1.6 micron (F160W) analysis the fluctuation power in the source subtracted F110W image is two orders of magnitude below the power in the image with all sources present. This leads to the conclusion that the 0.8--1.8 micron near infrared background is due to resolved galaxies in the redshift range z<8, with the majority of power in the redshift range of 0.5--1.5.

We study the relation between the star formation intensity of galaxies and the extinction by dust of their emitted light. We employ a detailed statistical analysis of Hubble Deep Field North data to show a clear positive correlation between the extinction and star formation intensity at all epochs from redshift 0.4 to 6.5. The extinction evidently increases with time for a given star formation intensity, consistent with the expected increase in the metallicity with time. Our observational results are well fitted at all epochs by a double power-law model with a fixed shape that simply shifts with redshift. The correlation between the extinction and the star formation intensity can be interpreted by combining two other trends: the correlation between the star formation rate and the gas content of galaxies, and the evolution of the dust-to-gas ratio in galaxies. If we assume that Kennicutt's observed relation for the former is valid at each redshift, then our findings imply an interesting variation in the dust-to-gas ratio in galaxies within each epoch and with time, and suggest new ways to investigate the cosmic evolution of this quantity.

Using the CHARA Array and the Palomar Testbed Interferometer, the chemically peculiar star $\lambda$ Bo\"{o}tis has been spatially resolved. We have measured the limb darkened angular diameter to be $\theta_{LD} = 0.533\pm0.029$ mas, corresponding to a linear radius of $R_{\star} = 1.70 \pm 0.10 R_\odot$. The measured angular diameter yields an effective temperature for $\lambda$ Boo of $T_{eff} = 8887 \pm 242$ K. Based upon literature surface gravity estimates spanning $\log{(g)} = 4.0-4.2$ $[\rm{cm s}^{-\rm{2}}]$, we have derived a stellar mass range of $M_{\star} = 1.1 - 1.7$ $M_\odot$. For a given surface gravity, the linear radius uncertainty contributes approximately $\sigma(M_\star) = 0.1-0.2 M_\odot$ to the total mass uncertainty. The uncertainty in the mass (i.e., the range of derived masses) is primarily a result of the uncertainty in the surface gravity. The upper bound of our derived mass range ($\log(g)=4.2, M_\star = 1.7\pm0.2 M_\odot$) is consistent with 100-300 MYr solar-metallicity evolutionary models. The mid-range of our derived masses ($\log(g)=4.1, M_\star = 1.3\pm0.2 M_\odot$) is consistent with 2-3 GYr metal-poor evolutionary models. A more definitive surface gravity determination is required to determine a more precise mass for $\lambda$ Boo.

We report observations of the recurrent nova RS Oph using long-baseline near-IR interferometry. We are able to resolve emission from the nova for several weeks after the February 2006 outburst. The near-IR source initially expands to a size of approximately 5 milli-arcseconds. However, beginning around day 10 the IR source appears to begin to shrink, reaching approximately 2 milli-arcseconds by day 100. We combine our measured angular diameters with previously available interferometric and photometric data to derive an emission measure for the source, and hence are able to determine the mass-loss rate of the nova in the days following the outburst.

NICMOS observations of the resolved object fluxes in the Hubble Deep Field North and the Hubble Ultra Deep Field are significantly below the fluxes attributed to a 1.4 - 1.8 microns Near InfraRed Background Excess (NIRBE) from previous low spatial resolution NIRS measurements. Tests placing sources in the NICMOS image with fluxes sufficient to account for the NIRBE indicate that the NIRBE flux must be either flat on scales greater than 100 arc second or clumped on scales of several arc minutes to avoid detection in the NICMOS image. A fluctuation analysis of the new NICMOS data shows a fluctuation spectrum consistent with that found at the same wavelength in deep 2MASS calibration images. The fluctuation analysis shows that the majority of the fluctuation power comes from resolved galaxies at redshifts of 1.5 and less and that the fluctuations observed in the earlier deep 2MASS observations can be completely accounted for with normal low redshift galaxies. Neither the NICMOS direct flux measurements nor the fluctuation analysis require an additional component of near infrared flux other than the flux from normal resolved galaxies in the redshift range between 0 and 7. The residual fluctuations in the angular range between 1 and 10 arc seconds is 1-2 nW m-2 sr-1 which is at or above several predictions of fluctuations from high redshift population III objects, but inconsistent with attributing the entire NIRBE to high redshift galaxies.

During the last two years we have used the Palomar Testbed Interferometer to observe several explosive variable stars, including V838 Monocerotis, V1663 Aquilae and recently RS Ophiuchi. We observed V838 Monocerotis approximately 34 months after its eruption, and were able to resolve the ejecta. Observations of V1663 Aql were obtained starting 9 days after peak brightness and continued for 10 days. We were able to resolve the milliarcsecond-scale emission and follow the expansion of the nova photosphere. When combined with radial-velocity information, these observations can be used to infer the distance to the nova. Finally we have resolved the recurrent nova RS Oph and can draw some preliminary conclusions regarding the emission morphology.

We use the NICMOS Treasury and ACS HUDF images to measure the extinction corrected star formation history for 4681 galaxies in the region common to both images utilizing the star formation rate distribution function and other techniques similar to those employed with the NICMOS and WFPC2 images in the HDFN. Unlike the HDFN the NICMOS region of the HUDF appears to lack highly luminous and high star formation rate galaxies at redshifts beyond 3. The HUDF provides a region that is completely uncorrelated to the HDFN and therefore provides and independent measure of the star formation history of the universe. The combined HUDF and HDFN star formation rates show an average star formation rate of 0.2 solar masses per yer per cubic megaparsec. The average SFR of the combined fields at z = 1-3 is 0.29 solar masses per year per cubic megaparsec while the average at z = 4-6 is 1.2 solar masses per year per cubic megaparsec. The SFRs at all redshifts are within 3 sigma of the average over all redshifts.

We present new K-band long baseline interferometer observations of three young stellar objects of the FU Orionis class, V1057 Cyg, V1515 Cyg and Z CMa-SE, obtained at the Keck Interferometer during its commissioning science period. The interferometer clearly resolves the source of near-infrared emission in all three objects. Using simple geometrical models we derive size scales (0.5-4.5 AU) for this emission. All three objects appear significantly more resolved than expected from simple models of accretion disks tuned to fit the broadband optical and infrared spectro-photometry. We explore variations in the key parameters that are able to lower the predicted visibility amplitudes to the measured levels, and conclude that accretion disks alone do not reproduce the spectral energy distributions and K-band visibilities simultaneously. We conclude that either disk models are inadequate to describe the near-infrared emission, or additional source components are needed. We hypothesize that large scale emission (10s of AU) in the interferometer field of view is responsible for the surprisingly low visibilities. This emission may arise in scattering by large envelopes believed to surround these objects.

SBS0335-052 is a well studied Blue Compact Dwarf galaxy with one of the lowest metallicities of any known galaxy. It also contains 6 previously identified Super Star Clusters. We combine archival HST NICMOS images in the Pa alpha line and the 1.6 micron continuum of the eastern component, SBS0335-052E, with other space and ground based data to perform a multi-wavelength analysis of the super star clusters. We concentrate on the southern most clusters, designated S1 and S2, which appear to be the youngest clusters and are the strongest emitters of Pa alpha, radio, and x-ray flux. Our analysis leads to a possible model for S1 and perhaps S2 as a cluster of very young, massive stars with strong stellar winds. The wind density can be high enough to absorb the majority of ionizing photons within less than 1000 AU of the stars, creating very compact HII regions that emit optically thick radiation at radio wavelengths. These winds would then effectively quench the photoionizing flux very close to the stars. This can delay the onset of negative feedback by photoionization and photodissociation on star formation in the clusters. This is significant since SBS0335-052E resembles the conditions that were probably common for high redshift star formation in galaxies near the epoch of reionization.

We present observations of the A7IV-V star Alderamin ($\alpha$ Cep, HR 8162, HD 203280) from the Georgia State University CHARA Array. These infrared interferometric angular size measurements indicate a non-circular projected disk brightness distribution for this known rapid rotator. The interferometric observations are modeled as arising from an elongated rigid atmosphere, with apparent polar and equatorial radii of $r_p=0.6753^{+0.0119}_{-0.0135}$ milliarcseconds (mas) and $r_e=0.8767^{+0.0293}_{-0.0183}$ mas, respectively, for a difference of $201\pm 32$ microarcseconds ($\mu$as), and with an axial ratio of $r_e/r_p =1.298\pm 0.051$. Using the Hipparcos distance of $14.96\pm 0.11$ pc, these angular measures translate to $2.18\pm0.05$ and $2.82\pm0.10$ $R_\odot$. The inclination of Alderamin to the line of sight indicated by this modeling is effectively edge-on ($i=88.2^{+1.8}_{-13.3}$). The star has a true rotational velocity of $283\pm 10$ km/s ($\sim83$% of breakup velocity), and a polar temperature of roughly 8400 K. Significantly, a necessary aspect of this modeling is a determination of the gravity darkening coefficient, which at a value of $\beta=0.084^{+0.026}_{-0.049}$ is consistent with a convective photosphere, as expected for an A7IV-V star. Our detailed characterization of this object allows us to investigate various scenarios for the angular momentum history of Alderamin and the appropriateness of certain stellar evolution models.

We present near infrared broad band and polarimetric images of the compact star forming cluster AFGL437 obtained with the NICMOS instrument aboard HST. Our high resolution images reveal a well collimated bipolar reflection nebulosity in the cluster and allow us to identify WK34 as the illuminating source. The scattered light in the bipolar nebulosity centered on this source is very highly polarized (up to 79%). Such high levels of polarization implies a distribution of dust grains lacking large grains, contrary to the usual dust models of dark clouds. We discuss the geometry of the dust distribution giving rise to the bipolar reflection nebulosity and make mass estimates for the underlying scattering material. We find that the most likely inclination of the bipolar nebulosity, south lobe inclined towards Earth, is consistent with the inclination of the large scale CO molecular outflow associated with the cluster, strengthening the identification of WK34 as the source powering it.

This paper describes the observations and data reduction techniques for the version 2.0 images and catalog of the NICMOS Ultra Deep Field Treasury program. All sources discussed in this paper are based on detections in the combined NICMOS F110W and F160W bands only. The NICMOS images are drizzled to 0.09 arc second pixels and aligned to the ACS UDF F850LP image which was rebinned to the same pixel scale. These form the NICMOS version 2.0 UDF images. The catalog sources are chosen with a conservative detection limit to avoid the inclusion of numerous spurious sources. The catalog contains 1293 objects in the 144 x 144 arc sececonds NICMOS subfield of the UDF. The 5 sigma signal to noise level is an average 0.6 arc second diameter aperture AB magnitude of ~27.7 at 1.1 and 1.6 microns. The catalog sources, listed in order of right ascension, satisfy a minimum signal to noise criterion of 1.4 sigma in at least 7 contiguous pixels of the combined F110W and F160W image

We take advantage of the Hubble Ultra Deep Field (UDF) data to study the restframe optical and ultra violet (UV) morphologies of the novel population of Distant Red Galaxies (DRGs). Six galaxies with J-Ks > 2.3 are found to Ks=21.5, five of which have photometric redshifts z_phot > 2, corresponding to a surface density of 0.9/arcmin^2. The surface brightness distributions of the z_phot > 2 galaxies are better represented by exponential disks than R^{1/4}-laws. Two of the z_phot > 2 galaxies are extended, while three have compact morphologies. The restframe optical morphology of the z_phot > 2 galaxies is quite different from the restframe UV morphology: all the galaxies have red central components which dominate in the NICMOS H_{160}-band images, and distinct off-center blue features which show up in (and often dominate) the ACS images. The mean measured effective radius of the z_phot > 2 galaxies is <r_e> =1.9+/-1.4 kpc, similar (within the errors) to the mean size of LBGs at similar redshifts. All the DRGs are resolved in the ACS images, while four are resolved in the NICMOS images. Two of the z_phot > 2 galaxies are bright X-ray sources and hence host AGN. The diverse restframe optical and UV morphological properties of DRGs derived here suggest that they have complex stellar populations, consisting of both evolved populations that dominate the mass and the restframe optical light, and younger populations, which show up as patches of star formation in the restframe UV light; in many ways resembling the properties of normal local galaxies. This interpretation is supported by fits to the broadband SEDs, which for all five z_phot > 2 are best represented by models with extended star formation histories and substantial amounts of dust.

We use all available fields with deep NICMOS imaging to search for J dropouts (H<28) at z~10. Our primary data set for this search were the two J+H NICMOS parallel fields taken with the ACS HUDF. The 5 sigma limiting mags were 28.6 in J and 28.5 in H. Several shallower fields were also used: J+H NICMOS frames available over the HDF North, the HDF South NICMOS parallel, and the ACS HUDF. The primary selection criterion was (J-H)>1.8. 11 such sources were found in all search fields using this criterion. 8 of these were clearly ruled out as credible z~10 sources, either as a result of detections (>2 sigma) blueward of J or their colors redward of the break (H-K~1.5). The nature of the 3 remaining sources could not be determined from the data. The number appears consistent with the expected contamination from low-z interlopers. Analysis of the stacked images for the 3 candidates also suggests contamination. Regardless of their true redshifts, the actual number of z~10 sources must be <=3. To assess the significance of these results, two lower redshift samples (a z~3.8 B-dropout and z~6 i-dropout sample) were projected to z~8-12 using a (1+z)^{-1} size scaling. They were added to the image frames, and the selection repeated, giving 15.6 and 4.8 J-dropouts, respectively. This suggests that to the limit of this probe (0.3 L*) there has been evolution from z~3.8 and possibly from z~6. This is consistent with the strong evolution already noted at z~6 and z~7.5 relative to z~3-4. Even assuming that 3 sources from this probe are at z~10, the rest-frame continuum UV (~1500 A) luminosity density at z~10 (integrated down to 0.3 L*) is just 0.19_{-0.09}^{+0.13}x that at z~3.8 (or 0.19_{-0.10}^{+0.15}x including cosmic variance). However, if none of our sources is at z~10, this ratio has a 1 sigma upper limit of 0.07. (abridged)

We present new optical, near infrared, and archival ultraviolet observations of XX Ophiuchi and AS 325, two proposed iron stars. These unusual stars have optical spectra dominated by emission lines arising from hydrogen as well as ionized metals such as iron, chromium, and titanium. Both stars have been classified as iron stars and a number of exotic models have been presented for their origin. Using two years of moderately high resolution optical spectroscopy, the first high signal-to-noise K-band spectroscopy of these sources (which reveals stellar photospheric absorption lines), and new near-infrared interferometric observations, we confirm that both systems are composed of two stars, likely binaries, containing a hot Be star with an evolved late-type secondary. The hydrogen emission features arise in the hot wind from the Be star while the corresponding P-Cygni absorption lines are produced from dense material in the expanding, radiation driven, wind around each system. The optical FeII emission lines are pumped by ultraviolet FeII absorption lines through fluorescence. Contrary to some claims in the literature, the spectral features of XX Oph and AS 325 are quite similar, evidence that they are comparable systems. We examine the variability of the spectral morphology and radial velocity motions of both sources. We also study the variability of XX Oph during a major photometric event and find that the spectral nature of the system varies during the event. A comparison of the velocity of the absorption line components in our new spectra with those in the literature show that the structure of the stellar wind from XX Oph has changed since the system was observed in 1951.

We have used long-baseline near-IR interferometry to resolve the peculiar eruptive variable V838 Mon and to provide the first direct measurement of its angular size. Assuming a uniform disk model for the emission we derive an apparent angular diameter at the time of observations (November-December 2004) of $1.83 \pm 0.06$ milli-arcseconds. For a nominal distance of $8\pm2$ kpc, this implies a linear radius of $1570 \pm 400 R_{\odot}$. However, the data are somewhat better fit by elliptical disk or binary component models, and we suggest that the emission may be strongly affected by ejecta from the outburst.

We report the results of a sensitive K-band survey of Herbig Ae/Be disk sizes using the 85-m baseline Keck Interferometer. Targets were chosen to span the maximum range of stellar properties to probe the disk size dependence on luminosity and effective temperature. For most targets, the measured near-infrared sizes (ranging from 0.2 to 4 AU) support a simple disk model possessing a central optically-thin (dust-free) cavity, ringed by hot dust emitting at the expected sublimation temperatures (T_sub~1000-1500K). Furthermore, we find a tight correlation of disk size with source luminosity R propto L^(1/2) for Ae and late Be systems (valid over more than 2 decades in luminosity), confirming earlier suggestions based on lower-quality data. Interestingly, the inferred dust-free inner cavities of the highest luminosity sources (Herbig B0-B3 stars) are under-sized compared to predictions of the optically-thin cavity model, likely due to optically-thick gas within the inner AU.

In the past year, the Multimission Archive at the Space Telescope Science Institute (MAST) has taken major steps in making MAST's holdings available using VO-defined protocols and standards, and in implementing VO-based tools. For example, MAST has implemented the Simple Cone Search protocol, and all MAST mission searches may be returned in the VOTable format, allowing other archives to use MAST data for their VO applications. We have made many of our popular High Level Science Products available through Simple Image Access Protocol (SIAP), and are implementing the VO Simple Spectral Access Protocol (SSAP). The cross correlation of VizieR catalogs with MAST missions is now possible, and illustrates the integration of VO services into MAST. The user can easily display the results from searches within MAST using the plotting tool VOPlot. MAST also participates in the NVO registry service. Thus, the user can harvest MAST holdings simultaneously with data from many other surveys and missions through the VO DataScope Data Inventory Service.

We have detected likely z~7-8 galaxies in the 144"x144" NICMOS observations of the Hubble Ultra Deep Field. Objects are required to be >3 sigma detections in both NICMOS bands, J and H. The selection criteria for this sample are (z-J)>0.8, (z-J) > 0.66 (J-H)+0.8, (J-H)<1.2, and no detection at <8500 Angstroms. The 5 selected sources have total magnitudes H~27. Four of the five sources are quite blue compared to typical lower-redshift dropout galaxies and are clustered within a 1 arcmin**2 region. Because all 5 sources are near the limit of the NICMOS data, we have carefully evaluated their reality. Each of the candidates is visible in different splits of the data and a median stack. We analyzed several noise images and estimate the number of spurious sources to be 1+/-1. A search using an independent reduction of this same data set clearly revealed 3 of the 5 candidates and weakly detected a 4th candidate, suggesting the contamination could be higher. For comparison with predictions from lower redshift samples we take a conservative approach and adopt four z~7-8 galaxies as our sample. With the same detection criteria on simulated datasets, assuming no-evolution from z~3.8, we predict 10 sources at z~7-8, or 14 if we use a more realistic (1+z)^{-1} size scaling. We estimate that the rest-frame continuum UV (~1800 Angstroms) luminosity density at z~7.5 (integrated down to 0.3L*) is just 0.20_{-0.08}^{+0.12}x that found at z~3.8. Effectively this sets an upper limit on the luminosity density down to 0.3 L*. This result is consistent with significant evolution at the bright end of the luminosity function from z~7.5 to z~3.8. Even with the lower UV luminosity density at z~7.5, it appears that galaxies could still play an important role in reionization at these redshifts, though definitive measurements remain to be made. [abridged]

We report on the i-dropouts detected in two exceptionally deep ACS fields (B_{435}, V_{606}, i_{775}, and z_{850} with 10 sigma limits of 28.8, 29.0, 28.5, and 27.8, respectively) taken in parallel with the UDF NICMOS observations. Using an i-z>1.4 cut, we find 30 i-dropouts over 21 arcmin^2 down to z_AB=28.1, or 1.4 i-dropouts arcmin^{-2}, with significant field-to-field variation (as expected from cosmic variance). This extends i-dropout searches some ~0.9^m further down the luminosity function than was possible in the GOODS field, netting a ~7x increase in surface density. An estimate of the size evolution for UV bright objects is obtained by comparing the composite radial flux profile of the bright i-dropouts (z<27.2) with scaled versions of the HDF-N + HDF-S U-dropouts. The best-fit is found with a (1+z)^{-1.57_{-0.53} ^{+0.50}} scaling in size (for fixed luminosity), extending lower redshift (1<z<5) trends to z~6. Adopting this scaling and the brighter i-dropouts from both GOODS fields, we make incompleteness estimates and construct a z~6 LF in the rest-frame continuum UV (~1350 A) over a 3.5 magnitude baseline, finding a shape consistent with that found at lower redshift. To evaluate the evolution in the LF from z~3.8, we make comparisons against different scalings of a lower redshift B-dropout sample. Though a strong degeneracy is found between luminosity and density evolution, our best-fit model scales as (1+z)^{-2.8} in number and (1+z)^0.1 in luminosity, suggesting a rest-frame continuum UV luminosity density at z~6 which is just 0.38_{-0.07} ^{+0.09}x that at z~3.8. Our inclusion of size evolution makes the present estimate lower than previous z~6 estimates.

We report novel, high-angular resolution interferometric measurements that imply the near-infrared nuclear emission in NGC 4151 is unexpectedly compact. We have observed the nucleus of NGC 4151 at 2.2 microns using the two 10-meter Keck telescopes as an interferometer and find a marginally resolved source ~0.1 pc in diameter. Our measurements rule out models in which a majority of the K band nuclear emission is produced on scales larger than this size. The interpretation of our measurement most consistent with other observations is that the emission mainly originates directly in the central accretion disk. This implies that AGN unification models invoking hot, optically thick dust may not be applicable to NGC 4151.

The original analysis of the star formation history in the NICMOS Deep images of the NHDF is extended to the entire NHDF utilizing NICMOS and WFPC2 archival data. The roughly constant star formation rate from redshifts 1 to 6 found in this study is consistent with the original results. Star formation rates from this study, Lyman break galaxies and sub-mm observations are now in concordance The spike of star formation at redshift 2 due to 2 ULIRGs in the small Deep NICMOS field is smoothed out in the larger area results presented here. The larger source base of this study allows comparison with predictions from hierarchical galaxy formation models. In general the observation are consistent with the predictions. The observed luminosity functions at redshifts 1-6 are presented for future comparisons with theoretical galaxy evolution calculations. Mid and far infrared properties of the sources are also calculated and compared with observations. A candidate for the VLA source VLA 3651+1221 is discussed.

We present results on two related topics: 1. A discussion of high redshift candidates (z>4.5), and 2. A study of very small galaxies at intermediate redshifts, both sets being detected in the region of the northern Hubble Deep Field covered by deep NICMOS observations at 1.6 and 1.1 microns. The high redshift candidates are just those with redshift z>4.5 as given in the recent catalog of Thompson, Weymann and Storrie-Lombardi, while the ``small galaxy'' sample is defined to be those objects with isophotal area <= 0.2 squ. arcsec and with photometric redshifts 1<z<4.5. Of the 19 possible high redshift candidates listed in the Thompson et al. catalog, 11 have (nominal) photometric redshifts less than 5.0. Of these, however, only 4 are ``robust'' in the sense of yielding high redshifts when the fluxes are randomly perturbed with errors comparable to the estimated measuring error in each wave band. For the 8 other objects with nominal photometric redshifts greater than 5.0, one (WFPC2 4--473) has a published spectroscopic redshift. Of the remaining 7, 4 are robust in the sense indicated above. Two of these form a close pair (NIC 586 and NIC 107). The redshift of the object having formally the highest redshift, at 6.56 (NIC118 = WFPC2 4--601), is problematic, since F606W and F814W flux are clearly present, and the nature of this object poses a dilemma. (abridged)

This letter matches the shape of the star formation intensity distribution function to empirical laws such as the Schmidt law. The shape of the distribution at a redshift of one is reproduced from the empirical Schmidt law with a critical density, a Schechter distribution of galaxy masses and the assumption that star formation occurs mainly in exponential disks. The shape of the distribution depends primarily on two values, the characteristic mass m* in the Schechter mass distribution and the characteristic radius re in the exponential disk. As these characteristic values evolve they will affect the shape of the distribution function. The expected direction of evolution of the parameters partially cancels each other leaving the distribution shape relatively invariant.

We report on an observing campaign in March 2001 to monitor the brightness of the later arriving Q0957+561 B image in order to compare with the previously published brightness observations of the (first arriving) A image. The 12 participating observatories provided 3543 image frames which we have analyzed for brightness fluctuations. From our classical methods for time delay determination, we find a 417.09 +/- 0.07 day time delay which should be free of effects due to incomplete sampling. During the campaign period, the quasar brightness was relatively constant and only small fluctuations were found; we compare the structure function for the new data with structure function estimates for the 1995--6 epoch, and show that the structure function is statistically non-stationary. We also examine the data for any evidence of correlated fluctuations at zero lag. We discuss the limits to our ability to measure the cosmological time delay if the quasar's emitting surface is time resolved, as seems likely.

We present angular size measurements of 22 oxygen-rich Mira variable stars. These data are part of a long term observational program using the Infrared Optical Telescope Array (IOTA) to characterize the observable behavior of these stars. Complementing the infrared angular size measurements, values for variable star phase, spectral type, bolometric flux and distance were established for stars in the sample; flux and distance led to values for effective temperature, and linear radius, respectively. Additionally, values for the K-[12] color excess were established for these stars, which is indicative of dusty mass loss. Stars with higher color excess are shown to be systematically 120 solar radii larger than their low color excess counterparts, regardless of period. This analysis appears to present a solution to a long-standing question presented by the evidence that some Mira angular diameters are indicative of first overtone pulsation, while other diameters are more consistent with fundamental pulsation. A simple examination of the resultant sizes of these stars in the context of pulsation mode is consistent with at least some of these objects pulsating in the fundamental mode.

To empirically calibrate the IR surface brightness fluctuation (SBF) distance scale and probe the properties of unresolved stellar populations, we measured fluctuations in 65 galaxies using NICMOS on the Hubble Space Telescope. The early-type galaxies in this sample include elliptical and S0 galaxies and spiral bulges in a variety of environments. Absolute fluctuation magnitudes in the F160W (1.6 micron) filter were derived for each galaxy using previously-measured I-band SBF and Cepheid variable star distances. F160W SBFs can be used to measure distances to early-type galaxies with a relative accuracy of ~10% provided that the galaxy color is known to ~0.035 mag or better. Near-IR fluctuations can also reveal the properties of the most luminous stellar populations in galaxies. Comparison of F160W fluctuation magnitudes and optical colors to stellar population model predictions suggests that bluer elliptical and S0 galaxies have significantly younger populations than redder ones, and may also be more metal-rich. There are no galaxies in this sample with fluctuation magnitudes consistent with old, metal-poor (t>5 Gyr, [Fe/H]<-0.7) stellar population models. Composite stellar population models imply that bright fluctuations in the bluer galaxies may be the result of an episode of recent star formation in a fraction of the total mass of a galaxy. Age estimates from the F160W fluctuation magnitudes are consistent with those measured using the H-beta Balmer line index. The two types of measurements make use of completely different techniques and are sensitive to stars in different evolutionary phases. Both techniques reveal the presence of intermediate-age stars in the early-type galaxies of this sample.

We present high-resolution imaging in the rest-frame optical of the weak radio source LBDS53W091. Previous optical spectroscopy has shown that this object has an evolved stellar population of age >3Gyr at z=1.55, determined from the amplitude of rest-frame UV spectral breaks. We have obtained deep Hubble Space Telescope imaging over 10 orbits with NICMOS camera 2, using the F160W H-band filter (1.6microns) which is a good approximation to the rest-frame R-band. Our observations reveal a radial light profile which is well fit by a de Vaucouleurs' r^{1/4} law, with a scale length of r_e=0.3" (3kpc/h_{50} for Omega_M=0.3 and Omega_{Lambda}=0). The elliptical morphology of the radio galaxy indicates a dynamically-evolved old system, consistent with the spectroscopic results. Some surrounding objects lie on the color:magnitude relation for a cluster at z=1.55, and are likely to be associated. The group of galaxies are somewhat more luminous than the fundamental plane of ellipticals at z=0 subject to (1+z)^4 surface brightness dimming, but are consistent with estimates of the fundamental plane at high redshift subject to passive luminosity evolution from a formation epoch of z>3.

M16=NGC 6611, the Eagle Nebula, is a well studied region of star formation and the source of a widely recognized Hubble Space Telescope (HST) image. High spatial resolution infrared observations with the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on HST reveal the detailed morphology of two embedded star formation regions that are heavily obscured at optical wavelengths. It is striking that only limited portions of the visually obscured areas are opaque at 2.2 microns. Although the optical images imply substantial columns of material, the infrared images show only isolated clumps of dense gas and dust. Rather than being an active factory of star production, only a few regions are capable of sustaining current star formation. Most of the volume in the columns may be molecular gas and dust, protected by capstones of dense dust. Two active regions of star formation are located at the tips of the optical northern and central large ``elephant trunk'' features shown in the WFPC2 images. They are embedded in two capstones of infrared opaque material that contains and trails behind the sources. Although the presence of these sources was evident in previous observations at the same and longer wavelengths, the NICMOS images provide a high resolution picture of their morphology. Two bright stars appear at the tip of the southern column and may be the result of recent star formation at the top of that column. These observations suggest that the epoch of star formation in M16 may be near its endpoint.

(first paragraph) We have obtained Halpha and Palpha emission line images covering the central 3 - 4 arcmin of M51 using the WFPC2 and NICMOS cameras on HST to study the high-mass stellar population. The 0.1 - 0.2 arcsec pixels provide 4.6 - 9 pc resolution in M51 and the Halpha/Palpha line ratios are used to obtain extinction estimates. A sample of 1373 Halpha emission regions is catalogued using an automated and uniform measurement algorithm. Their sizes are typically 10 - 100 pc. The luminosity function for the Halpha emission regions is obtained over the range L_{Halpha} = 10^{36} to 2 times 10^{39} erg s{-1}. The luminosity function is fit well by a power law with dN/dlnL proportional to L^{-1.01}). The power law is significantly truncated and no regions were found with observed L_{Halpha} above 2 times 10^{39} erg s^{-1} (uncorrected for extinction; the maximum seen in ground-based studies is approximately a factor of 5 higher, very likely due to blending of multiple regions). The extinctions derived here increase the maximum intrinsic luminosity to above 10^{40} erg s^{-1}). The logarithmically binned luminosity function is also somewhat steeper (alpha = -1.01) than that found ground-based imaging (alpha = -0.5 to -0.8) - probably also a result of our resolving regions which were blended in the ground-based images. The 2-point correlation function for the HII regions exhibits strong clustering on scales <= 2 arcsec or 96 pc.

This paper presents NICMOS images of the nucleus and emission line regions in NGC 1068. The location of the nucleus relative to the emission line features is established and the physics underlying the morphology is discussed.

We present infrared interferometric angular size measurements for the A7IV-V star Altair which indicate a non-circular projected disk brightness distribution. Given the known rapid rotation of this star, we model the data as arising from an elongated rigid atmosphere. To first order, an ellipse may be fit to our interferometric diameter measurements, with major and minor axes of 2a=3.461+-0.038 milliarcseconds (mas) and 2b=3.037+-0.069 mas, respectively, for a difference of 424+-79 microarcseconds between 2a and 2b, and with an axial ratio of a/b =1.140+-0.029. Assuming that the apparent oblateness of the photosphere is due to the star's rapid rotation, a more rigorous evaluation of the observation data in the context of a rigidly rotating Roche model shows that an estimate of v sin i = 210+-13 km/s can be derived that is independent of spectroscopic techniques. Also derived are values for the mean effective temperature, the mean linear radius, and an observational constraint upon the the relationship between rotation velocity and stellar inclination is established. Altair is the first main sequence star for which direct observations of an oblate photosphere have been reported, and the first star forwhich v sin i has been established from observations of the star's photospheric geometry.

Near-infrared (2.2 um) long baseline interferometric observations of Vega are presented. The stellar disk of the star has been resolved, and the data have been fitted with a limb darkened stellar disk of diameter Theta_{LD} = 3.28 +/- 0.01 mas. The derived effective temperature is T_eff = 9553 +/- 111 K. However, the residuals resulting from the stellar disk model appear to be significant and display organized structure. Instrumental artifacts, stellar surface structure, stellar atmosphere structure, and extended emission/scattering from the debris disk are discussed as possible sources of the residuals. While the current dataset cannot uniquely determine the origin of the residuals, the debris disk is found to be the most likely source. A simple debris disk model, with 3-6% of Vega's flux emanating from the disk at r < 4 AU, can explain the residuals.

We present photometric observations of an apparent Type Ia supernova (SN Ia) at a redshift of ~1.7, the farthest SN observed to date. SN 1997ff, was discovered in a repeat observation by the HST of the HDF-), and serendipitously monitored with NICMOS on HST throughout the GTO campaign. The SN type can be determined from the host galaxy type:an evolved, red elliptical lacking enough recent star formation to provide a significant population of core-collapse SNe. The class- ification is further supported by diagnostics available from the observed colors and temporal behavior of the SN, both of which match a typical SN Ia. The photo- metric record of the SN includes a dozen flux measurements in the I, J, and H bands spanning 35 days in the observed frame. The redshift derived from the SN photometry, z=1.7+/-0.1, is in excellent agreement with the redshift estimate of z=1.65+/-0.15 derived from the U_300,B_450,V_606,I_814,J_110,J_125,H_160, H_165,K_s photometry of the galaxy. Optical and near-infrared spectra of the host provide a very tentative spectroscopic redshift of 1.755. Fits to observations of the SN provide constraints for the redshift-distance relation of SNe~Ia and a powerful test of the current accelerating Universe hypothesis. The apparent SN brightness is consistent with that expected in the decelerating phase of the preferred cosmological model, Omega_M~1/3, Omega_Lambda~2/3. It is inconsistent with grey dust or simple luminosity evolution, candidate astro- physical effects which could mimic past evidence for an accelerating Universe from SNe Ia at z~0.5.We consider several sources of possible systematic error including lensing, SN misclassification, selection bias, and calibration errors. Currently, none of these effects appears likely to challenge our conclusions.

We present Hubble Space Telescope (HST) Near-Infrared Camera and Multi-object Spectrometer (NICMOS) observations of the reflection nebulosity associated with the T Tauri star HH 30. The images show the scattered light pattern characteristic of a highly inclined, optically thick disk with a prominent dustlane whose width decreases with increasing wavelength. The reflected nebulosity exhibits a lateral asymmetry in the upper lobe on the opposite side to that reported in previously published Wide Field Planetary Camera 2 (WFPC2) images. The radiation transfer model which most closely reproduces the data has a flared accretion disk with dust grains larger than standard interstellar medium grains by a factor of approximately 2.1. A single hotspot on the stellar surface provides the necessary asymmetry to fit the images and is consistent with previous modeling of the light curve and images. Photometric analysis results in an estimated extinction of Av>~80; however, since the photometry measures only scattered light rather than direct stellar flux, this a lower limit. The radiative transfer models require an extinction of Av = 7,900.

We use the ultraviolet and optical WFPC2 and near-infrared NICMOS images of the Hubble Deep Field North to measure and statistically compare an array of parameters for over 250 of the galaxies it contains. These parameters include redshift, rest-frame visible asymmetry and concentration, bolometric luminosity and extinction-corrected star formation rate. We find only one strong correlation, between bolometric luminosity and star formation rate, from which early-type galaxies noticeably deviate. When our asymmetry measurements are combined with those of a sample of nearby galaxies covering the full Hubble sequence, we find a weak correlation between redshift and rest-frame visible asymmetry, consistent with the qualitative evidence of galaxy morphological evolution from these and other deep Hubble Space Telescope images. The mean values of these asymmetry measurements show a monotonic increase with redshift interval over the range 0 < z < 2, increasing by a factor of approximately three. If this trend is real, it suggests that galaxy morphological evolution is a gradual process that is continuing through the present cosmological epoch. There is evidence that the dominant source of this evolution is the "minor" mergers of disk galaxies with smaller companions, which could also transform late-type spiral galaxies into early-type spiral galaxies. Interestingly, in contrast to local galaxies we find no correlations between star formation rate and either UV or visible asymmetry. This could arise if the star formation of high-redshift galaxies proceeds in episodes that are short (~ 100 Myr) relative to the time scales over which galaxy mergers produce strong asymmetries (~ 500 Myr), a result suggested by the high star formation rates of Lyman break galaxies.

We have carried out a high-resolution imaging search for the galaxy associated with the damped Lyman-alpha (DLA) absorber at z=1.859 toward the z_{em}=2.48 quasar QSO 1244+3443, using the HST and the NICMOS. Images were obtained in the broad filter F160W and the narrow filter F187N with camera 2 on NICMOS with the goal of detecting the rest-frame optical continuum and the H-alpha line emission from the DLA. After PSF subtraction, two weak features are seen at projected separations of 0.16-0.24" from the quasar. Parts of these features may be associated with the DLA absorber, although we cannot completely rule out that they could be artifacts of the point spread function (PSF). If associated with the DLA, the objects would be ~1-2 h_{70}^{-1} kpc in size with integrated flux densities of 2.5 and 3.3 mu Jy in the F160W filter, implying luminosities at lambda_{central}=5600 A in the DLA rest frame of 4.4-5.9 x 10^{9} h_{70}^{-2} L_{solar} at z=1.86, for q0=0.5. However, no significant H-alpha line emission is seen from these objects, suggesting low star formation rates (SFRs). Our 3 sigma upper limit on the SFR in the DLA is 1.3 h_{70}^{-2} M_{solar}/yr for q0 = 0.5 (2.4 h_{70}^{-2} M_{solar} yr^{-1} for q0 = 0.1). This together with our earlier result for LBQS 1210+1731 mark a significant improvement over previous constraints on the star formation rates of DLAs. A combination of low SFR and some dust extinction is likely to be responsible for the lack of H-alpha emission. Alternatively, the objects, may be associated with the quasar host galaxy. In any case, our observations suggest that the DLA is not a large bright proto-disk, but a compact object or a low-surface brightness galaxy. If the two features are PSF artifacts then the constraints on DLA properties are even more severe.

We measured infrared surface brightness fluctuation (SBF) distances to an isotropically-distributed sample of 16 distant galaxies with redshifts reaching 10,000 km/s using the near-IR camera and multi-object spectrometer (NICMOS) on the Hubble Space Telescope (HST). The excellent spatial resolution, very low background, and brightness of the IR fluctuations yielded the most distant SBF measurements to date. Twelve nearby galaxies were also observed and used to calibrate the F160W (1.6 micron) SBF distance scale. Of these, three have Cepheid variable star distances measured with HST and eleven have optical I-band SBF distance measurements. A distance modulus of 18.5 mag to the Large Magellanic Cloud was adopted for this calibration. We present the F160W SBF Hubble diagram and find a Hubble constant Ho=76 +/- 1.3 (1-sigma statistical) +/- 6 (systematic) km/s/Mpc. This result is insensitive to the velocity model used to correct for local bulk motions. Restricting the fit to the six most distant galaxies yields the smallest value of Ho=72 +/- 2.3 km/s/Mpc consistent with the data. This 6% decrease in the Hubble constant is consistent with the hypothesis that the Local Group inhabits an under-dense region of the universe, but is also consistent with the best-fit value of Ho=76 km/s/Mpc at the 1.5-sigma level.

This paper presents the star formation history in the NICMOS Northern Deep HDF. It uses the techniques of photometric redshifts and extinctions to correct for extinction of the ultra-violet flux. It presents a new method for correcting for surface brightness diming. It also predicts the 850 micron fluxes of the objects for comparison with SCUBA measurements

We present a photometric analysis of all high Galactic lattitude (|b| > 20d) broad-band parallel images taken by the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) instrument during its initial lifetime in HST Cycle 7. These images, taken through the F110W and F160W filters, reach a mean 3-sigma limiting magnitude of approximately 22 in both bands, and cover a total area of approximately 92 square arcminutes. The reddest of the 358 galaxies detected have F110W-F160W colors and F160W magnitudes consistent with being a combination of both dusty star-forming and evolved early-type galaxies at 1 < z < 2. The surface density of these galaxies is comparable to that of the population of Extremely Red Objects (EROs) discovered in ground-based surveys (~ 100 deg-2), suggesting that EROs also represent a combination of both galaxy types in this redshift range. Roughly 10% of the detected galaxies appear to be blue compact dwarf galaxies at z < 1, a result consistent with studies of the HST Medium Deep Survey Fields. The surface density of these objects down to magnitude of 22 in F160W is ~ 300 deg-2. None of the 681 point sources detected have F110W-F160W colors matching those expected for QSOs whose continua have been significantly reddened by internal dust. Our data limit the surface density of such QSOs to be < 50 deg-2 down the mean limiting magnitudes of the sample images. Since the surface density of QSOs detected on the basis of ultraviolet excess in optical surveys to comparable depth is ~ 100 deg-2, this argues against the suggestion that dust-reddened QSOs comprise the undetected majority of the QSO population. The F110W-F160W color can also be used to identify unreddened QSOs at z > 8, but we find no such candidates. This is consistent with the evidence that QSO space density declines sharply at z > 5.

We have obtained images of the Trapezium Cluster (140" x 140"; 0.3 pc x 0.3 pc) with the Hubble Space Telescope Near-Infrared Camera and Multi-Object Spectrometer (NICMOS). Combining these data with new ground-based K-band spectra (R=800) and existing spectral types and photometry and the models of D'Antona & Mazzitelli, we find that the distributions of ages of comparable samples of stars in the Trapezium, rho Oph, and IC 348 indicate median ages of \~0.4 Myr for the first two regions and ~1-2 Myr for the latter. The low-mass IMFs in these sites of clustered star formation are similar over a wide range of stellar densities and other environmental conditions. With current data, we cannot rule out modest variations in the substellar mass functions among these clusters. We then make the best estimate of the true form of the IMF in the Trapezium by using the evolutionary models of Baraffe et al. and an empirically adjusted temperature scale and compare this mass function to recent results for the Pleiades and the field. All of these data are consistent with an IMF that is flat or rises slowly from the substellar regime to about 0.6 Msun, and then rolls over into a power law that continues from about 1 Msun to higher masses with a slope similar to or somewhat larger than the Salpeter value of 1.35. For the Trapezium, this behavior holds from our completeness limit of ~0.02 Msun and probably, after a modest completeness correction, even from 0.01-0.02 Msun. These data include ~50 likely brown dwarfs. We test the predictions of theories of the IMF against various properties of the observed IMF.

In this article we investigate the morphology and stellar populations of high-redshift galaxies through multi-waveband HST imaging and ground-based spatially-resolved spectroscopy. We study the redshift evolution of galaxy morphology in the Hubble Deep Field, using the deep IDT-NICMOS near-infrared HST imaging coupled with spectroscopic and photometric redshifts. Using the multi-waveband data to compare the appearance of galaxies at the same rest-frame wavelengths reveals that "morphological k-corrections" (the change in appearance when viewing high-z objects at shorter rest-frame wavelengths) are only important in a minority of cases, and that galaxies were intrinsically more peculiar at high redshift. One example of significant morphological k-corrections is spiral galaxies, which often show more pronounced barred structure in the near-infrared than in the optical. Therefore, the apparent decline in the fraction of barred spirals at faint magnitudes in the optical HDF may be due to band-shifting effects at the higher redshifts, rather than intrinsic evolution. Using such features as the age-sensitive Balmer+4000Ang break, the spatially- resolved colours of distant galaxies in optical/near-infrared imaging can also be used to study their component stellar populations. We supplement this with deep Keck/LRIS spectroscopy of two extended sources: a chain galaxy at z=2.8 (HDF4-555.1, the "Hot Dog" - the brightest U-drop Lyman-break galaxy in the HDF) and a pair of z=4.04 gravitationally lensed arcs behind the cluster Abell 2390. The absence of measurable rotation across the z=2.8 chain galaxy implies that it is unikely to be a disk viewed edge on. With the resolution enhancement from lensing, we detect stellar populations of different ages in the z=4 arcs.

We report results of a high-resolution imaging search (in rest frame H-$\alpha$ and optical continuum) for the galaxy associated with the damped Lyman-$\alpha$ (DLA) absorber at $z=1.892$ toward the $z_{em}=2.543$ quasar LBQS 1210+1731, using HST/NICMOS. After PSF subtraction, a feature is seen in both the broad-band and narrow-band images, at a projected separation of 0.25$\arcsec$ from the quasar. If associated with the DLA, the object would be $\approx 2-3$ $h_{70}^{-1}$ kpc in size with a flux of $9.8 \pm 2.4$ $\mu$Jy in the F160W filter, implying a luminosity at $\lambda_{central}=5500$ {\AA} in the rest frame of $1.5 \times 10^{10}$ $h_{70}^{-2}$ L$_{\odot}$ at $z=1.89$, for $q_{0}=0.5$. However, no significant H-$\alpha$ emission is seen, suggesting a low star formation rate (SFR) (3 $\sigma$ upper limit of 4.0 $h_{70}^{-2}$ M$_{\odot}$ yr$^{-1}$), or very high dust obscuration. Alternatively, the object may be associated with the host galaxy of the quasar. H-band images obtained with the NICMOS camera 2 coronagraph show a much fainter structure $\approx 4-5$ $h_{70}^{-1}$ kpc in size and containing four knots of continuum emission, located 0.7$\arcsec$ away from the quasar. We have probed regions far closer to the quasar sight-line than in most previous studies of high-redshift intervening DLAs. The two objects we report mark the closest detected high-redshift DLA candidates yet to any quasar sight line. If the features in our images are associated with the DLA, they suggest faint, compact, somewhat clumpy objects rather than large, well-formed proto-galactic disks or spheroids.

We present near-infrared images obtained with the HST NICMOS camera for a sample of 9 luminous (LIGs: L_IR (8-1000 microns) >~ 10^11 L_sun) and 15 ultra-luminous (ULIGS: L_IR >~ 10^12 L_sun) infrared galaxies. The sample includes representative systems classified as warm (f_25/f_60 > 0.2) and cold (f_25/f_60 <~ 0.2) based on the mid-infrared colors and systems with nuclear emission lines classified as HII (i.e. starburst), QSO, Seyfert and LINER. The morphologies of the sample galaxies are diverse and provide further support for the idea that they are created by the collision or interactions of spiral galaxies. Although no new nuclei are seen in the NICMOS images, the NICMOS images do reveal new spiral structures, bridges, and circumnuclear star clusters...

We present positions, magnitudes, sizes and morphological classifications for 111 galaxies discovered in the Hubble Space Telescope NICMOS Camera 1 and Camera 2 parallel fields. We combine the magnitudes measured in the JHK-analog filters with those from deep ground-based images in V and/or R to measure photometric redshifts for 71 objects using Bruzual-Charlot population synthesis models. We find that these objects fall in the range z ~ 0.0 - 2.7, with a mean redshift of 0.8 and a mean luminosity of 1.6 L*. The NICMOS images reveal many of the objects to be ordered spirals and ellipticals similar to those in the local universe. However, we find a higher fraction (~ 14%) of morphologically peculiar and/or interacting galaxies than is observed among local galaxies (~ 3% - 4%). This is consistent with the results from other deep HST images including the Hubble Deep Field and Hubble Medium Deep Survey Field that the fraction of peculiar and interacting galaxies increases with redshift. As the NICMOS images of the sample galaxies cover their rest-frame near infrared and optical emission, this result increases confidence that such changes in morphology are genuine, as opposed to an effect produced by viewing galaxies in the rest-frame ultraviolet. We also find that at least 26 of the sample galaxies appear to be members of (non-interacting) pairs or groups, based on their proximity to one another and photometric redshifts. This is consistent with the results of recent ground-based optical surveys covering larger areas, and with the detection of galaxy groups and filaments at redshifts higher than those covered by our sample.

The gamma-ray blazar 3C 279 was monitored on a nearly daily basis with IUE, ROSAT and EGRET for three weeks between December 1992 and January 1993. During this period, the blazar was at a historical minimum at all wavelengths. Here we present the UV data obtained during the above multiwavelength campaign. A maximum UV variation of ~50% is detected, while during the same period the X-ray flux varied by no more than 13%. At the lowest UV flux level the average spectrum in the 1230-2700 A interval is unusually flat for this object (<alpha_UV>~1). The flattening could represent the lowest energy tail of the inverse Compton component responsible for the X-ray emission, or could be due to the presence of a thermal component at ~20000 K possibly associated with an accretion disk. The presence of an accretion disk in this blazar object, likely observable only in very low states and otherwise hidden by the beamed, variable synchrotron component, would be consistent with the scenario in which the seed photons for the inverse Compton mechanism producing the gamma-rays are external to the relativistic jet. We further discuss the long term correlation of the UV flux with the X-ray and gamma-ray fluxes obtained at various epochs. All UV archival data are included in the analysis. Both the X- and gamma-ray fluxes are generally well correlated with the UV flux, approximately with square root and quadratic dependences, respectively.

We present the first results of a survey of random fields with the slitless G141 ($\lambda_c = 1.5\mu, \Delta\lambda=0.8\mu$) grism on NICMOS. Approximately 64 square arcminutes have been observed at intermediate and high galactic latitudes. The 3$\sigma$ limiting line and continuum fluxes in each field vary from $7.5 \times 10^{-17}$ to $1 \times 10^{-17} erg/cm^2/sec$ and from H = 20 to 22, respectively. Our median and area weighted $3\sigma $ limiting line fluxes within a 4 pixel aperture are nearly identical at $4.1 \times 10^{-17} erg/cm^2/sec$ and are 60% deeper than the deepest narrow-band imaging surveys from the ground. We have identified 33 emission-line objects and derive their observed wavelengths, fluxes and equivalent widths. We argue that the most likely line identification is H$\alpha$ and that the redshift range probed is from 0.75 to 1.9. The 2$\sigma$ rest-frame equivalent width limits range from 9\AA to 130\AA with an average of 40\AA. The survey probes an effective co-moving volume of $10^5 h_{50}^{-3} Mpc^3$ for $q_0=0.5$. Our derived co-moving number density of emission line galaxies in the range $0.7 < z < 1.9$ is $3.3\times10^{-4} h_{50}^{3} Mpc^{-3}$, very similar to that of the bright Lyman break objects at $z \sim 3$. The objects with detected emission-lines have a median F160W magnitude of 20.4 (Vega scale) and a median H$\alpha$ luminosity of $2.7 \times 10^{42} erg/sec$. The implied star formation rates range from 1 to 324 M_{\odot}/yr, with an average [NII]6583,6548 corrected rate of 21 M_{\odot}/yr for H_0=50 km/s/Mpc and $q_0=0.5$ (34 M_{\odot}/yr for $q_0=0.1$).

We report the first near infrared (NIR) imaging of a circumstellar annular disk around the young (~8 Myr), Vega-like star, HR 4796A. NICMOS coronagraph observations at 1.1 and 1.6 microns reveal a ring-like symmetrical structure peaking in reflected intensity 1.05 arcsec +/- 0.02 arcsec (~ 70 AU) from the central A0V star. The ring geometry, with an inclination of 73.1 deg +/- 1.2 deg and a major axis PA of 26.8 deg +/- 0.6 deg, is in good agreement with recent 12.5 and 20.8 micron observations of a truncated disk (Koerner, et al. 1998). The ring is resolved with a characteristic width of less than 0.26 arcsec (17 AU) and appears abruptly truncated at both the inner and outer edges. The region of the disk-plane inward of ~60 AU appears to be relatively free of scattering material. The integrated flux density of the part of the disk that is visible (greater than 0.65 arcsec from the star) is found to be 7.5 +/- 0.5 mJy and 7.4 +/- 1.2 mJy at 1.1 and 1.6 microns, respectively. Correcting for the unseen area of the ring yields total flux densities of 12.8 +/- 1.0 mJy and 12.5 +/- 2.0 mJy, respectively (Vega magnitudes = 12.92 /+- 0.08 and 12.35 +/-0.18). The NIR luminosity ratio is evaluated from these results and ground-based photometry of the star. At these wavelengths Ldisk(lambda)/L*(lambda) = 1.4 +/- 0.2E-3 and 2.4 +/- 0.5E-3, giving reasonable agreement between the stellar flux scattered in the NIR and that which is absorbed in the visible and re-radiated in the thermal infrared. The somewhat red reflectance of the disk at these wavelengths implies mean particle sizes in excess of several microns, larger than typical interstellar grains. The confinement of material to a relatively narrow annular zone implies dynamical constraints on the disk particles by one or more as yet unseen bodies.

We present the discovery of a candidate substellar object in a survey of young stars in the solar vicinity using the sensitivity and spatial resolution afforded by the NICMOS coronagraph on the Hubble Space Telescope. The H=12.1 mag object was discovered approximately 2 arcsec from the TW Hydrae Association member CoD -33 7795 (TWA 5), and the photometry implies a spectral type M8-M8.5, with a temperature of approximately 2600K. We estimate that the probability of a chance alignment with a background objectof this nature is < 2 X 10^(-5), and therefore postulate the object (TWA 5B) is physically associated at a projected separation of 100 AU. Given the likely youth of the primary (~ 10 Myr), current brown dwarf cooling models predict a mass of approximately 20 Jupiter masses for TWA 5B.

This paper contains a catlog of sources observed in a portion of the Hubble Deep Field by the NICMOS instrument. Methods of observation and data analysis are discussed. The catalog has two parts. The first is a listing of all sources with signal to noise ratios greater than or equal to 2.5. The second part contains sources with signal to noise ratios less than 2.5. The second part is intended for researchers who wish to know the limits to which sources may appear at a location.

We present Keck LRIS spectroscopy along with NICMOS F110W (~J) and F160W (~H) images of the galaxy HDF4-473.0 (hereafter 4-473) in the Hubble Deep Field, with a detection of an emission line consistent with Ly-alpha at a redshift of z=5.60. Attention to this object as a high redshift galaxy was first drawn by Lanzetta, Yahil and Fernandez-Soto and appeared in their initial list of galaxies with redshifts estimated from the WFPC2 HDF photometry. It was selected by us for spectroscopic observation, along with others in the Hubble Deep Field, on the basis of the NICMOS F110W and F160W and WFPC2 photometry. For H_0 = 65 and q_0 = 0.125, use of simple evolutionary models along with the F814W (~I), F110W, and F160W magnitudes allow us to estimate the star formation rate (~13 M(solar)/yr). The colors suggest a reddening of E(B-V) ~ 0.06. The measured flux in the Ly-alpha line is approximately 1.0*10^(-17) ergs/cm/s and the restframe equivalent width, correcting for the absorption caused by intervening HI, is approximately 90AA. The galaxy is compact and regular, but resolved, with an observed FWHM of ~0.44". Simple evolutionary models can accurately reproduce the colors and these models predict the Ly-alpha flux to within a factor of 2. Using this object as a template shifted to higher redshifts, we calculate the magnitudes through the F814W and two NICMOS passbands for galaxies at redshifts 6 < z < 10.

We report high resolution imaging of the ultraluminous infrared galaxy Arp 220 at 1.1, 1.6, and 2.22 microns with NICMOS on the HST. The diffraction-limited images at 0.1--0.2 arcsecond resolution clearly resolve both nuclei of the merging galaxy system and reveal for the first time a number of luminous star clusters in the circumnuclear envelope. The morphologies of both nuclei are strongly affected by dust obscuration, even at 2.2 microns : the primary nucleus (west) presents a crescent shape, concave to the south and the secondary (eastern) nucleus is bifurcated by a dust lane with the southern component being very reddened. In the western nucleus, the morphology of the 2.2 micron emission is most likely the result of obscuration by an opaque disk embedded within the nuclear star cluster. The morphology of the central starburst-cluster in the western nucleus is consistent with either a circumnuclear ring of star formation or a spherical cluster with the bottom half obscured by the embedded dust disk. Comparison of cm-wave radio continuum maps with the near-infrared images suggests that the radio nuclei lie in the dust disk on the west and near the highly reddened southern component of the eastern complex. The radio nuclei are separated by 0.98 arcseconds (corresponding to 364 pc at 77 Mpc) and the half-widths of the infrared nuclei are approximately 0.2-0.5 arcseconds. At least 8, unresolved infrared sources -- probably globular clusters -- are also seen in the circumnuclear envelope at radii 2-7 arcseconds . Their near-infrared colors do not significantly constrain their ages.

We present the first high-resolution near-infrared images of the edge-on silhouette circumstellar disk, Orion 114-426, made using NICMOS on the Hubble Space Telescope. Images taken against the bright nebular background of the ionized hydrogen Pa$\alpha$ line at 1.87 micron show the major axis of the disk to be approximately 20% smaller than at 0.6 micron, from which we deduce the structure of the edge of the disk. Continuum images of diffuse polar lobes above and below the plane of the disk show a morphology and evolution with wavelength consistent with predictions for reflection nebulae in a diffuse envelope with large polar cavities, surrounding a thin, massless, Keplerian disk, centered on an otherwise hidden central star. We make use of our observations and reasonable assumptions about the underlying disk structure to show that the disk mass is at least 10 earth masses and plausibly $\geq 5\times 10^{-4}$ solar masses.