In performing cosmological N-body simulations, it is widely appreciated that the growth of structure on the largest scales within a simulation box will be inhibited by the finite size of the simulation volume. Following ideas set forth in Seto (1999), this paper shows that standard (a.k.a. 1-loop) cosmological perturbation theory (SPT) can be used to predict, in an approximate way, the deleterious effect of the box scale on the power spectrum of density fluctuations in simulation volumes. Alternatively, this approach can be used to quickly estimate post facto the effect of the box scale on power spectrum results from existing simulations. In this way SPT can help determine whether larger box sizes or other more-sophisticated methods are needed to achieve a particular level of precision for a given application (e.g. simulations to measure the non-linear evolution of baryon acoustic oscillations). I focus on SPT in this note and show that its predictions differ only by about a factor of two or less from the measured suppression inferred from both powerlaw and $\Lambda$CDM $N$-body simulations. It should be possible to improve the accuracy of these predictions through using more-sophisticated perturbation theory models. An appendix compares power spectrum measurements from the powerlaw simulations at outputs where box-scale effects are minimal to perturbation theory models and previously-published fitting functions. These power spectrum measurements are included with this paper to aid efforts to develop new perturbation theory models.

In setting up initial conditions for ensembles of cosmological N-body simulations there are, fundamentally, two choices: either maximizing the correspondence of the initial density field to the assumed fourier-space clustering or, instead, matching to real-space statistics and allowing the DC mode (i.e. overdensity) to vary from box to box as it would in the real universe. As a stringent test of both approaches, I perform ensembles of simulations using power law and a "powerlaw times a bump" model inspired by baryon acoustic oscillations (BAO), exploiting the self-similarity of these initial conditions to quantify the accuracy of the matter-matter two-point correlation results. The real-space method, which was originally proposed by Pen 1997 and implemented by Sirko 2005, performed well in producing the expected self-similar behavior and corroborated the non-linear evolution of the BAO feature observed in conventional simulations, even in the strongly-clustered regime (sigma8 >= 1). In revisiting the real-space method championed by Sirko 2005, it was also noticed that this earlier study overlooked an important integral constraint correction to the correlation function in results from the conventional approach that can be important in LambdaCDM simulations with Lbox =< 1 Gpc/h and on scales r >= Lbox / 10. Rectifying this shows that the fourier space and real space methods are about equally accurate and efficient for modeling the evolution and growth of the correlation function, contrary to previous claims. An appendix provides a useful independent-of-epoch analytic formula for estimating the importance of the integral constraint bias on correlation function measurements in LambdaCDM simulations.

Motivated by cosmological surveys that demand accurate theoretical modeling of the baryon acoustic oscillation (BAO) feature in galaxy clustering, we analyze N-body simulations in which a BAO-like gaussian bump modulates the linear theory correlation function \xi_L(r)=(r_0/r)^{n+3} of an underlying self-similar model with initial power spectrum P(k)=A k^n. These simulations test physical and analytic descriptions of BAO evolution far beyond the range of most studies, since we consider a range of underlying power spectra (n=-0.5, -1, -1.5) and evolve simulations to large effective correlation amplitudes (equivalent to \sigma_8=4-12 for r_bao = 100 Mpc/h). In all cases, non-linear evolution flattens and broadens the BAO bump in \xi(r) while approximately preserving its area. This evolution resembles a "diffusion" process in which the bump width \sigma_bao is the quadrature sum of the linear theory width and a length proportional to the rms relative displacement \Sigma_pair(r_bao}) of particle pairs separated by r_bao. For n=-0.5 and n=-1, we find no detectable shift of the location of the BAO peak, but the peak in the n=-1.5 model shifts steadily to smaller scales, following r_peak}/r_bao = 1-1.08(r_0/r_bao)^{1.5}. The "SimpleRG" perturbation theory scheme and, to a lesser extent, standard 1-loop perturbation theory are fairly successful at explaining the non-linear evolution of the fourier power spectrum of our models. Analytic models also explain why the \xi(r) peak shifts much more for n=-1.5 than for n >= -1, though no ab initio model we have examined reproduces all of our numerical results. Simulations with L_box = 10 r_bao} and L_box = 20 r_bao yield consistent results for \xi(r) at the BAO scale, provided one corrects for the integral constraint imposed by the uniform density box.

The paucity of observed dwarf galaxies in the Local Group relative to the abundance of predicted dark matter halos remains one of the greatest puzzles of the LCDM paradigm. Solving this puzzle now requires not only matching the numbers of objects but also understanding the details of their star formation histories. We present a summary of such histories derived from the HST data using the color-magnitude diagram fitting method. To reduce observational uncertainties, we condense the data into five cumulative parameters -- the fractions of stellar mass formed in the last 1, 2, 5, and 10 Gyr, and the mean stellar age. We interpret the new data with a phenomenological model based on the mass assembly histories of dark matter halos and the Schmidt law of star formation. The model correctly predicts the radial distribution of the dwarfs and the fractions of stars formed in the last 5 and 10 Gyr. However, in order to be consistent with the observations, the model requires a significant amount of recent star formation in the last 2 Gyr. Within the framework of our model, this prolonged star formation can be achieved by adding a stochastic variation in the density threshold of the star formation law. The model results are not sensitive to late gas accretion, the slope of the Schmidt law, or the details of cosmic reionization. A few discrepancies still remain: our model typically predicts too large stellar masses, only a modest population of ultra-faint dwarfs, and a small number of dwarfs with anomalously young stellar populations. Nevertheless, the observed star formation histories of Local Group dwarfs are generally consistent the expected star formation in cold dark matter halos.

We present two-dimensional inviscid hydrodynamic simulations of a protoplanetary disk with an embedded planet, emphasizing the evolution of potential vorticity (the ratio of vorticity to density) and its dependence on numerical resolutions. By analyzing the structure of spiral shocks made by the planet, we show that progressive changes of the potential vorticity caused by spiral shocks ultimately lead to the excitation of a secondary instability. We also demonstrate that very high numerical resolution is required to both follow the potential vorticity changes and identify the location where the secondary instability is first excited. Low-resolution results are shown to give the wrong location. We establish the robustness of a secondary instability and its impact on the torque onto the planet. After the saturation of the instability, the disk shows large-scale non-axisymmetry, causing the torque on the planet to oscillate with large amplitude. The impact of the oscillating torque on the protoplanet's migration remains to be investigated.