As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

We investigate planetesimal accretion via a direct N-body simulation of an annulus at 1 AU orbiting a 1 $M_{\odot}$ star. The planetesimal ring, which initially contains N = $10^6$ bodies is evolved into the oligarchic growth phase. Unlike previous lower resolution studies, we find that the mass distribution of planetesimals develops a bump at intermediate mass after the oligarchs form. This feature marks a boundary between growth modes. The smallest planetesimals are packed tightly enough together to populate mean motion resonances with the oligarchs, which heats the small bodies, enhancing their growth. If we depopulate most of the resonances by decreasing the width of the annulus, this effect becomes weaker. To clearly demonstrate the dynamics driving these growth modes, we also examine the evolution of a planetary embryo embedded in an annulus of collisionless planetesimals. In this case, we find that the resonances push planetesimals away from the embryo, decreasing the surface density of the bodies adjacent to the embryo. This effect only occurs when the annulus is wide enough and the mass resolution of the planetesimals is fine enough to populate the resonances. The bump we observe in the mass distribution resembles the 100 km power law break seen in the size distribution of asteroid belt objects. Although the bump produced in our simulations occurs at a size larger than 100 km, we show that the bump location is sensitive to the initial planetesimal mass, which implies that this feature is potentially useful for constraining planetesimal formation models.

Over the solar-activity cycle, there are extended periods where significant discrepancies occur between the spacecraft-observed total (unsigned) open magnetic flux and that determined from coronal models. In this article, the total open heliospheric magnetic flux is computed using two different methods and then compared with results obtained from in-situ interplanetary magnetic-field observations. The first method uses two different types of photospheric magnetic-field maps as input to the Wang Sheeley Arge (WSA) model: i) traditional Carrington or diachronic maps, and ii) Air Force Data Assimilative Photospheric Flux Transport model synchronic maps. The second method uses observationally derived helium and extreme-ultraviolet coronal-hole maps overlaid on the same magnetic-field maps in order to compute total open magnetic flux. The diachronic and synchronic maps are both constructed using magnetograms from the same source, namely the National Solar Observatory Kitt Peak Vacuum Telescope and Vector Spectromagnetograph. The results of this work show that the total open flux obtained from observationally derived coronal holes agrees remarkably well with that derived from WSA, especially near solar minimum. This suggests that, on average, coronal models capture well the observed large-scale coronal-hole structure over most of the solar cycle. Both methods show considerable deviations from total open flux deduced from spacecraft data, especially near solar maximum, pointing to something other than poorly determined coronal-hole area specification as the source of these discrepancies.

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.

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.

We present the results of a campaign to search for subdwarf O (sdO) star pulsators among bright field stars. The motivation for this project is the recent discovery by Randall et al. (2011), of four rapidly pulsating sdO stars in the globular cluster Omega Cen, with Teff near 50,000 K, 5.4 < log g < 6.0, and hydrogen-rich atmospheres. The only previously known sdO pulsator is significantly hotter at 68,500 K and log g = 6.1. All of the sdO pulsators identified so far are fainter than 17.4 in the V band and, thus, are poor candidates for an in-depth follow-up with asteroseismology. We therefore obtained high S/N light curves and spectroscopy for a number of field sdO stars to attempt to discover bright counterparts to these stars, particularly the Omega Cen pulsators. Our primary sample consisted of 19 sdO stars with hydrogen-rich atmospheres, log N(He)/N(H) < -1.0, effective temperatures in the range 40,000 K < Teff < 67,000 K, and surface gravities 5.3 < log g < 6.1. We also observed 17 additional helium-rich sdO stars with log N(He)/N(H) > -0.1 and similar temperatures and gravities. To date, we have found no detectable pulsations at amplitudes above 0.08% (4 times the mean noise level) in any of the 36 field sdO stars that we observed. The presence of pulsations in Omega Cen sdO stars and their apparent absence in seemingly comparable field sdO stars is perplexing. While very suggestive, the significance of this result is difficult to assess more completely right now due to remaining uncertainties about the temperature width and purity of the Omega Cen instability strip and the existence of any sdO pulsators with weaker amplitudes than the current detection limit in globular clusters.