Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravitohydrodynamics 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-2 J. KIM ET AL. FOR THE AGORA COLLABORATION 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.

Dwarf spheroidal galaxies (dSphs) are extremely gas-poor, dark matter-dominated galaxies, which make them ideal to test the predictions of the cold dark matter (CDM) model. We argue that the removal of the baryonic component from gas-rich dwarf irregular galaxies, the progenitors of dSphs, can substantially reduce their central density. Thus, it may play an important role in alleviating one of the problems of the CDM model related to the structure of relatively massive satellite galaxies of the Milky Way (MW). Traditionally, collisionless cosmological N-body simulations are used when confronting theoretical predictions with observations. However, these simulations assume that the baryon fraction everywhere in the Universe is equal to the cosmic mean, an assumption which is incorrect for dSphs. We find that the combination of (i) the lower baryon fraction in dSphs compared to the cosmic mean and (ii) the concentration of baryons in the inner part of the MW halo can go a long way towards explaining the observed circular velocity profiles of dSphs. We perform controlled numerical simulations that mimic the effects of baryons. From these we find that the blowing away of baryons by ram pressure, when the dwarfs fall into larger galaxies, decreases the circular velocity profile of the satellite and reduces the density in the central ∼200-500 pc by a factor of (1 − f b ) 4 ≈ 0.5, where f b is the cosmological fraction of baryons. Additionally, the enhanced baryonic mass in the central regions of the parent galaxy generates tidal forces, which are larger than those experienced by subhaloes in traditional N-body simulations. Increased tidal forces substantially alter circular velocity profiles for satellites with pericentres less than 50 kpc. We show that these two effects are strong enough to bring the predictions of subhaloes from CDM simulations into agreement with the observed structure of MW dSphs, regardless of the details of the baryonic processes.

Despite recent success in forming realistic present-day galaxies, simulations still form the bulk of their stars earlier than observations indicate. We investigate the process of stellar mass assembly in low-mass field galaxies, a dwarf and a typical spiral, focusing on the effects of radiation from young stellar clusters on the star formation histories. We implement a novel model of star formation (SF) with a deterministic low efficiency per free-fall time, as observed in molecular clouds. Stellar feedback is based on observations of star-forming regions, and includes radiation pressure from massive stars, photoheating in H II regions, supernovae, and stellar winds. We find that stellar radiation has a strong effect on the formation of lowmass galaxies, especially at z > 1, where it efficiently suppresses SF by dispersing cold and dense gas, preventing runaway growth of the stellar component. This behaviour is evident in a variety of observations but had so far eluded analytical and numerical models without radiation feedback. Compared to supernovae alone, radiation feedback reduces the SF rate by a factor of ∼ 100 at z 2, yielding rising SF histories which reproduce recent observations of Local Group dwarfs. Stellar radiation also produces bulgeless spiral galaxies and may be responsible for excess thickening of the stellar disc. The galaxies also feature rotation curves and baryon fractions in excellent agreement with current data. Lastly, the dwarf galaxy shows a very slow reduction of the central dark matter density caused by radiation feedback over the last ∼ 7 Gyr of cosmic evolution.

PDS 144 is a pair of Herbig Ae stars that are separated by 5.35 ′′ on the sky. It has previously been shown to have an A2Ve Herbig Ae star viewed at 83 ○ inclination as its northern member and an A5Ve Herbig Ae star as its southern member. Direct imagery revealed a disk occulting PDS 144 N -the first edge-on disk observed around a Herbig Ae star. The lack of an obvious disk in direct imagery suggested PDS 144 S might be viewed face-on or not physically associated with PDS 144 N. Multi-epoch HST imagery of PDS 144 with a 5 yr baseline demonstrates PDS 144 N & S are comoving and have a common proper motion with TYC 6782-878-1. TYC 6782-878-1 has previously been identified as a member of Upper Sco sub-association A at d = 145 ± 2 pc with an age of 5 -10 Myr. Ground-based imagery reveals jets and a string of HH knots extending 13 ′ (possibly further) which are aligned to within 7°± 6°on the sky. By combining proper motion data and the absence of a dark mid-plane with radial velocity data, we measure the inclination of PDS 144 S to be i = 73°± 7°. The radial velocity of the jets from PDS 144 N & S indicates they, and therefore their disks, are misaligned by 25°± 9°. This degree of misalignment is similar to that seen in T-Tauri wide binaries.

We report the discovery and orbital determination of 14 trans-Neptunian objects (TNOs) from the ESSENCE Supernova Survey difference imaging data set. Two additional objects discovered in a similar search of the SDSS-II Supernova Survey database were recovered in this effort. ESSENCE repeatedly observed fields far from the solar system ecliptic (Ϫ21Њ ! b ! Ϫ5Њ), reaching limiting magnitudes per observation of and I ≈ 23.1 R ≈ . We examine several of the newly detected objects in detail, including 2003 UC 414 , which orbits entirely 23.7 between Uranus and Neptune and lies very close to a dynamical region that would make it stable for the lifetime of the solar system. 2003 SS 422 and 2007 TA 418 have high eccentricities and large perihelia, making them candidate members of an outer class of TNOs. We also report a new member of the "extended" or "detached" scattered disk, 2004 VN 112 , and verify the stability of its orbit using numerical simulations. This object would have been visible to ESSENCE for only ∼2% of its orbit, suggesting a vast number of similar objects across the sky. We emphasize that off-ecliptic surveys are optimal for uncovering the diversity of such objects, which in turn will constrain the history of gravitational influences that shaped our early solar system.