We make detailed theoretical predictions for the assembly properties of the Local Group (LG) in the standard LambdaCDM cosmological model. We use three cosmological N-body dark matter simulations from the CLUES project, which are designed to reproduce the main dynamical features of the matter distribution down to the scale of a few Mpc around the LG. Additionally, we use the results of an unconstrained simulation with a sixty times larger volume to calibrate the influence of cosmic variance. We characterize the Mass Aggregation History (MAH) for each halo by three characteristic times, the formation, assembly and last major merger times. A major merger is defined by a minimal mass ratio of 10:1. We find that the three LGs share a similar MAH with formation and last major merger epochs placed on average \approx 10 - 12 Gyr ago. Between 12% and 17% of the halos in the mass range 5 x 10^11 Msol/h < M_h < 5 x 10^12 Msol/h have a similar MAH. In a set of pairs of halos within the same mass range, a fraction of 1% to 3% share similar formation properties as both halos in the simulated LG. An unsolved question posed by our results is the dynamical origin of the MAH of the LGs. The isolation criteria commonly used to define LG-like halos in unconstrained simulations do not narrow down the halo population into a set with quiet MAHs, nor does a further constraint to reside in a low density environment. The quiet MAH of the LGs provides a favorable environment for the formation of disk galaxies like the Milky Way and M31. The timing for the beginning of the last major merger in the Milky Way dark matter halo matches with the gas rich merger origin for the thick component in the galactic disk. Our results support the view that the specific large and mid scale environment around the Local Group play a critical role in shaping its MAH and hence its baryonic structure at present.

Contradicting results have been reported in the literature with respect to the performance of the numerical techniques employed for the study of supersonic turbulence. We aim at characterising the performance of different particle-based and grid-based techniques on the modelling of decaying supersonic turbulence. Four different grid codes (ENZO, FLASH, TVD, ZEUS) and three different SPH codes (GADGET, PHANTOM, VINE) are compared. We additionally analysed two calculations denoted as PHANTOM A and PHANTOM B using two different implementations of artificial viscosity. Our analysis indicates that grid codes tend to be less dissipative than SPH codes, though details of the techniques used can make large differences in both cases. For example, the Morris & Monaghan viscosity implementation for SPH results in less dissipation (PHANTOM B and VINE versus GADGET and PHANTOM A). For grid codes, using a smaller diffusion parameter leads to less dissipation, but results in a larger bottleneck effect (our ENZO versus FLASH runs). As a general result, we find that by using a similar number of resolution elements N for each spatial direction means that all codes (both grid-based and particle-based) show encouraging similarity of all statistical quantities for isotropic supersonic turbulence on spatial scales k<N/32 (all scales resolved by more than 32 grid cells), while scales smaller than that are significantly affected by the specific implementation of the algorithm for solving the equations of hydrodynamics. At comparable numerical resolution, the SPH runs were on average about ten times more computationally intensive than the grid runs, although with variations of up to a factor of ten between the different SPH runs and between the different grid runs. (abridged)

We have performed time-dependent numerical simulations of the interstellar medium (ISM) which account for galactic shear and magnetic fields, vertical gravity, and a radiative cooling function for atomic gas. This allows us to study the magnetorotational instability (MRI) in cloudy, vertically-stratified disks. As in previous unstratified models, we find that thermal instability interacts with MRI-driven turbulence and galactic shear to produce a network of cold, dense, filamentary clouds embedded in a warm diffuse ambient medium. This structure strongly resembles the morphology of HI gas observed in the 21 cm line. There is significant thermally-unstable gas, but the density and temperature distributions retain the twin peaks of the classical two-phase ISM. We analyze the vertical distributions of density and various pressure terms, and address what supports the ISM vertically. Turbulent velocities of the cold gas increase as the cold mass fraction decreases, but are generally low ~1-3 km/s near the midplane; they increase to > 5 km/s at high z. Finally, we argue that in the outer parts of galactic disks, MRI is likely able to prevent the development of self-gravitating instabilities and hence suppress star formation, even if cold gas is present.

This paper reports on 3D numerical simulations of dynamics and thermodynamics in the diffuse ISM. Our models are local, account for sheared galactic rotation, magnetic fields, and realistic cooling, and resolve scales ~1-200 pc. This combination permits the study of quasi-steady turbulence in a cloudy medium representing the warm/cold atomic ISM. Turbulence is driven by the magnetorotational instability (MRI); our models are the first to study the saturated state of MRI under strongly inhomogeneous conditions, with cloud/intercloud contrasts of about 100. For volume-averaged densities n=0.25 - 4 cm^3, the saturated-state velocity dispersion ranges from 8-1 km/s, scaling as n^-0.77. The MRI is therefore likely quite important in driving turbulence in low-density ISM regions, both at high latitudes in the inner Galaxy,and throughout the far outer Galaxy. The MRI may be key to suppressing star formation at large radii in spiral galaxies, where without MRI-driven turbulence, a thin, cold layer could form and become gravitationally instable. Our presentation also includes analysis of thermal structure (temperature and density PDFs are broadened and include extreme departures from equilibrium, but retain the classical two-phase bimodal character); results on the distribution of clump masses (the mass spectrum peaks at about 100 solar masses); comparisons of saturated-state MRI scalings with single-phase simulation results (we find B^2 is independent of n); and examples of synthetic HI line profile maps (showing that physical clumps are not easily distinguished in velocity components, and vice versa).

The structure and dynamics of diffuse gas in the Milky Way and other disk galaxies may be strongly influenced by thermal and magnetorotational instabilities (TI and MRI) on scales of about 1-100 pc. We initiate a study of these processes, using two-dimensional numerical hydrodynamic and magnetohydrodynamic (MHD) simulations with conditions appropriate for the atomic interstellar medium (ISM). We demonstrate, consistent with previous work, that nonlinear development of ``pure TI'' produces a network of filaments that condense into cold clouds at their intersections, yielding a distinct two-phase warm/cold medium within about 20 Myr. TI-driven turbulent motions of the clouds saturate at subsonic amplitudes for uniform initial P/k=2000 K cm^-3. MRI has previously been studied in near-uniform media; our simulations include both TI+MRI models, which begin from uniform-density conditions, and cloud+MRI models, which begin with a two-phase cloudy medium. Both the TI+MRI and cloud+MRI models show that MRI develops within a few galactic orbital times, just as for a uniform medium. The mean separation between clouds can affect which MRI mode dominates the evolution. Provided intercloud separations do not exceed half the MRI wavelength, we find the MRI growth rates are similar to those for the corresponding uniform medium. This opens the possibility, if low cloud volume filling factors increase MRI dissipation times compared to those in a uniform medium, that MRI-driven motions in the ISM could reach amplitudes comparable to observed HI turbulent linewidths.

Spectroscopic measurements of Mira variable stars, as a function of phase, probe the stellar atmospheres and underlying pulsation mechanisms. For example, measuring variations in TiO, VO, and ZrO with phase can be used to help determine whether these molecular species are produced in an extended region above the layers where Balmer line emission occurs or below this shocked region. Using the same methods, the Balmer-line increment, where the strongest Balmer line at phase zero is H-delta and not H-alpha can be measured and explanations tested, along with another peculiarity, the absence of the H-epsilon line in the spectra of Miras when other Balmer lines are strong. We present new spectra covering the spectral range from 6200 Angstroms to 9000 Angstroms of 20 Mira variables. A relationship between variations in the CaII IR triplet and H-alpha as a function of phase support the hypothesis that H-epsilon's observational characteristics result from an interaction of H-epsilon photons with the CaII H line. New periods and epochs of variability are also presented for each star.